
We investigate the constraints on the reionization history of the Universe
from a joint analysis of the cosmic microwave background and neutral hydrogen
fraction data. The $\tanh$ parameterization and principal component analysis
methods are applied to the reionization history respectively. The commonly used
$\tanh$ parameterization is over simplistic when the neutral hydrogen fraction
data is taken into account. Using the principal component analysis method, the
reconstructed reionization history is consistent with the neutral hydrogen
fraction data. With PCA method, we reconstruct the neutral hydrogen fraction at
$z=9.75$ as $x_\mathrm{HI}=0.69^{+0.30}_{0.32}$ for $6<z<20$ range
reconstruction, and $x_\mathrm{HI}=0.76^{+0.22}_{0.27}$ for $6<z<30$ range
reconstruction. These results suggest that the Universe began to reionize at
redshift no later than $z=10$ at $95\%$ confidence level.

The gravitational wave (GW) as a standard siren directly determines the
luminosity distance from the gravitational waveform without reference to the
specific cosmological model, of which the redshift can be obtained separately
by means of the electromagnetic counterpart like GW events from binary neutron
stars and massive black hole binaries (MBHBs). To see to what extent the
standard siren can reproduce the presumed dipole anisotropy written in the
simulated data of standard siren events from typical configurations of GW
detectors, we find that (1) for the Laser Interferometer Space Antenna with
different MBHB models during fiveyear observations, the cosmic isotropy can be
ruled out at $3\sigma$ confidence level (C.L.) and the dipole direction can be
constrained roughly around $20\%$ at $2\sigma$ C.L., as long as the dipole
amplitude is larger than $0.03$, $0.06$ and $0.025$ for MBHB models Q3d, pop
III and Q3nod with increasing constraining ability, respectively; (2) for
Einstein Telescope with no less than $200$ standard siren events, the cosmic
isotropy can be ruled out at $3\sigma$ C.L. if the dipole amplitude is larger
than $0.06$, and the dipole direction can be constrained within $20\%$ at
$3\sigma$ C.L. if the dipole amplitude is near $0.1$; (3) for the DeciHertz
Interferometer Gravitational wave Observatory with no less than $100$ standard
siren events, the cosmic isotropy can be ruled out at $3\sigma$ C.L. for dipole
amplitude larger than $0.03$ , and the dipole direction can even be constrained
within $10\%$ at $3\sigma$ C.L. if the dipole amplitude is larger than $0.07$.
Our work manifests the promising perspective of the constraint ability on the
cosmic anisotropy from the standard siren approach.

We study the hydrodynamics of bubble expansion in cosmological firstorder
phase transition in the FriedmannLema\^{\i}treRobertsonWalker (FLRW)
background with probe limit. Different from previous studies for fast
firstorder phase transition in flat background, we find that, for slow
firstorder phase transition in FLRW background with a given peculiar velocity
of the bubble wall, the efficiency factor of energy transfer into bulk motion
of thermal fluid is significantly reduced, thus decreasing the
previouslythought dominated contribution from sound wave to the stochastic
gravitationalwave background.

Hydrodynamics with both vector and axial currents is under study within a
simple holographic model, consisting of canonical $U(1)_V\times U(1)_A$ gauge
fields in asymptotically AdS$_5$ black brane. When gravitational backreaction
is taken into account, the chiral electric separation effect (CESE), namely the
generation of an axial current as the response to an external electric field,
is realized naturally. Via fluid/gravity calculation, firstorder transports in
the hydrodynamic constitutive relations are evaluated analytically in terms of
vector chemical potential $\mu$ and axial chemical potential $\mu_5$, as well
as the temperature $T$. Apart from the proportionality factor $\mu\mu_5$, the
CESE conductivity is found to depend on the dimensionless quantities $\mu/T$
and $\mu_5/T$ nontrivially. As a complementary study, frequencydependence of
the transports within this holographic model is revealed through the linear
response analysis, demonstrating agreement with the results obtained from
fluid/gravity calculation.

In recent years, superfluid dark matter (SfDM) has become a competitive model
of emergent modified Newtonian dynamics (MOND) scenario: MOND phenomenons
naturally emerge as a derived concept due to an extra force mediated between
baryons by phonons as a result of axionlike particles condensed as superfluid
at galactic scales; Beyond galactic scales, these axionlike particles behave as
normal fluid without phononmediated MONDlike force between baryons, therefore
SfDM also maintains the usual success of $\Lambda$CDM at cosmological scales.
In this paper, we use gravitational waves (GWs) to probe the relevant parameter
space of SfDM. GWs through BoseEinstein condensate (BEC) could propagate with
a speed slightly deviation from the speedoflight due to the change in the
effective refractive index, which depends on the SfDM parameters and GWsource
properties. We find that Five hundred meter Aperture Spherical Telescope
(FAST), Square Kilometre Array (SKA) and International Pulsar Timing Array
(IPTA) are the most promising means as GW probe of relevant parameter space of
SfDM. Future spacebased GW detectors are also capable of probing SfDM if a
multimessenger approach is adopted.

