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 brief report, we investigate the existence of 4-dimensional static spherically symmetric black holes (BHs) in the Einstein-complex-scalar-Gauss-Bonnet (EcsGB) gravity with an arbitrary potential $V(\phi)$ and a coupling $f(\phi)$ between the scalar field $\phi$ and the Gauss-Bonnet (GB) term. We find that static regular BH solutions with complex scalar hairs do not exist. This conclusion does not depend on the coupling between the GB term and the scalar field, nor on the scalar potential $V(\phi)$ and the presence of a cosmological constant $\Lambda$ (which can be either positive or negative), as longer as the scalar field remains complex and is regular across the horizon.

With the arrival of the era of gravitational wave astronomy, the strong gravitational field regime will be explored soon in various aspects. In this article, we provide a general review over cylindrical systems in Einstein's theory of general relativity. In particular, we first review the general properties, both local and global, of several important solutions of Einstein's field equations, including the Levi-Civita and Lewis solutions and their extensions to include the cosmological constant and matter fields, and pay particular attention to properties that represent the generic features of the theory, such as the formation of the observed extragalactic jets and gravitational Faraday rotation. We also review studies of cylindrical wormholes, gravitational collapse and Hoop conjecture, and polarizations of gravitational waves. In addition, by rigorously defining cylindrically symmetric spacetimes, we clarify various (incorrect) claims existing in the literature, regarding to the generality of such spacetimes.

The physics of the dark sector has remained one of the controversial areas of modern cosmology at present and hence it naturally attracts massive attention to the scientific community. With the developments of the astronomical data, the physics of the dark sector is becoming much more transparent than it was some twenty years back. The detection of gravitational waves (GWs) has now opened a cluster of possibilities in the cosmological regime. Being motivated by the detection of GWs and its possible impact on the physics of dark matter and dark energy, in this work we focus on the interacting dark energy models. Assuming the simplest possibility in which the vacuum energy with equation-of-state $w_x =-1$ is allowed to interact with the pressureless dark matter, we have extracted the constraints of the cosmological parameters. We find that the addition of the GWs data to the CMB measurements significantly improves up to a factor 4 of the parameters space, and up to a factor 2 for the full combination of current cosmological datasets, namely CMB+BAO+Pantheon+RSD+R16+CC+WL. The most affected parameters by the inclusion of the GWs are $\Omega_ch^2$, $\theta_{MC}$, $\xi$, and the derived parameters $\Omega_{m0}$, $\sigma_8$ and $H_0$.

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.

In this paper, we construct an effective rotating loop quantum black hole (LQBH) solution, starting from the spherical symmetric LQBH by applying the Newman-Janis algorithm modified by Azreg-A\"{i}nou's non-complexification procedure, and study the effects of loop quantum gravity (LQG) on its shadow. Given the rotating LQBH, we discuss its horizon, ergosurface, and regularity as $r \to 0$. Depending on the values of the specific angular momentum $a$ and the polymeric function $P$ arising from LQG, we find that the rotating solution we obtained can represent a regular black hole, a regular extreme black hole, or a regular spacetime without horizon (a non-black-hole solution). We also study the effects of LQG and rotation and show that, in addition to the specific angular momentum, the polymeric function also causes deformations in the size and shape of the black hole shadow. Interestingly, for a given value of $a$ and inclination angle $\theta_0$, the apparent size of the shadow monotonically decreases, and the shadow gets more distorted with increasing $P$. We also consider the effects of $P$ on the deviations from the circularity of the shadow and find that the deviation from circularity increases with increasing $P$ for fixed values of $a$ and $\theta_0$. Additionally, we explore the observational implications of $P$ in comparing with the latest Event Horizon Telescope (EHT) observation of the supermassive black hole, M87*. The connection between the shadow radius and quasinormal modes in the eikonal limit as well as the deflection of massive particles are also considered.

We investigate the consequences of different regularizations and ambiguities in loop cosmological models on the predictions in the scalar and tensor primordial spectrum of the cosmic microwave background using the dressed metric approach. Three models, standard loop quantum cosmology (LQC), and two modified loop quantum cosmologies (mLQC-I and mLQC-II) arising from different regularizations of the Lorentzian term in the classical Hamiltonian constraint are explored for chaotic inflation in spatially-flat Friedmann-Lema\^itre-Robertson-Walker (FLRW) universe. In each model, two different treatments of the conjugate momentum of the scale factor are considered. The first one corresponds to the conventional treatment in dressed metric approach, and the second one is inspired from the hybrid approach which uses the effective Hamiltonian constraint. For these two choices, we find the power spectrum to be scale-invariant in the ultra-violet regime for all three models, but there is at least a $10\%$ relative difference in amplitude in the infra-red and intermediate regimes. All three models result in significant differences in the latter regimes. In mLQC-I, the magnitude of the power spectrum in the infra-red regime is of the order of Planck scale irrespective of the ambiguity in conjugate momentum of the scale factor. The relative difference in the amplitude of the power spectrum between LQC and mLQC-II can be as large as $50\%$ throughout the infra-red and intermediate regimes. Differences in amplitude due to regularizations and ambiguities turn out to be small in the ultra-violet regime.

In this work, we revisit the dynamics of pre-inflationary universe with a family of $\alpha-$attractor potentials, in the framework of loop quantum cosmology, in which the big bang singularity is generically resolved purely with quantum geometric effects, and replaced by a quantum bounce. At the bounce, the background evolution is divided into two distinct classes, the first is dominated by the kinetic energy of the inflaton field and the second by the potential energy. In both classes, we find the physically viable initial conditions numerically that provide not only the slow-roll inflation, but also sufficient e-folds to be compatible with observations. In the entire range of kinetic energy dominated initial conditions (except some subsets of Models 2 and 4), the background evolution prior to reheating is always split into three different phases: bouncing, transition and slow-roll inflation. In the bouncing phase, the numerical evolution of the scale factor is independent not only of the initial data, but also the inflationary potentials, as long as it is dominated by the kinetic energy, and can be well approximated by an analytical solution, whereas in the potential energy dominated case, such approximated results do not exist. Moreover, we study the phase space analysis for a class of $\alpha-$attractor potentials, and discuss the phase space trajectories for physically viable initial conditions of the inflaton field.

Similar to the Schwarzschild coordinates for spherical black holes, the Baldwin, Jeffery and Rosen (BJR) coordinates for plane gravitational waves are often singular, and extensions beyond such singularities are necessary, before studying asymptotic properties of such spacetimes at the null infinity of the plane, on which the gravitational waves propagate. The latter is closely related to the studies of memory effects and soft graviton theorems. In this paper, we point out that in the BJR coordinates all the spacetimes are singular physically at the focused point $u = u_s$, except for the two cases: (1) $\alpha =1/2, \; \forall \; \chi_n$; and (2) $\alpha =1, \; \chi_i = 0\; (i = 1, 2, 3)$, where $\chi_n$ are the coefficients in the expansion $\chi \equiv \left[{\mbox{det}}\left(g_{ab}\right) \right]^{1/4} = \left(u - u_s\right)^{\alpha}\sum_{n = 0}^{\infty}\chi_n \left(u - u_s\right)^n$ with $\chi_0 \not= 0$, the constant $\alpha \in (0, 1]$ characterizes the strength of the singularities, and $g_{ab}$ denotes the reduced metric on the two-dimensional plane orthogonal to the propagation direction of the wave. Therefore, the hypersurfaces $u= u_s$ already represent the boundaries of such spacetimes, and the null infinity does not belong to them. As a result, they cannot be used to study properties of plane gravitational waves at null infinities, including memory effects and soft graviton theorems.

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.

Recently the background independent nonperturbative quantization technique has been extended to theories of gravity and the corresponding quantum effective cosmology has been derived, which provides us with necessary avenue to explore the pre-inflationary dynamics. Brans-Dicke loop quantum cosmology (BDLQC) is one such theory whose effective background dynamics is considered in this article. Starting with a quantum bounce, we explore the pre-inflationary dynamics of a universe sourced by a scalar field with the Starobinsky potential in BDLQC. Our study is based on the idea that though Einstein's and Jordan's frames are classically equivalent up to a conformal transformation in Brans-Dick theory, it is not true after quantization. Taking the Jordan frame as the physical one we explore in detail the bouncing scenario which is followed by a phase of slow roll inflation. The three phases of the evolution of the universe namely, bouncing, transition from quantum bounce to classical universe and the slow roll inflation, are noted for an initially kinetic energy dominated bounce. In addition, to be consistent with observational data, we also find out the allowed data space of initial conditions that would produce 60 e-folds of expansion during the slow roll inflation.

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 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.

The effective field theory (EFT) of inflation provides an essential picture to explore the effects of the unknown high energy physics in the single scalar field inflation models. For a generic EFT of inflation, possible high energy corrections to simple slow-roll inflation can modify both the propagating speed and dispersion relations of the cosmological scalar and tensor perturbations. With the arrival of the era of precision cosmology, it is expected that these high energy corrections become more important and have to be taken into account in the analysis with future precise observational data. In this paper we study the observational predictions of the EFT of inflation by using the third-order uniform asymptotic approximation method. We calculate explicitly the primordial power spectra, spectral indices, running of the spectral indices for both scalar and tensor perturbations, and the ratio between tensor and scalar spectra. These expressions are all written in terms of the Hubble flow parameters and the flow of four new slow-roll parameters and expanded up to the next-to-leading order in the slow-roll expansions so they represent the most accurate results obtained so far in the literature. The flow of the four new slow-roll parameters, which arise from the four new operators introduced in the action of the EFT of inflation, can affect the primordial perturbation spectra at the leading-order and the corresponding spectral indices at the next-to-leading order.

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 effective field theory (EFT) of inflation provides a natural framework to study the new physical effects on primordial perturbations. Recently a healthy extension of the EFT of inflation with high-order operators has been proposed, which avoids ghosts and meanwhile leads to a nonlinear dispersion relation of the scalar perturbations. This paper is devoted to studying the effects of these high-order operators by using the uniform asymptotic approximation method. In particular, we first construct the approximate analytical solution to the mode function of the scalar perturbations. Because of the presence of the high-order operators, the perturbation modes usually experience a period of non-adiabatic evolution before they cross the Hubble radius, which could lead to the production of excited states and modifications of the primordial perturbation spectrum. However, we show that the modified power spectrum is still nearly scale-invariant and the presence of the high-order operators can only affect the overall amplitude of the spectrum. In particular, after showing explicitly the impact of these new effects on particle production rate and perturbation spectrum, we explore their origin in detail.

We study the evolution of spatially flat Friedmann-Lema\^{i}tre-Robertson-Walker universe for chaotic and Starobinsky potentials in the framework of modified loop quantum cosmologies. These models result in a non-singular bounce as in loop quantum cosmology, but with far more complex modified Friedmann dynamics with higher order than quadratic terms in energy density. For the kinetic energy dominated bounce, we obtain analytical solutions using different approximations and compare with numerical evolution for various physical variables. The relative error turns out to be less than $0.3\%$ in the bounce regime for both of the potentials. Generic features of dynamics, shared with loop quantum cosmology, are established using analytical and numerical solutions. Detailed properties of three distinct phases in dynamics separating bounce regime, transition stage and inflationary phase are studied. For the potential energy dominated bounce, we qualitatively describe its generic features and confirm by simulations that they all lead to the desired slow-roll phase in the chaotic inflation. However, in the Starobinsky potential, the potential energy dominated bounce cannot give rise to any inflationary phase. Finally, we compute the probability for the desired slow-roll inflation to occur in the chaotic inflation and as in loop quantum cosmology, find a very large probability for the universe to undergo inflation.

We consider a varieties of quintessence scalar field models in a homogeneous and isotropic geometry of the universe with zero spatial curvature aiming to provide stringent constraints using a series of cosmological data sets, namely, the cosmic microwave background observations (CMB), baryon acoustic oscillations (BAO), joint light curve analysis (JLA) from supernovae type Ia, redshift space distortions (RSD), and the cosmic chronometers (CC). From the qualitative evolution of the models, we find all of them are able to execute a fine transition from the past decelerating phase to the presently accelerating expansion where in addition, the equation of state of the scalar field (also the effective equation of state) might be close to that of the $\Lambda$CDM cosmology depending on its free parameters. From the observational analyses, we find that the scalar field parameters are unconstrained irrespective of all the observational datasets. In fact, we find that the quintessence scalar field models are pretty much determined by the CMB observations since any of the external datasets such as BAO, JLA, RSD, CC does not add any constraining power to CMB. Additionally, we observe a strong negative correlation between the parameters $H_0$ (present value of the Hubble parameter), $\Omega_{m0}$ (density parameter for the matter sector, i.e., cold dark matter plus baryons) exists, while no correlation between $H_0$, and $\sigma_8$ (amplitude of the matter fluctuation) are not correlated. We also comment that the present models are unable to reconcile the tension on $H_0$. Finally, we conclude our work with the Bayesian analyses which report that the non-interacting $\Lambda$CDM model is preferred over all the quintessence scalar field models.

