An idealized "test" object in general relativity moves along a geodesic. However, if the object has a finite mass, this will create additional curvature in the spacetime, causing it to deviate from geodesic motion. If the mass is nonetheless sufficiently small, such an effect is usually treated perturbatively and is known as the gravitational self-force due to the object. This issue is still an open problem in gravitational physics today, motivated not only by basic foundational interest, but also by the need for its direct application in gravitational-wave astronomy. In particular, the observation of extreme-mass-ratio inspirals by the future space-based detector LISA will rely crucially on an accurate modeling of the self-force driving the orbital evolution and gravitational wave emission of such systems. In this paper, we present a novel derivation, based on conservation laws, of the basic equations of motion for this problem. They are formulated with the use of a quasilocal (rather than matter) stress-energy-momentum tensor---in particular, the Brown-York tensor---so as to capture gravitational effects in the momentum flux of the object, including the self-force. Our formulation and resulting equations of motion are independent of the choice of the perturbative gauge. We show that, in addition to the usual gravitational self-force term, they also lead to an additional "self-pressure" force not found in previous analyses, and also that our results correctly recover known formulas under appropriate conditions. Our approach thus offers a fresh geometrical picture from which to understand the self-force fundamentally, and potentially useful new avenues for computing it practically.

We revisit for 714 SNeIa the discrepancy between the red-shift associated to the distance modulus $\mu$ and the spectroscopic red-shift. Previous work has shown that the total red-shift $z$ might be a combination of the expansion red-shift $z_{\rm C}$ and of a static, blue or red shift $z_{\rm LSV}(r)$, $r$ being the comoving distance. The latter is due to the energy non-conservation of the photon propagating through Electro-Magnetic (EM) background fields (host galaxy, intergalactic and Milky Way), under Lorentz(-Poincar\'e) Symmetry Violation (LSV), associated to the Standard-Model Extension (SME). The non-conservation stems from the vacuum expectation value of the vector and tensor LSV fields. For zero radiation $\Omega_{\rm rad}$ and curvature $\Omega_{\rm k}$ densities, and matter density $\Omega_m = 0.28$, the SN1a positions in the ($\mu$, z) plan are recovered according to the different strengths, orientations, alignments and space-time dependencies of the EM fields and LSV components. The LSV vacuum energy may be thus tantamount to $\Omega_\Lambda \simeq 0.7$, but unrelated to an accelerated expansion. We present models with red or blue-shifts $z_{\rm LSV}$, below $10\%$ of $z$. The $\nu$ frequency variation is below $10^{-19} \Delta \nu/\nu$ per m.

We propose a space-based interferometer surveying the gravitational wave (GW) sky in the milli-Hz to $\mu$-Hz frequency range. By the 2040s', the $\mu$-Hz frequency band, bracketed in between the Laser Interferometer Space Antenna (LISA) and pulsar timing arrays, will constitute the largest gap in the coverage of the astrophysically relevant GW spectrum. Yet many outstanding questions related to astrophysics and cosmology are best answered by GW observations in this band. We show that a $\mu$-Hz GW detector will be a truly overarching observatory for the scientific community at large, greatly extending the potential of LISA. Conceived to detect massive black hole binaries from their early inspiral with high signal-to-noise ratio, and low-frequency stellar binaries in the Galaxy, this instrument will be a cornerstone for multimessenger astronomy from the solar neighbourhood to the high-redshift Universe.

In the presence of Lorentz Symmetry Violation (LSV) associated with the Standard-Model Extension (SME), we have recently shown the non-conservation of the energy-momentum tensor of a light-wave crossing an Electro-Magnetic (EM) background field even when the latter and the LSV are constant. Incidentally, for a space-time dependent LSV, the presence of an EM field is not necessary. Herein, we infer that in a particle description, the energy non-conservation for a photon implies violation of frequency invariance in vacuo giving rise to a red or blue shift. We discuss the potential consequences on cosmology.

Within the Standard Model Extension (SME), we expand our previous findings on four classes of violations of Super-Symmetry (SuSy) and Lorentz Symmetry (LoSy), differing in the handedness of the Charge conjugation-Parity-Time reversal (CPT) symmetry and in whether considering the impact of photinos on photon propagation. The violations, occurring at the early universe high energies, show visible traces at present in the Dispersion Relations (DRs). For the CPT-odd classes ($V_{\mu}$ breaking vector) associated with the Carroll-Field-Jackiw (CFJ) model, the DRs and the Lagrangian show for the photon an effective mass, gauge invariant, proportional to $|{\vec V}|$. The group velocity exhibits a classic dependency on the inverse of the frequency squared. For the CPT-even classes ($k_{F}$ breaking tensor), when the photino is considered, the DRs display also a massive behaviour inversely proportional to a coefficient in the Lagrangian and to a term linearly dependent on $k_{F}$. All DRs display an angular dependence and lack LoSy invariance. In describing our results, we also point out the following properties: i) the appearance of complex or simply imaginary frequencies and super-luminal speeds and ii) the emergence of bi-refringence. Finally, we point out the circumstances for which SuSy and LoSy breakings, possibly in presence of an external field, lead to the non-conservation of the photon energy-momentum tensor. We do so for both CPT sectors.

