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})$.

We report results from searches for new physics with low-energy electronic recoil data recorded with the XENON1T detector, with an exposure of 0.65 t-y and an unprecedentedly low background rate of $76 \pm 2_{stat}$ events/(t y keV) between 1-30 keV. An excess over known backgrounds is observed below 7 keV, rising towards lower energies and prominent between 2-3 keV. The solar axion model has a 3.5$\sigma$ significance, and a three-dimensional 90% confidence surface is reported for axion couplings to electrons, photons, and nucleons. This surface is inscribed in the cuboid defined by $g_{ae} < 3.7 \times 10^{-12}$, $g_{ae}g_{an}^{eff} < 4.6 \times 10^{-18}$, and $g_{ae}g_{a\gamma} < 7.6\times10^{-22}~{GeV}^{-1}$, and excludes either $g_{ae}=0$ or $g_{ae}g_{a\gamma}=g_{ae}g_{an}^{eff}=0$. The neutrino magnetic moment signal is similarly favored over background at 3.2$\sigma$ and a confidence interval of $\mu_{\nu} \in (1.4,~2.9)\times10^{-11}\mu_B$ (90% C.L.) is reported. Both results are in strong tension with stellar constraints. The excess can also be explained by $\beta$ decays of tritium at 3.2$\sigma$ significance with a corresponding tritium concentration in xenon of $(6.2 \pm 2.0) \times 10^{-25}$ mol/mol. Such a trace amount can be neither confirmed nor excluded with current knowledge of production and reduction mechanisms. The significances of the solar axion and neutrino magnetic moment hypotheses are decreased to 2.1$\sigma$ and 0.9$\sigma$, respectively, if an unconstrained tritium component is included in the fitting. With respect to bosonic dark matter, the excess favors a monoenergetic peak at ($2.3 \pm 0.2$) keV (68% C.L.) with a 3.0$\sigma$ global (4.0$\sigma$ local) significance over background. This analysis sets the most restrictive direct constraints to date on pseudoscalar and vector bosonic dark matter for most masses between 1 and 210 keV/c$^2$.

The transonic phenomenon of black hole accretion and the existence of the photon sphere are the characteristics of strong gravitational fields near a black hole horizon. In this work, we study spherical accretion flow onto a general parametrized spherically symmetric black hole spacetimes. For this purpose, we analyze the accretion process of various perfect fluids, such as the isothermal fluid of ultra-stiff, ultra-relativistic, and sub-relativistic types and polytropic fluid, respectively. The influences of extra parameters beyond the Schwarzschild black hole in the general parameterized spherically symmetric black hole on the flow behaviors of the above-mentioned test fluids are studied in detail. In addition, by studying the accretion of ideal photon gas, we further discuss the correspondence between the sonic radius of accreting photon gas and the photon sphere for the general parameterized spherically symmetric black hole. Some possible future extensions of our analysis are also discussed.

We detail the sensitivity of the liquid xenon (LXe) DARWIN observatory to solar neutrinos via elastic electron scattering. We find that DARWIN will have the potential to measure the fluxes of five solar neutrino components: $pp$, $^7$Be, $^{13}$N, $^{15}$O and $pep$. The precision of the $^{13}$N, $^{15}$O and $pep$ components is hindered by the double-beta decay of $^{136}$Xe and, thus, would benefit from a depleted target. A high-statistics observation of $pp$ neutrinos would allow us to infer the values of the weak mixing angle, $\sin^2\theta_w$, and the electron-type neutrino survival probability, $P_e$, in the electron recoil energy region from a few keV up to 200 keV for the first time, with relative precision of 5% and 4%, respectively, at an exposure of 300 ty. An observation of $pp$ and $^7$Be neutrinos would constrain the neutrino-inferred solar luminosity down to 0.2%. A combination of all flux measurements would distinguish between the high (GS98) and low metallicity (AGS09) solar models with 2.1-2.5$\sigma$ significance, independent of external measurements from other experiments or a measurement of $^8$B neutrinos through coherent elastic neutrino-nucleus scattering in DARWIN. Finally, we demonstrate that with a depleted target DARWIN may be sensitive to the neutrino capture process of $^{131}$Xe.