In this paper, we investigate a scenario in which late time cosmic acceleration might arise due to coupling between dark matter and baryonic matter without resorting to dark energy or large scale modification of gravity associated with extra degrees of freedom. The scenario can give rise to late time acceleration in Jordan frame and no acceleration in Einstein frame - \textit{generic modification of gravity} caused by disformal coupling. Using a simple parametrization of the coupling function, in maximally disformal case, we constrain the model parameters by using the age constraints due to globular cluster data. We also obtain observational constraints on the parameters using $H(z)+SNIa+BAO$ data sets. In this case, we distinguish between phantom and non phantom acceleration and show that the model can give rise to phantom behavior in a narrow region of parameter space.

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.

While loop quantum cosmology (LQC) predicts a robust quantum bounce of the background evolution of a Friedmann-Robertson-Walker (FRW) spacetime prior to the standard slow-roll inflation, whereby the big bang singularity is resolved, there are several different quantization procedures to cosmological perturbations, for instance, {\em the deformed algebra, dressed metric, and hybrid quantizations}. This paper devotes to study the quantum bounce effects of primordial perturbations in the hybrid approach. The main discrepancy of this approach is the effective positive mass at the quantum bounce for the evolution of the background that is dominated by the kinetic energy of the inflaton field at the bounce, while this mass is always nonpositive in the dressed metric approach. It is this positivity of the effective mass that violates the adiabatic evolution of primordial perturbations at the initial moments of the quantum bounce. With the assumption that the evolution of the background is dominated by the kinetic energy of the inflaton at the bounce, we find that the effective potentials for both scalar and tensor perturbations can be well approximately described by a P\"{o}schl-Teller (PT) potential, which allows us to find analytical solutions of perturbations, and from these analytical expressions we are able to study the non-adiabatic evolution of primordial perturbations in details. In particular, we derive their quantum bounce effects and investigate their observational constraints. In addition, the impacts of quantum bounce effects on the non-Gaussinity and their implication on the explanations of observed power asymmetry in CMB have also been explored.

Highly precise pulsar timing is very important for understanding the nature of a neutron star, and it can even be used to detect gravitational waves. Unfortunately, the accuracy of the pulsar timing is seriously affected by the spin-down irregularities of pulsars.

Highly precise pulsar timing is very important for understanding the nature of neutron stars and the implied physics, even being used to detect directly gravitational waves. Unfortunately, the accuracy of the pulsar timing is seriously affected by the spin-down irregularities, such as spin fluctuations with a manifestation of low frequency structures (the so-called red noise processes), and various activities of magnetospheres. Except from the random timing noise, significant periodic or quasi-periodic components could also be found in some timing residuals of pulsars, indicating the presence of unmodelled deterministic effects. The physical origins of these effects still remain unexplained. In this paper, we suggest a new mechanism involving the de Haas-van Alphen magnetic oscillation, which could trigger the observed low frequency structures. We find that the proposed magnetic oscillation period is about 1-$10^{2}$ yr, which is about $10^{-4}$ times as long as the classical characteristic time scale of interior magnetic field evolution for a normal neutron star. Due to the magnetic oscillation, we estimate both the braking index range between $10^{-5}$ and $10^{5}$ and the residuals range between 4 to 906 ms for some specific samples of quasi-periodic pulsars. Those are consistent with the pulsar timing observations. Nonlinear phenomena of de Haas-van Alphen effect and the Condon domain structure are considered. Avalanche-like magnetized process of domain wall motion (rapidly magnetic energy release) could associated with some special emitting manifestations of neuron star radiation.

In this paper, we systematically study spacetimes of gravitational plane waves in Einstein-aether theory. Due to the presence of the timelike aether vector field, now the problem in general becomes overdetermined. In particular, for the linearly polarized plane waves, there are five independent vacuum Einstein-aether field equations for three unknown functions. Therefore, solutions exist only for particular choices of the four free parameters $c_{i}$'s of the theory. We find that there exist eight cases, in two of which any form of gravitational plane waves can exist, similar to that in general relativity, while in the other six cases, gravitational plane waves exist only in particular forms. Beyond these eight cases, solutions either do not exist or are trivial (simply representing a Minkowski spacetime with a constant or dynamical aether field.).

We numerically study the gravitational collapse of a massless scalar field with spherical symmetry in Einstein-aether theory, and show that apparent, spin-0 and dynamical universal horizons (dUHs) can be all formed. The spacetime and the aether field are well-behaved and regular, including regions nearby these horizons (but away from the center of spherical symmetry). The spacetime outside the apparent and spin-0 horizons settles down to a static configuration, and some of such resulting static black holes were already found numerically in the literature. On the other hand, the proper distance of the outermost dUH from the apparent (or spin-0) horizon keeps increasing on $\mbox{\ae}$ther-orthogonal time slices. This indicates that the outermost dUH is evolving into the causal boundary, even for excitations with large speeds of propagation.

Qualitative dynamics of three different loop quantizations of spatially flat isotropic and homogeneous models is studied using effective spacetime description of the underlying quantum geometry. These include the standard loop quantum cosmology (LQC), its recently revived modification (referred to as mLQC-I), and another related modification of LQC (mLQC-II) whose dynamics is studied in detail for the first time. Various features of LQC, including quantum bounce and pre-inflationary dynamics, are found to be shared with the mLQC-I and mLQC-II models. We study universal properties of dynamics for chaotic inflation, fractional monodromy inflation, Starobinsky potential, non-minimal Higgs inflation, and an exponential potential. We find various critical points and study their stability, which reveal various qualitative similarities in the post-bounce phase for all these models. The pre-bounce qualitative dynamics of LQC and mLQC-II turns out to be very similar, but is strikingly different from that of mLQC-I. In the dynamical analysis, some of the fixed points turn out to be degenerate for which center manifold theory is used. For all these potentials, non-perturbative quantum gravitational effects always result in a non-singular inflationary scenario with a phase of super-inflation succeeded by the conventional inflation. We show the existence of inflationary attractors, and obtain scaling solutions in the case of the exponential potential. Since all of the models agree with general relativity at late times, our results are also of use in classical theory where qualitative dynamics of some of the potentials has not been studied earlier.

We systematically study the preinflationary dynamics of the spatially flat Friedmann-Lemaitre-Robertson-Walker universe filled with a single scalar field that has the generalized $\alpha-$attractor potentials, in the framework of loop quantum cosmology, in which the big bang singularity is replaced generically by a non-singular quantum bounce due to purely quantum geometric effects. The evolution can be divided into two different classes, one is dominated initially (at the quantum bounce) by the kinetic energy of the scalar field, and one is not. In both cases, we identify numerically the physically viable initial conditions that lead to not only a slow-roll inflationary phase, but also enough $e$-folds to be consistent with observations, and find that the output of such a viable slow-roll inflationary phase is generic. In addition, we also show that in the case when the evolution of the universe is dominated initially by the kinetic energy of the scalar field (except for a very small set in the phase space), the evolution before reheating is aways divided into three different phases: {\em bouncing, transition and slow-roll inflation}. This universal feature does not depend on the initial conditions of the system nor on the specific potentials of the scalar field, as long as it is dominated initially by the kinetic energy of the scalar field at the bounce. Moreover, we carry out phase space analyses for the models under consideration and compare our results with the power-law and Starobinsky potentials.

We demonstrate a compact 29.3 GHz visible astro-comb covering the spectrum from 560nm to 700nm. A 837 MHz Yb:fiber laser frequency comb phase locked to a Rb clock served as the seed comb to ensure the frequency stability and high side mode suppression ratio. After the visible super-continuum generation, a cavity-length-fixed Fabry-Perot cavity made by ultra-low expansion glass was utilized to filter the comb teeth for eliminating the rapid active dithering. The mirrors were home-made complementary chirped mirrors pair with zero dispersion and high reflection to guarantee no mode skipping. These filtered comb teeth were clearly resolved in an astronomical spectrograph of 49,000 resolution, exhibiting sharp linetype, zero noise floor, and uniform exposure amplitude.

Using the Wald formalism, we investigate the thermodynamics of charged black holes in D-dimensional stationary spacetimes with or without rotations in Einstein-\ae ther-Maxwell theory. In particular, assuming the existence of a scaling symmetry of the action, we obtain the Smarr integral formula, which can be applied to both Killing and universal horizons. When restricted to 4-dimensional spherically symmetric spacetimes, previous results obtained by a different method are re-derived.

We constrain cosmological parameters using only Hubble parameter data and quantify the impact of future Hubble parameter measurements on parameter estimation for the most typical dark energy models. We first constrain cosmological parameters using 52 current Hubble parameter data including the Hubble constant measurement from the Hubble Space Telescope. Then we simulate the baryon acoustic oscillation signals from WFIRST (Wide-Field Infrared Survey Telescope) covering the redshift range of $z\in [0.5, 2]$ and the redshift drift data from E-ELT (European Extremely Large Telescope) in the redshift range of $z\in [2,5]$. It is shown that solely using the current Hubble parameter data could give fairly good constraints on cosmological parameters. Compared to the current Hubble parameter data, with the WFIRST observation the $H(z)$ constraints on dark energy would be improved slightly, while with the E-ELT observation the $H(z)$ constraints on dark energy is enormously improved.

By a method of population synthesis we construct a model of dark matter (DM) accretion onto binary black holes (BHs), and investigate the merger rate of the binary BHs. We find that the merger rate can weakly increase (less than 10\%). However, the DM accretion can efficiently enhance the masses of binary BHs. In our model, the result for $Z = 0.01$ without the DM accretion can not explain the GW170104, GW170814, and GW150914, while with the DM accretion, it can cover all observations well. For the higher metallicity ($Z = 0.02$), our model can not explain the mergers of high-mass binary BHs like GW170104, GW170814, and GW150914. We estimate that the merger rate of binary BHs lies between $\rm55Gpc^{-3}yr^{-1}$ to $\rm197Gpc^{-3}yr^{-1}$.

In this paper, we carry out a systematic analysis of the theoretical and observational constraints on the dimensionless coupling constants $c_i$ ($i=1,2,3,4$) of the Einstein-aether theory, taking into account the events GW170817 and GRB 170817A. The combination of these events restricts the deviation of the speed $c_T$ of the spin-2 graviton to the range, $- 3\times 10^{-15} < c_T -1 < 7\times 10^{-16}$, which for the Einstein-aether theory implies $\left|c_{13}\right| \le 10^{-15}$ with $c_{ij} \equiv c_{i} + c_{j}$. The rest of the constraints are divided into two groups: those on the ($c_1, c_{14}$)-plane and those on the ($c_2, c_{14}$)-plane, except the strong-field constraints. The latter depend on the sensitivities $\sigma_{\ae}$ of neutron stars, which are not known at present in the new ranges of the parameters found in this paper.

We present a systematic study of the cosmological dynamics resulting from an effective Hamiltonian, recently derived in loop quantum gravity using Thiemann's regularization and earlier obtained in loop quantum cosmology (LQC) by keeping the Lorentzian term explicit in the Hamiltonian constraint. We show that quantum geometric effects result in higher than quadratic corrections in energy density in comparison to LQC causing a non-singular bounce. Dynamics can be described by the Hamilton's or the Friedmann-Raychaudhuri equations, but the map between the two descriptions is not one-to-one. A careful analysis resolves the tension on symmetric versus asymmetric bounce in this model, showing that the bounce must be asymmetric and symmetric bounce is physically inconsistent, in contrast to the standard LQC. In addition, the current observations only allow a scenario where the pre-bounce branch is asymptotically de Sitter, similar to a quantization of the Schwarzschild interior in LQC, and the post-bounce branch yields the classical general relativity. For a quadratic potential, we find that a slow-roll inflation generically happens after the bounce, which is quite similar to what happens in LQC.