Comments are due on a recent paper by McGruder III (2017) in which the author deals with the concept of gravitational repulsion in the context of the Schwarzschild-Droste solution. Repulsion (deceleration) for ingoing particles into a black hole is a concept proposed several times starting from Droste himself in 1916. It is a coordinate effect appearing to an observer at a remote distance from the black hole and when coordinate time is employed. Repulsion has no bearing and relation to the local physics of the black hole, and moreover it cannot be held responsible for accelerating outgoing particles. Thereby, the energy boost of cosmic rays cannot be produced by repulsion.

The photon mass, $m_\gamma$, can in principle be constrained using measurements of the dispersion measures (DMs) of fast radio bursts (FRBs), once the FRB redshifts are known. The DM of the repeating FRB 121102 is known to $< 1$\%, a host galaxy has now been identified with high confidence,and its redshift, $z$, has now been determined with high accuracy: $z = 0.19273(8)$. Taking into account the plasma contributions to the DM from the Intergalactic medium (IGM) and the Milky Way, we use the data on FRB 121102 to derive the constraint $m_\gamma \lesssim 2.2 \times 10^{-14}$ eV c$^{-2}$ ($3.9 \times 10^{-50}$ kg). Since the plasma and photon mass contributions to DMs have different redshift dependences, they could in principle be distinguished by measurements of more FRB redshifts, enabling the sensitivity to $m_\gamma$ to be improved.

The gravitational waves emitted by binary systems with extreme mass ratios carry unique astrophysical information expected to be probed by the next generation of gravitational wave detectors such as LISA. The detection of these binaries rely on an accurate modeling of the gravitational self-force that drives their orbital evolution. Although the theoretical formalism to compute the self-force has been largely established, the mathematical tools needed to implement it are still under development, and the self-force computation remains an open problem. We present here a frequency-domain implementation of the particle-without-particle (PwP) technique previously developed for the computation of the scalar self-force -- a helpful testbed for the gravitational self-force.

Our understanding of the universe relies mostly on electromagnetism. As photons are the messengers, fundamental physics is concerned in testing their properties. Photon mass upper limits have been earlier set through pulsar observations, but new investigations are offered by the excess of dispersion measure (DM) sometimes observed with pulsar and magnetar data at low frequencies, or with the fast radio bursts (FRBs), of yet unknown origin. Arguments for the excess of DM do not reach a consensus, but are not mutually exclusive. Thus, we remind that for massive electromagnetism, dispersion goes as the inverse of the frequency squared. Thereby, new avenues are offered also by the recently operating ground observatories in 10-80 MHz domain and by the proposed Orbiting Low Frequency Antennas for Radio astronomy (OLFAR). The latter acts as a large aperture dish by employing a swarm of nano-satellites observing the sky for the first time in the 0.1 - 15 MHz spectrum. The swarm must be deployed sufficiently away from the ionosphere to avoid distortions especially during the solar maxima, terrestrial interference and offer stable conditions for calibration during observations.