We study the production of gravitational waves during oscillations of the
inflaton around the minimum of a cuspy potential after inflation. We find that
a cusp in the potential can trigger copious oscillon formation, which sources a
characteristic energy spectrum of gravitational waves with double peaks. The
discovery of such a doublepeak spectrum could test the underlying inflationary
physics.

In this paper, we simulate the SuperEddington accreting massive black holes
(SEAMBHs) as the candles to probe the cosmology for the first time. SEAMBHs
have been demonstrated that they can provide a new tool for estimating
cosmological distance. Thus we create a series of mock data sets of SEAMBHs
specially in hight redshift region to check the abilities of them on probing
the cosmology. To fulfill the potential of the SEAMBHs on the cosmology, we
apply the simulated data to three projects. The first is exploring their
abilities to constrain the cosmological parameters, in which we combine
different data sets of current observations such as the cosmic microwave
background (CMB) from {\it Planck} and type Ia supernovae (SNe Ia) from Joint
Lightcurve Analysis (JLA). We find that the high redshift SEAMBHs can help to
break the degeneracies of the background cosmological parameters constrained by
{\it Planck} and JLA, thus give a much tighter constraints of the cosmological
parameters. The second is using the high redshift SEAMBHs as the complements of
the low redshift JLA to constrain the early expansion rate. Our results show
that these high redshift SEAMBHs are very powerful on constraining the early
Hubble rate, thus can give us more information of the expansion history of our
Universe, which is also crucial for the test of $\Lambda$CDM model in high
redshift region. Finally, we check the SEAMBHs candles' abilities to
reconstruct the equation of state of dark energy at high redshift. In summary,
our results show that the SEAMBHs, as the rare candles in high redshift region,
can provide us a new and independent observation to probe the evolution of the
universe in the future.

We discuss a scenario that the dark matter in late universe emerges as the
holographic stress energy tensor on the hypersurface in higher dimensional flat
bulk. Firstly we construct a toy model with a de Sitter hypersurface as the
holographic screen in flat spacetime. After adding the baryonic matter on the
screen, both of the dark matter and dark energy can be described by the
BrownYork stress energy tensor. From the Hamiltonian constraint equation in
higher dimensional spacetime, we find an interesting relation between the dark
matter and baryonic matter density parameters, by using the Lambda cold dark
matter parameterization. We further combine this holographic embedding of
emergent dark matter with brane world scenario and present a new
parameterization for the Friedmann equation, which can be reduced to our toy
constraint in the current universe. We also comment on the connection with the
Verlinde's emergent gravity, where the dark matter is regarded as the elastic
response of the baryonic matter on the de Sitter spacetime background. We show
that from the holographic de Sitter model with elasticity, the TullyFisher
relation of the dark matter distribution in the galaxy scale can be derived.

We present a fast numerical screened halo model algorithm (CHAM) for modeling
nonlinear power spectrum for the alternative models to LCDM. This method has
three obvious advantages. First of all, it is not being restricted to a
specific dark energy/modified gravity model. In principle, all of the screened
scalartensor theories can be applied. Second, the least assumptions are made
in the calculation. Hence, the physical picture is very easily understandable.
Third, it is very predictable and does not rely on the calibration from Nbody
simulation. As an example, we show the case of HuSawicki f(R) gravity. In this
case, the typical CPU time with the current parallel Python script (8 threads)
is roughly within $10$ minutes. The resulting spectra are in a good agreement
with Nbody data within a few percentage accuracy up to k~1 h/Mpc.

We present a minimal bottomup extension of the ChernSimons bulk action for
holographic translational symmetry breaking that naturally gives rise to pair
density waves. We construct stationary inhomogeneous black hole solutions in
which both the U(1) symmetry and spatially translational symmetry are
spontaneously broken at finite temperature and charge density. This novel
solution provides a dual description of a superconducting phase intertwined
with charge, current and parity orders.