We analytically investigate the Schwinger pair production in the de Sitter background by using {\em the uniform asymptotic approximation method}, and show that the equation of motion in general has two turning points, and the nature of these points could be single, double, real or complex, depending on the choice of the free parameters involved in the theory. Different natures of these points lead to different electric currents. In particular, when $\beta \equiv m^2/H^2-9/4$ is positive, both turning points are complex, and the electric current due to the Schwinger process is highly suppressed, where $m$ and $H$ denote, respectively, the mass of the particle and the Hubble parameter. For the turning points to be real, it is necessary to have $\beta < 0$, and the more negative of $\beta$, the easier to produce particles. In addition, when $\beta < 0$, we also study the particle production when the electric field $E$ is very weak. We find that the electric current in this case is proportional to $E^{1/2 - \sqrt{|\beta|}}$, which is strongly enhanced in the weak electric field limit when $m < \sqrt{2} H$.

In this paper, we study the pre-inflationary dynamics for the power-law potential $(V(\phi) \propto \phi^n)$ with $n<2$ in the framework of loop quantum cosmology. In the case where the kinetic energy of the inflaton dominates at the initial, the evolution of the universe can always be divided into three different phases prior to preheating: {\em bouncing, transition and slow-roll inflation}. During the bouncing phase, the evolution of the expansion factor is independent not only on the initial conditions but also the inflationary potentials, and is given explicitly by an analytical solution. In contrast, for the potential energy dominated initial conditions, this universality is lost. We also obtain total number of e-folds during the slow-roll inflation, whereby physically viable models are identified. In addition, we present phase space analysis for the inflationary potentials under consideration and compare our results with the ones obtained previously for different potentials.

We study the generalized $\alpha$ attractor model in the context of the late time cosmic acceleration. The model interpolates between the scaling freezing and thawing dark energy models. In the slow roll region, the original potential is modified whereas the modification ceases in the asymptotic region and the effective potential behaves as the quadratic one. In our setting, the field rolls slowly around the present epoch and mimics the scaling behavior in the future. We obtain observational constraints on the model parameters by using an integrated data base (SN+Hubble+BAO+CMB).

The evolutions of the flat FLRW universe and its linear perturbations are studied systematically in the dressed metric approach of LQC. When it is dominated by the kinetic energy of the inflaton at the quantum bounce, the evolution of the background can be divided into three different phases prior to the preheating, {\em bouncing, transition and slow-roll inflation}. During the bouncing phase, the evolution is independent of not only the initial conditions, but also the inflationary potentials. In particular, the expansion factor can be well described by the same exact solution in all the cases considered. In contrast, in the potential dominated case such a universality is lost. It is because of this universality that the linear perturbations are also independent of the inflationary models and obtained exactly. During the transition phase, the evolutions of the background and its linear perturbations are found explicitly, and then matched to the ones given in the other two phases. Hence, once the initial conditions are imposed, the linear scalar and tensor perturbations will be uniquely determined. Considering two different sets of initial conditions, one imposed during the contracting phase and the other at the bounce, we calculate the Bogoliubov coefficients and find that the two sets yield the same results and all lead to particle creations at the onset of the inflation. Due to the pre-inflationary dynamics, the scalar and tensor power spectra become scale-dependent. Comparing with the Planck 2015 data, we find constraints on the total e-folds that the universe must have expanded since the bounce, in order to be consistent with current observations.

In this paper, we study the non-projectable 2d Ho\v{r}ava gravity coupled with a non-relativistic scalar field, where the coupling is in general non-minimal and of the form $f(\phi)R$, where $f(\phi)$ is an arbitrary function of the scalar field $\phi$, and $R$ denotes the 2d Ricci scalar. In particular, we first investigate the Hamiltonian structure, and show that there are two-first and two-second class constraints, similar to the pure gravity case, but now the local degree of freedom is one, due to the presence of the scalar field. Then, we present various exact stationary solutions of this coupled system, and find that some of them represent black holes but now with universal horizons as their boundaries. At these horizons, the Hawking radiations are thermal with temperatures proportional to their surface gravities, which normally depend on the non-linear dispersion relations of the particles radiated, similar to the (3+1)-dimensional case.

Loop quantum cosmology (LQC) provides a resolution of the classical big bang singularity in the deep Planck era. The evolution, prior to the usual slow-roll inflation, naturally generates excited states at the onset of the slow-roll inflation. It is expected that these quantum gravitational effects could leave its fingerprints on the primordial perturbation spectrum and non-Gaussianity, and lead to some observational evidences in the cosmic microwave background (CMB). While the impact of the quantum effects on the primordial perturbation spectrum has been already studied and constrained by current data, in this paper we continue studying such effects on the non-Gaussianity of the primordial curvature perturbations. In this paper, we present detailed and analytical calculations of the non-Gaussianity and show explicitly that the corrections due to quantum effects are in the same magnitude of the slow-roll parameters in the observable scales and thus are well within current observational constraints. Despite this, we show that the non-Gaussianity in the squeezed limit can be enhanced at superhorizon scales and further, these effects may yield a large statistical anisotropy on the power spectrum through the Erickcek-Kamionkowski-Carroll mechanism.

In this Letter, we study analytically the evolutions of the flat Friedmann-Lemaitre-Robertson-Walker (FLRW) universe and its linear perturbations in the framework of {\em the dressed metric approach} in loop quantum cosmology (LQC). Assuming that the evolution of the background is dominated by the kinetic energy of the inflaton at the quantum bounce, we find that both evolutions of the background and its perturbations are independent of the inflationary potentials during the pre-inflationary phase. During this period the effective potentials of the perturbations can be well approximated by a P\"oschl-Teller (PT) potential, from which we find analytically the mode functions and then calculate the corresponding Bogoliubov coefficients at the onset of the slow-roll inflation, valid for any inflationary model with a single scalar field. Imposing the Bunch-Davies (BD) vacuum in the contracting phase prior to the bounce when the modes are all inside the Hubble horizon, we show that particles are generically created due to the pre-inflation dynamics. Matching them to those obtained in the slow-roll inflationary phase, we investigate the effects of the pre-inflation dynamics on the scalar and tensor power spectra and find features that can be tested by current and forthcoming observations. In particular, to be consistent with the Planck 2015 data, we find that the universe must have expanded at least $141$ e-folds since the bounce.

In the framework of a model based on the gravitational field of the Kerr black hole, we turn to investigate the kinematic behavior of extragalactic jets. We analytically calculate the observable velocities and accelerations along any geodesic. Then, by numerical calculations, we apply our results to a geodesic, typical of the M87 jet, and probe our results by confrontation to recent observations. A transition from non-relativistic to ultrarelativistic speeds at subparsec scale is highlighted. This transition comes sooner and more abruptly than in models based on magnetic paradigm, which indicates that we need a weaker magnetic field to explain observed synchrotron radiation. We attribute the ejection phenomenon to the repulsive effect of the gravitomagnetic Kerr field.

The slow-roll inflation for a single scalar field that couples to the Gauss-Bonnet (GB) term represents an important higher-order curvature correction inspired by string theory. With the arrival of the era of precision cosmology, it is expected that the high-order corrections become more and more important. In this paper we study the observational predictions of the slow-roll inflation with the GB term by using the third-order uniform asymptotic approximation method. We calculate explicitly the primordial power spectra, spectral indices, running of the spectral indices for both scalar and tensor perturbations, and the ratio between tensor and scalar spectra. These expressions are all written in terms of the Hubble and GB coupling flow parameters and expanded up to the next-to-leading order in the slow-roll expansions. The upper bounds of errors of the approximations at the third-order are $0.15\%$, so they represent the most accurate results obtained so far in the literature. We expect that the understanding of the GB corrections in the primordial spectra and their constraints by forthcoming observational data will provide clues for the UV complete theory of quantum gravity, such as the string/M-theory.

In this paper, we revisit the issue of static hairs of black holes in gravitational theories with broken Lorentz invariance in the case that the speed $c_{\phi}$ of the khronon field becomes infinitely large, $c_{\phi} = \infty$, for which the sound horizon of the khronon field coincides with the universal horizon, and the boundary conditions at the sound horizon reduce to those given normally at the universal horizons. As a result, less boundary conditions are present in this extreme case in comparison with the case $c_{\phi} = $ finite. Then, it would be expected that static hairs might exist. However, we show analytically that even in this case static hairs still cannot exist, based on a decoupling limit analysis. We also consider the cases in which $c_{\phi}$ is finite but with $c_{\phi} \gg 1$, and obtain the same conclusion.

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.

Ho\v{r}ava gravity at a Lifshitz point is a theory intended to quantize gravity by using techniques of traditional quantum field theories. To avoid Ostrogradsky's ghosts, a problem that has been plaguing quantization of general relativity since the middle of 1970's, Ho\v{r}ava chose to break the Lorentz invariance by a Lifshitz-type of anisotropic scaling between space and time at the ultra-high energy, while recovering (approximately) the invariance at low energies. With the stringent observational constraints and self-consistency, it turns out that this is not an easy task, and various modifications have been proposed, since the first incarnation of the theory in 2009. In this review, we shall provide a progress report on the recent developments of Ho\v{r}ava gravity. In particular, we first present four most-studied versions of Ho\v{r}ava gravity, by focusing first on their self-consistency and then their consistency with experiments, including the solar system tests and cosmological observations. Then, we provide a general review on the recent developments of the theory in three different but also related areas: (i) universal horizons, black holes and their thermodynamics; (ii) non-relativistic gauge/gravity duality; and (iii) quantization of the theory. The studies in these areas can be generalized to other gravitational theories with broken Lorentz invariance.

We study analytically quantum tunneling of relativistic and non-relativistic particles at both Killing and universal horizons of Einstein-Maxwell-aether black holes, after high-order curvature corrections are taken into account, for which the dispersion relation of the particles becomes nonlinear. Our results at the Killing horizons confirm the previous ones, i.e., at high frequencies the corresponding radiation remains thermal and the nonlinearity of the dispersion does not alter the Hawking radiation significantly. In contrary, non-relativistic particles are created at universal horizons and are radiated out to infinity. The radiation also has a thermal spectrum, and the corresponding temperature takes the form, $T^{z}_{UH} = 2\kappa_{UH} (z-1)/(2\pi z)$, where $z$ denotes the power of the leading term in the nonlinear dispersion relation, $\kappa_{UH}$ is the surface gravity of the universal horizon, defined by peering behavior of ray trajectories at the universal horizon. We also study the Smarr formula by assuming that: (a) the entropy is proportional to the area of the universal horizon, and (b) the first law of black hole thermodynamics holds, whereby we derive the Smarr mass, which in general is different from the total mass obtained at infinity. This indicates that one or both of these assumptions must be modified.

Presolar grains constitute remnants of stars that existed before the formation of the solar system. In addition to providing direct information on the materials from which the solar system formed, these grains provide ground-truth information for models of stellar evolution and nucleosynthesis. Here we report the in-situ identification of two unique presolar graphite grains from the primitive meteorite LaPaz Icefield 031117. Based on these two graphite grains, we estimate a bulk presolar graphite abundance of 5(-3)(+7) ppm in this meteorite. One of the grains (LAP-141) is characterized by an enrichment in 12C and depletions in 33,34S, and contains a small iron sulfide subgrain, representing the first unambiguous identification of presolar iron sulfide. The other grain (LAP-149) is extremely 13C-rich and 15N-poor, with one of the lowest 12C/13C ratios observed among presolar grains. Comparison of its isotopic compositions with new stellar nucleosynthesis and dust condensation models indicates an origin in the ejecta of a low-mass CO nova. Grain LAP-149 is the first putative nova grain that quantitatively best matches nova model predictions, providing the first strong evidence for graphite condensation in nova ejecta. Our discovery confirms that CO nova graphite and presolar iron sulfide contributed to the original building blocks of the solar system.

In this paper, we provide a systematic investigation of high-order primordial perturbations with nonlinear dispersion relations due to quantum gravitational effects in the framework of {\em uniform asymptotic approximations}. Because of these effects, the equation of motion of the mode function in general has multiple-turning points. After obtaining analytically approximated solutions to any order in different regions, associated with different types of turning points, we match them to the third one. To this order the errors are less than $0.15\%$. General expressions of the power spectra of the primordial tensor and scalar perturbations are derived explicitly. We also investigate effects of back-reactions of the quantum gravitational corrections, and make sure that inflation lasts long enough in order to solve the underlying problems, such as flatness, horizon and monopole. Then, we study various features of the spectra that are observationally relevant. In particular, under a moderate assumption about the energy scale of the underlying theory of quantum gravity, we have shown that the quantum gravitational effects may alter significantly the ratio between the tensor and scalar power spectra, thereby providing a natural mechanism to alleviate the tension between observations and certain inflationary models, including the one with a quadratic potential.

