
We present our current best estimate of the plausible observing scenarios for
the Advanced LIGO, Advanced Virgo and KAGRA gravitationalwave detectors over
the next several years, with the intention of providing information to
facilitate planning for multimessenger astronomy with gravitational waves. We
estimate the sensitivity of the network to transient gravitationalwave signals
for the third (O3), fourth (O4) and fifth observing (O5) runs, including the
planned upgrades of the Advanced LIGO and Advanced Virgo detectors. We study
the capability of the network to determine the sky location of the source for
gravitationalwave signals from the inspiral of binary systems of compact
objects, that is BNS, NSBH, and BBH systems. The ability to localize the
sources is given as a skyarea probability, luminosity distance, and comoving
volume. The median sky localization area (90\% credible region) is expected to
be a few hundreds of square degrees for all types of binary systems during O3
with the Advanced LIGO and Virgo (HLV) network. The median sky localization
area will improve to a few tens of square degrees during O4 with the Advanced
LIGO, Virgo, and KAGRA (HLVK) network. We evaluate sensitivity and localization
expectations for unmodeled signal searches, including the search for
intermediate mass black hole binary mergers.

Accounting for stellar activity is a crucial component of the search for
eversmaller planets orbiting stars of all spectral types. We use Doppler
imaging methods to demonstrate that starspot induced radial velocity
variability can be effectively reduced for moderately rotating, fully
convective stars. Using starspot distributions extrapolated from sunspot
observations, we adopt typical M dwarf starspot distributions with low contrast
spots to synthesise line profile distortions. The distortions are recovered
using maximum entropy regularised fitting and the corresponding stellar radial
velocities are measured. The procedure is demonstrated for a lateM star
harbouring an orbiting planet in the habitable zone. The technique is effective
for stars with vsini = 110 km/s, reducing the stellar noise contribution by
factors of nearly an order of magnitude. With a carefully chosen observing
strategy, the technique can be used to determine the stellar rotation period
and is robust to uncertainties such as unknown stellar inclination. While
demonstrated for latetype M stars, the procedure is applicable to all spectral
types.

The terrestrial fossil record shows a significant variation in the extinction
and origination rates of species during the past half billion years. Numerous
studies have claimed an association between this variation and the motion of
the Sun around the Galaxy, invoking the modulation of cosmic rays, gamma rays
and comet impact frequency as a cause of this biodiversity variation. However,
some of these studies exhibit methodological problems, or were based on coarse
assumptions (such as a strict periodicity of the solar orbit). Here we
investigate this link in more detail, using a model of the Galaxy to
reconstruct the solar orbit and thus a predictive model of the temporal
variation of the extinction rate due to astronomical mechanisms. We compare
these predictions as well as those of various reference models with
paleontological data. Our approach involves Bayesian model comparison, which
takes into account the uncertainties in the paleontological data as well as the
distribution of solar orbits consistent with the uncertainties in the
astronomical data. We find that various versions of the orbital model are not
favored beyond simpler reference models. In particular, the distribution of
mass extinction events can be explained just as well by a uniform random
distribution as by any other model tested. Although our negative results on the
orbital model are robust to changes in the Galaxy model, the Sun's coordinates
and the errors in the data, we also find that it would be very difficult to
positively identify the orbital model even if it were the true one. (In
contrast, we do find evidence against simpler periodic models.) Thus while we
cannot rule out there being some connection between solar motion and
biodiversity variations on the Earth, we conclude that it is difficult to give
convincing positive conclusions of such a connection using current data.

Nearby gammaray bursts (GRBs) are likely to have represented a significant
threat to life on the Earth. Recent observations suggest that a significant
source of such bursts is compact binary mergers in globular clusters. This link
between globular clusters and GRBs offers the possibility to find time
intervals in the past with higher probabilities of a nearby burst, by tracing
globular cluster orbits back in time. Here we show that the expected flux from
such bursts is not flat over the past 550 Myr but rather exhibits three broad
peaks, at 70, 180 and 340 Myr ago. The main source for nearby GRBs for all
three time intervals is the globular cluster 47 Tuc, a consequence of its large
mass and high stellar encounter rate, as well as the fact that it is one of the
globular clusters which comes quite close to the Sun. Mass extinction events
indeed coincide with all three time intervals found in this study, although a
chance coincidence is quite likely. Nevertheless, the identified time intervals
can be used as a guide to search for specific signatures of GRBs in the
geological record around these times.