Recently a novel four-dimensional Einstein-Gauss-Bonnet (4EGB) theory of gravity was proposed by Glavan and Lin [D. Glavan and C. Lin, Phys. Rev. Lett. 124, 081301 (2020)] which bypasses the Lovelock' s theorem with the re-scaling of the Gauss-Bonnet coupling constant $\alpha \to \alpha/(D-4)$ in the limit $D\to 4$. Here we study the physical properties of the electromagnetic radiation emitted from a thin accretion disk around the static spherically symmetric black hole in the 4EGB gravity. For this purpose, we assume the disk is in a steady-state and in hydrodynamic and thermodynamic equilibrium so that the emitted electromagnetic radiation is a black body spectrum. We study in detail the effects of the Gauss-Bonnet coupling constant $\alpha$ in 4EGB gravity on the energy flux, temperature distribution, and electromagnetic spectrum of the disk. It is shown that with the increases of the parameter $\alpha$, the energy flux, temperature distribution, and electromagnetic spectrum of the accretion disk all increase. Besides, we also show that the accretion efficiency increases as the growth of the parameter $\alpha$. Our results indicate that the thin accretion disk around the static spherically symmetric black hole in the 4EGB gravity is hotter, more luminosity, and more efficient than that around a Schwarzschild black hole with the same mass for a positive $\alpha$, while it is cooler, less luminosity, and less efficient for a negative $\alpha$.

In this paper, we study the motion of test particles around two exact charged black-hole solutions in Einstein-{\AE}ther theory. Specifically, we first consider the quasi-periodic oscillations (QPOs) and their resonances generated by the particle moving in the Einstein-{\AE}ther black hole and then turn to study the periodic orbits of the massive particles. For QPOs, we drop the usually adopted assumptions $\nu_U=\nu_\theta$, $\nu_L=\nu_r$, and $\nu_U/\nu_L=3/2$ with $\nu_U$ ($\nu_L$) and $\nu_r$ ($\nu_\theta$) being the upper (lower) frequency of QPOs and radial (vertical) epicyclic frequency of the orbiting particles, respectively. Instead, we put-forward a new working ansatz for which the Keplerian radius is much closer to that of the innermost stable circular orbit and explore in details the effects of the {\ae}ther field on the frequencies of QPOs. We then realize good curves for the frequencies of QPOs, which fit to data of three microquasars very well by ignoring any effects of rotation and magnetic fields. For the study of the periodic orbits of massive particle moving in the two Einstein-{\AE}ther black hole solutions, we first consider the geodesic equations in the two black holes and solve the geodesic equation for one of black hole (the black hole with \ae{}ther parameter $c_{14}=0$ but $c_{123}\neq 0$) analytically. The innermost stable circular orbits are also analyzed and we find the isco radius increases with increasing $c_{13}$ for the first type black hole while decreases with increasing $c_{14}$ for the second one. We also obtain several periodic orbits and find that they share similar taxonomy schemes as the periodic equatorial orbits in the Schwarzschild/Kerr metrics. The equations for null geodesics are also briefly considered, where we study circular photon orbits and bending angles for gravitational lensing.

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.

Xenon dual-phase time projection chambers designed to search for Weakly Interacting Massive Particles have so far shown a relative energy resolution which degrades with energy above $\sim$200 keV. This has limited their sensitivity in the search for rare events like the neutrinoless double-beta decay of $^{136}$Xe at its $Q$-value, $Q_{\beta\beta}\simeq$ 2.46 MeV. For the XENON1T dual-phase time projection chamber, we demonstrate that the relative energy resolution at 1 $\sigma/\mu$ is as low as (0.79$\pm$0.02) % in its one-ton fiducial mass, and for single-site interactions at $Q_{\beta\beta}$. We achieve this by using a signal correction method to rectify the saturation effects of the signal readout system. The very good energy resolution from keV to MeV energies demonstrated in XENON1T opens up new windows for the xenon dual-phase dark matter detectors to simultaneously search for other rare events.

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.

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

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

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