We present here the general expressions for the acceleration of massive test particles along the symmetry axis of the Kerr metric, and then study the main properties of this acceleration in different regions of the spacetime. In particular, we show that there exists a region near the black hole in which the gravitational field is repulsive. We provide possible physical interpretations about the role of this effect in terms of the different conserved parameters. The studies of these geodesics are important not only to understand better the structure of the Kerr spacetime but also to its use as a possible mechanism for the production of extragalactic jets. Our results are obtained with the help of expressing the geodesics of the Kerr spacetime in terms of the Weyl coordinates.

We first derive the primordial power spectra, spectral indices and runnings of both scalar and tensor perturbations of a flat inflationary universe to the second-order approximations of the slow-roll parameters, in the framework of loop quantum cosmology with the inverse-volume quantum corrections. This represents an extension of our previous work in which the parameter $\sigma$ was assumed to be an integer, where $\sigma$ characterizes the quantum corrections and in general can take any of values from the range $\sigma \in (0, 6]$. Restricting to the first-order approximations of the slow-roll parameters, we find corrections to the results obtained previously in the literature, and point out the causes for such errors. To our best knowledge, these represent the most accurate calculations of scalar and tensor perturbations given so far in the literature. Then, fitting the perturbations to the recently released data by Planck (2015), we obtain the most severe constraints for various values of $\sigma$. Using these constraints as our referring point, we discuss whether these quantum gravitational corrections can lead to measurable signatures in the future cosmological observations. We show that, depending on the value of $\sigma$, the scale-dependent contributions to the relativistic inflationary spectra due to the inverse-volume corrections could be well within the range of the detectability of the forthcoming generations of experiments, such as the Stage IV experiments.

Loop quantum cosmology (LQC) provides promising resolutions to the trans-Planckian issue and initial singularity arising in the inflationary models of general relativity. In general, due to different quantization approaches, LQC involves two types of quantum corrections, the holonomy and inverse-volume, to both of the cosmological background evolution and perturbations. In this paper, using {\em the third-order uniform asymptotic approximations}, we derive explicitly the observational quantities of the slow-roll inflation in the framework of LQC with these quantum corrections. We calculate the power spectra, spectral indices, and running of the spectral indices for both scalar and tensor perturbations, whereby the tensor-to-scalar ratio is obtained. We expand all the observables at the time when the inflationary mode crosses the Hubble horizon. As the upper error bounds for the uniform asymptotic approximation at the third-order are $\lesssim 0.15\%$, these results represent the most accurate results obtained so far in the literature. It is also shown that with the inverse-volume corrections, both scalar and tensor spectra exhibit a deviation from the usual shape at large scales. Then, using the Planck, BAO and SN data we obtain new constraints on quantum gravitational effects from LQC corrections, and find that such effects could be within the detection of the forthcoming experiments.

In the framework of the Einstein-Maxwell-aether theory, we present two new classes of exact charged black hole solutions, which are asymptotically flat and possess the universal as well as Killing horizons. We also construct the Smarr formulas, and calculate the temperatures of the horizons, using the Smarr mass-area relation. We find that, in contrast to the neutral case, such obtained temperature is not proportional to its surface gravity at any of the two kinds of the horizons. Einstein-Maxwell-aether black holes with the cosmological constant and their topological cousins are also given.

We derive the primordial power spectra and spectral indexes of the density fluctuations and gravitational waves in the framework of loop quantum cosmology (LQC) with holonomy and inverse-volume corrections, by using the uniform asymptotic approximation method to its third-order, at which the upper error bounds are $\lesssim 0.15\%$, and accurate enough for the current and forthcoming cosmological observations. Then, using the Planck, BAO and SN data we obtain the tightest constraints on quantum gravitational effects from LQC corrections, and find that such effects could be well within the detection of the current and forthcoming cosmological observations.

In this paper, we first generalize the definition of stationary universal horizons to dynamical ones, and then show that (dynamical) universal horizons can be formed from realistic gravitational collapse. This is done by constructing analytical models of a collapsing spherically symmetric star with finite thickness in Einstein-aether theory.

$k$-inflation represents the most general single-field inflation, in which the perturbations usually obey an equation of motion with a time-dependent sound speed. In this paper, we study the observational predictions of the $k$-inflation by using the high-order uniform asymptotic approximation method. We calculate explicitly the slow-roll expressions of the power spectra, spectral indices, and running of the spectral indices for both the scalar and tensor perturbations. These expressions are all written in terms of the Hubble and sound speed flow parameters. It is shown that the previous results obtained by using the first-order uniform asymptotic approximation have been significantly improved by the high-order corrections of the uniform asymptotic approximations. Furthermore, we also check our results by comparing them with the ones obtained by other approximation methods, including the Green's function method, WKB approximation, and improved WKB approximation, and find the relative errors.

Effective data visualization is a key part of the discovery process in the era of big data. It is the bridge between the quantitative content of the data and human intuition, and thus an essential component of the scientific path from data into knowledge and understanding. Visualization is also essential in the data mining process, directing the choice of the applicable algorithms, and in helping to identify and remove bad data from the analysis. However, a high complexity or a high dimensionality of modern data sets represents a critical obstacle. How do we visualize interesting structures and patterns that may exist in hyper-dimensional data spaces? A better understanding of how we can perceive and interact with multi dimensional information poses some deep questions in the field of cognition technology and human computer interaction. To this effect, we are exploring the use of immersive virtual reality platforms for scientific data visualization, both as software and inexpensive commodity hardware. These potentially powerful and innovative tools for multi dimensional data visualization can also provide an easy and natural path to a collaborative data visualization and exploration, where scientists can interact with their data and their colleagues in the same visual space. Immersion provides benefits beyond the traditional desktop visualization tools: it leads to a demonstrably better perception of a datascape geometry, more intuitive data understanding, and a better retention of the perceived relationships in the data.

The uniform asymptotic approximation method provides a powerful, systematically-improved, and error-controlled approach to construct accurate analytical approximate solutions of mode functions of perturbations of the Friedmann-Robertson-Walker universe, designed especially for the cases where the relativistic linear dispersion relation is modified after gravitational quantum effects are taken into account. These include models from string/M-Theory, loop quantum cosmology and Ho\v{r}ava-Lifshitz quantum gravity. In this paper, we extend our previous studies of the first-order approximations to high orders for the cases where the modified dispersion relation (linear or nonlinear) has only one-turning point (or zero). We obtain the general expressions for the power spectra and spectral indices of both scalar and tensor perturbations up to the third-order, at which the error bounds are $\lesssim 0.15\%$. As an application of these formulas, we calculate the power spectra and spectral indices in the slow-roll inflation with a nonlinear power-law dispersion relation. To check the consistency of our formulas, we further restrict ourselves to the relativistic case, and calculate the corresponding power spectra, spectral indices and runnings to the second-order. Then, we compare our results with the ones obtained by the Green function method, and show explicitly that the results obtained by these two different methods are consistent within the allowed errors.

We develop a technique to construct analytical solutions of the linear perturbations of inflation with a nonlinear dispersion relation, due to quantum effects of the early universe. Error bounds are given and studied in detail. The analytical solutions describe the exact evolution of the perturbations extremely well even when only the first-order approximations is used.

Recently BICEP2 found that the vanishing of the tensor-to-scalar ratio $r$ is excluded at $7\sigma$ level, and its most likely value is $r=0.2^{+0.07}_{-0.05}$ at $1 \sigma$ level. This immediately causes a tension with the Planck constraint $r < 0.11$. In addition, it also implies that the inflaton (in single field slow-roll inflation) experienced a Planck excursion during inflation $\Delta\phi /M_{Pl} \geq {\cal{O}}(1)$, whereby the effective theory of inflation becomes questionable. In this brief report, we show that the inflationary paradigm is still robust, even after the quantum effects are taken into account. Moreover, these effects also help to relax the tension on the different values of $r$ given by BICEP2 and Planck.

In this paper, we study static vacuum solutions of quantum gravity at a fixed Lifshitz point in (2+1) dimensions, and present all the diagonal solutions in closed forms in the infrared limit. The exact solutions represent spacetimes with very rich structures: they can represent generalized BTZ black holes, Lifshitz space-times or Lifshitz solitons, in which the spacetimes are free of any kind of space-time singularities, depending on the choices of the free parameters of the solutions. We also find several classes of exact static non-diagonal solutions, which represent similar space-time structures as those given in the diagonal case. The relevance of these solutions to the non-relativistic Lifshitz-type gauge/gravity duality is discussed.

In this paper, we first propose a universal coupling between the gravity and matter in the framework of the Ho\v{r}ava-Lifshitz theory of gravity with an extra U(1) symmetry for both the projectable and non-projectable cases. Then, using this universal coupling we study the post-Newtonian approximations and obtain the parameterized post-Newtonian (PPN) parameters in terms of the coupling constants of the theory. Contrary to the previous works in which only two PPN parameters were calculated, we obtain {\it all} PPN parameters. We then, for the first time in either projectable or non-projectable case, find that all the solar system tests carried out so far are satisfied in a large region of the parameters space. In particular, the same results obtained in general relativity can be easily realized here. A remarkable feature is that the solar system tests impose no constraint on the parameter $\lambda$ appearing in the kinetic part of the action. As a result, the solar system tests, when combined with the condition for avoidance of strong coupling, do not lead to an upper bound on the energy scale $M_{*}$ that suppresses higher dimensional operators in the theory. This is in sharp contrast to other versions of the HL theory.

In this paper, we study 3-point correlation function of primordial gravitational waves generated in the de Sitter background in the framework of the general covariant Ho\v{r}ava-Lifshitz gravity with an arbitrary coupling constant $\lambda$. We find that, at cubic order, the interaction Hamiltonian receives contributions from four terms built of the 3-dimensional Ricci tensor $R_{ij}$ of the leaves $t = $ constant. In particular, the 3D Ricci scalar $R$ yields the same $k$-dependence as that in general relativity, but with different magnitude due to coupling with the $U(1)$ field $A$ and a UV history. Interestingly, the two terms $R_{ij}R^{ij}$ and $\left(\nabla^{i}R^{jk}\right)\left(\nabla_{i}R_{jk}\right)$ exhibit peaks at the squeezed limit. We show that this is due to the effects of the polarization tensors. The signal generated by the fourth term, $R^i_j R^j_k R^k_i$, favors the equilateral shape when spins of the three tensor fields are the same, but peaks in between the equilateral and squeezed limits when spins are mixed. The consistency with the recently-released Planck observations on non-Gaussianity is also discussed and is found that $\left(H/M_*\right)^2\left(H/M_{pl}\right) \le 10^{-8}$, where $M_{*}$ denotes the suppression energy of high-order operators, $M_{pl}$ the Planck mass, and $H$ the energy of inflation.

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.

We present a technique, {\em the uniform asymptotic approximation}, to construct accurate analytical solutions of the linear perturbations of inflation after quantum effects of the early universe are taken into account, for which the dispersion relations generically become nonlinear. We construct explicitly the error bounds associated with the approximations and then study them in detail. With the understanding of the errors and the proper choice of the Liouville transformations of the differential equations of the perturbations, we show that the analytical solutions describe the exact evolution of the linear perturbations extremely well even only in the first-order approximations. As an application of the approximate analytical solutions, we calculate the power spectra and indices of scalar and tensor perturbations in the slow-roll inflation, and find that the amplitudes of the power spectra get modified due to the quantum effects, while the power spectrum indices remain the same as in the linear case.

In this paper, we study inflation in the general covariant Ho\v{r}ava-Lifshitz gravity without the projectability condition. We write down explicitly the equations of the linear scalar perturbations of the FRW universe for a single scalar field without specifying to any gauge. Applying these equations to a particular gauge, we are able to obtain a master equation of the perturbations, in contrast to all the other versions of the theory without the projectability condition. This is because in the current version of the theory it has the same degree of freedom of general relativity. To calculate the power spectrum and index, we first define the initial conditions as the ones that minimize the energy of the ground state. Then, we obtain the full solutions of the equation of motion, by using the WKB approximations. From these solutions, we calculate the power spectrum and spectrum index of the comoving curvature perturbations and find the corrections due to the high order spatial derivative terms of the theory to those standard results obtained in general relativity. Remarkably, partly of corrections are a direct consequence of the non-projectability condition. It is also shown that the perturbations are still of scale-invariance, and the results obtained in the general covariant Ho\v{r}ava-Lifshitz gravity without the projectability condition are consistent with all current cosmological observations.

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.