We show that photons may be redshifted or blueshifted when interacting with the field of an overcritical dipole, which incorporates the one-loop QED corrections coming from vacuum polarization. Using the effective metric, it follows that such effect depends on the polarization of the photon. The shifts, plotted against the azimuthal angle for various values of the magnetic field, may show an intensity comparable to the gravitational redshift for a magnetar. To obtain these results, we have corrected previous literature.

In the context of Standard Model Extensions (SMEs), we analyse four general classes of Super Symmetry (SuSy) and Lorentz Symmetry (LoSy) breaking, leading to {observable} imprints at our energy scales. The photon dispersion relations show a non-Maxwellian behaviour for the CPT (Charge-Parity-Time reversal symmetry) odd and even sectors. The group velocities exhibit also a directional dependence with respect to the breaking background vector (odd CPT) or tensor (even CPT). In the former sector, the group velocity may decay following an inverse squared frequency behaviour. Thus, we extract a massive and gauge invariant Carroll-Field-Jackiw photon term in the Lagrangian and show that the mass is proportional to the breaking vector. The latter is estimated by ground measurements and leads to a photon mass upper limit of $10^{-19}$ eV or $2 \times 10^{-55}$ kg and thereby to a potentially measurable delay at low radio frequencies.

Our understanding of the universe at large and small scales relies largely on electromagnetic observations. As photons are the messengers, fundamental physics has a concern in testing their properties, including the absence of mass. We use Cluster four spacecraft data in the solar wind at 1 AU to estimate the mass upper limit for the photon. We look for deviations from Amp\`ere's law, through the curlometer technique for the computation of the magnetic field, and through the measurements of ion and electron velocities for the computation of the current. We show that the upper bound for $m_\gamma$ lies between $1.4 \times 10^{-49}$ and $3.4 \times 10^{-51}$ kg, and thereby discuss the currently accepted lower limits in the solar wind.

The frequency-dependent time delays in fast radio bursts (FRBs) can be used to constrain the photon mass, if the FRB redshifts are known, but the similarity between the frequency dependences of dispersion due to plasma effects and a photon mass complicates the derivation of a limit on $m_\gamma$. The dispersion measure (DM) of FRB 150418 is known to $\sim 0.1$%, and there is a claim to have measured its redshift with an accuracy of $\sim 2$%, but the strength of the constraint on $m_\gamma$ is limited by uncertainties in the modelling of the host galaxy and the Milky Way, as well as possible inhomogeneities in the intergalactic medium (IGM). Allowing for these uncertainties, the recent data on FRB 150418 indicate that $m_\gamma \lesssim 1.8 \times 10^{-14}$eV c$^{-2}$ ($3.2 \times 10^{-50}$kg), if FRB 150418 indeed has a redshift $z = 0.492$ as initially reported. In the future, the different redshift dependences of the plasma and photon mass contributions to DM can be used to improve the sensitivity to $m_\gamma$ if more FRB redshifts are measured. For a fixed fractional uncertainty in the extra-galactic contribution to the DM of an FRB, one with a lower redshift would provide greater sensitivity to $m_\gamma$.

We apply our method of indirect integration, described in Part I, at fourth order, to the radial fall affected by the self-force. The Mode-Sum regularisation is performed in the Regge-Wheeler gauge using the equivalence with the harmonic gauge for this orbit. We consider also the motion subjected to a self-consistent and iterative correction determined by the self-force through osculating stretches of geodesics. The convergence of the results confirms the validity of the integration method. This work complements and justifies the analysis and the results appeared in Int. J. Geom. Meth. Mod. Phys., 11, 1450090 (2014).

The Regge-Wheeler-Zerilli (RWZ) wave-equation describes Schwarzschild-Droste black hole perturbations. The source term contains a Dirac distribution and its derivative. We have previously designed a method of integration in time domain. It consists of a finite difference scheme where analytic expressions, dealing with the wave-function discontinuity through the jump conditions, replace the direct integration of the source and the potential. Herein, we successfully apply the same method to the geodesic generic orbits of EMRI (Extreme Mass Ratio Inspiral) sources, at second order. An EMRI is a Compact Star (CS) captured by a Super Massive Black Hole (SMBH). These are considered the best probes for testing gravitation in strong regime. The gravitational wave-forms, the radiated energy and angular momentum at infinity are computed and extensively compared with other methods, for different orbits (circular, elliptic, parabolic, including zoom-whirl).

Radial fall has historically played a momentous role. It is one of the most classical problems, the solutions of which represent the level of understanding of gravitation in a given epoch. A {\it gedankenexperiment} in a modern frame is given by a small body, like a compact star or a solar mass black hole, captured by a supermassive black hole. The mass of the small body itself and the emission of gravitational radiation cause the departure from the geodesic path due to the back-action, that is the self-force. For radial fall, as any other non-adiabatic motion, the instantaneous identity of the radiated energy and the loss of orbital energy cannot be imposed and provide the perturbed trajectory. In the first part of this letter, we present the effects due to the self-force computed on the geodesic trajectory in the background field. Compared to the latter trajectory, in the Regge-Wheeler, harmonic and all others smoothly related gauges, a far observer concludes that the self-force pushes inward (not outward) the falling body, with a strength proportional to the mass of the small body for a given large mass; further, the same observer notes an higher value of the maximal coordinate velocity, this value being reached earlier on during infall. In the second part of this letter, we implement a self-consistent approach for which the trajectory is iteratively corrected by the self-force, this time computed on osculating geodesics. Finally, we compare the motion driven by the self-force without and with self-consistent orbital evolution. Subtle differences are noticeable, even if self-force effects have hardly the time to accumulate in such a short orbit.

We adopt the Dirac-Detweiler-Whiting radiative and regular effective field in curved spacetime. Thereby, we derive straightforwardly the first order perturbative correction to the geodesic of the background in a covariant form, for the extreme mass ratio two-body problem. The correction contains the self-force contribution and a background metric dependent term.

Gravitational waves coming from Super Massive Black Hole Binaries (SMBHBs) are targeted by both Pulsar Timing Array (PTA) and Space Laser Interferometry (SLI). The possibility of a single SMBHB being tracked first by PTA, through inspiral, and later by SLI, up to merger and ring down, has been previously suggested. Although the bounding parameters are drawn by the current PTA or the upcoming Square Kilometer Array (SKA), and by the New Gravitational Observatory (NGO), derived from the Laser Interferometer Space Antenna (LISA), this paper also addresses sequential detection beyond specific project constraints. We consider PTA-SKA, which is sensitive from 10^(-9) to p x 10^(-7) Hz (p=4, 8), and SLI, which operates from s x 10^(-5) up to 1 Hz (s = 1, 3). A SMBHB in the range 2x 10^(8) - 2 x 10^(9) solar masses (the masses are normalised to a (1+z) factor, the red shift lying between z = 0.2 and z=1.5) moves from the PTA-SKA to the SLI band over a period ranging from two months to fifty years. By combining three Super Massive Black Hole (SMBH)-host relations with three accretion prescriptions, nine astrophysical scenarios are formed. They are then related to three levels of pulsar timing residuals (50, 5, 1 ns), generating twenty-seven cases. For residuals of 1 ns, sequential detection probability will never be better than 4.7 x 10^(-4) y^(-2) or 3.3 x 10^(-6) y^(-2) (per year to merger and per year of survey), according to the best and worst astrophysical scenarios, respectively; put differently this means one sequential detection every 46 or 550 years for an equivalent maximum time to merger and duration of the survey. The chances of sequential detection are further reduced by increasing values of the s parameter (they vanish for s = 10) and of the SLI noise, and by decreasing values of the remnant spin. REST OF THE ABSTRACT IN THE PDF FILE.

We intend to develop part of the theoretical tools needed for the detection of gravitational waves coming from the capture of a compact object, 1-100 solar masses, by a Supermassive Black Hole, up to a 10 billion solar masses, located at the centre of most galaxies. The analysis of the accretion activity unveils the star population around the galactic nuclei, and tests the physics of black holes and general relativity. The captured small mass is considered a probe of the gravitational field of the massive body, allowing a precise measurement of the particle motion up to the final absorption. The knowledge of the gravitational signal, strongly affected by the self-force - the orbital displacement due to the captured mass and the emitted radiation - is imperative for a successful detection. The results include a strategy for wave equations with a singular source term for all type of orbits. We are now tackling the evolution problem, first for radial fall in Regge- Wheeler gauge, and later for generic orbits in the harmonic or de Donder gauge for Schwarzschild-Droste black holes. In the Extreme Mass Ratio Inspiral, the determination of the orbital evolution demands that the motion of the small mass be continuously corrected by the self-force, i.e. the self-consistent evolution. At each of the integration steps, the self-force must be computed over an adequate number of modes; further, a differential-integral system of general relativistic equations is to be solved and the outputs regularised for suppressing divergences. Finally, for the provision of the computational power, parallelisation is under examination.

On the basis of a recently proposed strategy of finite element integration in time domain for partial differential equations with a singular source term, we present a fourth order algorithm for non-rotating black hole perturbations in the Regge-Wheeler gauge. Herein, we address even perturbations induced by a particle plunging in. The forward time value at the upper node of the $(r^*,t)$ grid cell is obtained by an algebraic sum of i) the preceding node values of the same cell, ii) analytic expressions, related to the jump conditions on the wave function and its derivatives, iii) the values of the wave function at adjacent cells. In this approach, the numerical integration does not deal with the source and potential terms directly, for cells crossed by the particle world line. This scheme has also been applied to circular and eccentric orbits and it will be object of a forthcoming publication.

Perturbations of Schwarzschild-Droste black holes in the Regge-Wheeler gauge benefit from the availability of a wave equation and from the gauge invariance of the wave function, but lack smoothness. Nevertheless, the even perturbations belong to the C\textsuperscript{0} continuity class, if the wave function and its derivatives satisfy specific conditions on the discontinuities, known as jump conditions, at the particle position. These conditions suggest a new way for dealing with finite element integration in time domain. The forward time value in the upper node of the $(t, r^*$) grid cell is obtained by the linear combination of the three preceding node values and of analytic expressions based on the jump conditions. The numerical integration does not deal directly with the source term, the associated singularities and the potential. This amounts to an indirect integration of the wave equation. The known wave forms at infinity are recovered and the wave function at the particle position is shown. In this series of papers, the radial trajectory is dealt with first, being this method of integration applicable to generic orbits of EMRI (Extreme Mass Ratio Inspiral).

We establish the jump conditions for the wavefunction and its derivatives through the formal solutions of the wave equation. These conditions respond to the requirement of continuity of the perturbations at the position of the particle and they are given for any mode at first order. Using these jump conditions, we then propose a new method for computing the radiated waveform without direct integration of the source term. We consider this approach potentially applicable to generic orbits.

Free fall has signed the greatest markings in the history of physics through the leaning Pisa tower, the Cambridge apple tree and the Einstein lift. The perspectives offered by the capture of stars by supermassive black holes are to be cherished, because the study of the motion of falling stars will constitute a giant step forward in the understanding of gravitation in the regime of strong field. After an account on the perception of free fall in ancient times and on the behaviour of a gravitating mass in Newtonian physics, this chapter deals with last century debate on the repulsion for a Schwarzschild black hole and mentions the issue of an infalling particle velocity at the horizon. Further, black hole perturbations and numerical methods are presented, paving the way to the introduction of the self-force and other back-action related methods. The impact of the perturbations on the motion of the falling particle is computed via the tail, the back-scattered part of the perturbations, or via a radiative Green function. In the former approach, the self-force acts upon the background geodesic; in the latter, the geodesic is conceived in the total (background plus perturbations) field. Regularisation techniques (mode-sum and Riemann-Hurwitz $z$ function) intervene to cancel divergencies coming from the infinitesimal size of the particle. An account is given on the state of the art, including the last results obtained in this most classical problem, together with a perspective encompassing future space gravitational wave interferometry and head-on particle physics experiments. As free fall is patently non-adiabatic, it requires the most sophisticated techniques for studying the evolution of the motion. In this scenario, the potential of the self-consistent approach, by means of which the background geodesic is continuously corrected by the self-force contribution, is examined.

Through detection by low gravitational wave space interferometers, the capture of stars by supermassive black holes will constitute a giant step forward in the understanding of gravitation in strong field. The impact of the perturbations on the motion of the star is computed via the tail, the back-scattered part of the perturbations, or via a radiative Green function. In the former approach, the self-force acts upon the background geodesic, while in the latter, the geodesic is conceived in the total (background plus perturbations) field. Regularisations (mode-sum and Riemann-Hurwitz $\zeta$ function) intervene to cancel divergencies coming from the infinitesimal size of the particle. The non-adiabatic trajectories require the most sophisticated techniques for studying the evolution of the motion, like the self-consistent approach.

The coalescence rate of two neutron stars (NS) is revisited. For estimation of the number of bound NS-NS and the probability of their coalescence in a timescale $\tau$, the galactic star formation history, directly derived from observations, and the evolution of massive stars are considered. The newly established galactic merging rate is $(1.7\pm 1.0) \times 10^{-5} yr^{-1}$, while the local merging rate, including the contribution of elliptical galaxies, is about a factor of two higher, $3.4 \times 10^{-5} yr^{-1}$. Using the present data basis on galaxy distribution in the local universe and the expected sensitivity of the first generation of laser beam interferometers, we estimate that one event should occur every 125 years for LIGO and one event each 148 years for VIRGO. The situation is considerably improved for advanced-LIGO since we predict that 6 events per year should be detected whereas for a recently proposed VIRGO new configuration, the event rate might increase up to 3 events every two years.

In this paper we present new estimates of the coalescence rate of neutron star binaries in the local universe and we discuss its consequences for the first generations of ground based interferometers. Our approach based on both evolutionary and statistical methods gives a galactic merging rate of 1.7 10$^{-5}$ yr$^{-1}$, in the range of previous estimates 10$^{-6}$ - 10$^{-4}$ yr$^{-1}$. The local rate which includes the contribution of elliptical galaxies is two times higher, in the order of 3.4 10$^{-5}$ yr$^{-1}$. We predict one detection every 148 and 125 years with initial VIRGO and LIGO, and up to 6 events per year with their advanced configuration. Our recent detection rate estimates from investigations on VIRGO future improvements are quoted.

The concept of a measurement of the yet unevaluated Hannay angle, by means of an Earth-bound satellite, adiabatically driven by the Moon, is shown herein. Numerical estimates are given for the angles, the orbital displacements, the shortening of the orbital periods, for different altitudes. It is concluded that the Hannay effect is measurable in high Earth orbits, by means of atomic clocks, accurate Time & Frequency transfer system and precise positioning.

The Hannay angle has been previously studied for a celestial circular restricted three-body system by means of an adiabatic approach. In the present work, three main results are obtained. Firstly, a formal connection between perturbation theory and the Hamiltonian adiabatic approach shows that both lead to the Hannay angle; it is thus emphasised that this effect is already contained in classical celestial mechanics, although not yet defined nor evaluated separately. Secondly, a more general expression of the Hannay angle, valid for an action-dependent potential is given; such a generalised expression takes into account that the restricted three-body problem is a time-dependent, two degrees of freedom problem even when restricted to the circular motion of the test body. Consequently, (some of) the eccentricity terms cannot be neglected {\it a priori}. Thirdly, we present a new numerical estimate for the Earth adiabatically driven by Jupiter. We also point out errors in a previous derivation of the Hannay angle for the circular restricted three-body problem, with an action-independent potential.

This work deals with the motion of a radially falling star in Schwarzschild geometry and correctly identifies radiation reaction terms by the perturbative method. The results are: i) identification of all terms up to first order in perturbations, second in trajectory deviation, and mixed terms including lowest order radiation reaction terms; ii) renormalisation of all divergent terms by the $\zeta$ Riemann and Hurwitz functions. The work implements a method previously identified by one of the authors and corrects some current misconceptions and results.