We investigate the impact of different redshift distributions of random
samples on the baryon acoustic oscillations (BAO) measurements of
$D_V(z)r_\mathrm{d}^\mathrm{fid}/r_\mathrm{d}$ from the twopoint correlation
functions (2PCF) of galaxies in the Data Release 12 (DR12) of the Baryon
Oscillation Spectroscopic Survey (BOSS). Big surveys, such as BOSS, usually
assign redshifts to the random samples by randomly drawing values from the
measured redshift distributions of the data, which would necessarily introduce
fiducial signals of fluctuations into the random samples, weakening the signals
of BAO, if the cosmic variance cannot be ignored. We propose a smooth function
of redshift distribution that fits the data well to populate the random galaxy
samples. The resulting cosmological parameters match the input parameters of
the mock catalogue very well. The significance of BAO signals has been improved
by $0.33\sigma$ for a lowredshift (LOWZ) sample and by $0.03\sigma$ for a
constantstellarmass (CMASS) sample, though the absolute values do not change
significantly. Given the precision of the measurements of current cosmological
parameters, it would be appreciated for the future improvements on the
measurements of galaxy clustering.

Recent studies of the fluctuations of an open string in AdS space show some
pieces of evidence that the string with a worldsheet horizon could be a dual
description of SYK model, as they saturate universal chaos bound and share the
same symmetry. An open string hangs from the AdS boundary to the horizon of
black brane could be dual to a 0+1 dimensional boundary state. To be specific,
we find that the fluctuation of the string in charged BTZ black hole has an
asymptotic scaling symmetry, and its Euclidean IR fixed point is governed by
the quadratic order of Schwarzian action, which is just the low energy
effective theory of the SYK model. Considering the open string worldsheet also
has natural reparametrization symmetry, we conjecture that the action of the
string worldsheet is a dual description of SYK state.

The gravitational waves from compact binary systems are viewed as a standard
siren to probe the evolution of the universe. This paper summarizes the
potential and ability to use the gravitational waves to constrain the
cosmological parameters and the dark sector interaction in the Gaussian process
methodology. After briefly introducing the method to reconstruct the dark
sector interaction by the Gaussian process, the concept of standard sirens and
the analysis of reconstructing the dark sector interaction with LISA are
outlined. Furthermore, we estimate the constraint ability of the gravitational
waves on cosmological parameters with ET. The numerical methods we use are
Gaussian process and the MarkovChain MonteCarlo. Finally, we also forecast
the improvements of the abilities to constrain the cosmological parameters with
ET and LISA combined with the \it Planck.

We investigate in details the gravitational wave (GW) from the firstorder
phase transition (PT) in the extended standard model of particle physics with a
dimensionsix operator, which is capable of exhibiting the recently discovered
slow firstorder PT in addition to the usually studied fast firstorder PT. To
simplify the discussion, it is sufficient to work with an example of a toy
model with the sextic term, and we propose an unified description for both slow
and fast firstorder PTs. We next study the full oneloop effective potential
of the model with fixed/running renormalizationgroup (RG) scales. Compared to
the prediction of GW energy density spectrum from the fixed RG scale, we find
that the presence of running RG scale could amplify the peak amplitude by
amount of one order of magnitude while shift the peak frequency to the lower
frequency regime, and the promising regime of detection within the sensitivity
ranges of various spacebased GW detectors shrinks down to a lower cutoff
value of the sextic term rather than the previous expectation.

We investigate constraints on Lorentz invariance violation in the neutrino
sector from a joint analysis of big bang nucleosynthesis and the cosmic
microwave background. The effect of Lorentz invariance violation during the
epoch of big bang nucleosynthesis changes the predicted helium4 abundance,
which influences the power spectrum of the cosmic microwave background at the
recombination epoch. In combination with the latest measurement of the
primordial helium4 abundance, the Planck 2015 data of the cosmic microwave
background anisotropies give a strong constraint on the deformation parameter
since adding the primordial helium measurement breaks the degeneracy between
the deformation parameter and the physical dark matter density.

We study the effect of inhomogeneity, which is induced by the graviton mass
in massive gravity, on the mutual information and the chaotic behavior of a
2+1dimensional field theory from the gauge/gravity duality. When the system is
nearhomogeneous, the mutual information increases as the graviton mass grows.
However, when the system is far from homogeneity, the mutual information
decreases as the graviton mass increases. By adding the perturbations of energy
into the system, we investigate the dynamical mutual information in the shock
wave geometry. We find that the greater perturbations disrupt the mutual
information more rapidly, which resembles the butterfly effect in chaos theory.
Besides, the greater inhomogeneity reduces the dynamical mutual information
more quickly just as in the static case.