In this paper, we study non-Gaussianity generated by a single scalar field in slow-roll inflation in the framework of the non-relativistic general covariant Ho\v{r}ava-Lifshitz theory of gravity with the projectability condition and an arbitrary coupling constant $\lambda$, where $\lambda$ characterizes the deviation of the theory from general relativity (GR) in the infrared. We find that the leading effect of self-interaction, in contrary to the case of minimal scenario of GR, is in general of the order $\hat{\alpha}_{n} \epsilon^{3/2}$, where $\epsilon$ is a slow-roll parameter, and $\hat{\alpha}_{n} (n = 3, 5)$ are the dimensionless coupling coefficients of the six-order operators of the Lifshitz scalar, and have no contributions to power spectra and indices of both scalar and tensor. The bispectrum, comparing with the standard one given in GR, is enhanced, and gives rise to a large value of the nonlinearity parameter $f_{\text{NL}}$.We study how the modified dispersion relation with high order moment terms affects the evaluation of the mode function and in turn the bispectrum, and show explicitly that the mode function takes various asymptotic forms during different periods of its evolution. In particular, we find that it is in general of superpositions of oscillatory functions, instead of plane waves like in the minimal scenario of GR. This results in a large enhancement of the folded shape in the bispectrum.

We study spherically symmetric, stationary vacuum configurations in general covariant theory (U(1) extension) of Ho\v{r}ava-Lifshitz gravity with the projectability condition and an arbitrary coupling constant $\lambda$, and obtain all the solutions in closed forms. If the gauge field $A$ and the Newtonian prepotential $\varphi$ do not directly couple to matter fields, the theory is inconsistent with solar system tests for $\lambda\not=1$, no matter how small $|\lambda-1|$ is. This is shown to be true also with the most general ansatz of spherical (but not necessarily stationary) configurations. Therefore, to be consistent with observations, one needs either to find a mechanism to restrict $\lambda$ precisely to $\lambda_{GR}=1$, or to consider $A$ and/or $\varphi$ as parts of the 4-dimensional metric on which matter fields propagate. In the latter, requiring that the line element be invariant not only under the foliation-preserving diffeomorphism but also under the local U(1) transformations, we propose the replacements, $N \rightarrow N - \upsilon(A - {\cal{A}})/c^2$ and $N^i \rightarrow N^i+N\nabla^{i}\varphi$, where $\upsilon$ is a dimensionless coupling constant to be constrained by observations, $N$ and $N^i$ are, respectively, the lapse function and shift vector, and ${\cal{A}} \equiv - \dot{\varphi} + N^i\nabla_{i}\varphi + N(\nabla_{i}\varphi)^2/2$. With this prescription, we show explicitly that the aforementioned solutions are consistent with solar system tests for both $\lambda=1$ and $\lambda\not=1$, provided that $|\upsilon-1|<10^{-5}$. From this result, the physical and geometrical interpretations of the fields $A$ and $\varphi$ become clear. However, it still remains to be understood how to obtain such a prescription from the action principle.

In this paper, we show that the spin-0 gravitons appearing in Horava-Lifshitz gravity without the projectability condition can be eliminated by extending the gauge symmetries of the foliation-preserving diffeomorphisms to include a local U(1) symmetry. As a result, the problems of stability, ghost, strong coupling, and different speeds in the gravitational sector are automatically resolved. In addition, with the detailed balance condition softly breaking, the number of independent coupling constants can be significantly reduced (from more than 70 down to 15), while the theory is still UV complete and possesses a healthy IR limit, whereby the prediction powers of the theory are considerably improved. The strong coupling problem in the matter sector can be cured by introducing an energy scale $M_{*}$, so that $M_{*} < \Lambda_{\omega}$, where $M_{*}$ denotes the suppression energy of high order derivative terms, and $\Lambda_{\omega}$ the would-be strong coupling energy scale.

In this paper, we study inflation in the framework of the nonrelativistic general covariant theory of the Ho\v{r}ava-Lifshitz gravity with the projectability condition and an arbitrary coupling constant $\lambda$. We find that the Friedmann-Robterson-Walker (FRW) universe is necessarily flat in such a setup. We work out explicitly the linear perturbations of the flat FRW universe without specifying to a particular gauge, and find that the perturbations are different from those obtained in general relativity, because of the presence of the high-order spatial derivative terms. Applied the general formulas to a single scalar field, we show that in the sub-horizon regions, the metric and scalar field are tightly coupled and have the same oscillating frequencies. In the super-horizon regions, the perturbations become adiabatic, and the comoving curvature perturbation is constant. We also calculate the power spectra and indices of both the scalar and tensor perturbations, and express them explicitly in terms of the slow roll parameters and the coupling constants of the high-order spatial derivative terms. In particular, we find that the perturbations, of both scalar and tensor, are almost scale-invariant, and the spectrum indices are the same as those given in GR, but the ratio of the scalar and tensor power spectra depends on the high-order spatial derivative terms, and can be different from that of GR significantly.

We systematically study the evolution of the Friedmann-Robertson-Walker (FRW) universe coupled with a cosmological constant $\Lambda$ and a perfect fluid that has the equation of state $p=w\rho$, where $p$ and $\rho$ denote, respectively, the pressure and energy density of the fluid, and $w$ is an arbitrary real constant. Depending on the specific values of $w,\; \Lambda$, and the curvature $k$ of 3-dimensional space, we separate all of the solutions into various cases. In each case the main properties of the evolution are given in detail, including the periods of deceleration and/or acceleration, and the existence of big bang, big crunch, and big rip singularities. In some cases, errors in classification and interpretation appearing in standard textbooks have been corrected.

We present a fully nonlinear study of long wavelength cosmological perturbations within the framework of the projectable Horava-Lifshitz gravity, coupled to a single scalar field. Adopting the gradient expansion technique, we explicitly integrate the dynamical equations up to any order of the expansion, then restrict the integration constants by imposing the momentum constraint. While the gradient expansion relies on the long wavelength approximation, amplitudes of perturbations do not have to be small. When the $\lambda\to 1$ limit is taken, the obtained nonlinear solutions exhibit a continuous behavior at any order of the gradient expansion, recovering general relativity in the presence of a scalar field and the "dark matter as an integration constant". This is in sharp contrast to the results in the literature based on the "standard" (and naive) perturbative approach where in the same limit, the perturbative expansion of the action breaks down and the scalar graviton mode appears to be strongly coupled. We carry out a detailed analysis on the source of these apparent pathologies and determine that they originate from an improper application of the perturbative approximation in the momentum constraint. We also show that there is a new branch of solutions, valid in the regime where $|\lambda-1|$ is smaller than the order of perturbations. In the limit $\lambda\to1$, this new branch allows the theory to be continuously connected to general relativity (plus "dark matter").

We consider an extended theory of Horava-Lifshitz gravity with the detailed balance condition softly breaking, but without the projectability condition. With the former, the number of independent coupling constants is significantly reduced. With the latter and by extending the original foliation-preserving diffeomorphism symmetry $ {{Diff}}(M, {\cal{F}})$ to include a local U(1) symmetry, the spin-0 gravitons are eliminated. Thus, all the problems related to them disappear, including the instability, strong coupling, and different speeds in the gravitational sector. When the theory couples to a scalar field, we find that the scalar field is not only stable in both the ultraviolet (UV) and infrared (IR), but also free of the strong coupling problem, because of the presence of high-order spatial derivative terms of the scalar field. Furthermore, applying the theory to cosmology, we find that due to the additional U(1) symmetry, the Friedmann-Robertson-Walker (FRW) universe is necessarily flat. We also investigate the scalar, vector, and tensor perturbations of the flat FRW universe, and derive the general linearized field equations for each kind of the perturbations.

We systematically study black holes in the Horava-Lifshitz (HL) theory by following the kinematic approach, in which a horizon is defined as the surface at which massless test particles are infinitely redshifted. Because of the nonrelativistic dispersion relations, the speed of light is unlimited, and test particles do not follow geodesics. As a result, there are significant differences in causal structures and black holes between general relativity (GR) and the HL theory. In particular, the horizon radii generically depend on the energies of test particles. Applying them to the spherical static vacuum solutions found recently in the nonrelativistic general covariant theory of gravity, we find that, for test particles with sufficiently high energy, the radius of the horizon can be made as small as desired, although the singularities can be seen in principle only by observers with infinitely high energy. In these studies, we pay particular attention to the global structure of the solutions, and find that, because of the foliation-preserving-diffeomorphism symmetry, ${Diff}(M,{\cal{F}})$, they are quite different from the corresponding ones given in GR, even though the solutions are the same. In particular, the ${Diff}(M,{\cal{F}})$ does not allow Penrose diagrams. Among the vacuum solutions, some give rise to the structure of the Einstein-Rosen bridge, in which two asymptotically flat regions are connected by a throat with a finite non-zero radius. We also study slowly rotating solutions in such a setup, and obtain all the solutions characterized by an arbitrary function $A_{0}(r)$. The case $A_{0} = 0$ reduces to the slowly rotating Kerr solution obtained in GR.

We study the strong coupling problem in the Horava-Melby-Thompson setup of the Horava-Lifshitz gravity with an arbitrary coupling constant $\lambda$, generalized recently by da Silva, where $\lambda$ describes the deviation of the theory in the infrared from general relativity that has $\lambda_{GR} = 1$. We find that a scalar field in the Minkowski background becomes strong coupling for processes with energy higher than $\Lambda_{\omega} [\equiv (M_{pl}/c_1)^{3/2} M_{pl}|\lambda - 1|^{5/4}]$, where generically $c_1 \ll M_{pl}$. However, this problem can be cured by introducing a new energy scale $M_{*}$, so that $M_{*} < \Lambda_{\omega}$, where $M_{*}$ denotes the suppression energy of high order derivative terms of the theory.

Detailed balance and projectability conditions are two main assumptions when Horava recently formulated his theory of quantum gravity - the Horava-Lifshitz (HL) theory. While the latter represents an important ingredient, the former often believed needs to be abandoned, in order to obtain an ultraviolet stable scalar field, among other things. In this paper, because of several attractive features of this condition, we revisit it, and show that the scalar field can be stabilized, if the detailed balance condition is allowed to be softly broken. Although this is done explicitly in the non-relativistic general covariant setup of Horava-Melby-Thompson with an arbitrary coupling constant $\lambda$, generalized lately by da Silva, it is also true in other versions of the HL theory. With the detailed balance condition softly breaking, the number of independent coupling constants can be still significantly reduced. It is remarkable to note that, unlike other setups, in this da Silva generalization, there exists a master equation for the linear perturbations of the scalar field in the flat Friedmann-Robertson-Walker background.

In this paper, we investigate three important issues: stability, ghost and strong coupling, in the Horava-Melby-Thompson setup of the Horava-Lifshitz theory with $\lambda \not= 1$, generalized recently by da Silva. We first develop the general linear scalar perturbations of the Friedmann-Robertson-Walker (FRW) universe with arbitrary spatial curvature, and find that an immediate by-product of the setup is that, in all the inflationary models described by a scalar field, the FRW universe is necessarily flat. Applying them to the case of the Minkowski background, we find that it is stable, and, similar to the case $\lambda = 1$, the spin-0 graviton is eliminated. The vector perturbations vanish identically in the Minkowski background. Thus, similar to general relativity, a free gravitational field in this setup is completely described by a spin-2 massless graviton even with $\lambda \not= 1$. We also study the ghost problem in the FRW background, and find explicitly the ghost-free conditions. To study the strong coupling problem, we consider two different kinds of spacetimes all with the presence of matter, one is cosmological and the one is static. We find that the coupling becomes strong for a process with energy higher than $M_{pl} |c_{\psi}|^{5/2}$ in the flat FRW background, and $M_{pl}|c_{\psi}|^{3}$ in a static weak gravitational field, where $|c_{\psi}| \equiv |(1-\lambda)/(3 \lambda -1)|^{1/2}$.

Inclusion of $f(R)$ term in the action of Horava-Lifshitz quantum gravity with projectability but without detailed balance condition is investigated, where $R$ denotes the 3-spatial dimensional Ricci scalar. Conditions for the spin-0 graviton to be free of ghosts and instability are studied. The requirement that the theory reduce to general relativity in the IR makes the scalar mode unstable in the Minkowski background but stable in the de Sitter. It is remarkable that the dark sector, dark matter and dark energy, of the universe has a naturally geometric origin in such a setup. Bouncing universes can also be constructed. Scalar perturbations in the FRW backgrounds with non-zero curvature are presented.

