The Weak Equivalence Principle (WEP), stating that two bodies of different compositions and/or mass fall at the same rate in a gravitational field (universality of free fall), is at the very foundation of General Relativity. The MICROSCOPE mission aims to test its validity to a precision of $10^{-15}$, two orders of magnitude better than current on-ground tests, by using two masses of different compositions (titanium and platinum alloys) on a quasi-circular trajectory around the Earth. This is realised by measuring the accelerations inferred from the forces required to maintain the two masses exactly in the same orbit. Any significant difference between the measured accelerations, occurring at a defined frequency, would correspond to the detection of a violation of the WEP, or to the discovery of a tiny new type of force added to gravity. MICROSCOPE's first results show no hint for such a difference, expressed in terms of E\"otv\"os parameter $\delta(Ti,Pt)=[-1\pm{}9{\rm (stat)}\pm{}9{\rm (syst)}] \times{}10^{-15}$ (both 1$\sigma$ uncertainties) for a titanium and platinum pair of materials. This result was obtained on a session with 120 orbital revolutions representing 7\% of the current available data acquired during the whole mission. The quadratic combination of 1$\sigma$ uncertainties leads to a current limit on $\delta$ of about $1.3\times{}10^{-14}$.

For a Schwarzchild black hole of mass $M$, we consider a test particle falling from rest at infinity and becoming trapped, at late time, on the unstable circular orbit of radius $r=4GM/c^2$. When the particle is endowed with a small mass, $\mu\ll M$, it experiences an effective gravitational self-force, whose conservative piece shifts the critical value of the angular momentum and the frequency of the asymptotic circular orbit away from their geodesic values. By directly integrating the self-force along the orbit (ignoring radiative dissipation), we numerically calculate these shifts to $O(\mu/M)$. Our numerical values are found to be in agreement with estimates first made within the Effective One Body formalism, and with predictions of the first law of black-hole-binary mechanics (as applied to the asymptotic circular orbit). Our calculation is based on a time-domain integration of the Lorenz-gauge perturbation equations, and it is a first such calculation for an unbound orbit. We tackle several technical difficulties specific to unbound orbits, illustrating how these may be handled in more general cases of unbound motion. Our method paves the way to calculations of the self-force along hyperbolic-type scattering orbits. Such orbits can probe the two-body potential down to the "light ring", and could thus supply strong-field calibration data for eccentricity-dependent terms in the Effective One Body model of merging binaries.

According to the Weak Equivalence Principle, all bodies should fall at the same rate in a gravitational field. The MICROSCOPE satellite, launched in April 2016, aims to test its validity at the $10^{-15}$ precision level, by measuring the force required to maintain two test masses (of titanium and platinum alloys) exactly in the same orbit. A non-vanishing result would correspond to a violation of the Equivalence Principle, or to the discovery of a new long-range force. Analysis of the first data gives $\delta\rm{(Ti,Pt)}= [-1 \pm 9 (\mathrm{stat}) \pm 9 (\mathrm{syst})] \times 10^{-15}$ (1$\sigma$ statistical uncertainty) for the titanium-platinum E\"otv\"os parameter characterizing the relative difference in their free-fall accelerations.

The data analysis of the gravitational wave signals emitted by coalescing neutron star binaries requires the availability of an accurate analytical representation of the dynamics and waveforms of these systems. We propose an effective-one-body (EOB) model that describes the general relativistic dynamics of neutron star binaries from the early inspiral up to merger. Our EOB model incorporates an enhanced attractive tidal potential motivated by recent analytical advances in the post-Newtonian and gravitational self-force description of relativistic tidal interactions. No fitting parameters are introduced for the description of tidal interaction in the late, strong-field dynamics. We compare the model energetics and the gravitational wave phasing with new high-resolution multi-orbit numerical relativity simulations of equal-mass configurations with different equations of state. We find agreement within the uncertainty of the numerical data for all configurations. Our model is the first semi-analytical model which captures the tidal amplification effects close to merger. It thereby provides the most accurate analytical representation of binary neutron star dynamics and waveforms currently available.