Direct dark matter detection experiments based on a liquid xenon target are leading the search for dark matter particles with masses above $\sim$ 5 GeV/c$^2$, but have limited sensitivity to lighter masses because of the small momentum transfer in dark matter-nucleus elastic scattering. However, there is an irreducible contribution from inelastic processes accompanying the elastic scattering, which leads to the excitation and ionization of the recoiling atom (the Migdal effect) or the emission of a Bremsstrahlung photon. In this letter, we report on a probe of low-mass dark matter with masses down to about 85 MeV/c$^2$ by looking for electronic recoils induced by the Migdal effect and Bremsstrahlung, using data from the XENON1T experiment. Besides the approach of detecting both scintillation and ionization signals, we exploit an approach that uses ionization signals only, which allows for a lower detection threshold. This analysis significantly enhances the sensitivity of XENON1T to light dark matter previously beyond its reach.

We report constraints on light dark matter (DM) models using ionization signals in the XENON1T experiment. We mitigate backgrounds with strong event selections, rather than requiring a scintillation signal, leaving an effective exposure of $(22 \pm 3)$ tonne-days. Above $\sim\!0.4\,\mathrm{keV}_\mathrm{ee}$, we observe $<1 \, \text{event}/(\text{tonne} \times \text{day} \times \text{keV}_\text{ee})$, which is more than one thousand times lower than in similar searches with other detectors. Despite observing a higher rate at lower energies, no DM or CEvNS detection may be claimed because we cannot model all of our backgrounds. We thus exclude new regions in the parameter spaces for DM-nucleus scattering for DM masses $m_\chi$ within $3-6\,\mathrm{GeV}/\mathrm{c}^2$, DM-electron scattering for $m_\chi > 30\,\mathrm{MeV}/\mathrm{c}^2$, and absorption of dark photons and axion-like particles for $m_\chi$ within $0.186 - 1 \, \mathrm{keV}/\mathrm{c}^2$.

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.

In this paper we present a new black hole solution surrounded by dark matter halo in the galactic center using the mass model of M87 and that coming from the Universal Rotation Curve (URC) dark matter profile representing family of spiral galaxies. In both cases the DM halo density is cored with a size $r_0$ and a central density $\rho_0$: $\rho(r)= \rho_0/(1+r/r_0)(1+(r/r_0)^2)$. Since $r_0\rho_0=120 M_{\odot}$/pc$^2$ [Donato et al., MNRAS, 397, 1169, 2009], then by varying the central density one can reproduce the DM profile in any spiral. Using the Newman-Jains method we extend our solution to obtain a rotating black hole surrounded by dark matter halo. We find that, the apparent shape of the shadow beside the black hole spin $a$, it also depends on the central density of the surrounded dark matter $\rho_0$. As a specific example we consider the galaxy M87, with a central density $\rho_0=6.9\,\times 10^{6} M_{\odot}$/kpc$^3$ and a core radius $r_0=91.2$ kpc. In the case of M87, our analyses show that the effect of dark matter on the size of the black hole shadow is almost negligible compared to the shadow size of the Kerr vacuum solution hence the angular diameter $ 42$ $\mu$as remains almost unaltered when the dark matter is considered. For a small totally dark matter dominated spiral such as UGC 7232, we find similar effect of dark matter on the shadow images compared to the M87. However, in specific conditions having a core radius comparable to the black hole mass and dark matter with very high density, we show that the shadow images decreases compared to the Kerr vacuum black hole. The effect of dark matter on the apparent shadow shape can shed some light in future observations as an indirect way to detect dark matter using the shadow images.

We report the first experimental results on spin-dependent elastic weakly interacting massive particle (WIMP) nucleon scattering from the XENON1T dark matter search experiment. The analysis uses the full ton year exposure of XENON1T to constrain the spin-dependent proton-only and neutron-only cases. No significant signal excess is observed, and a profile likelihood ratio analysis is used to set exclusion limits on the WIMP-nucleon interactions. This includes the most stringent constraint to date on the WIMP-neutron cross section, with a minimum of $6.3\times10^{-42}$ cm$^2$ at 30 GeV/c${}^2$ and 90% confidence level. The results are compared with those from collider searches and used to exclude new parameter space in an isoscalar theory with an axial-vector mediator.