The stability of de Sitter spacetime in Horava-Lifshitz theory of gravity with projectability but without detailed balance condition is studied. It is found that, in contrast to the case of the Minkowski background, the spin-0 graviton now is stable for any given $\xi$, and free of ghost for $\xi \le 0$ in the infrared limit, where $\xi$ is the dynamical coupling constant.

Horava and Melby-Thompson recently proposed a new version of the Horava-Lifshitz theory of gravity, in which the spin-0 graviton is eliminated by introducing a Newtonian pre-potential $\phi$ and a local U(1) gauge field $A$. In this paper, we first derive the corresponding Hamiltonian, super-momentum constraints, the dynamical equations, and the equations for $\phi$ and $A$, in the presence of matter fields. Then, we apply the theory to cosmology, and obtain the modified Friedmann equation and the conservation law of energy, in addition to the equations for $\phi$ and $A$. When the spatial curvature is different from zero, terms behaving like dark radiation and stiff-fluid exist, from which, among other possibilities, bouncing universe can be constructed. We also study linear perturbations of the FRW universe with any given spatial curvature $k$, and derive the most general formulas for scalar perturbations. The vector and tensor perturbations are the same as those recently given by one of the present authors [A. Wang, Phys. Rev. D{\bf 82}, 124063 (2010)] in the setup of Sotiriou, Visser and Weinfurtner. Applying these formulas to the Minkowski background, we have shown explicitly that the scalar and vector perturbations of the metric indeed vanish, and the only remaining modes are the massless spin-2 gravitons.

In this paper, we consider two different issues, stability and strong coupling, raised lately in the newly-proposed Horava-Lifshitz (HL) theory of quantum gravity with projectability condition. We find that all the scalar modes are stable in the de Sitter background, due to two different kinds of effects, one from high-order derivatives of the spacetime curvature, and the other from the exponential expansion of the de Sitter space. Combining these effects properly, one can make the instability found in the Minkowski background never appear even for small-scale modes, provided that the IR limit is sufficiently closed to the relativistic fixed point. At the fixed point, all the modes become stabilized. We also show that the instability of Minkowski spacetime can be cured by introducing mass to the spin-0 graviton. The strong coupling problem is investigated following the effective field theory approach, and found that it cannot be cured by the Blas-Pujolas-Sibiryakov mechanism, initially designed for the case without projectability condition, but might be circumvented by the Vainshtein mechanism, due to the non-linear effects. In fact, we construct a class of exact solutions, and show explicitly that it reduces smoothly to the de Sitter spacetime in the relativistic limit.

Spacetime singularities are studied in both the $D+d$-dimensional string theory and its $D$-dimensional effective theory, obtained by the Kaluza-Klein compactification. It is found that spacetime singularities in the low dimensional effective theory may or may not remain after lifted to the $D+d$-dimensional string theory, depending on particular solutions. It is also found that there exist cases in which spacetime singularities appearing in high/low dimensional spacetimes do not necessarily happen on the same surfaces.

We study cosmological vector and tensor perturbations in Horava-Lifshitz gravity, adopting the most general Sotiriou-Visser-Weinfurtner generalization without the detailed balance but with projectability condition. After deriving the general formulas in a flat FRW background, we find that the vector perturbations are identical to those given in general relativity. This is true also in the non-flat cases. For the tensor perturbations, high order derivatives of the curvatures produce effectively an anisotropic stress, which could have significant efforts on the high-frequency modes of gravitational waves, while for the low-frenquency modes, the efforts are negligible. The power spectrum is scale-invariant in the UV regime, because of the particular dispersion relations. But, due to lower-order corrections, it will eventually reduce to that given in GR in the IR limit. Applying the general formulas to the de Sitter and power-law backgrounds, we calculate the power spectrum and index, using the uniform approximations, and obtain their analytical expressions in both cases.

We study spherically symmetric static spacetimes generally filled with an anisotropic fluid in the nonrelativistic general covariant theory of gravity. In particular, we find that the vacuum solutions are not unique, and can be expressed in terms of the $U(1)$ gauge field $A$. When solar system tests are considered, severe constraints on $A$ are obtained, which seemingly pick up the Schwarzschild solution uniquely. In contrast to other versions of the Horava-Lifshitz theory, non-singular static stars made of a perfect fluid without heat flow can be constructed, due to the coupling of the fluid with the gauge field. These include the solutions with a constant pressure. We also study the general junction conditions across the surface of a star. In general, the conditions allow the existence of a thin matter shell on the surface. When applying these conditions to the perfect fluid solutions with the vacuum ones as describing their external spacetimes, we find explicitly the matching conditions in terms of the parameters appearing in the solutions. Such matching is possible even without the presence of a thin matter shell.

Five-dimensional spacetimes of two orbifold 3-branes are studied, by assuming that {\em the two 3-branes are spatially homogeneous, isotropic, and independent of time}, following the so-called "bulk-based" approach. The most general form of the metric is obtained, and the corresponding field equations are divided into three groups, one is valid on each of the two 3-branes, and the third is valid in the bulk. The Einstein tensor on the 3-branes is expressed in terms of the discontinuities of the first-order derivatives of the metric coefficients. Thus, once the metric is known in the bulk, the distribution of the Einstein tensor on the two 3-branes is uniquely determined. As applications, we consider two different cases, one is in which the bulk is locally $AdS_{5}$, and the other is where it is vacuum. In some cases, it is shown that the universe is first decelerating and then accelerating. The global structure of the bulk as well as the 3-branes is also studied, and found that in some cases the solutions may represent the collision of two orbifold 3-branes. The applications of the formulas to the studies of the cyclic universe and the cosmological constant problem are also pointed out.

We study the radion stability and radion mass in the framework of the Horava-Witten (HW) heterotic M-Theory on $S^{1}/Z_{2}$, and find that the radion is stable and its mass can be of the order of GeV. The gravity is localized on the visible brane, and the spectrum of the gravitational Kaluza-Klein (KK) modes is discrete and can have a mass gap of TeV. The corrections to the 4D Newtonian potential from the higher order gravitational KK modes are exponentially suppressed. Applying such a setup to cosmology, we find the generalized Friedmann-like equations on each of the two orbifold branes.

Dynamical models of prototype gravastars were constructed in order to study their stability. The models are the Visser-Wiltshire three-layer gravastars, in which an infinitely thin spherical shell of stiff fluid divides the whole spacetime into two regions, where the internal region is de Sitter, and the external is Schwarzschild. It is found that in some cases the models represent the "bounded excursion" stable gravastars, in which the thin shell is oscillating between two finite radii, while in other cases they collapse until the formation of black holes. In the phase space, the region for the ``bounded excursion" gravastars is very small in comparison to that of black holes, but not empty. Therefore, although the existence of gravastars cannot be excluded from such dynamical models, our results do indicate that, even if gravastars indeed exist, they do not exclude the existence of black holes.

We present a systematical study of brane worlds in string theory on $S^{1}/Z_{2}$. Starting with the toroidal compactification of the NS/NS sector in (D+d) dimensions, we first obtain an effective $D$-dimensional action, and then compactify one of the $(D-1)$ spatial dimensions by introducing two orbifold branes as its boundaries. By combining the Gauss-Codacci and Lanczos equations, we write down explicitly the general gravitational field equations on each of the two branes, while using distribution theory we express the matter field equations on the branes in terms of the discontinuities of the first derivatives of the matter fields. Afterwards, we address three important issues: (i) the hierarchy problem; (ii) the radion mass; and (iii) the localization of gravity, the 4-dimensional Newtonian effective potential and the Yukawa corrections due to the gravitational high-order Kaluza-Klein (KK) modes. With a very conservative estimation, we find that the radion mass is of the order of $10^{-2} GeV$. The gravity is localized on the visible brane, and the spectrum of the gravitational KK modes is discrete and can be of the order of TeV. The corrections to the 4-dimensional Newtonian potential from the higher order of gravitational KK modes are exponentially suppressed and can be safely neglected in current experiments. In an appendix, we also present a systematical and pedagogical study of the Gauss-Codacci equations and Israel's junction conditions across a (D-1)-dimensional hypersurface, which can be either spacelike or timelike.

We study perturbations of a scalar field cosmology in Horava-Lifshitz gravity, adopting the most general setup without detailed balance but with the projectability condition. We derive the generalized Klein-Gordon equation, which is sixth-order in spatial derivatives. Then we investigate scalar field perturbations coupled to gravity in a flat Friedmann-Robertson-Walker background. In the sub-horizon regime, the metric and scalar field modes have independent oscillations with different frequencies and phases except in particular cases. On super-horizon scales, the perturbations become adiabatic during slow-roll inflation driven by a single field, and the comoving curvature perturbation is constant.

We use current and future simulated data of the growth rate of large scale structure in combination with data from supernova, BAO, and CMB surface measurements, in order to put constraints on the growth index parameters. We use a recently proposed parameterization of the growth index that interpolates between a constant value at high redshifts and a form that accounts for redshift dependencies at small redshifts. We also suggest here another exponential parameterization with a similar behaviour. The redshift dependent parametrizations provide a sub-percent precision level to the numerical growth function, for the full redshift range. Using these redshift parameterizations or a constant growth index, we find that current available data from galaxy redshift distortions and Lyman-alpha forests is unable to put significant constraints on any of the growth parameters. For example both $\Lambda$CDM and flat DGP are allowed by current growth data. We use an MCMC analysis to study constraints from future growth data, and simulate pessimistic and moderate scenarios for the uncertainties. In both scenarios, the redshift parameterizations discussed are able to provide significant constraints and rule out models when incorrectly assumed in the analysis. The values taken by the constant part of the parameterizations as well as the redshift slopes are all found to significantly rule out an incorrect background. We also find that, for our pessimistic scenario, an assumed constant growth index over the full redshift range is unable to rule out incorrect models in all cases. This is due to the fact that the slope acts as a second discriminator at smaller redshifts and therefore provide a significant test to identify the underlying gravity theory.

We systematically study spherically symmetric static spacetimes filled with a fluid in the Horava-Lifshitz theory of gravity with the projectability condition, but without the detailed balance. We establish that when the spacetime is spatially Ricci flat the unique vacuum solution is the de Sitter Schwarzshcild solution, while when the spacetime has a nonzero constant curvature, there exist two different vacuum solutions; one is an (Einstein) static universe, and the other is a new spacetime. This latter spacetime is maximally symmetric and not flat. We find all the perfect fluid solutions for such spacetimes, in addition to a class of anisotropic fluid solutions of the spatially Ricci flat spacetimes. To construct spacetimes that represent stars, we investigate junction conditions across the surfaces of stars and obtain the general matching conditions with or without the presence of infinitely thin shells. It is remarkable that, in contrast to general relativity, the radial pressure of a star does not necessarily vanish on its surface even without the presence of a thin shell, due to the presence of high order derivative terms. Applying the junction conditions to our explicit solutions, we show that it is possible to match smoothly these solutions (all with nonzero radial pressures) to vacuum spacetimes without the presence of thin matter shells on the surfaces of stars.

In this brief report, we summarize our recent studies in brane cosmology in both string theory and M-Theory on $S^{1}/Z_{2}$. In such setups, we find that the radion is stable and its mass, with a very conservative estimation, can be of the order of $0. 1 \sim 0.01$ GeV. The hierarchy problem can be addressed by combining the large extra dimension, warped factor, and tension coupling mechanisms. Gravity is localized on the visible brane, and the spectrum of the gravitational Kaluza-Klein (KK) modes is discrete and can have a mass gap of TeV. The corrections to the 4D Newtonian potential from the higher order gravitational KK modes are exponentially suppressed. Applying such setups to cosmology, we find that a late transient acceleration of the universe seems to be the generic feature of the theory, due to the interaction between branes and bulk. A bouncing early universe is also rather easily realized.