We show in detail that the $k_T$factorization for exclusive processes is
gaugedependent and inconsistent.

In the collinear factorization of the form factor for the transition
$\gamma^* \pi^0 \to \pi$ the hard part contains double log terms. These terms
will spoil the perturbative expansion of the hard part. A simple exponentiation
for resummation leads to divergent results. We study the resummation of these
double log's. We make an analysis to show the origin of the double log's. With
the understanding of the origin one can introduce soft factors and Nonstandard
LightCone Wave Functions(NLCWF) to derive a factorized form for the form
factor, where the hard part does not contain the double log's. There is a
perturbative relation between NLCWF and the standard LightCone Wave
Function(LCWF). Beside the renormalization scale $\mu$ the introduced NLCWF's
and soft factors have extra scales to characterize the double log's. Using the
evolutions of the extra scales and the relation we can do the resummation of
the double log's perturbatively in sense that LCWF's are the only
nonpertubative objects in the resumed formula. Our results with some models of
LCWF show that there is a significant difference between numerical predictions
with or without the resummation.

A new method is proposed to calculate wave functions in $k_T$factorization
in \cite{LiMi} as a comment about our paper \cite{FMW}. We point out that the
results obtained with the method are in conflict with the translation
invariance and depend on the chosen contours for loopintegrals. Therefore, the
method is in principle unacceptable and the results with the method cannot be
correct.

In the $k_T$factorization for exclusive processes, the nontrivial
$k_T$dependence of perturbative coefficients, or hard parts, is obtained by
taking offshell partons. This brings up the question of whether the
$k_T$factorization is gauge invariant. We study the $k_T$factorization for
the case $\pi \gamma^* \to \gamma$ at oneloop in a general covariant gauge.
Our results show that the hard part contains a lightcone singularity that is
absent in the Feynman gauge, which indicates that the $k_T$factorization is
{\it not} gauge invariant. These divergent contributions come from the
$k_T$dependent wave function of $\pi$ and are not related to a special
process. Because of this fact the $k_T$factorization for any process is not
gauge invariant and is violated. Our study also indicates that the
$k_T$factorization used widely for exclusive Bdecays is not gauge invariant
and is violated.

The k_Tfactorization has been widely used for exclusive decays of Bmesons.
In this factorization the pertubative coefficients are extracted from
scattering of offshell partons. Because the offshellness of partons the
extracted pertubative coefficients in general are not gaugeinvariant. We show
that these perturbative coefficients contain gaugedependent singularities
beyond treelevel. This leads to that the k_Tfactorization is gaugedependent
and violated in the general covariant gauge. This is not consistent with
expectations of physics. In the factorization there exist two different wave
functions for a Bmeson in general. We show that the two wave functions satisfy
different evolution equations whose solutions are used to resum large log's.
Based on this fact, the rseummation of large log's at leading log approximation
in the widely used k_Tfactorization with only one wave function is
inconsistent.

In the collinear factorization of the form factor for the transition
$\gamma^* \pi^0 \to \gamma$ the hard part contains double log terms as $\ln^2
x$ with $x$ as the momentum fraction of partons from 0 to 1. A simple
exponentiation for resummation leads to divergent results. We study the
resummation of these $\ln^2 x$ terms. We show that the $\ln^2 x$ terms come
partly from the lightcone wave function(LCWF) and partly from the form factor.
We introduce a jet factor to factorize the $\ln^2 x$ term in the form factor.
To handel the $\ln^2 x$ terms from the LCWF we introduce a nonstandard
lightcone wave function(NLCWF) with the gauge links off the lightcone
direction. An interesting relation between two wave function is found. With the
introduced NLCWF and the jet factor we can refactorize the form factor and
obtain a new hard part which does not contain terms with $\ln^2 x$. Beside the
renormalization scale $\mu$ the introduce NLCWF and jet factor have extra
scales to characterize their $x$behaviors. Using the evolutions of the extra
scales and the relation we can do the resummation perturbatively in sense that
the LCWF is the only nonpertubative object in the resumed formula. Our results
with some models of LCWF show that there is a significant difference between
numerical predictions with the resummation and that without the resummation,
and the resummed predictions can describe the experimental data.