We present first results on the scalar WIMP-pion coupling from 1 t$\times$yr of exposure with the XENON1T experiment. This interaction is generated when the WIMP couples to a virtual pion exchanged between the nucleons in a nucleus. In contrast to most non-relativistic operators, these pion-exchange currents can be coherently enhanced by the total number of nucleons, and therefore may dominate in scenarios where spin-independent WIMP-nucleon interactions are suppressed. Moreover, for natural values of the couplings, they dominate over the spin-dependent channel due to their coherence in the nucleus. Using the signal model of this new WIMP-pion channel, no significant excess is found, leading to an upper limit cross section of $6.4\times10^{-46}$ cm$^2$ (90 % confidence level) at 30 GeV/c$^2$ WIMP mass.

Recently, the nonperturbative quantization scheme of loop quantum gravity has been extended to the Brans-Dicke theory and the corresponding loop quantum Brans-Dicke cosmology has been derived, which provides an essential platform to explore inflationary models in this framework. In this paper, we consider two inflation models, the Starobinsky and $\alpha$-attractor inflation whose cosmological predictions are in excellent agreement with Planck data, and study systematically their pre-inflationary dynamics as well as the slow-roll inflation. We show that for both models, the background evolution of a flat Friedmann-Lema\^{i}tre-Robertson-Walker universe in general can be divided into three different phases: the pre-inflationary quantum phase, quantum-to-classical transition, and the slow-roll inflation. The pre-inflationary dynamics are dominated by the quantum geometry effects of loop quantum Brans-Dicke cosmology and the corresponding Universe could be either initially expanding or contracting, depending on the initial velocity of inflaton field. It is shown that the detailed evolution of pre-inflationary quantum phase also depend on specific inflation models. After the pre-inflationary quantum phase, the universe gradually evolves into the slow-roll inflation with some of initial conditions for Starobinsky and $\alpha$-attractor potentials. In addition, to be consistent with observational data, we also find the restricted parameter space of initial conditions that could produce at least $60$ $e$-folds during the slow-roll inflation.

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.

We report on a search for Weakly Interacting Massive Particles (WIMPs) using 278.8 days of data collected with the XENON1T experiment at LNGS. XENON1T utilizes a liquid xenon time projection chamber with a fiducial mass of $(1.30 \pm 0.01)$ t, resulting in a 1.0 t$\times$yr exposure. The energy region of interest, [1.4, 10.6] $\mathrm{keV_{ee}}$ ([4.9, 40.9] $\mathrm{keV_{nr}}$), exhibits an ultra-low electron recoil background rate of $(82\substack{+5 \\ -3}\textrm{ (sys)}\pm3\textrm{ (stat)})$ events/$(\mathrm{t}\times\mathrm{yr}\times\mathrm{keV_{ee}})$. No significant excess over background is found and a profile likelihood analysis parameterized in spatial and energy dimensions excludes new parameter space for the WIMP-nucleon spin-independent elastic scatter cross-section for WIMP masses above 6 GeV/c${}^2$, with a minimum of $4.1\times10^{-47}$ cm$^2$ at 30 GeV/c${}^2$ and 90% confidence level.

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

We report the first dark matter search results from XENON1T, a $\sim$2000-kg-target-mass dual-phase (liquid-gas) xenon time projection chamber in operation at the Laboratori Nazionali del Gran Sasso in Italy and the first ton-scale detector of this kind. The blinded search used 34.2 live days of data acquired between November 2016 and January 2017. Inside the (1042$\pm$12) kg fiducial mass and in the [5, 40] $\mathrm{keV}_{\mathrm{nr}}$ energy range of interest for WIMP dark matter searches, the electronic recoil background was $(1.93 \pm 0.25) \times 10^{-4}$ events/(kg $\times$ day $\times \mathrm{keV}_{\mathrm{ee}}$), the lowest ever achieved in a dark matter detector. A profile likelihood analysis shows that the data is consistent with the background-only hypothesis. We derive the most stringent exclusion limits on the spin-independent WIMP-nucleon interaction cross section for WIMP masses above 10 GeV/c${}^2$, with a minimum of 7.7 $\times 10^{-47}$ cm${}^2$ for 35-GeV/c${}^2$ WIMPs at 90% confidence level.

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.

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.