Singularities in $(3+1)$-dimensional Horava-Lifshitz (HL) theory of gravity are studied. These singularities can be divided into scalar, non-scalar curvature, and coordinate singularities. Because of the foliation-preserving diffeomorphisms of the theory, the number of scalars that can be constructed from the extrinsic curvature tensor $K_{ij}$, the 3-dimensional Riemann tensor and their derivatives is much large than that constructed from the 4-dimesnional Riemann tensor and its derivatives in general relativity (GR). As a result, even for the same spacetime, it may be singular in the HL theory but not in GR. Two representative families of solutions with projectability condition are studied, one is the (anti-) de Sitter Schwarzschild solutions, and the other is the Lu-Mei-Pope (LMP) solutions written in a form satisfying the projectability condition - the generalized LMP solutions. The (anti-) de Sitter Schwarzschild solutions are vacuum solutions of both HL theory and GR, while the LMP solutions with projectability condition satisfy the HL equations coupled with an anisotropic fluid with heat flow. It is found that the scalars $K$ and $K_{ij}K^{ij}$ are singular only at the center for the de Sitter Schwarzschild solution, but singular at both the center and $ r = (3M/|\Lambda|)^{1/3}$ for the anti-de Sitter Schwarzschild solution. The singularity at $ r = (3M/|\Lambda|)^{1/3}$ is absent in GR. In addition, all the generalized LMP solutions have two scalar curvature singularities, located at either $r = 0$ and $r=r_{s} > 0$, or $r=r_{1}$ and $r= r_{2}$ with $r_{2} > r_{1} > 0$, or $r=r_{s} > 0$ and $r = \infty$, depending on the choice of the free parameter $\lambda$.

In the Horava-Lifshitz theory of quantum gravity, two conditions -- detailed balance and projectability -- are usually assumed. The breaking of projectability simplifies the theory, but it leads to serious problems with the theory. The breaking of detailed balance leads to a more complicated form of the theory, but it appears to resolve some of the problems. Sotiriou, Visser and Weinfurtner formulated the most general theory of Horava-Lifshitz type without detailed balance. We compute the linear scalar perturbations of the FRW model in this form of HL theory. We show that the higher-order curvature terms in the action lead to a gravitational effective anisotropic stress on small scales. Specializing to a Minkowski background, we study the spin-0 scalar mode of the graviton, using a gauge-invariant analysis, and find that it is stable in both the infrared and ultraviolet regimes for $0 \le \xi \le 2/3$. However, in this parameter range the scalar mode is a ghost.

In this paper, we study the implementation of brane worlds in type II string theory. Starting with the NS/NS sector of type II string, we first compactify the $(D+d_{+} + d_{-})$-dimensional spacetime, and reduce the corresponding action to a D-dimensional effective action, where the topologies of $M_{d_{+}}$ and $M_{d_{-}}$ are arbitrary. We further compactify one of the $(D-1)$ spatial dimensions on an $S^{1}/Z_{2}$ orbifold, and derive the gravitational and matter field equations both in the bulk and on the branes. Then, we investigate two key issues in such a setup: (i) the radion stability and radion mass; and (ii) the localization of gravity, and the corresponding Kaluza-Klein (KK) modes. We show explicitly that the radion is stable and its mass can be in the order of $GeV$. In addition, the gravity is localized on the visible brane, and its spectrum of the gravitational KK towers is discrete and can have a mass gap of $TeV$, too. The high order Yukawa corrections to the 4-dimensional Newtonian potential is exponentially suppressed, and can be negligible. Applying such a setup to cosmology, we obtain explicitly the field equations in the bulk and the generalized Friedmann equations on the branes.

The growth rate of matter perturbation and the expansion rate of the Universe can be used to distinguish modified gravity and dark energy models in explaining cosmic acceleration. We explore here the inclusion of spatial curvature into the growth factor. We expand previous results using the approximation $\Omega_{m}^\gamma$ and then suggest a new form, $f_a=\Omega_m^\gamma+(\gamma-4/7)\Omega_k$, as an approximation for the growth factor when the curvature $\Omega_k$ is not negligible, and where the growth index $\gamma$ is usually model dependent. The expression recovers the standard results for the curved and flat $\Lambda$CDM and Dvali-Gabadadze-Porrati models. Using the best fit values of $\Omega_{m0}$ and $\Omega_{k0}$ to the expansion/distance measurements from Type Ia supernovae, baryon acoustic oscillation, WMAP5, and $H(z)$ data, we fit the growth index parameter to current growth factor data and obtain $\gamma_{\Lambda}(\Omega_{k} \not= 0) = 0.65^{+0.17}_{-0.15}$ and $\gamma_{DGP}(\Omega_{k} \not= 0) = 0.53^{+0.14}_{-0.12}$. For the $\Lambda$CDM model, the 1-$\sigma$ observational bounds are found consistent with theoretical value, unlike the case for the Dvali-Gabadadze-Porrati model. We also find that the current data we used is not enough to put significant constraints when the 3 parameters in $f_a$ are fit simultaneously. Importantly, we find that, in the presence of curvature, the analytical expression proposed for $f_a$ provides a better fit to the growth factor than other forms and should be useful for future high precision missions and studies.

We study thermodynamics of cosmological models in the Horava-Lifshitz theory of gravity, and systematically investigate the evolution of the universe filled with a perfect fluid that has the equation of state $p=w\rho$, where $p$ and $\rho$ denote, respectively, the pressure and energy density of the fluid, and $w$ is an arbitrary real constant. Depending on specific values of the free parameters involved in the models, we classify all of them into various cases. In each case the main properties of the evolution are studied in detail, including the periods of deceleration and/or acceleration, and the existence of big bang, big crunch, and big rip singularities. We pay particular attention on models that may give rise to a bouncing universe.

Dynamical models of prototype gravastars are constructed and studied. The models are the Visser-Wiltshire three-layer gravastars, in which an infinitely thin spherical shell of a perfect fluid with the equation of state $p = (1-\gamma)\sigma$ divides the whole spacetime into two regions, where the internal region is de Sitter, and the external is Schwarzschild. When $\gamma < 1$ and $\Lambda \not= 0$, it is found that in some cases the models represent stable gravastars, and in some cases they represent "bounded excursion" stable gravastars, where the thin shell is oscillating between two finite radii, while in some other cases they collapse until the formation of black holes. However, when $\gamma \ge 1$, even with $\Lambda \not= 0$, only black holes are found. In the phase space, the region for both stable gravastars and "bounded excursion" gravastars is very small in comparison to that of black holes, although it is not completely empty.

Spherically symmetric expansionfree distributions are systematically studied. The whole set of field equations and junction conditions are presented for a general distribution of dissipative anisotropic fluid (principal stresses unequal), and the expansionfree condition is integrated. In order to understand the physical meaning of expansionfree motion, two different definitions for the radial velocity of a fluid element are discussed. It is shown that the appearance of a cavity is inevitable in the expansionfree evolution. The nondissipative case is considered in detail and the Skripkin model is recovered.

Orbifold branes in string theory are investigated, and the general field equations on the branes are given explicitly for type II and heterotic string. It is shown that the effective cosmological constant on each of the two branes can be easily lowered to its current observational value, using large extra dimensions. This is also true for type I string. The radion stability is studied by using the Goldberger-Wise mechanism, and shown explicitly that it is stable. Therefore, brane world of string theory provides a viable and built-in mechanism for solving the long-standing cosmological constant problem. Applying the formulas to cosmology, we obtain the generalized Friedmann equations on each of the two branes.

Recently, in Gong {\em et al} \cite{GWW07} and Wang and Santos \cite{WS07} it was shown that the effective cosmological constant on each of the two orbifold branes can be easily lowered to its current observational value, by using the large extra dimensions in the framework of both M-Theory and string theory on $S^{1}/Z_{2}$. In this paper, we study the current acceleration of the universe, using the formulas developed in \cite{WS07}. We first construct explicitly time-dependent solution to the 10-dimensional bulk of the Neveu-Schwarz/Neveu-Schwarz sector, compactified on a 5-dimensional torus. Then, we write down the generalized Friedmann equations on each of the two dynamical branes, and fit the models to the 182 gold supernova Ia data and the BAO parameter from SDSS, using both of our MINUIT and Monte-Carlo Markov Chain (MCMC) codes. With the best fitting values of the parameters involved as initial conditions, we integrate the generalized Friedmann equations numerically and find the future evolution of the universe. We find that it depends on the choice of the radion potentials $V_{4}^{(I)} (I = 1, 2)$ of the branes. In particular, when choosing them to be the Goldberger-Wise potentials, $V_{4}^{(I)} = \lambda_{4}^{(I)} (\psi^{2} - {v_{I}}^{2})^{2}$, we find that the current acceleration of the universe driven by the effective cosmological constant is only temporary. Due to the effects of the potentials, the universe will be finally in its decelerating expansion phase again. We also study the proper distance between the two branes, and find that it remains almost constant during the whole future evolution of the universe in all the models considered.

We use the Monte-Carlo Markov Chain method to explore the dark energy property and the cosmic curvature by fitting two popular dark energy parameterizations to the observational data. The new 182 gold supernova Ia data and the ESSENCE data both give good constraint on the DE parameters and the cosmic curvature for the dark energy model $w_0+w_a z/(1+z)$. The cosmic curvature is found to be $|\Omega_k|\la 0.03$. For the dark energy model $w_0+w_a z/(1+z)^2$, the ESSENCE data gives better constraint on the cosmic curvature and we get $|\Omega_k|\leq 0.02$.

Orbifold branes are studied in the framework of the 11-dimensional Horava-Witten heterotic M-Theory. It is found that the effective cosmological constant can be easily lowered to its current observational value by the mechanism of large extra dimensions. The domination of this constant over the evolution of the universe is only temporarily. Due to the interaction of the bulk and the branes, the universe will be in its decelerating expansion phase again in the future, whereby all problems connected with a far future de Sitter universe are resolved.

Although there is mounting observational evidence that the cosmic expansion is undergoing a late-time acceleration, the physical mechanism behind such a phenomenon is yet unknown. In this paper, we investigate a holographic dark energy (HDE) model with interaction between the components of the dark sector in the light of current cosmological observations. We use both the new \emph{gold} sample of 182 type Ia supernovae (SNe Ia) and the 192 SNe Ia ESSENCE data, the baryon acoustic oscillation measurement from the Sloan Digital Sky Survey and the shift parameter from the three-year Wilkinson Microwave Anisotropy Probe data. In agreement with previous results, we show that these observations suggest a very weak coupling or even a noninteracting HDE. The phantom crossing behavior in the context of these scenarios is also briefly discussed.

With the help of a masslike function which has dimension of energy and equals to the Misner-Sharp mass at the apparent horizon, we show that the first law of thermodynamics of the apparent horizon $dE=T_AdS_A$ can be derived from the Friedmann equation in various theories of gravity, including the Einstein, Lovelock, nonlinear, and scalar-tensor theories. This result strongly suggests that the relationship between the first law of thermodynamics of the apparent horizon and the Friedmann equation is not just a simple coincidence, but rather a more profound physical connection.

The energy conditions provide a very promising model-independent study of the current acceleration of the universe. However, in order to connect these conditions with observations, one often needs first to integrate them, and then find the corresponding constraints on some observational variables, such as the distance modulus. Those integral forms can be misleading, and great caution is needed when one interprets them physically. A typical example is that the transition point of the deceleration parameter $q(z)$ is at about $z \simeq 0.76$ in the $\Lambda$CDM model. However, with the same model when we consider the dimensionless Hubble parameter $E(z)$, which involves the integration of $q(z)$, we find that $E(z)$ does not cross the line of $q(z) = 0$ before $z = 2$. Therefore, to get the correct result, we cannot use the latter to determine the transition point. With these in mind, we carefully study the constraints from the energy conditions, and find that, among other things, the current observational data indeed strongly indicate that our universe has ocne experienced an accelerating expansion phase between the epoch of galaxy formation and the present.

The energy conditions give upper bounds on the luminosity distance. We apply these upper bounds to the 192 essence supernova Ia data to show that the Universe had experienced accelerated expansion. This conclusion is drawn directly from the distance modulus-reshift graph. In addition to be a very simple method, this method is also totally independent of any cosmological model. From the degeneracy of the distance modulus at low redshift, we argue that the choice of $w_0$ for probing the property of dark energy is misleading. One explicit example is used to support this argument.

We have investigated the thermodynamical properties of dark energy. Assuming that the dark energy temperature $T\sim a^{-n}$ and considering that the volume of the Universe enveloped by the apparent horizon relates to the temperature, we have derived the dark energy entropy. For dark energy with constant equation of state $w>-1$ and the generalized Chaplygin gas, the derived entropy can be positive and satisfy the entropy bound. The total entropy, including those of dark energy, the thermal radiation and the apparent horizon, satisfies the generalized second law of thermodynamics. However, for the phantom with constant equation of state, the positivity of entropy, the entropy bound, and the generalized second law cannot be satisfied simultaneously.