According to the conjecture "complexity equals action," the complexity of a
holographic state is equal to the action of a WheelerDeWitt (WDW) patch of
black holes in antide Sitter space. In this paper we calculate the action
growth of charged black holes with a single horizon, paying attention to the
contribution from a spacelike singularity inside the horizon. We consider two
kinds of such charged black holes: one is a charged dilaton black hole, and the
other is a BornInfeld black hole with $\beta^2 Q^2<1/4$. In both cases,
although an electric charge appears in the black hole solutions, the inner
horizon is absent, instead a spacelike singularity appears inside the horizon.
We find that the action growth of the WDW patch of the charged black hole is
finite and satisfies the Lloyd bound. As a check, we also calculate the action
growth of a charged black hole with a phantom Maxwell field. In this case,
although the contributions from the bulk integral and the spacelike singularity
are individually divergent, these two divergences just cancel each other and a
finite action growth is obtained. But in this case, the Lloyd bound is violated
as expected.

Environment, such as the accretion disk, could modify the signal of the
gravitational wave from the astrophysical black hole binaries. In this article,
we model the matter field around the intermediatemass binary black holes by
means of an axionlike scalar field and investigate their joint evolution. In
details, we consider the equal mass binary black holes surrounded by a shell of
axionlike scalar field both in spherical symmetric and nonspherical symmetric
cases, and with different strength of the scalar field. Our result shows that
the environmental scalar field could essentially modify the dynamics. Firstly,
in the spherical symmetric case, with increasing of the scalar field strength,
the number of circular orbit of the binary black hole is reduced. It means that
the scalar field could significantly accelerate the merger process. Secondly,
once the scalar field strength exceeds certain critical value, the scalar field
could collapse into a third black hole with its mass being larger than the
binary. Consequently, the new black hole collapsed from the environmental
scalar field could accrete the binary promptly and the binary collides headon
between each other. In this process, there is almost no any quadrupole signal
produced, namely the gravitational wave is greatly suppressed. Thirdly, when
the scalar field strength is relatively smaller than the critical value, the
black hole orbit could develop eccentricity through the accretion of the scalar
field. Fourthly, during the initial stage of the inspire, the gravitational
attractive force from the axionlike scalar field could induce a sudden turn in
the binary orbits, hence result in a transient wiggle in the gravitational
waveform. Finally, in the nonspherical case, the scalar field could
gravitationally attract the binary moving toward the mass center of the scalar
field and slow down the merger process.

The direct detection of gravitational wave by Laser Interferometer
GravitationalWave Observatory indicates the coming of the era of
gravitationalwave astronomy and gravitationalwave cosmology. It is expected
that more and more gravitationalwave events will be detected by currently
existing and planned gravitationalwave detectors. The gravitational waves open
a new window to explore the Universe and various mysteries will be disclosed
through the gravitationalwave detection, combined with other cosmological
probes. The gravitationalwave physics is not only related to gravitation
theory, but also is closely tied to fundamental physics, cosmology and
astrophysics. In this review article, three kinds of sources of gravitational
waves and relevant physics will be discussed, namely gravitational waves
produced during the inflation and preheating phases of the Universe, the
gravitational waves produced during the firstorder phase transition as the
Universe cools down and the gravitational waves from the three phases:
inspiral, merger and ringdown of a compact binary system, respectively. We will
also discuss the gravitational waves as a standard siren to explore the
evolution of the Universe.

We perform a forecast analysis of the ability of the LISA spacebased
interferometer to reconstruct the dark sector interaction using gravitational
wave (GW) standard sirens at high redshift. We employ Gaussian process methods
to reconstruct the distanceredshift relation in a model independent way. We
adopt simulated catalogues of standard sirens given by merging massive black
hole binaries (MBHBs) visible by LISA, with an electromagnetic (EM) counterpart
detectable by future telescopes. The catalogues are constructed considering
three different astrophysical scenarios for the evolution of MBHB mergers based
on the semianalytic model of E. Barausse (2012). We first use these standard
siren datasets to assess the potential of LISA in reconstructing a possible
interaction between vacuum dark energy and dark matter. Then we combine the
LISA cosmological data with supernovae data simulated for the Dark Energy
Survey (DES). We consider two scenarios distinguished by the time duration of
the LISA mission: 5 and 10 years. Using only LISA standard siren data, the dark
sector interaction can be well reconstructed from redshift $z\sim1$ to $z\sim3$
(5 yr) and $z\sim1$ up to $z\sim5$ (10 yr), though the reconstruction is
inefficient at lower redshift. When combined with the DES datasets, the
interaction is well reconstructed in the whole redshift region from $z\sim0$ to
$z\sim3$ (5 yr) and $z\sim0$ to $z\sim5$ (10 yr). MBHB standard sirens can thus
be used to constrain the dark sector interaction at redshift ranges not
reachable by usual supernovae datasets which probe only the $z\lesssim 1.5$
range. GW standard sirens will not only constitute a complementary and
alternative way, with respect to familiar EM observations, to probe the cosmic
expansion, but will also provide new tests to constrain possible deviations
from the standard $\Lambda$CDM dynamics, especially at high redshift.