We present results of searches for vector and pseudo-scalar bosonic super-WIMPs, which are dark matter candidates with masses at the keV-scale, with the XENON100 experiment. XENON100 is a dual-phase xenon time projection chamber operated at the Laboratori Nazionali del Gran Sasso. A profile likelihood analysis of data with an exposure of 224.6 live days $\times$ 34\,kg showed no evidence for a signal above the expected background. We thus obtain new and stringent upper limits in the $(8-125)$\,keV/c$^2$ mass range, excluding couplings to electrons with coupling constants of $g_{ae} > 3\times10^{-13}$ for pseudo-scalar and $\alpha'/\alpha > 2\times10^{-28}$ for vector super-WIMPs, respectively. These limits are derived under the assumption that super-WIMPs constitute all of the dark matter in our galaxy.

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.

This paper aims to obtain an analytic approximation to the evolution of circular orbits governed by the Earth's $J_2$ and the luni-solar gravitational perturbations. Assuming that the lunar orbital plane coincides with the ecliptic plane, Allan and Cook (1964) derived an analytic solution to the orbital plane evolution of circular orbits. Using their result as an intermediate solution, we establish an approximate analytic model with lunar orbital inclination and its node regression be taken into account. Finally, an approximate analytic expression is derived, which is accurate compared to the numerical results except for the resonant cases when the period of the reference orbit approximately equals the integer multiples (especially 1 or 2 times) of lunar node regression period.

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.

DARk matter WImp search with liquid xenoN (DARWIN) will be an experiment for the direct detection of dark matter using a multi-ton liquid xenon time projection chamber at its core. Its primary goal will be to explore the experimentally accessible parameter space for Weakly Interacting Massive Particles (WIMPs) in a wide mass-range, until neutrino interactions with the target become an irreducible background. The prompt scintillation light and the charge signals induced by particle interactions in the xenon will be observed by VUV sensitive, ultra-low background photosensors. Besides its excellent sensitivity to WIMPs above a mass of 5 GeV/c2, such a detector with its large mass, low-energy threshold and ultra-low background level will also be sensitive to other rare interactions. It will search for solar axions, galactic axion-like particles and the neutrinoless double-beta decay of 136-Xe, as well as measure the low-energy solar neutrino flux with <1% precision, observe coherent neutrino-nucleus interactions, and detect galactic supernovae. We present the concept of the DARWIN detector and discuss its physics reach, the main sources of backgrounds and the ongoing detector design and R&D efforts.

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.

Astrophotonics is the next-generation approach that provides the means to miniaturize near-infrared (NIR) spectrometers for upcoming large telescopes and make them more robust and inexpensive. The target requirements for our spectrograph are: a resolving power of about 3000, wide spectral range (J and H bands), free spectral range of about 30 nm, high on-chip throughput of about 80% (-1dB) and low crosstalk (high contrast ratio) between adjacent on-chip wavelength channels of less than 1% (-20dB). A promising photonic technology to achieve these requirements is Arrayed Waveguide Gratings (AWGs). We have developed our first generation of AWG devices using a silica-on-silicon substrate with a very thin layer of silicon-nitride in the core of our waveguides. The waveguide bending losses are minimized by optimizing the geometry of the waveguides. Our first generation of AWG devices are designed for H band and have a resolving power of around 1500 and free spectral range of about 10 nm around a central wavelength of 1600 nm. The devices have a footprint of only 12 mm x 6 mm. They are broadband (1450-1650 nm), have a peak on-chip throughput of about 80% (-1 dB) and contrast ratio of about 1.5% (-18 dB). These results confirm the robustness of our design, fabrication and simulation methods. Currently, the devices are designed for Transverse Electric (TE) polarization and all the results are for TE mode. We are developing separate J- and H-band AWGs with higher resolving power, higher throughput and lower crosstalk over a wider free spectral range to make them better suited for astronomical applications.

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.

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.

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

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

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

We study Hawking-like radiation in a Friedmann-Robertson-Walker (FRW) universe using the quasi-classical WKB/tunneling method which pictures this process as a "tunneling" of particles from behind the apparent horizon. The correct temperature of the Hawking-like radiation from the FRW spacetime is obtained using a canonical invariant tunneling amplitude. In contrast to the usual quantum mechanical WKB/tunneling problem where the tunneling amplitude has only a spatial contribution, we find that the tunneling amplitude for FRW spacetime (i.e. the imaginary part of the action) has both spatial and temporal contributions. In addition we study back reaction and energy conservation of the radiated particles and find that the tunneling probability and change in entropy, ${\cal S}$ are related by the relationship: $\Gamma\propto\exp[-\Delta {\cal S}]$ which differs from the standard result $\Gamma\propto\exp[\Delta {\cal S}]$. By regarding the whole FRW universe as an isolated adiabatic system the change in the total entropy is zero. Then splitting the entropy between interior and exterior parts of the horizon ($\Delta {\cal S}_{total}=\Delta {\cal S}_{int} + \Delta {\cal S}_{ext}=0$), we can explain the origin of the minus sign difference with the usual result: our $\Delta {\cal S}$ is for the interior region while the standard result from black hole physics is for the exterior region.