We complete the analytical determination, at the 4th post-Newtonian (4PN) approximation, of the conservative dynamics of gravitationally interacting two-point-mass systems. This completion is obtained by resolving the infra-red ambiguity which had blocked a previous 4PN calculation [P.Jaranowski and G.Sch\"afer, Phys. Rev. D 87, 081503(R) (2013)] by taking into account the 4PN breakdown of the usual near-zone expansion due to infinite-range tail-transported temporal correlations found long ago [L.Blanchet and T.Damour, Phys. Rev. D 37, 1410 (1988)]. This leads to a Poincar\'e-invariant 4PN-accurate effective action for two masses, which mixes instantaneous interaction terms (described by a usual Hamiltonian) with a (time-symmetric) nonlocal-in-time interaction.

We compute the conservative piece of the gravitational self-force (GSF) acting on a particle of mass m_1 as it moves along an (unstable) circular geodesic orbit between the innermost stable circular orbit (ISCO) and the light ring of a Schwarzschild black hole of mass m_2>> m_1. More precisely, we construct the function h_{uu}(x) = h_{\mu\nu} u^{\mu} u^{\nu} (related to Detweiler's gauge-invariant "redshift" variable), where h_{\mu\nu} is the regularized metric perturbation in the Lorenz gauge, u^{\mu} is the four-velocity of m_1, and x= [Gc^{-3}(m_1+m_2)\Omega]^{2/3} is an invariant coordinate constructed from the orbital frequency \Omega. In particular, we explore the behavior of h_{uu} just outside the "light ring" at x=1/3, where the circular orbit becomes null. Using the recently discovered link between h_{uu} and the piece a(u), linear in the symmetric mass ratio \nu, of the main radial potential A(u,\nu) of the Effective One Body (EOB) formalism, we compute a(u) over the entire domain 0<u<1/3. We find that a(u) diverges at the light-ring as ~0.25 (1-3u)^{-1/2}, explain the physical origin of this divergence, and discuss its consequences for the EOB formalism. We construct accurate global analytic fits for a(u), valid on the entire domain 0<u<1/3 (and possibly beyond), and give accurate numerical estimates of the values of a(u) and its first 3 derivatives at the ISCO, as well as the O(\nu) shift in the ISCO frequency. In previous work we used GSF data on slightly eccentric orbits to compute a certain linear combination of a(u) and its first two derivatives, involving also the O(\nu) piece \bar d(u) of a second EOB radial potential {\bar D}(u,\nu). Combining these results with our present global analytic representation of a(u), we numerically compute {\bar d}(u)$ on the interval 0<u\leq 1/6.

The gravitational wave signal from a binary neutron star inspiral contains information on the nuclear equation of state. This information is contained in a combination of the tidal polarizability parameters of the two neutron stars and is clearest in the late inspiral, just before merger. We use the recently defined tidal extension of the effective one-body formalism to construct a controlled analytical description of the frequency-domain phasing of neutron star inspirals up to merger. Exploiting this analytical description we find that the tidal polarizability parameters of neutron stars can be measured by the advanced LIGO-Virgo detector network from gravitational wave signals having a reasonable signal-to-noise ratio of $\rho=16$. This measurability result seems to hold for all the nuclear equations of state leading to a maximum mass larger than $1.97M_\odot$. We also propose a promising new way of extracting information on the nuclear equation of state from a coherent analysis of an ensemble of gravitational wave observations of separate binary merger events.

Using accurate numerical relativity simulations of (nonspinning) black-hole binaries with mass ratios 1:1, 2:1 and 3:1 we compute the gauge invariant relation between the (reduced) binding energy $E$ and the (reduced) angular momentum $j$ of the system. We show that the relation $E(j)$ is an accurate diagnostic of the dynamics of a black-hole binary in a highly relativistic regime. By comparing the numerical-relativity $E^{\rm NR} (j)$ curve with the predictions of several analytic approximation schemes, we find that, while the usual, non-resummed post-Newtonian-expanded $E^{\rm PN} (j)$ relation exhibits large and growing deviations from $E^{\rm NR} (j)$, the prediction of the effective one-body formalism, based purely on known analytical results (without any calibration to numerical relativity), agrees strikingly well with the numerical-relativity results.

If the fundamental constants of nature have a cosmic spatial variation, there will in general be extra forces with a preferred direction in space which violate the equivalence principle. We show that the millimeter-precision Apache Point Observatory Lunar Laser-ranging Operation provides a very sensitive probe of such variation that has the capability of detecting a cosmic gradient of the ratio between the quark masses and the strong interaction scale at the level with the gradient of ln (m_quark/Lambda_QCD) ~ 2.6 x 10^-6 Glyr^-1, which is comparable to the cosmic gradients suggested by the recently reported measurements of Webb et al. We also point out the capability of presently planned improved equivalence principle tests, at the Delta g/g < 10^-17 level, to probe similar cosmic gradients.

The aim of these lectures is to give an introduction to several topics which lie at the intersection of string theory, gravity theory and gravity phenomenology. One successively reviews: (i) the "membrane" approach to the dissipative dynamics of classical black holes, (ii) the current experimental tests of gravity, and their theoretical interpretation, (iii) some aspects of the string-inspired phenomenology of the gravitational sector, and (iv) some possibilities for observing string-related signals in cosmology (including a discussion of gravitational wave signals from cosmic superstrings).

We study to what extent wormholes can mimic the observational features of black holes. It is surprisingly found that many features that could be thought of as ``characteristic'' of a black hole (endowed with an event horizon) can be closely mimicked by a globally static wormhole, having no event horizon. This is the case for: the apparently irreversible accretion of matter down a hole, no-hair properties, quasi-normal-mode ringing, and even the dissipative properties of black hole horizons, such as a finite surface resistivity equal to 377 Ohms. The only way to distinguish the two geometries on an observationally reasonable time scale would be through the detection of Hawking's radiation, which is, however, too weak to be of practical relevance for astrophysical black holes. We point out the existence of an interesting spectrum of quantum microstates trapped in the throat of a wormhole which could be relevant for storing the information ``lost'' during a gravitational collapse.

We study the gravitational emission, in Superstring Theory, from fundamental strings exhibiting cusps. The classical computation of the gravitational radiation signal from cuspy strings features strong bursts in the special null directions associated to the cusps. We perform a quantum computation of the gravitational radiation signal from a cuspy string, as measured in a gravitational wave detector using matched filtering and located in the special null direction associated to the cusp. We study the quantum statistics (expectation value and variance) of the measured filtered signal and find that it is very sharply peaked around the classical prediction. Ultimately, this result follows from the fact that the detector is a low-pass filter which is blind to the violent high-frequency quantum fluctuations of both the string worldsheet, and the incoming gravitational field.

Using the Effective One Body approach, that includes nonperturbative resummed estimates for the damping and conservative parts of the compact binary dynamics, we compute the recoil during the late inspiral and the subsequent plunge of non-spinning black holes of comparable masses moving in quasi-circular orbits. Further, using a prescription that smoothly connects the plunge phase to a perturbed single black hole, we obtain an estimate for the total recoil associated with the binary black hole coalescence. We show that the crucial physical feature which determines the magnitude of the terminal recoil is the presence of a ``burst'' of linear momentum flux emitted slightly before coalescence. When using the most natural expression for the linear momentum flux during the plunge, together with a Taylor-expanded $(v/c)^4$ correction factor, we find that the maximum value of the terminal recoil is $\sim 74$ km/s and occurs for a mass ratio $m_2/m_1 \simeq 0.38$. We comment, however, on the fact that the above `best bet estimate' is subject to strong uncertainties because the location and amplitude of the crucial peak of linear momentum flux happens at a moment during the plunge where most of the simplifying analytical assumptions underlying the Effective One Body approach are no longer justified. Changing the analytical way of estimating the linear momentum flux, we find maximum recoils that range between 49 and 172 km/s. (Abridged)

The gravitational wave (GW) signals emitted by a network of cosmic strings are reexamined in view of the possible formation of a network of cosmic superstrings at the end of brane inflation. The reconnection probability $p$ of intersecting fundamental or Dirichlet strings might be much smaller than 1, and the properties of the resulting string network may differ significantly from those of ordinary strings (which have $p=1$). In addition, it has been recently suggested that the typical length of newly formed loops may differ by a factor $\epsilon \ll 1$ from its standard estimate. Here, we analyze the effects of the two parameters $p$ and $\epsilon$ on the GW signatures of strings. We consider both the GW bursts emitted from cusps of oscillating string loops, which have been suggested as candidate sources for the LIGO/VIRGO and LISA interferometers, and the stochastic GW background, which may be detectable by pulsar timing observations. In both cases we find that previously obtained results are \textit{quite robust}, at least when the loop sizes are not suppressed by many orders of magnitude relative to the standard scenario. We urge pulsar observers to reanalyze a recently obtained 17-year combined data set to see whether the large scatter exhibited by a fraction of the data might be due to a transient GW burst activity of some sort, e.g. to a near cusp event.

This review is made of two parts which are related to Poincar\'e in different ways. The first part reviews the work of Poincar\'e on the Theory of (Special) Relativity. One emphasizes both the remarkable achievements of Poincar\'e, and the fact that he never came close to what is the essential conceptual achievement of Einstein: changing the concept of time. The second part reviews a topic which probably would have appealed to Poincar\'e because it involves several mathematical structures he worked on: chaotic dynamics, discrete reflection groups, and Lobachevskii space. This topic is the hidden role of Kac-Moody algebras in the billiard description of the asymptotic behaviour of certain Einstein-matter systems near a cosmological singularity. Of particular interest are the Einstein-matter systems arising in the low-energy limit of superstring theory. These systems seem to exhibit the highest-rank hyperbolic Kac-Moody algebras, and notably E(10), as hidden symmetries.

The structure of the general, inhomogeneous solution of (bosonic) Einstein-matter systems in the vicinity of a cosmological singularity is considered. We review the proof (based on ideas of Belinskii-Khalatnikov-Lifshitz and technically simplified by the use of the Arnowitt-Deser-Misner Hamiltonian formalism) that the asymptotic behaviour, as one approaches the singularity, of the general solution is describable, at each (generic) spatial point, as a billiard motion in an auxiliary Lorentzian space. For certain Einstein-matter systems, notably for pure Einstein gravity in any spacetime dimension D and for the particular Einstein-matter systems arising in String theory, the billiard tables describing asymptotic cosmological behaviour are found to be identical to the Weyl chambers of some Lorentzian Kac-Moody algebras. In the case of the bosonic sector of supergravity in 11 dimensional spacetime the underlying Lorentzian algebra is that of the hyperbolic Kac-Moody group E(10), and there exists some evidence of a correspondence between the general solution of the Einstein-three-form system and a null geodesic in the infinite dimensional coset space E(10)/K (E(10)), where K (E(10)) is the maximal compact subgroup of E(10).

The structure of the general, inhomogeneous solution of (bosonic) Einstein-matter systems in the vicinity of a cosmological singularity is considered. We review the proof (based on ideas of Belinskii-Khalatnikov-Lifshitz and technically simplified by the use of the Arnowitt-Deser-Misner Hamiltonian formalism) that the asymptotic behaviour, as one approaches the singularity, of the general solution is describable, at each (generic) spatial point, as a billiard motion in an auxiliary Lorentzian space. For certain Einstein-matter systems, notably for pure Einstein gravity in any spacetime dimension D and for the particular Einstein-matter systems arising in String theory, the billiard tables describing asymptotic cosmological behaviour are found to be identical to the Weyl chambers of some Lorentzian Kac-Moody algebras. In the case of the bosonic sector of supergravity in 11 dimensional spacetime the underlying Lorentzian algebra is that of the hyperbolic Kac-Moody group E_{10}, and there exists some evidence of a correspondence between the general solution of the Einstein-three-form system and a null geodesic in the infinite dimensional coset space E_{10} / K (E_{10}), where K (E_{10}) is the maximal compact subgroup of E_{10}.

We provide a method for analytically constructing high-accuracy templates for the gravitational wave signals emitted by compact binaries moving in inspiralling eccentric orbits. By contrast to the simpler problem of modeling the gravitational wave signals emitted by inspiralling {\it circular} orbits, which contain only two different time scales, namely those associated with the orbital motion and the radiation reaction, the case of {\it inspiralling eccentric} orbits involves {\it three different time scales}: orbital period, periastron precession and radiation-reaction time scales. By using an improved `method of variation of constants', we show how to combine these three time scales, without making the usual approximation of treating the radiative time scale as an adiabatic process. We explicitly implement our method at the 2.5PN post-Newtonian accuracy. Our final results can be viewed as computing new `post-adiabatic' short period contributions to the orbital phasing, or equivalently, new short-period contributions to the gravitational wave polarizations, $h_{+,\times}$, that should be explicitly added to the `post-Newtonian' expansion for $h_{+,\times}$, if one treats radiative effects on the orbital phasing of the latter in the usual adiabatic approximation. Our results should be of importance both for the LIGO/VIRGO/GEO network of ground based interferometric gravitational wave detectors (especially if Kozai oscillations turn out to be significant in globular cluster triplets), and for the future space-based interferometer LISA.

After recalling the definition of black holes, and reviewing their energetics and their classical thermodynamics, one expounds the conjecture of Bekenstein, attributing an entropy to black holes, and the calculation by Hawking of the semi-classical radiation spectrum of a black hole, involving a thermal (Planckian) factor. One then discusses the attempts to interpret the black-hole entropy as the logarithm of the number of quantum micro-states of a macroscopic black hole, with particular emphasis on results obtained within string theory. After mentioning the (technically cleaner, but conceptually more intricate) case of supersymmetric (BPS) black holes and the corresponding counting of the degeneracy of Dirichlet-brane systems, one discusses in some detail the ``correspondence'' between massive string states and non-supersymmetric Schwarzschild black holes.

We review some string-inspired theoretical models which incorporate a correlated spacetime variation of coupling constants while remaining naturally compatible both with phenomenological constraints coming from geochemical data (Oklo; Rhenium decay) and with present equivalence principle tests. Barring unnatural fine-tunings of parameters, a variation of the fine-structure constant as large as that recently ``observed'' by Webb et al. in quasar absorption spectra appears to be incompatible with these phenomenological constraints. Independently of any model, it is emphasized that the best experimental probe of varying constants are high-precision tests of the universality of free fall, such as MICROSCOPE and STEP. Recent claims by Bekenstein that fine-structure-constant variability does not imply detectable violations of the equivalence principle are shown to be untenable.

This paper discusses the IAU Resolutions B1.3, B1.4, B1.5 and B1.9 (2000) that were adopted during the 24th General Assembly in Manchester, 2000 and provides details and explanations for these Resolutions. It is explained why they present significant progress over the corresponding IAU 1991 Resolutions and why they are necessary in the light of present accuracies in astrometry, celestial mechanics and metrology. In fact most of these Resolutions are consistent with astronomical models and software already in use. The metric tensors and gravitational potentials of both the Barycentric Celestial Reference System and Geocentric Celestial Reference System are defined and discussed. The necessity and relevance of the two celestial reference systems are explained. The transformations of coordinates and gravitational potentials are discussed. Potential coefficients parameterizing the post-Newtonian gravitational potentials are expounded. Simplified versions of the time transformations suitable for modern clock accuracies are elucidated. Various approximations used in the Resolutions are explicated and justified. Some models (e.g. for higher spin moments) that serve the purpose for estimating orders of magnitude have actually never been published before.

We explore spherically symmetric stationary solutions, generated by ``stars'' with regular interiors, in purely massive gravity. We reexamine the claim that the resummation of non-linear effects can cure, in a domain near the source, the discontinuity exhibited by the linearized theory as the mass m of the graviton tends to zero. First, we find analytical difficulties with this claim, which appears not to be robust under slight changes in the form of the mass term. Second, by numerically exploring the inward continuation of the class of asymptotically flat solutions, we find that, when m is ``small'', they all end up in a singularity at a finite radius, well outside the source, instead of joining some conjectured ``continuous'' solution near the source. We reopen, however, the possibility of reconciling massive gravity with phenomenology by exhibiting a special class of solutions, with ``spontaneous symmetry breaking'' features, which are close, near the source, to general relativistic solutions and asymptote, for large radii, a de Sitter solution of curvature ~m^2.

It is shown in detail that the dynamics of the Einstein-dilaton-p-form system in the vicinity of a spacelike singularity can be asymptotically described, at a generic spatial point, as a billiard motion in a region of Lobachevskii space (realized as an hyperboloid in the space of logarithmic scale factors). This is done within the Hamiltonian formalism, and for an arbitrary number of spacetime dimensions $D \geq 4$. A key role in the derivation is played by the Iwasawa decomposition of the spatial metric, and by the fact that the off-diagonal degrees of freedom, as well as the p-form degrees of freedom, get ``asymptotically frozen'' in this description. For those models admitting a Kac-Moody theoretic interpretation of the billiard dynamics we outline how to set up an asymptotically equivalent description in terms of a one-dimensional non-linear sigma-model formally invariant under the corresponding Kac-Moody group.

We explore the cosmological solutions of classes of non-linear bigravity theories. These theories are defined by effective four-dimensional Lagrangians describing the coupled dynamics of two metric tensors, and containing, in the linearized limit, both a massless graviton and an ultralight one. We focus on two paradigmatic cases: the case where the coupling between the two metrics is given by a Pauli-Fierz-type mass potential, and the case where this coupling derives from five-dimensional brane constructions. We find that cosmological evolutions in bigravity theories can be described in terms of the dynamics of two ``relativistic particles'', moving in a curved Lorenzian space, and connected by some type of nonlinear ``spring''. Classes of bigravity cosmological evolutions exhibit a ``locking'' mechanism under which the two metrics ultimately stabilize in a bi-de-Sitter configuration, with relative (constant) expansion rates. In the absence of matter, we find that a generic feature of bigravity cosmologies is to exhibit a period of cosmic acceleration. This leads us to propose bigravity as a source of a new type of dark energy (``tensor quintessence''), exhibiting specific anisotropic features. Bigravity could also have been the source of primordial inflation.

We discuss the fully non-linear formulation of multigravity. The concept of universality classes of effective Lagrangians describing bigravity, which is the simplest form of multigravity, is introduced. We show that non-linear multigravity theories can naturally arise in several different physical contexts: brane configurations, certain Kaluza-Klein reductions and some non-commutative geometry models. The formal and phenomenological aspects of multigravity (including the problems linked to the linearized theory of massive gravitons) are briefly discussed.

We explore a version of the cosmological dilaton-fixing and decoupling mechanism in which the dilaton-dependence of the low-energy effective action is extremized for infinitely large values of the bare string coupling $g_s^2 = e^{\phi}$. We study the efficiency with which the dilaton $\phi$ runs away towards its ``fixed point'' at infinity during a primordial inflationary stage, and thereby approximately decouples from matter. The residual dilaton couplings are found to be related to the amplitude of the density fluctuations generated during inflation. For the simplest inflationary potential, $V (\chi) = {1/2} m_{\chi}^2 (\phi) \chi^2$, the residual dilaton couplings are shown to predict violations of the universality of gravitational acceleration near the $\Delta a / a \sim 10^{-12}$ level. This suggests that a modest improvement in the precision of equivalence principle tests might be able to detect the effect of such a runaway dilaton. Under some assumptions about the coupling of the dilaton to dark matter and/or dark energy, the expected time-variation of natural ``constants'' (in particular of the fine-structure constant) might also be large enough to be within reach of improved experimental or observational data.

In a recently proposed scenario, where the dilaton decouples while cosmologically attracted towards infinite bare string coupling, its residual interactions can be related to the amplitude of density fluctuations generated during inflation, and are large enough to be detectable through a modest improvement on present tests of free-fall universality. Provided it has significant couplings to either dark matter or dark energy, a runaway dilaton can also induce time-variations of the natural "constants" within the reach of near-future experiments.

We compare recent numerical results, obtained within a ``helical Killing vector'' (HKV) approach, on circular orbits of corotating binary black holes to the analytical predictions made by the effective one body (EOB) method (which has been recently extended to the case of spinning bodies). On the scale of the differences between the results obtained by different numerical methods, we find good agreement between numerical data and analytical predictions for several invariant functions describing the dynamical properties of circular orbits. This agreement is robust against the post-Newtonian accuracy used for the analytical estimates, as well as under choices of resummation method for the EOB ``effective potential'', and gets better as one uses a higher post-Newtonian accuracy. These findings open the way to a significant ``merging'' of analytical and numerical methods, i.e. to matching an EOB-based analytical description of the (early and late) inspiral, up to the beginning of the plunge, to a numerical description of the plunge and merger. We illustrate also the ``flexibility'' of the EOB approach, i.e. the possibility of determining some ``best fit'' values for the analytical parameters by comparison with numerical data.

We briefly review three aspects of string cosmology: (1) the ``stochastic'' approach to the pre-big bang scenario, (2) the presence of chaos in the generic cosmological solutions of the tree-level low-energy effective actions coming out of string theory, and (3) the remarkable link between the latter chaos and the Weyl groups of some hyperbolic Kac-Moody algebras. Talk given at the Francqui Colloquium ``Strings and Gravity: Tying the Forces Together'' (Brussels, October 2001).

It has been recently pointed out that any primary galactic population of Weakly Interacting Massive Particles (WIMP) generates, through collisions with solar matter, a secondary population of ``slow'' WIMPs trapped in the inner solar system. We show that taking into account this ``slow'' solar-system population dramatically enhances the possibility to probe the existence of stable massive neutrinos (of a 4th generation) in underground neutrino experiments. Our work suggests that a reanalysis of existing underground neutrino data should be able to bring extremely tight constraints on the possible existence of a stable massive 4th neutrino.

It is shown that the never ending oscillatory behaviour of the generic solution, near a cosmological singularity, of the massless bosonic sector of superstring theory can be described as a billiard motion within a simplex in 9-dimensional hyperbolic space. The Coxeter group of reflections of this billiard is discrete and is the Weyl group of the hyperbolic Kac-Moody algebra E$_{10}$ (for type II) or BE$_{10}$ (for type I or heterotic), which are both arithmetic. These results lead to a proof of the chaotic (``Anosov'') nature of the classical cosmological oscillations, and suggest a ``chaotic quantum billiard'' scenario of vacuum selection in string theory.

The strong beams of high-frequency gravitational waves (GW) emitted by cusps and kinks of cosmic strings are studied in detail. As a consequence of these beams, the stochastic ensemble of GW's generated by a cosmological network of oscillating loops is strongly non Gaussian, and includes occasional sharp bursts that stand above the ``confusion'' GW noise made of many smaller overlapping bursts. Even if only 10% of all string loops have cusps these bursts might be detectable by the planned GW detectors LIGO/VIRGO and LISA for string tensions as small as $G \mu \sim 10^{-13}$. In the implausible case where the average cusp number per loop oscillation is extremely small, the smaller bursts emitted by the ubiquitous kinks will be detectable by LISA for string tensions as small as $G \mu \sim 10^{-12}$. We show that the strongly non Gaussian nature of the stochastic GW's generated by strings modifies the usual derivation of constraints on $G \mu$ from pulsar timing experiments. In particular the usually considered ``rms GW background'' is, when $G \mu \gaq 10^{-7}$, an overestimate of the more relevant confusion GW noise because it includes rare, intense bursts. The consideration of the confusion GW noise suggests that a Grand Unified Theory (GUT) value $ G \mu \sim 10^{-6}$ is compatible with existing pulsar data, and that a modest improvement in pulsar timing accuracy could detect the confusion noise coming from a network of cuspy string loops down to $ G \mu \sim 10^{-11}$. The GW bursts discussed here might be accompanied by Gamma Ray Bursts.

Some time ago, it was found that the never-ending oscillatory chaotic behaviour discovered by Belinsky, Khalatnikov and Lifshitz (BKL) for the generic solution of the vacuum Einstein equations in the vicinity of a spacelike ("cosmological") singularity disappears in spacetime dimensions $D= d+1>10$. Recently, a study of the generalization of the BKL chaotic behaviour to the superstring effective Lagrangians has revealed that this chaos is rooted in the structure of the fundamental Weyl chamber of some underlying hyperbolic Kac-Moody algebra. In this letter, we show that the same connection applies to pure gravity in any spacetime dimension $\geq 4$, where the relevant algebras are $AE_d$. In this way the disappearance of chaos in pure gravity models in $D > 10$ dimensions becomes linked to the fact that the Kac-Moody algebras $AE_d$ are no longer hyperbolic for $d > 9$.

A new Solar System population of Weakly Interacting Massive Particle (WIMP) dark matter has been proposed to exist. We investigate the implications of this population on indirect signals in neutrino telescopes (due to WIMP annihilations in the Earth) for the case when the WIMP is the lightest neutralino of the MSSM, the minimal supersymmetric extension of the standard model. The velocity distribution and capture rate of this new population is evaluated and the flux of neutrino-induced muons from the center of the Earth in neutrino telescopes is calculated. We show that the effects of the new population can be crucial for masses around 60-120 GeV, where enhancements of the predicted muon flux from the center of the Earth by up to a factor of 100 compared to previously published estimates occur. As a result of the new WIMP population, neutrino telescopes should be able to probe a much larger region of parameter space in this mass range.

Some spatially homogeneous Bianchi type I cosmological models filled with homogeneous ``electric'' $p$-form fields are shown to mimic the never-ending oscillatory behaviour of generic string cosmologies established recently. The validity of the ``Kasner-free-flights plus collisions-on-potential-walls'' picture is also illustrated in the case of known, non-chaotic, superstring solutions.

Cusps of cosmic strings emit strong beams of high-frequency gravitational waves (GW). As a consequence of these beams, the stochastic ensemble of gravitational waves generated by a cosmological network of oscillating loops is strongly non Gaussian, and includes occasional sharp bursts that stand above the rms GW background. These bursts might be detectable by the planned GW detectors LIGO/VIRGO and LISA for string tensions as small as $G \mu \sim 10^{-13}$. The GW bursts discussed here might be accompanied by Gamma Ray Bursts.

It is shown that the general solution near a spacelike singularity of the Einstein-dilaton-p-form field equations relevant to superstring theories and M-theory exhibits an oscillatory behaviour of the Belinskii-Khalatnikov-Lifshitz type. String dualities play a significant role in the analysis.

Recently, a new Solar System population of weakly interacting massive particle (WIMP) dark matter has been proposed to exist. We investigate the implications of this population on indirect signals in neutrino telescopes (due to WIMP annihilations in the Earth) for the case when the WIMP is the lightest neutralino of the MSSM, the minimal supersymmetric extension of the standard model. The velocity distribution and capture rate of this new population is evaluated and the flux of neutrino-induced muons from the center of the Earth in neutrino telescopes is calculated. The strength of the signal is very sensitive to the velocity distribution of the new population. We analytically estimate this distribution using the approximate conservation of the component of the WIMP angular momentum orthogonal to the ecliptic plane. The non-linear problem of combining a fixed capture rate from the standard galactic WIMP population with one rising linearly with time from the new population to obtain the present-day annihilation rate in the Earth is also solved analytically. We show that the effects of the new population can be crucial for masses below around 150 GeV, where enhancements of the predicted muon flux from the center of the Earth by up to a factor of 100 compared to previously published estimates occur. As a result of the new WIMP population, the next generation of neutrino telescopes should be able to probe a much larger region of parameter space in the mass range 60-130 GeV.

We study the radiation reaction on cosmic strings due to the emission of dilatonic, gravitational and axionic waves. After verifying the (on average) conservative nature of the time-symmetric self-interactions, we concentrate on the finite radiation damping force associated with the half-retarded minus half-advanced ``reactive'' fields. We revisit a recent proposal of using a ``local back reaction approximation'' for the reactive fields. Using dimensional continuation as convenient technical tool, we find, contrary to previous claims, that this proposal leads to antidamping in the case of the axionic field, and to zero (integrated) damping in the case of the gravitational field. One gets normal positive damping only in the case of the dilatonic field. We propose to use a suitably modified version of the local dilatonic radiation reaction as a substitute for the exact (non-local) gravitational radiation reaction. The incorporation of such a local approximation to gravitational radiation reaction should allow one to complete, in a computationally non-intensive way, string network simulations and to give better estimates of the amount and spectrum of gravitational radiation emitted by a cosmologically evolving network of massive strings.

Big Bang Nucleosynthesis (BBN) is studied within the framework of a two-parameter family of tensor-scalar theories of gravitation, with nonlinear scalar-matter coupling function a(phi). We run a BBN code modified by tensor-scalar gravity, and impose that the theoretically predicted BBN yields of Deuterium, Helium and Lithium lie within some conservative observational ranges. It is found that large initial values of a(phi) (corresponding to initial cosmological expansion rates much larger than standard) are compatible with observed BBN yields. However, the BBN-inferred upper bound on the cosmological baryon density is insignificantly modified by considering tensor-scalar gravity. Taking into account the effect of e^+ e^- annihilation together with the subsequent effect of the matter-dominated era (which both tend to decouple phi from matter), we find that the present value of the scalar coupling, i.e. the present level of deviation from Einstein's theory, must be, for compatibility with BBN, smaller than alpha_0^2 < 10^{-6.5} beta^{-1} (Omega_{matter} h^2 / 0.15)^{-3/2} when beta > 0.5.

It is shown that a large class of higher-order (i.e. non-quadratic) scalar kinetic terms can, without the help of potential terms, drive an inflationary evolution starting from rather generic initial conditions. In many models, this kinetically driven inflation (or "k-inflation" for short) rolls slowly from a high-curvature initial phase, down to a low-curvature phase and can exit inflation to end up being radiation-dominated, in a naturally graceful manner. We hope that this novel inflation mechanism might be useful in suggesting new ways of reconciling the string dilaton with inflation.

We describe in detail how perturbations due to the planets can cause a sub-population of WIMPs captured by scattering in surface layers of the Sun to evolve to have orbits which no longer intersect the Sun. We argue that such WIMPs, if their orbit has a semi-major axis less than 1/2 of Jupiter's, can persist in the solar system for cosmological timescales. This leads to a new, previously unanticipated WIMP population intersecting the Earth's orbit. The WIMP-nucleon cross sections required for this population to be significant are precisely those in the range predicted for SUSY dark matter, lying near the present limits obtained by direct underground dark matter searches using cyrogenic detectors. Thus, if a WIMP signal is observed in the next generation of detectors, a potentially measurable signal due to this new population must exist. This signal, lying in the keV range for Germanium detectors, would be complementary to that of galactic halo WIMPs. A comparison of event rates, anisotropies, and annual modulations would not only yield additional confirmation that any claimed signal is indeed WIMP-based, but would also allow one to gain information on the nature of the underlying dark matter model.

Perturbations due to the planets combined with the non-Coulomb nature of the gravitational potential in the Sun imply that WIMPs that are gravitationally captured by scattering in surface layers of the Sun can evolve into orbits that no longer intersect the Sun. For orbits with a semi-major axis $ < 1/2$ of Jupiter's orbit, such WIMPs can persist in the solar system for $ > 10^9$ years, leading to a previously unanticipated population intersecting Earth's orbit. For WIMPs detectable in the next generation of detectors, this population can provide a complementary signal, in the keV range, to that of galactic halo dark matter.

We map the general relativistic two-body problem onto that of a test particle moving in an effective external metric. This effective-one-body approach defines, in a non-perturbative manner, the late dynamical evolution of a coalescing binary system of compact objects. The transition from the adiabatic inspiral, driven by gravitational radiation damping, to an unstable plunge, induced by strong spacetime curvature, is predicted to occur for orbits more tightly bound than the innermost stable circular orbit in a Schwarzschild metric of mass M = m1 + m2. The binding energy, angular momentum and orbital frequency of the innermost stable circular orbit for the time-symmetric two-body problem are determined as a function of the mass ratio.

We formulate the basic postulate of pre-big bang cosmology as one of ``asymptotic past triviality'', by which we mean that the initial state is a generic perturbative solution of the tree-level low-energy effective action. Such a past-trivial ``string vacuum'' is made of an arbitrary ensemble of incoming gravitational and dilatonic waves, and is generically prone to gravitational instability, leading to the possible formation of many black holes hiding singular space-like hypersurfaces. Each such singular space-like hypersurface of gravitational collapse becomes, in the string-frame metric, the usual big-bang t=0 hypersurface, i.e. the place of birth of a baby Friedmann universe after a period of dilaton-driven inflation. Specializing to the spherically-symmetric case, we review and reinterpret previous work on the subject, and propose a simple, scale-invariant criterion for collapse/inflation in terms of asymptotic data at past null infinity. Those data should determine whether, when, and where collapse/inflation occurs, and, when it does, fix its characteristics, including anisotropies on the big bang hypersurface whose imprint could have survived till now. Using Bayesian probability concepts, we finally attempt to answer some fine-tuning objections recently moved to the pre-big bang scenario.

Binary systems comprising at least one neutron star contain strong gravitational field regions and thereby provide a testing ground for strong-field gravity. Two types of data can be used to test the law of gravity in compact binaries: binary pulsar observations, or forthcoming gravitational-wave observations of inspiralling binaries. We compare the probing power of these two types of observations within a generic two-parameter family of tensor-scalar gravitational theories. Our analysis generalizes previous work (by us) on binary-pulsar tests by using a sample of realistic equations of state for nuclear matter (instead of a polytrope), and goes beyond a previous study (by C.M. Will) of gravitational-wave tests by considering more general tensor-scalar theories than the one-parameter Jordan-Fierz-Brans-Dicke one. Finite-size effects in tensor-scalar gravity are also discussed.

We draw attention to the possibility that inflation (i.e. accelerated expansion) might continue after the end of slow roll, during a period of fast oscillations of the inflaton field \phi . This phenomenon takes place when a mild non-convexity inequality is satisfied by the potential V(\phi). The presence of such a period of \phi-oscillation-driven inflation can substantially modify reheating scenarios. In some models the effect of these fast oscillations might be imprinted on the primordial perturbation spectrum at cosmological scales.

We study the deflection of light (and the redshift, or integrated time delay) caused by the time-dependent gravitational field generated by a localized material source lying close to the line of sight. Our calculation explicitly takes into account the full, near-zone, plus intermediate-zone, plus wave-zone, retarded gravitational field. Contrary to several recent claims in the literature, we find that the deflections due to both the wave-zone 1/r gravitational wave and the intermediate-zone 1/r^2 retarded fields vanish exactly. The leading total time-dependent deflection caused by a localized material source, such as a binary system, is proven to be given by the quasi-static, near-zone quadrupolar piece of the gravitational field, and therefore to fall off as the inverse cube of the impact parameter.

The existence of a dilaton (or moduli) with gravitational-strength coupling to matter imposes stringent constraints on the allowed energy scale of cosmic strings, $\eta$. In particular, superheavy gauge strings with $\eta \sim 10^{16} GeV$ are ruled out unless the dilaton mass $m_{\phi} \gsim 100 TeV$, while the currently popular value $m_{\phi} \sim 1 TeV$ imposes the bound $\eta \lsim 3 \times 10^{11} GeV$. Similar constraints are obtained for global topological defects. Some non-standard cosmological scenarios which can avoid these constraints are pointed out.

General relativity predicts that energy and momentum conservation laws hold and that preferred frames do not exist. The parametrised post-Newtonian formalism (PPN) phenomenologically quantifies possible deviations from general relativity. The PPN parameter alpha_3 (which identically vanishes in general relativity) plays a dual role in that it is associated both with a violation of the momentum conservation law, and with the existence of a preferred frame. By considering the effects of alpha_3 neq 0 in certain binary pulsar systems, it is shown that alpha_3 < 2.2 x 10^-20 (90% CL). This limit improves on previous results by several orders of magnitude, and shows that pulsar tests of alpha_3 rank (together with Hughes-Drever-type tests of local Lorentz invariance) among the most precise null experiments of physics.

Some recently discovered nonperturbative strong-field effects in tensor-scalar theories of gravitation are interpreted as a scalar analog of ferromagnetism: "spontaneous scalarization". This phenomenon leads to very significant deviations from general relativity in conditions involving strong gravitational fields, notably binary-pulsar experiments. Contrary to solar-system experiments, these deviations do not necessarily vanish when the weak-field scalar coupling tends to zero. We compute the scalar "form factors" measuring these deviations, and notably a parameter entering the pulsar timing observable gamma through scalar-field-induced variations of the inertia moment of the pulsar. An exploratory investigation of the confrontation between tensor-scalar theories and binary-pulsar experiments shows that nonperturbative scalar field effects are already very tightly constrained by published data on three binary-pulsar systems. We contrast the probing power of pulsar experiments with that of solar-system ones by plotting the regions they exclude in a generic two-dimensional plane of tensor-scalar theories.

Gravity being a long-range force, one might {\it a priori} expect the Universe's global matter distribution to select a preferred rest frame for local gravitational physics. The phenomenology of preferred-frame effects, in the strong-gravitational field context of binary pulsars, is described by two parameters $\hat{\alpha}_1$ and $\hat{\alpha}_2$. These parameters vanish identically in general relativity, and reduce, in the weak-field limit, to the two parametrized post-Newtonian (PPN) parameters ${\alpha}_1$ and ${\alpha}_2$. We derive a limit of $|\hat{\alpha}_1| < 1.7\times 10^{-4}$ (90\%~C.L.) using the very low eccentricity binary pulsar PSR J2317+1439, improving by a factor of 3 on previous limits.