We construct a class of analytic solutions with two free parameters to the five-dimensional Einstein field equations, which represents the collision of two timelike 3-branes. We study the local and global properties of the spacetime, and find that spacelike singularities generically develop after the collision, due to the mutual focus of the two branes. Non-singular spacetime can be constructed only in the case where both of the two branes violate the energy conditions.

The new 182 gold supernova Ia data, the baryon acoustic oscillation measurement and the shift parameter determined from the Sloan Digital Sky Survey and the three-year Wilkinson Microwave Anisotropy Probe data are combined to reconstruct the dark energy equation of state parameter $w(z)$ and the deceleration parameter $q(z)$. We find that the strongest evidence of acceleration happens around the redshift $z\sim 0.2$ and the stringent constraints on $w(z)$ lie in the redshift range $z\sim 0.2-0.5$. At the sweet spot, $-1.2<w(z)<-0.6$ for the dark energy parametrization $w(z)=w_0+w_a z/(1+z)^2$ at the $3\sigma$ confidence level. The transition redshift $z_t$ when the Universe underwent the transition from deceleration to acceleration is derived to be $z_t=0.36^{+0.23}_{-0.08}$. The combined data is also applied to find out the geometry of the Universe, and we find that at the $3\sigma$ confidence level, $|\Omega_k|\alt 0.05$ for the simple one parameter dark energy model, and $-0.064<\Omega_k<0.028$ for the $\Lambda$CDM model.

We have investigated the thermodynamical properties of the Universe with dark energy. Adopting the usual assumption in deriving the constant co-moving entropy density that the physical volume and the temperature are independent, we observed some strange thermodynamical behaviors. However, these strange behaviors disappeared if we consider the realistic situation that the physical volume and the temperature of the Universe are related. Based on the well known correspondence between the Friedmann equation and the first law of thermodynamics of the apparent horizon, we argued that the apparent horizon is the physical horizon in dealing with thermodynamics problems. We have concentrated on the volume of the Universe within the apparent horizon and considered that the Universe is in thermal equilibrium with the Hawking temperature on the apparent horizon. For dark energy with $w\ge -1$, the holographic principle and the generalized second law are always respected.

The Taiwan-America Occultation Survey (TAOS) aims to determine the number of small icy bodies in the outer reach of the Solar System by means of stellar occultation. An array of 4 robotic small (D=0.5 m), wide-field (f/1.9) telescopes have been installed at Lulin Observatory in Taiwan to simultaneously monitor some thousand of stars for such rare occultation events. Because a typical occultation event by a TNO a few km across will last for only a fraction of a second, fast photometry is necessary. A special CCD readout scheme has been devised to allow for stellar photometry taken a few times per second. Effective analysis pipelines have been developed to process stellar light curves and to correlate any possible flux changes among all telescopes. A few billion photometric measurements have been collected since the routine survey began in early 2005. Our preliminary result of a very low detection rate suggests a deficit of small TNOs down to a few km size, consistent with the extrapolation of some recent studies of larger (30--100 km) TNOs.

We study the existence of black holes in a homogeneous and isotropic expanding Friedmann-Robertson-Walker (FRW) universe dominated by dark energy. We show that black holes can exist in such a universe by considering some specific McVittie solutions. Although these solutions violate all three energy conditions, the FRW background does satisfy the weak energy condition.

We discuss the phenomenological model in which the potential energy of the quintessence field depends linearly on the energy density of the spatial curvature. We find that the pressure of the scalar field takes a different form when the potential of the scalar field also depends on the scale factor and the energy momentum tensor of the scalar field can be expressed as the form of a perfect fluid. A general coupling was proposed to explain the current accelerating expansion of the Universe and solve the fine-tuning problem.

The compactification of M theory with time dependent hyperbolic internal space gives an effective scalar field with exponential potential which provides a transient acceleration in Einstein frame in four dimensions. Ordinary matter and radiation are present in addition to the scalar field coming from compactification. We find that we have to fine-tune the initial conditions of the scalar field so that our Universe experiences acceleration now. During the evolution of our Universe, the volume of the internal space increases about 12 times. The time variation of the internal space results in a large time variation of the fine structure constant which violates the observational constraint on the variation of the fine structure constant. The large variation of the fine structure constant is a generic feature of transient acceleration models.

We show that the most general dark energy model that possesses a scaling solution $\rho_\phi\propto a^n$ is the k-essence model, which includes both of the quintessence and tachyon models. The exact scaling solutions are then derived. The potential that gives the tracking solution in which dark energy exactly tracks the background matter field is the inverse squared potential. The quintessence field with exponential potential can be obtained from the k-essence field with the inverse squared potential. We also find the fixed points and study their main properties, whereby the scalar field dominant fixed point is identified.

We propose a new parametrization of the deceleration parameter to study its time-variation behavior. The advantage of parameterizing the deceleration parameter is that we do not need to assume any underlying theory of gravity. By fitting the model to the 157 gold sample supernova Ia data, we find strong evidence that the Universe is currently accelerating and it accelerated in the past. By fitting the model to the 115 nearby and Supernova Legacy Survey supernova Ia data, the evidence that the Universe is currently accelerating is weak, although there is still a strong evidence that the Universe once accelerated in the past. The results obtained from the 157 gold sample supernova Ia data and those from the 115 supernova Ia data are not directly comparable because the two different data sets measure the luminosity distance up to different redshifts. We then use the Friedmann equation and a dark energy parametrization to discuss the same problem. When we fit the model to the supernova Ia data alone, we find weak evidence that the Universe is accelerating and the current matter density is higher than that measured from other experiments. After we add the Sloan Digital Sky Survey data to constrain the dark energy model, we find that the behavior of the deceleration parameter is almost the same as that obtained from parameterizing the deceleration parameter.

The gravitational collapse of a spherically symmetric star, made of a dust fluid, $\rho_{DM}$, in a background of dark energy, $p = w\rho,\; (w < -1/3)$, is studied. It is found that when only dark energy is present, black holes are {\em never} formed. When both of them are present, black holes can be formed, due to the condensation of the dust fluid. Initially the dust fluid may not play an important role, but, as time increases, it will dominate the collapse and finally leads to formation of black holes. This result remains true even when the interaction between the dust fluid and dark energy does not vanish. When $w < -1$ (phantoms), some models can also be interpreted as representing the death of a white hole that ejects both dust and phantoms. The ejected matter re-collapses to form a black hole.

All homothetic self-similar solutions of the Brans-Dicke scalar field in three-dimensional spacetime with circular symmetry are found in closed form.

Gauge choice for a spherically symmetric 3-brane embedded in a D-dimensional bulk with arbitrary matter fields on and off the brane is studied. It is shown that Israel's junction conditions across the brane restrict severely the dependence of the matter fields on the spacetime coordinates. As examples, a scalar field or a Yang-Mills potential can be only either time-dependent or radial-coordinate dependent for the chosen gauge, while for a perfect fluid it must be co-moving.

Recently, the gravitational collapse of an infinite cylindrical thin shell of matter in an otherwise empty spacetime with two hypersurface orthogonal Killing vectors was studied by Gon\c{c}alves [Phys. Rev. {\bf D65}, 084045 (2002).]. By using three "alternative" criteria for trapped surfaces, the author claimed to have shown that {\em they can never form either outside or on the shell, regardingless of the matter content for the shell, except at asymptotical future null infinite}. Following Penrose's original idea, we first define trapped surfaces in cylindrical spacetimes in terms of the expansions of null directions orthogonal to the surfaces, and then show that the first criterion used by Gon\c{c}alves is incorrect. We also show that his analysis of non-existence of trapped surfaces in vacuum is incomplete. To confirm our claim, we present an example that is a solution to the vacuum Einstein field equations and satisfies all the regular conditions imposed by Gon\c{c}alves. After extending the solution to the whole spacetime, we show explicitly that trapped surfaces exist in the extended region.

We analyze an oscillating universe model in brane world scenario. The oscillating universe cycles through a series of expansions and contractions and its energy density is dominated by dust matter at early-time expansion phase and by phantom dark energy at late-time expansion phase. We find that the period of the oscillating universe is not sensitive to the tension of the brane, but sensitive to the equation-of-state parameter $w$ of the phantom dark energy, and the ratio of the period to the current Hubble age approximately varies from 3 to 9 when the parameter $w$ changes from -1.4 to -1.1. The fraction of time that the oscillating universe spends in the coincidence state is also comparable to the period of the oscillating universe. This result indicates that the coincidence problem can be significantly ameliorated in the oscillating universe without singularity.

Here we study some general properties of spherical shear-free collapse. Its general solution when imposing conformal flatness is reobtained [1,2] and matched to the outgoing Vaidya spacetime. We propose a simple model satisfying these conditions and study its physical consequences. Special attention deserve, the role played by relaxational processes and the conspicuous link betweeen dissipation and density inhomogeneity.

Recent type Ia supernovas data seemingly favor a dark energy model whose equation of state $w(z)$ crosses -1 very recently, which is a much more amazing problem than the acceleration of the universe. In this paper we show that it is possible to realize such a crossing without introducing any phantom component in a Gauss-Bonnet brane world with induced gravity, where a four dimensional curvature scalar on the brane and a five dimensional Gauss-Bonnet term in the bulk are present. In this realization, the Gauss-Bonnet term and the mass parameter in the bulk play a crucial role.

The problem of kink stability of isothermal spherical self-similar flow in newtonian gravity is revisited. Using distribution theory we first develop a general formula of perturbations, linear or non-linear, which consists of three sets of differential equations, one in each side of the sonic line and the other along it. By solving the equations along the sonic line we find explicitly the spectrum, $k$, of the perturbations, whereby we obtain the stability criterion for the self-similar solutions. When the solutions are smoothly across the sonic line, our results reduce to those of Ori and Piran. To show such obtained perturbations can be matched to the ones in the regions outside the sonic line, we study the linear perturbations in the external region of the sonic line (the ones in the internal region are identically zero), by taking the solutions obtained along the line as the boundary conditions. After properly imposing other boundary conditions at spatial infinity, we are able to show that linear perturbations, satisfying all the boundary conditions, exist and do not impose any additional conditions on $k$. As a result, the complete treatment of perturbations in the whole spacetime does not alter the spectrum obtained by considering the perturbations only along the sonic line.

We study a cosmological model in which phantom dark energy is coupled to dark matter by phenomenologically introducing a coupled term to the equations of motion of dark energy and dark matter. This term is parameterized by a dimensionless coupling function $\delta$, Hubble parameter and the energy density of dark matter, and it describes an energy flow between the dark energy and dark matter. We discuss two cases: one is the case where the equation-of-state $\omega_e$ of the dark energy is a constant; the other is that the dimensionless coupling function $\delta$ is a constant. We investigate the effect of the interaction on the evolution of the universe, the total lifetime of the universe, and the ratio of the period when the universe is in the coincidence state to its total lifetime. It turns out that the interaction will produce significant deviation from the case without the interaction.

We extend Green's function method developed by Stewart and Gong to calculate the power spectrum of density perturbation in the case with a time-dependent sound speed, and explicitly give the power spectrum and spectral index up to second-order corrections in the slow-roll expansion. The case of tachyon inflation is included as a special case.

The LASCO C2 and C3 coronagraphs recorded a unique coronal mass ejection on April 2, 1999. The event did not have the typical three-part CME structure and involved a small filament eruption without any visibile overlying streamer ejecta. The event exhibited an unusually clear signature of a wave propagating at the CME flanks. The speed and density of the CME front and flanks were consistent with the existence of a shock. To better establish the nature of the white light wave signature, we employed a simple MHD simulation using the LASCO measurements as constraints. Both the measurements and the simulation strongly suggest that the white light feature is the density enhancement from a fast-mode MHD shock. In addition, the LASCO images clearly show streamers being deflected when the shock impinges on them. It is the first direct imaging of this interaction.

The Taiwanese-American Occultation Survey (TAOS) will detect objects in the Kuiper Belt, by measuring the rate of occultations of stars by these objects, using an array of three to four 50cm wide-field robotic telescopes. Thousands of stars will be monitored, resulting in hundreds of millions of photometric measurements per night. To optimize the success of TAOS, we have investigated various methods of gathering and processing the data and developed statistical methods for detecting occultations. In this paper we discuss these methods. The resulting estimated detection efficiencies will be used to guide the choice of various operational parameters determining the mode of actual observation when the telescopes come on line and begin routine observations. In particular we show how real-time detection algorithms may be constructed, taking advantage of having multiple telescopes. We also discuss a retrospective method for estimating the rate at which occultations occur.

A solution to the Einstein field equations that represents a rigidly rotating dust accompanied by a thin matter shell of the same type is found.