Recently, a Hamilton-Jacobi method beyond semiclassical approximation in black hole physics was developed by \emph{Banerjee} and \emph{Majhi}\cite{beyond0}. In this paper, we generalize their analysis of black holes to the case of Friedmann-Robertson-Walker (FRW) universe. It is shown that all the higher order quantum corrections in the single particle action are proportional to the usual semiclassical contribution. The corrections to the Hawking-like temperature and entropy of apparent horizon for FRW universe are also obtained. In the corrected entropy, the area law involves logarithmic area correction together with the standard inverse power of area term.

The influence of the generalized uncertainty principle (GUP) and extended uncertainty principle (EUP) on the thermodynamics of the Friedmann-Robertson-Walker (FRW) universe has been investigated. It is shown that the entropy of the apparent horizon of the FRW universe gets a correction if one considers the effect of the GUP or EUP. Moreover, starting with the modified entropy and applying the first law of thermodynamics, $dE=TdS$, to the apparent horizon of the FRW universe, we obtain the modified Friedmann equations. The influence of the GUP or EUP on the thermodynamics of the FRW universe provides a deep insight into the understanding of the quantum gravity or large length scale corrections to the dynamics of the FRW universe.

A model is proposed for the observed combination of power-law and thermal emission of keV X-rays from rotationally powered pulsars. For gamma-ray pulsars with accelerators very many stellar radii above the neutron star surface, 100 MeV curvature gamma-rays from $e^{-}$ or $e^{+}$ flowing starward out of such accelerators are converted to electron-positron pairs on closed field lines all around the star. These pairs strongly affect X-ray emission from near the star in two ways. (1) The pairs are a source of synchrotron emission immediately following their creation in regions where $B \sim 10^{10}$ G. This emission, in the photon energy range 0.1 keV $\lesssim E_{X} \lesssim$ 5 MeV, has a power-law spectrum with energy index 0.5 and X-ray luminosity that depends on the backflow current, and is typically $\sim 10^{33}$ \lum. (2) The pairs ultimately form a cyclotron resonance ``blanket'' surrounding the star except for two holes along the open field line bundles which pass through it. In such a blanket the gravitational pull on electron-positron pairs toward the star is balanced by the hugely amplified push of outflowing surface emitted X-rays wherever cyclotron resonance occurs. Because of it the neutron star is surrounded by a leaky ``hohlraum'' of hot blackbody radiation with two small holes, which prevents direct X-ray observation of a heated polar cap of a gamma-ray pulsar. Weakly spin-modulated radiation from the blanket together with more strongly spin-modulated radiation from the holes through it would then dominate observed low energy (0.1--10 keV) emission. For non-gamma-ray pulsars, in which no such accelerators with their accompanying extreme relativistic backflow toward the star are ...

Spinning superfluid neutrons in the core of a neutron star interact strongly with co-existing superconducting protons. One consequence is that the outward(inward) motion of core superfluid neutron vortices during spin-down(up) of a neutron star may alter the core's magnetic field. Such core field changes are expected to result in movements of the stellar crust and changes in the star's surface magnetic field which reflect those in the core below. Observed magnitudes and evolution of the spin-down indices of canonical pulsars are understood as a consequence of such surface field changes. If the growing crustal strains caused by the changing core magnetic field configuration in canonical spinning-down pulsars are relaxed by large scale crust-cracking events, special properties are predicted for the resulting changes in spin-period. These agree with various glitch observations, including glitch activity, permanent shifts in spin-down rates after glitches in young pulsars, the intervals between glitches, families of glitches with different magnitudes in the same pulsar, the sharp drop in glitch intervals and magnitudes as pulsar spin-periods approach 0.7s, and the general absence of glitching beyond this period.