We investigate the constraint ability of the gravitational wave (GW) as the
standard siren on the cosmological parameters by using the thirdgeneration
gravitational wave detector: the Einstein Telescope. We simulate the luminosity
distances and redshift measurements from 100 to 1000 GW events. We use two
different algorithms to constrain the cosmological parameters. For the Hubble
constant $H_0$ and dark matter density parameter $\Omega_m$, we adopt the
Markov chain Monte Carlo approach. We find that with about 500600 GW events we
can constrain the Hubble constant with an accuracy comparable to
\textit{Planck} temperature data and \textit{Planck} lensing combined results,
while for the dark matter density, GWs alone seem not able to provide the
constraints as good as for the Hubble constant; the sensitivity of 1000 GW
events is a little lower than that of \textit{Planck} data. It should require
more than 1000 events to match the \textit{Planck} sensitivity. Yet, for
analyzing the more complex dynamical property of dark energy, i.e., the
equation of state $w$, we adopt a new powerful nonparametric method: the
Gaussian process. We can reconstruct $w$ directly from the observational
luminosity distance at every redshift. In the low redshift region, we find that
about 700 GW events can give the constraints of $w(z)$ comparable to the
constraints of a constant $w$ by \textit{Planck} data with type Ia supernovae.
Those results show that GWs as the standard sirens to probe the cosmological
parameters can provide an independent and complementary alternative to current
experiments.

In this note we first review the degenerate vacua arising from the BMS
symmetries. According to the discussion in [1] one can define BMSanalogous
supertranslation and superrotation for spacetime with black hole in Gaussian
null coordinates. In the leading and subleading orders of near horizon
approximation, the infinitely degenerate black hole solutions are derived by
considering Einstein equations with or without cosmological constant, and they
are related to each other by the diffeomorphism generated by horizon
supertranslation. Higher order results and degenerate Rindler horizon solutions
also are given in appendices.

In this article, we give a definition of apparent horizon in a two
dimensional general dilaton gravity theory. With this definition, we construct
the mechanics of the apparent horizon by introducing a quasilocal energy of
the theory. Our discussion generalizes the apparent horizons mechanics in
general spherically symmetric spactimes in four or higher dimensions to the two
dimensional dilaton gravity case.

Recently a ComplexityAction (CA) duality conjecture has been proposed, which
relates the quantum complexity of a holographic boundary state to the action of
a WheelerDeWitt (WDW) patch in the antide Sitter (AdS) bulk. In this paper we
further investigate the duality conjecture for stationary AdS black holes and
derive some exact results for the growth rate of action within the
WheelerDeWitt (WDW) patch at late time approximation, which is supposed to be
dual to the growth rate of quantum complexity of holographic state. Based on
the results from the general $D$dimensional ReissnerNordstr\"{o}m (RN)AdS
black hole, rotating/charged Ba\~{n}adosTeitelboimZanelli (BTZ) black hole,
KerrAdS black hole and charged GaussBonnetAdS black hole, we present a
universal formula for the action growth expressed in terms of some
thermodynamical quantities associated with the outer and inner horizons of the
AdS black holes. And we leave the conjecture unchanged that the stationary AdS
black hole in Einstein gravity is the fastest computer in nature.

The gravitational collapse of a massless scalar field enclosed with a
perfectly reflecting wall in a spacetime with a cosmological constant $\Lambda$
is investigated. The mass scaling for the gapped collapse $ M_{AH}M_g \propto
(\epsilon_c\epsilon)^\xi$ is confirmed and a new time scaling for the gapped
collapse $T_{AH}T_g\propto(\epsilon_c\epsilon)^\zeta$ is found. We find that
both of these two critical exponents depend on the combination $\Lambda R^2$,
where $R$ is the radial position of the reflecting wall. Especially, we find an
evolution of the critical exponent $\xi$ from $0.37$ in the confined asymptotic
dS case with $\Lambda R^2=1.5$ to $0.7$ in asymptotic AdS case ($\Lambda
R^2\rightarrow\infty$), while the critical exponent $\zeta$ varies from $0.10$
to $0.26$, which shows the new critical behavior for the gapped collapse is
essentially different from the one in the Choptuik's case.