# Pinpointing gravitational waves via astrometric gravitational wave antennas

Mariateresa Crosta, Mario Gilberto Lattanzi, Christophe Le Poncin-Lafitte, Mario Gai, Qi Zhaoxiang, Alberto Vecchiato

► To cite this version:

Mariateresa Crosta, Mario Gilberto Lattanzi, Christophe Le Poncin-Lafitte, Mario Gai, Qi Zhaoxiang, et al. Pinpointing gravitational waves via astrometric gravitational wave antennas. Scientific Reports, 2024, 14 (1), pp.5074. 10.1038/s41598-024-55671-9. hal-04487697

HAL Id: hal-04487697

https://hal.science/hal-04487697v1

Submitted on 6 Mar 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_2.png" /></figure>

Check for updates

OPEN

# Pinpointing gravitational waves via astrometric gravitational wave antennas

Mariateresa Crosta<fnref n="1" />✉, Mario Gilberto Lattanzi<fnref n="1" />, Christophe Le Poncin-Lafitte², Mario Gai<fnref n="1" />, Qi Zhaoxiang³ & Alberto Vecchiato<fnref n="1" />

The direct detection of gravitational waves by ground-based optical interferometers has opened a new window in astronomy. Besides, the sensitivity of these linear detectors to the direction of arrival of an incoming gravitational wave is limited compared to current prospects of high-precision, space-based, astrometry. Indeed, advanced methods of differential relativistic astrometry offer a unique opportunity to overcome that situation. Here, we present a novel concept for a gravitational wave antenna that uses angles between close pairs of point-like sources as natural (angular) "arms" to characterise the very tiny variations in angular separations induced by a passing gravitational wave. The proposed new astrometric gravitational wave observable proves to be a powerful tool to substantially enhance the effect of gravitational waves of different strengths by exploiting optical resolution to the fullest. Then, by optically multiplexing three (or more) of such astrometric "arms", it would be also possible to pinpoint source directions to unprecedented levels.

Experimental confirmation of gravitational waves (GWs) through the LIGO<fnref n="1" />,² and VIRGO³ antennas has gained great impulse to the search and characterization of candidate GW sources. New ground-based experiments worldwide are already⁴, or in the process of⁵⁻⁷, joining the effort, and the LISA mission⁸ will soon implement similar concepts in space. The prime objective of this effort is the complete characterization of GWs, i.e., to determine their amplitude and frequency spectra, and then pinpoint their directions for multi-wavelength "profiling". A review on gravitational-wave physics and astronomy in the present decade can be found in Ref.⁹.

Astrometric (angular) observations are generally targeted for the accurate determination of directions to the incoming photons and their change with time. Therefore they collect photons that have interacted with different time-dependent gravitational fields along their path to the observer. With the advent of highly accurate astrometric and radial velocity measurements in space, observation reduction models compliant with General Relativity (GR) have become a necessity. As, when compared to the targeted measurement precision levels, the relatively small amount of space-time curvature due to the Sun and all of the other relevant Solar System (SS) masses (including the Earth-Moon system) affecting incoming photons and space-borne "observers" alike can no longer be ignored. The weak gravity regime influences electromagnetic propagation on a much wider domain than its strong counterpart.

Gravitational wave detection via astrometry was explored by some authors¹⁰⁻¹⁵. These authors consider the extra-shifts on a single light direction induced by passing GWs, i.e., on photons propagating to an observer from within the SS. Such a detection is not considered promising, since it requires at least nanoarcsecond (nas) accuracy¹⁶⁻¹⁸, and demands knowledge of satellite attitude to similar levels. Also, initial investigations on the potential of the Gaia astrometry in GW research focused on periodic GW signals with period shorter than the Gaia operating time¹⁹,²⁰, or on secular effects²¹ on QSO proper motions for longer periods, i.e. for ultra-low GW frequencies of cosmological nature.

As a matter of fact, the astrometric observable considered so far, like in the case of Gaia, is expressed as the direction cosine between incoming stellar light and the observer attitude-tetrad²². Its strict application to a passing GW, being the tiny angular variations to stellar directions in the argument of the direction cosines, imposes a very strong requirement on the knowledge of satellite attitude at, or below, the microarcsecond (μas) accuracy, beyond current feasibility.

<footnote n="1">¹Astrophysical Observatory of Turin (OATo)-INAF, via Osservatorio 20, 10025 Pino, Torinese, Italy. ²SYRTE, Observatoire de Paris, Université PSL, CNRS, Sorbonne Université, LNE, 61 avenue de l'Observatoire, 75014 Paris, France. ³Chinese Academy of Sciences (CAS), Shanghai Astronomical Observatory (SHAO), 80 Nandan Road, Shanghai 200030, China. ✉email: mariateresa.crosta@inaf.it</footnote>
Disentangling GW signals also requires that the SS background metric is consistently developed to account for possible background (natural) systematics, thus preventing unwanted effects in modelling our astrometric observable. With $v/c \sim 10^{-4}$ rad, $v$ being the typical velocity of the relevant SS metric sources, terms in the null geodesic solution should be retained to $v/c \sim 10^{-4}$ at least, i.e., below the nano-arcsec (nas) level¹⁸.

To overcome the above limitations, here we present a novel idea of an Astrometric GW Antenna that fully exploits a GR formulation for differential astrometry, i.e., uses as antenna arms angles between close stellar pairs as measured at the observer's location.

We first introduce a different fundamental observation equation for the astrometric GW detection based on sufficiently narrow angles between double stellar-like sources each materializing pairs of local lines of sight (LOS's). We then give a first evaluation of its potential impact on GW science and its practical feasibility, in the light of new, specialized and much improved concepts for space astrometry missions²³⁻²⁵.

# Results

## The fundamental observation equation for the astrometric gravitational wave antenna

Let us assume the global metric due to both SS sources and passing or standing GW perturbations at the observer's location in the form

$$g_{\alpha\beta} = g^{SS}_{\alpha\beta} + h^{GW}_{\alpha\beta} = \eta_{\alpha\beta} + \sum_{(a)} h^{SS}_{(a)\alpha\beta} + h^{GW}_{\alpha\beta} + O(h^2), \quad (1)$$

where $\eta_{\alpha\beta}$ is the flat Minkowskian metric with signature $(-, +, +, +)$ and the subscript (a) stands for the a-source. The cosine of the angle between two observed (obs) light directions $\ell_1^\alpha$ and $\ell_2^\alpha$ writes (see Methods for the derivation)

$$\cos \psi_{1,2} = g_{\alpha\beta} (\ell_1^\alpha \ell_2^\beta)_{\text{obs}}, \quad (2)$$

where $\ell^\alpha$ is the null tangent unit four-vector projected on the rest space of the local barycentric observer, namely, for our SS, the observer at rest relative to the barycentric celestial reference system (BCRS) with coordinates $(t, x^i)$. Definition (2) guarantees that $\psi_{1,2}$ is an observed quantity; its differential nature greatly relaxes precision requirements on the knowledge of satellite attitude (see comments in Methods) and payload (thermal and mechanical) stability.

Similarly to what just done for the metric, stellar light directions can be represented as the SS part (due to the background metric) plus a perturbation shift, i.e. $\delta\ell^\alpha$, attributed purely to the passing GW:

$$\ell_{\text{obs}}^{\alpha} = \ell^{\alpha(\text{SS})} + \delta\ell^{\alpha} + O(\delta\ell^2), \quad (3)$$

in particular for the SS part

$$\ell_{\text{obs}}^{\alpha(\text{SS})} = \ell_0^{\alpha} + \epsilon \ell_{(1)}^{\alpha} + \epsilon^2 \ell_{(2)}^{\alpha} + \epsilon^3 \ell_{(3)}^{\alpha} + \epsilon^4 \ell_{(4)}^{\alpha} + O(\epsilon^5), \quad (4)$$

where $\epsilon = v/c$, being $v$ the typical velocity of each relevant SS metric source, and the subscripts in parenthesis indicate the order of approximation. Then, the right-hand side of (2) can be further simplified, assuming $\delta\ell$ starts at $\sim \epsilon^4$, as

$$\cos \psi_{1,2} = \cos \psi_{1,2}^{\text{SS}} + \eta_{\alpha\beta} (\ell_{10}^{\alpha} \delta\ell_{2}^{\beta} + \ell_{20}^{\alpha} \delta\ell_{1}^{\beta}) + h_{\alpha\beta}^{\text{GW}} \ell_{10}^{\alpha} \ell_{20}^{\beta} + O(\epsilon^5) + O(h^2), \quad (5)$$

$\ell_{10}^\alpha$ representing the unperturbed (Minkowskian) light direction to star i=1,2. Alternatively, (5) can be reformulated as

$$\cos \psi_{1,2} = \cos \psi_{1,2}^{\text{SS}} + F_{1,2}^{\text{GW}}, \quad (6)$$

where we collect all the GW terms related to the theoretical modelling in the expression $F_{1,2}^{\text{GW}} = \eta_{\alpha\beta} (\ell_{10}^\alpha \delta\ell_{2}^\beta + \ell_{20}^\alpha \delta\ell_{1}^\beta) + h_{\alpha\beta}^{\text{GW}} \ell_{10}^\alpha \ell_{20}^\beta + O(\epsilon^5) + O(h^2)$. Note that the right-hand side of Eqs. (5) and (6) is a spacetime invariant that can be evaluated in any coordinate system.

On the other hand, since the passing GW produces an extra shift on light deflection, we expect a perturbation $\delta\psi^{\text{GW}}$ to the undisturbed angle between two light directions.

Then, we can expand the left-hand side of our observation Eq. (2) (i.e., its observational part) as follows:

$$\cos(\psi_{1,2}^{\text{SS}} + \delta\psi_{1,2}^{\text{GW}}) = \cos(\psi_{1,2}^{\text{SS}}) \cos(\delta\psi_{1,2}^{\text{GW}}) - \sin(\psi_{1,2}^{\text{SS}}) \sin(\delta\psi_{1,2}^{\text{GW}}). \quad (7)$$

Considering, now, both Eqs. (6) and (7), we obtain

$$\cos(\psi_{1,2}^{\text{SS}}) \cos(\delta\psi_{1,2}^{\text{GW}}) - \sin(\psi_{1,2}^{\text{SS}}) \sin(\delta\psi_{1,2}^{\text{GW}}) = \cos \psi_{1,2}^{\text{SS}} + F_{1,2}^{\text{GW}}. \quad (8)$$

Since $\delta\psi^{\text{GW}} \ll 1$ and for $\delta\psi^{\text{GW}} \ll \psi^{\text{SS}}$, the previous equation reduces to consider the unknown angle due to the passing GW, i.e.

$$\delta\psi_{1,2}^{\text{GW}} = -\frac{F_{1,2}^{\text{GW}}}{\sin(\psi_{1,2}^{\text{SS}})}. \quad (9)$$
The new expression (9), deduced from observation Eq. (2), provides the observed perturbation to the $\psi_{1,2}^{SS}$ angle caused by the passing GW in the $F_{1,2}^{GW}$ term. Angular perturbations can be made larger through the factor $(\sin(\psi_{1,2}^{SS}))^{-1}$ that acts as a “signal amplifier” for the GW detection. That this amplification factor is not arbitrary is clearly a consequence of the condition $\delta\psi^{GW} \ll \psi^{SS}$ and the fact that (see below for further clarification on this point) an optical device, a telescope, implementing Eq. (9) has to be able to measure the angular distance of the 1,2 star-like pair materializing one arm of the astrometric antenna, and this is mainly limited by its resolving power.

Terms like the $1/\sin(\psi)$ factor are rather typical when the effects of gravitational perturbations need to be put to test as this is accomplished by utilizing measurements of angles between directions to pairs of stellar-like objects. An example is that in, e.g.,²⁶ where, after modeling the deflection of light from a target source due to a spherical mass along the line-of-sight, the Authors turn the discussion into the actual way of measuring the effect by using the angle $\theta$ between the target and a reference source, a direct analogue of our angle $\psi$ between two sources forming one arm of the antenna. As, in close analogy with linear GW antennas, the angle $\psi$ represents one angular arm of our proposed astrometric GW antenna. It is by increasing the length L of the arms of linear antennas that the effect of the GW perturbation becomes easier to measure through the relation $\delta L \sim h^{GW} L$; analogously, it is by being able to measure smaller and smaller separation angles (i.e., increasing the resolution power of the optics used) that we can increase the measurability on the GW-induced effect on an (angular) arm of the astrometric antenna via our relation (9), i.e., $\delta\psi \sim h^{GW}/\sin(\psi)$.

As we are putting forth a novel operational principle for measuring GWs that takes great advantage from the quantity $\psi_{1,2}^{SS}$, here it will suffice to indicate how that can be estimated/derived from actual measurements of the pair separation. What we measure is actually the angle $\psi_{1,2}(t_i)$ between point-like sources '1' and '2' at time $t_i$. The measurements $\psi_{1,2}(t_i)$ are taken with high cadence, i.e., with frequency $\omega_S \gg \omega_{GW}$, the oscillating frequency of the GW (assumed 'monochromatic') we seek to unveil. This last condition ensures that the Nyquist-Shannon (sampling) theorem is satisfied and, at the same time, sufficient statistics is built to beat (single) measurement noise.

Considering the average $\langle \psi_{1,2}(t_i) \rangle_N$ over the N separations $\psi_{1,2}(t_i)$ taken over a measurement session, with $N \gg 1$, we have: $\langle \psi_{1,2}(t_i) \rangle_N = \langle \psi_{1,2}^{SS} \rangle_N + \langle \delta\psi_{1,2}^{GW}(t_i) \rangle_N \simeq \hat{\psi}_{1,2}^{SS}$, as we can make the average $\langle \delta\psi_{1,2}^{GW}(t_i) \rangle_N$ as small as needed, much like it is done in signal spectral analysis. And, with this evaluation of the unperturbed angular separation of our antenna arm (the point-like pair) directly from the observations it follows that: $\delta\hat{\psi}_{1,2}^{GW}(t_i) \equiv \psi_{1,2}(t_i) - \hat{\psi}_{1,2}^{SS}$.

This is, in principle, how the GW perturbation to the antenna angle can be estimated directly from the measurements, along with $\hat{\psi}_{1,2}^{SS}$, and then used in Eq. (9) to build the observation equations from which strength and direction to the GW source can be recovered.

## The operating principle of the GW astrometric antenna

The measurement and initial data processing protocol just sketched supports, at least in principle, the practicability of our concept for an astrometric GW antenna.

To provide quantitative examples suggestive of the potential of Eq. (9), it will suffice to consider a plane GW, of a given frequency $\omega_{GW}$, linearly polarized (i.e., $A_\times = 0$, see below). Indeed, the actual application of Eq. (9) requires the definition of the following ingredients: (1) the $h^{GW}$ at the observer, (2) a pair of unperturbed local lines-of-sight, $\ell_{i0}^\alpha$ (see Eq. 3), and (3), the corresponding shift $\delta\ell$ due to the passing GW from the geodesic equation with metric (1).

The general form of the GW perturbation $\delta\ell$ can be expressed as a function of argument $\tilde{k}_\alpha x^\alpha$, namely $h_{ij}^{GW}(\tilde{k}_\alpha x^\alpha)$, with tangent vector $\tilde{k}^\alpha = \tilde{k}^0 \partial_0^\alpha + \tilde{k}^i \partial_i^\alpha$, where $p^i = \tilde{k}^i/\tilde{k}^0$ is the direction of the GW propagation. Then, in the linearized regime, the gravitational wave shift is recovered within the suitable astrometric models, via the geodesic integration with respect to a suitable parameter $\sigma$ (see “Methods” for details), namely:

$$ \int_{\sigma_{\text{obs}}}^{\sigma_{*}} d\delta\ell^i = \frac{\ell_0^i + p^i}{2(1 + p \cdot \ell_0)} \ell_0^j \ell_0^k \int_{\sigma_{\text{obs}}}^{\sigma_{*}} dh_{jk}^{GW}(\sigma) - \frac{\ell_0^j}{2} \int_{\sigma_{\text{obs}}}^{\sigma_{*}} dh_{ij}^{GW}(\sigma) + O(\epsilon^5), \quad (10) $$

that, in the far away zone, reduces to

$$ \delta\ell^i = \frac{\ell_0^i + p^i}{2(1 + p \cdot \ell_0)} \ell_0^j \ell_0^k h_{jk}^{GW}(\sigma_{\text{obs}}) - \frac{\ell_0^j}{2} h_{ij}^{GW}(\sigma_{\text{obs}}) + O(\epsilon^5), \quad (11) $$

explicitly showing that $\delta\ell^i \propto h^{GW}$, and in agreement with what found in Ref.¹¹,¹⁴. In such a case, our expression (9) becomes

$$ \delta\psi_{1,2}^{GW} = -\frac{h_{ij}^{GW}(\sigma_{\text{obs}})}{2 \sin(\psi_{1,2}^{SS})} \left\{ \frac{[(\ell_{10} \cdot \ell_{20}) + (\ell_{10} \cdot p)]\ell_{20}^i \ell_{20}^j}{(1 + p \cdot \ell_{20})} + \frac{[(\ell_{20} \cdot \ell_{10}) + (\ell_{20} \cdot p)]\ell_{10}^i \ell_{10}^j}{(1 + p \cdot \ell_{10})} \right\} + O(\epsilon^5), \quad (12) $$

which shows also the dependence on the scalar products between the SS Minkowskian directions of photon propagation from stars 1 and 2 (i.e., alternatively, the angle $\psi_{10,20}$) and on the scalar products of each star direction to the GW source (angles $\psi_{10,p}$ or $\psi_{20,p}$).

Lastly, it is worth stressing once more that Eq. (9) holds only for $\delta\psi^{GW} \ll \psi^{SS}$. Increasing the optical resolution power to infinity, i.e., for the pair separation angle $\psi$ going to zero, would make our observable
<figure><img src="image_3.png" /><figcaption>**Figure 1.** In (a) a possible configuration of a 3-LOS telescope, the (x,y,z) triad, with respect to the chosen directions: the local (SS) vectors $\ell_{20}$, $\ell_{50}$ and $\ell_{10}$ to stars 3,5 and 1, respectively. The angles $\psi_{ij}$, representing the instantaneous angular distances of the stellar pairs are all assumed small at ~0.01" for (b) and ~0.001" for (c). With the same colour coding, the other two panels show, fixed the resolution limit, the astrometric signals $\delta\psi_{ij}^{GW}$ (from Eq. (12)) caused by a GW passing along: (b) the z direction (green line) with $A_+ \sim 10^{-18}$ and frequency $10^3$Hz; (c) the x direction (red line) with amplitude $A_+ \sim 10^{-19}$ and frequency 0.1 Hz. The other curves show the perfect anti-correlation of the GW-induced signals along the other two LOS's of the astrometric antenna.</figcaption></figure>

$\delta\psi^{GW}$ degenerate. This can be easily seen from Eq. (6) where, with $\psi \to 0$, the two cosine terms on both sides of this equation $\to 1$, leaving $F^{GW} = 0$. As $F^{GW}$ does not depend on $\psi$ (Eq. 12), this can only mean that our observation Eq. (6) through Eq. (9) (or Eq. 12) degenerate and would no longer be useful to measure the effect of a passing $h^{GW}$.

Equation (12) above governs the operating principle of (one angular arm) of the astrometric antenna, and clearly shows its direct relation with the direction to an incoming GW. For, the GW term $h_{ij}$ (and its time variation) will mostly characterize the detection (amplitude and phase term), while the factor within the curly brackets will assume specific patterns according to the direction of the incoming GW relative to the observer (spatial) orientation.

## Pinpointing sources of gravitational waves

The spatial orientation of the observer, i.e., an optical system endowed with three viewing directions (or lines of sight-LOS) multiplexed onto a common focal plane, is idealized in Fig. 1a as a **x**, **y**, **z** triad and schematically illustrated in Fig. 2.

The local (SS) unperturbed direction to star 1 is along the same direction at which LOS **z** is pointing. Also, $\ell_{20}$, the local direction to star 2, is at a very small angular separation to $\ell_{10}$, i.e., $\cos(\psi_{10,20}) \sim 1$.

Similarly to the LOS in the **z** direction of our schematic 3-way telescope, one can imagine to have the other two viewing directions, **x** and **y**, aligned along the local directions $\ell_{30}$ and $\ell_{50}$ to stars 3 and 5, respectively (Fig. 1a), while directions $\ell_{40}$ and $\ell_{60}$, to stars 4 and 6, are the corresponding close optical companions. Close stellar (or stellar-like) pairs here mean that we are always concerned with angles $\psi_{i0,j0} \lesssim 0.01$ arcsecond, a number representative of the operational resolution limit reached with telescopes already operating (e.g. HST, ESA's Euclid) or that will soon operate in space (like, e.g., NASA's NextGen Space Telescope, CNSA's CSST or see further below for more) at optical wavelengths, i.e. $\ge 550$ nm, including the near-IR (to ~2 micron).

In the “Transverse and Traceless” (TT) standard gauge the GW components are $h_{0i} = 0$, $\delta^{ij}h_{ij} = 0$, $\delta^{ij}h_{jk,i} = 0$, with only two independent degrees of freedom, the two amplitudes $A_+$ and $A_\times$. Taking advantage of the property that in the TT gauge only the components perpendicular to the direction of propagation survive, we proceed to show, without losing too much to generality, the principle of astrometrically measuring a GW reaching a space-born observer from within the SS, including its direction.

We utilize representative cases of GW strains known from the literature. The first column of Table 1 presents the amplitude, $A_+$, of metric perturbations from possible GW sources as described in Refs.²⁷,²⁸. The last three columns provide, for each $A_+$, the maximum (angular) perturbations resulting on the local unperturbed (angular)
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Schematic illustration of a 3-LOS telescope endowed with a beam combination module on a common focal plane. Differently from Fig. 1a the GW direction Z is 70° from the z axis of the (x,y,z) triad. The output contemplates the case reported in the fourth row of Table 1 with $\omega \sim 10^{-3}$Hz and $A \sim 10^{-19}$Hz for $\psi_{i,j} \sim 0.01''$. The minimum signal (the blue dots) points towards the GW direction.</figcaption></figure>

separations of the pairs (the “arms” of the astrometric antenna) along the three telescope axes (Fig. 1a). The perturbation signals can be calculated from Eq. (12) for a plane GW traveling along the z direction and $\psi_{i_0j_0}$ set to the values in Table 1. Gravitational strains as large as $h \sim 10^{-18}$ are associated with SN core collapse events, thus, they are short lived and at frequencies $\omega \sim 10^3$ Hz. To the range of high frequency sources, $\omega=[10-1,000]$ Hz, belong also the cases of coalescing compact binary systems (NS-NS, pairs of stellar black holes, BH\*-BH\*, or NS-BH\*). Besides, at $\omega < 1$ Hz and in a range of characteristic strain amplitudes spanning 3 orders of magnitude, from $\sim 10^{-18}$ to $10^{-21}$, one finds not only coalescing super-massive BH's, but also the significantly more numerous population of resolved and unresolved Milky Way binaries, with at least one degenerate companion. The fourth row in Table 1 reports, as an example, the output of a linearized GW signal produced by an hypothetical binary system, at the distance of 100pc, with masses $m_1 = 20M_\odot$ and $m_2 = 15M_\odot$, an orbital separation of $1 R_\odot$, and frequency $\omega \sim 10^{-3}$Hz.

Figure 1b shows 5 ms (i.e. 5 times the simulated period) of the angular perturbation $\delta\psi_{i,j}^{GW}$ experienced by the three “arms” (the angles $\psi_{i_0j_0}$) under the strain of the high frequency-high amplitude case mentioned before. As expected, with a GW propagating in the direction of the z-axis, the 1–2 stellar pair is practically unperturbed, while all the action is with the pairs along the x and y axes.

If statistically meaningful measurements of $\sim 5\mu\text{as}$ variations in angular separations, as in Fig. 1b, although possible, are yet to be proven, variations 10 times larger are already within the capabilities of present-day, or soon-to-be, space astrometry missions, like those mentioned earlier (see Ref. 23 and references therein). This is possible also because actual (angular) resolution limits, governing the amplifications of the perturbation signal associated with a GW (Eq. 12), can be made significantly higher by utilizing optimal calibrations procedures$^{29–31}$. Signal frequencies of 1 kHz, as in Fig. 1b are hard to achieve with high signal-to-noise ratio, as sampling above the Nyquist-Shannon frequency would call for integration times shorter than 0.5 ms. “Fast” exposure times are at 25 ms in the case of the astrometer FGS aboard HST$^{32}$ (and references therein), and only 100 ms on Euclid$^{33}$ (and references therein), clearly not enough. And, simply allowing for much shorter exposure times might not improve things as, given the size of the telescope mirrors on these missions, it might become difficult to reach sufficiently faint magnitudes, thus forcing the 3-LOS telescope to limit its orientations to directions with sufficiently bright stellar pairs.

Besides, as mentioned earlier, GW sources with similar amplitudes, or even larger, are expected at more comfortable frequencies, making this scenario a clear case for GW science with an astrometric antenna.

Figure 1c illustrates the progress of $\delta\psi_{i,j}^{GW}$ with time for the GW strain amplitude in row 2 of Table 1; $\omega = 0.1$ Hz, while the separations of the angular arms are here set to 0.001", instead of 0.01", thus generating variations of a few $\mu\text{as}$. Therefore, this figure and the considerations on the cases presented in Table 1 confirm that an astrometric antenna, capable of monitoring periodic signals of amplitude $> 1 \mu\text{as}$, would ideally possess the ability of measuring GW's associated with a range of coalescing massive BH's events.

Finally, the last two rows of Table 1 refers to smaller amplitudes as expected for GW's from core collapse events (at $\omega \sim 100$ Hz), or, at much shorter frequencies, from resolved and unresolved MW binaries. These cases appear beyond today's technology on current missions, including payloads that will be flying into orbit within the next few years. Therefore, further expanding the access to the physics of GW's for an astrometric antenna would require dedicated technological developments (see e.g., Refs. 23,34).
<table><caption>**Table 1.** Angular perturbations, max(δψ<sup>GW</sup><sub>i,j</sub>) (μas), for different linear strains of amplitude A<sub>+</sub> (radians) propagating along z. The values derive from Eq. (12) after setting the three antenna arms ψ<sub>i0,j0</sub> to the values in the second column (in arcsec). As expected, decreasing ψ<sub>i0,j0</sub> by one order of magnitude, the corresponding δψ<sup>GW</sup><sub>i,j</sub> is amplified by the same amount. This is the case for the values reported in the second row compared to its direct analogue in first. Note that along the GW direction of propagation δψ<sup>GW</sup><sub>1,2</sub> is not null, as the direction to star 2, forming the 1-2 pair, although quite small, is not coincident with the LOS to star 1.</caption><thead><tr><th>A<sub>+</sub> (radians)</th><th>ψ<sub>i0,i0</sub> (arcsec)</th><th>max(δψ<sup>GW</sup><sub>1,2</sub>) (μas)</th><th>max(δψ<sup>GW</sup><sub>3,4</sub>) (μas)</th><th>max(δψ<sup>GW</sup><sub>5,6</sub>) (μas)</th></tr></thead><tbody><tr><td>10<sup>-18</sup></td><td>0.01</td><td>4.85 × 10<sup>-15</sup></td><td>4.85</td><td>4.85</td></tr><tr><td>10<sup>-18</sup></td><td>0.001</td><td>4.85 × 10<sup>-16</sup></td><td>48.57</td><td>48.57</td></tr><tr><td>10<sup>-19</sup></td><td>0.01</td><td>4.85 × 10<sup>-16</sup></td><td>0.49</td><td>0.49</td></tr><tr><td>6 × 10<sup>-19</sup></td><td>0.01</td><td>3.02 × 10<sup>-15</sup></td><td>2.57</td><td>2.57</td></tr><tr><td>10<sup>-20</sup></td><td>0.01</td><td>4.85×10<sup>-17</sup></td><td>4.85 × 10<sup>-2</sup></td><td>4.85 × 10<sup>-2</sup></td></tr><tr><td>10<sup>-21</sup></td><td>0.01</td><td>0</td><td>4.85 × 10<sup>-3</sup></td><td>4.85 × 10<sup>-3</sup></td></tr></tbody></table>

The priority of this article was to prove the concept of a novel idea that might develop into new space-borne astrometric instrumentation for the study of GW's. From this perspective, it is evident that the angular arms of a 3-LOS GW antenna, as schematically depicted in Fig. 1a, would register perfectly correlated astrometric signals with amplitudes and phases depending on actual orientations of incoming GW's. For the ideal examples in Fig. 1b,c, these data would immediately tell the direction of arrival of the gravitation strain causing the coherent "trembling" in the antenna arms (the chosen stellar pairs) as that of minimum signal. Therefore, the ultimate precision with which a direction can be recovered is sub-mas, independently from when an actual implementation of our astrometric antenna concept might be sent into orbit, as the sources will be in the Gaia Catalogs (see, e.g., Ref. 35), and therefore their absolute coordinates known to high accuracy. Or, "worst case", one could say the actual directional uncertainty is comparable to the angles separating the stellar pairs used as arms, i.e., ~10 mas or less. Either case, these are unprecedented numbers and pinpointing source directions to such accuracies would tremendously help multi-wavelength and multi-messenger follow-up investigations.

As for the stability of the LOS's (often referred to as telescope pointing stability) during elementary photon integration times (individual exposures), this enters the error budget at the level of the elementary angular separations measurements of the pairs, each materializing an astrometric antenna arm. Assuming the pointing stability as a random process, its effect adds in quadrature to the angular size of an (aberration free) optical diffraction pattern of the point-like sources forming each pair. Pointing stabilities of 1 mas can be routinely achieved on modern space observatories (see, e.g., 30, and references therein, for the HST), and this must be compared to diffractions patterns of ~60 mas for 2m circular apertures at λ~600 nm, therefore contributing a negligible effect. Besides, even in less ideal conditions, the differential nature of our fundamental measurements is such that pointing errors are minimized as common mode effect on angular separations, especially at the resolutions of interest here.

## Discussion

Immediately after the initial development at the idea of an antenna for measuring GW's using astrometry from a telescope in space (see sec. 6.4 in Ref. 36), we began efforts to prove, with simulations and laboratory tests, the feasibility of an astrometric antenna based on the precepts described in this article (Refs. 24,25,37 and references therein). We refer to those studies for results and discussions on what can be currently said on implementation issues like: (1) the very possibility to build a 3-LOS multiplexing telescope (Ref. 38 and references therein), (2) the limit of centering accuracies of star-like images on digital detectors, (3) actual (beyond Rayleigh's) resolution limits for the antenna arms (depending not only on optics and detection system, but also on magnitude and color of the stellar pairs), (4) other natural (intrinsic or cosmic) causes of astrometric noise as stellar activity, and (5) identify (via spatial laser metrology of critical degrees-of-freedom) and deal with instrumental noise mimicking unwanted variations of the antenna arms ψ<sub>i0,j0</sub>.

To date, before any attempt at extracting GW signals from Gaia-like astrometric data, one would generally consider the end of the global reduction process in order to obtain the best possible knowledge of satellite attitude (orientation) and instrumental behaviour at all of relevant time scales. In such a context, looking for variations in the direction to a single source on the sky at the nas (~5 × 10<sup>-15</sup> rad) level implies the knowledge of an "absolute" reference, e.g. the telescope LOS, at comparable precision ~10<sup>-15</sup>. This is an impossible requirement on the reconstructed attitude of a science satellite. The differential technique proposed in Refs. 24,25 aims at nas measurement over a distance of the order of 1", thus implicitly reducing the relative precision requirement to ~10<sup>-9</sup>, with an improvement of six orders of magnitude. Also, instrument calibration requirements are strongly alleviated. In Gaia, we have a variation of the electro-optical response of hundreds of mas over the 0.5° field, calibrated to the μas level.

Assuming a linear model, the corresponding electro-optical response variation for an astrometric GW telescope over 1" would be in the range of hundreds of μas, i.e., a comparable calibration "power" would scale the measurement reliability to the nas regime. Actually, optimal optics design, exploiting the higher-than-linear decrease of many aberrations close to the optical axis, would reduce instrumental contribution even further 37.

In addition to the above, we will have to simulate much more realistic scenarios (i.e., more general forms of GW's and use real-sky pairs) and conditions (realistic noise levels) to investigate viable strategies for the actual
retrieval of amplitude and phase (carrying the direction information, Figs. 1b,c and 2). However, especially for these aspects, help would certainly come from the large amount of work done, and proven on real data by the LIGO and VIRGO collaborations.

Nevertheless, as we show in this work, the astrometric observable could amplify the GW-induced signal if one takes into account the angle between two space-like directions of light³⁶ in the framework of general relativistic astrometry (see e.g. Ref.²²). In such a case, one is in principle exempted from dealing with satellite's attitude and the GW astrometric measurement can be translated into an observation equation accounting for a wide range of frequencies.

The diversity of GW frequencies that can potentially be treated with Eq. (12) and appropriate modeling of gravitational shifts could improve PTA and/or LISA observations⁸,³⁹, and the low frequency domain due to periodic sources (e.g., Galactic binary WDs identified by Gaia), thus further enhancing the mapping of the Milky Way substructures⁴⁰. Moreover, the advantage of Eq. (9) is the possibility to exploit a large number of null geodesics, so to better scrutinize the GW direction, a critical aspect of the GW detection and multiwavelength characterization. This same feature would also enable tests on GW polarization modes by combining different telescope orientations.

Finally, the same principle of a space-born multi-LOS telescope discussed here can push investigations on possible ground-based realisations of an astrometric GW antenna. Such a development could support present and future linear interferometer by working as signal sentinel and by spotting almost immediately directions to incoming signals, thanks to the fully correlated nature of the GW-induced astrometric signals expected from the combined LOS's (Figs. 1b,c, and 2).

In conclusion, the potential of the relativistic astrometric observable advocated here for the astrometric detection and precise identification of gravitational waves, by using pairs of natural stars, will significantly add to the best GW detection procedures in use. In fact, the availabiltiy of close stellar pairs opens to the exploitation of a large number of configurations, and therefore to monitoring gravitational waves coming from any direction; thus providing, at the same time, extensive statistics to uncover the properties of a GW source. This helps in a truly complementary and independent way all of the efforts dedicated to multiband GW searches bridging low and high frequencies at different redshifts (see, for example, Refs.⁴¹⁻⁴³). And, in case the two unperturbed LOS's of a single astrometric antenna arm, although angularly very close, are actually related to two stars at different distances, it would be also possible to investigate time retarded effects and test GW speed. Last, a suitable choice of the strain $h^{GW}$ could pave the way for new GW tests on gravitation interaction with photons: it would suffice modeling the $F^{GW}$ function appropriately.

# Methods

## The cosine expression for the GW observation equation

The cosine of two light directions (i,j) is defined in General Relativity as:

$$ \cos \psi_{i,j} = \frac{P(u)_{\alpha\beta} k_i^{\alpha} k_j^{\beta}}{(\sqrt{P(u)_{\alpha\beta} k_i^{\alpha} k_j^{\beta}})(\sqrt{P(u)_{\alpha\beta} k_j^{\alpha} k_i^{\beta}})} \quad (13) $$

where $P(u)_{\alpha\beta} = g_{\alpha\beta} + u_{\alpha}u_{\beta}$ is the operator that projects with respect to the local barycentric observer $u^{\alpha} = 1/\sqrt{-g_{00}}$. Then the photon 4-momentum can be decomposed as

$$ k^{\alpha} = -(u|k)u^{\alpha} + l^{\alpha}, \quad (14) $$

where $(u|k) = g_{\alpha\beta}u^{\alpha}k^{\beta}$ and $l^{\alpha}$ is the spatial null vector projected on the rest space of $u^{\alpha}$, Defining

$$ \bar{k}^{\alpha} = -\frac{k^{\alpha}}{(u|k)}, \quad l^{\alpha} = -\frac{l^{\alpha}}{(u|k)} = \bar{k}^{\alpha} - u^{\alpha}. \quad (15) $$

the cosine for two directions i,j simplifies as

$$ \cos \psi_{i,j} = g_{\alpha\beta} (\ell_i^{\alpha} \ell_j^{\beta})_{\text{obs}}. \quad (16) $$

The BCRS metric is defined by IAU resolutions as a post-Newtonian (pN) solution of the Einstein field equations. Thus, one has to take into account terms of this metric accurate to the order of the GW perturbations sought for. Dropping the sum symbol, let us express the metric including both sources as:

$$ g_{\alpha\beta} = \eta_{\alpha\beta} + \epsilon h_{(1)\alpha\beta}^{SS} + \epsilon^2 h_{(2)\alpha\beta}^{SS} + \epsilon^3 h_{(3)\alpha\beta}^{SS} + \epsilon^4 h_{(4)\alpha\beta}^{SS} + h_{\alpha\beta}^{GW} + O(\epsilon^5) \quad (17) $$

where $\epsilon$ is of the order of $v/c$, being $v$ the typical velocity of each relevant SS metric source, and the subscripts in parenthesis indicate the order of approximation in $\epsilon$. We assume that the GW perturbations are of order $\epsilon^4$ at best, i.e., at the nanoarcsecond level. Like for the SS metric, the SS contribution to light direction can be approximated as:

$$ \ell_{\text{obs}}^{\alpha(\text{SS})} = \ell_0^{\alpha} + \epsilon \ell_{(1)}^{\alpha} + \epsilon^2 \ell_{(2)}^{\alpha} + \epsilon^3 \ell_{(3)}^{\alpha} + \epsilon^4 \ell_{(4)}^{\alpha} + O(\epsilon^5). \quad (18) $$

From the assumptions above, one finally finds:
$$ \cos \psi_{i,j} = \cos \psi_{i,j}^{SS} + \eta_{\alpha\beta}(\ell_{i0}^{\alpha} \delta \ell_{j}^{\beta} + \ell_{j0}^{\alpha} \delta \ell_{i}^{\beta})_{\text{obs}} + h_{\alpha\beta}^{\text{GW}} \ell_{i0}^{\alpha} \ell_{j0}^{\beta} + O(\epsilon^5). \quad (19) $$

## The total cosine versus an additive GW direction cosine

In this section we clarify the consequences, when in the presence of a passing GW, of utilizing the cosine as observable by simply extending what it is done in the context of the Gaia mission, where the direction cosine refers to the angle of the incoming light to the observer (satellite) attitude-tetrad $E_{\hat{a}}^{\alpha}$.

Let us consider a priori assumption that the effect of a GW is that of adding a cosine term to the direction cosine associated with the SS metric, i.e., denoting with $\cos(\hat{\psi})$ the cosine of the angle $\hat{\psi}$ of a light direction to the tetrad, the total (tot) cosine is given by:

$$ \cos(\hat{\psi})_{\text{tot}} = \cos(\hat{\psi})^{\text{SS}} + \cos(\hat{\psi})^{\text{GW}}. \quad (20) $$

With the extra assumption that the correction $\delta\ell_{(GW)}$, induced by the GW to the SS line-of-sight at the observer $\ell_{(SS)}$, is (formally) known, Eq. (20) would read:

$$ (\eta_{\mu\nu} + h_{\mu\nu}^{\text{SS}})(\ell_{(SS)}^{\mu} E_{\hat{a}}^{\nu}) + (\eta_{\mu\nu} + h_{\mu\nu}^{\text{GW}})(\delta\ell_{(GW)}^{\mu} E_{\hat{a}}^{\nu}) \approx \cos(\hat{\psi})^{\text{SS}} + \eta_{\mu\nu} \delta\ell_{(GW)}^{\mu} E_{\hat{a}}^{\nu}, \quad (21) $$

implying no explicit dependence on the strain $h_{GW}$.

If, on the other hand, the GW perturbation is not a known part of the observed direction at the observer, one would have:

$$ \begin{aligned} & (\eta_{\mu\nu} + h_{\mu\nu}^{\text{SS}})[(\ell^{\mu} + \delta\ell^{\mu})E_{\hat{a}}^{\nu}] + (\eta_{\mu\nu} + h_{\mu\nu}^{\text{GW}})[(\ell^{\mu} + \delta\ell^{\mu})E_{\hat{a}}^{\nu}] \\ & \approx \cos(\hat{\psi})^{\text{SS}} + 2\eta_{\mu\nu}\delta\ell^{\mu}E_{\hat{a}}^{\nu} + \eta_{\mu\nu}\ell^{\mu}E_{\hat{a}}^{\nu} + h_{\mu\nu}^{\text{GW}}\ell_{0}^{\mu}E_{\hat{a}}^{\nu}, \end{aligned} \quad (22) $$

with the flat Minkowskian metric term entering twice in the Equation for the $\ell^{\alpha}$ component. Then, the simple addition to the “background” direction cosine in the SS metric (i.e., BCRS) would introduce the flat Minkowskian contribution twice without a priori disentanglement from the local-line-of-sight of the GW component or, if the GW shift is considered as a separate part, it would imply to discard a priori the $h^{GW}$ strain in the observation equation.

It is only when we drop the assumption made with Eq. (20), i.e. when working directly with the total cosine, that we finally recover an expression similar to observation equation in the Main article as:

$$ \cos(\hat{\psi})_{\text{tot}} = (\eta_{\mu\nu} + h_{\mu\nu}^{\text{SS}} + h_{\mu\nu}^{\text{GW}})[(\ell^{\mu} + \delta\ell^{\mu})E_{\hat{a}}^{\nu}] \approx \cos(\hat{\psi})^{\text{SS}} + \eta_{\mu\nu}\delta\ell^{\mu}E_{\hat{a}}^{\nu} + h_{\mu\nu}^{\text{GW}}\ell^{\mu}E_{\hat{a}}^{\nu}. \quad (23) $$

However, this form of the observation equation explicitly depends on the observer orientation, i.e., on satellite attitude, with the drawbacks in relation to attitude errors heavily affecting the error budget of the measurements.

## The astrometric gravitational wave shift

The four tangent vector to a null geodesic satisfies the well known conditions:

$$ k^{\alpha} \nabla_{\alpha} k^{\beta} = 0, \quad k^{\alpha} k_{\alpha} = 0, \quad (24) $$

$\nabla_{\alpha}$ being the covariant derivative associated with the spacetime metric. The decomposition of the photon 4-momentum with respect to an observer $u^{\alpha}$ implies that the trajectory is parametrized by $\sigma$ such that

$$ \bar{k}^{\alpha} = -\frac{k^{\alpha}}{(u|k)} = \frac{dx^{\alpha}}{d\sigma}, \quad (25) $$

and Eq. (24) becomes

$$ \bar{k}^{\alpha} \nabla_{\alpha} \bar{k}^{\beta} = - \frac{d \ln [-(u|k)]}{d\sigma} \bar{k}^{\beta}, \quad (26) $$

which is related to the affine parameter $\lambda$ by $d\sigma = -(u|k)d\lambda$.

It is easy to check that in the case of a static observer

$$ \frac{d \ln [-(u|k)]}{d\sigma} = \ell^{\alpha} \ell^{\beta} k_{\alpha\beta} - \ell^{\alpha} a_{\alpha} = -\ell^{\alpha} \ell^{\beta} \theta_{\alpha\beta} - \ell^{\alpha} a_{\alpha}, \quad (27) $$

where the two spatial fields coming from the splitting of the covariant derivative of $u$, i.e., $\nabla_{\beta} u^{\alpha} = -a^{\alpha} u_{\beta} - k^{\alpha} \beta$, the acceleration vector $a^{\alpha}$ and the kinematical tensor $k^{\alpha}_{\beta} = \omega^{\alpha}_{\beta} - \theta^{\alpha}_{\beta}$ are expressed as a combination of the vorticity and expansions of the congruence of curves related to fiducial observers $u^{\alpha}$. Thus, the geodesic equation transforms into

$$ \frac{d\bar{k}^{\alpha}}{d\sigma} + \Gamma^{\alpha}_{\mu\nu} \bar{k}^{\mu} \bar{k}^{\nu} - [\ell^{\mu} \ell^{\nu} \theta_{\mu\nu} + \ell^{\mu} a_{\mu}] \bar{k}^{\alpha} = 0, \quad (28) $$

or
$$ \frac{d\ell^\alpha}{d\sigma} + \Gamma^\alpha_{\mu\nu} \ell^\mu (\ell^\nu + u^\nu) + a^\alpha - k^\alpha_\sigma \ell^\sigma - [\ell^\mu \ell^\nu \theta_{\mu\nu} + \ell^\mu a_\mu] (\ell^\alpha + u^\alpha) = 0, \quad (29) $$

for the unknown local line-of-sight $\ell^\alpha$. In case of static observers, the kinematical fields reduce to

$$ a^i = \partial_0 h_{0i} - \frac{1}{2} \partial_i h_{00} = h_{0i,0} - \frac{1}{2} h_{00,i}, \quad \theta_{ij} = \frac{1}{2} h_{ij,0}, \quad \omega_{ij} = -h_{0[i,j]}. \quad (30) $$

Note that considering the metric as an approximate solution of the Einstein field equation composed of a background part plus a GW perturbation, also the affine coefficients can be split respectively into two parts

$$ \Gamma^\alpha_{\mu\nu} = \frac{1}{2} g^{\alpha\rho} (g_{\rho\beta,\gamma} + g_{\rho\gamma,\beta} - g_{\beta\gamma,\rho}) \approx \Gamma^{\alpha(\text{SS})}_{\beta\gamma} + \Gamma^{\alpha(\text{GW})}_{\beta\gamma}, \quad (31) $$

at the first order of the perturbation. The same for the parameter sigma

$$ d\sigma = -(u|k)d\lambda = -[(g^{(\text{SS})}_{\alpha\beta} + h^{(\text{GW})}_{\alpha\beta}) u^\alpha k^\beta]d\lambda. \quad (32) $$

The TT gauge choice implies that the second term of equation (32) does not contribute.

All of the above implies the possibility again to split Eq. (29) into the Solar System part plus the GW one, thus it allows to integrate separately each term. As a matter of fact, at the order of $\epsilon^4$, it is possible to isolate from Eq. (29) the contribution of the GW part and obtain

$$ \frac{d\delta\ell^i}{d\sigma} \approx -\frac{1}{2} \ell_0^j \ell_0^k (2h_{ij,k}^{\text{GW}} - h_{jk,i}^{\text{GW}}) - \ell_0^j h_{ij,0}^{\text{GW}} + \frac{\ell_0^j \ell_0^k \ell_0^i}{2} h_{jk,0}^{\text{GW}}. \quad (33) $$

Since $d\sigma = d\lambda + O(\epsilon^2)$ and assuming for the photon trajectory $x^0(\sigma) = x^0_{\text{obs}} + \sigma + O(\epsilon^2)$ and $x^i(\sigma) = x^i_{\text{obs}} + \ell_0^i \sigma + O(\epsilon^2)$, the argument of $h^{\text{GW}}$ becomes:

$$ \tilde{k}_{\alpha} x^{\alpha} = \tilde{k}_{0} \sigma (1 + p \cdot \ell_{0}) + \tilde{\psi} + O(\epsilon^{2}) \quad (34) $$

where $\tilde{\psi} = \tilde{k}_\alpha x^\alpha_{\text{obs}}$ can be considered a phase term. Then,

$$ h_{ij,0}^{\text{GW}} = \frac{1}{(1 + p \cdot \ell_0)} \frac{dh^{\text{GW}}}{d\sigma}, \quad h_{ij,k}^{\text{GW}} = \frac{p^k}{(1 + p \cdot \ell_0)} \frac{dh^{\text{GW}}}{d\sigma}. \quad (35) $$

Via a direct integration of Eq. (33) we easily obtain the gravitational shift of the local direction:

$$ \int_{\sigma_{\text{obs}}}^{\sigma_{*}} d\delta\ell^i = \frac{\ell_0^i + p^i}{2(1 + p \cdot \ell_0)} \ell_0^j \ell_0^k \int_{\sigma_{\text{obs}}}^{\sigma_{*}} dh_{jk}^{(\text{GW})}(\sigma) - \frac{\ell_0^j}{2} \int_{\sigma_{\text{obs}}}^{\sigma_{*}} dh_{ij}^{(\text{GW})}(\sigma) + O(\epsilon^5), \quad (36) $$

which coincides with the result in Ref.¹⁴ when the distance to the stellar source is many gravitational waves away, namely the detection occurs in the far-away wave zone.

Any variation induced by the GW on the orthonormal basis is in principle absorbed as second order effects, or in the calculations if proportional to the GW strain. In conclusion, the derivation of Eqs. (9) and (12) were obtained by solving the geodesic for the null four vector with respect to its affine parameter and taking into account the kinematical geometrical proprieties of the congruence of the fiducial observers and their variations, as expressions presented here show, consistently with the GR theory of measurement and the results of the cited Authors.

## Data availability

All data generated or analysed during this study are included in this published article.

Received: 9 September 2023; Accepted: 26 February 2024

Published online: 01 March 2024

## References

1. Abbott, B. P. et al. Observation of gravitational waves from a binary black hole merger. Phys. Rev. Lett. 116, 061102. https://doi.org/10.1103/PhysRevLett.116.061102 (2016).

2. LIGO Scientific Collaboration et al. Advanced LIGO. Classical Quant. Gravity. 32, 074001, https://doi.org/10.1088/0264-9381/32/7/074001 (2015).

3. Acernese, F. et al. Increasing the astrophysical reach of the advanced virgo detector via the application of squeezed vacuum states of light. Phys. Rev. Lett. 123, 231108. https://doi.org/10.1103/PhysRevLett.123.231108 (2019).

4. Collaboration Kagra. et al. KAGRA: 2.5 generation interferometric gravitational wave detector. Nat. Astron. 3, 35–40. https://doi.org/10.1038/s41550-018-0658-y (2019).

5. Punturo, M. et al. The Einstein Telescope: A third-generation gravitational wave observatory. Classical Quant. Gravity 27, 194002. https://doi.org/10.1088/0264-9381/27/19/194002 (2010).

6. Abbott, B. P. et al. Exploring the sensitivity of next generation gravitational wave detectors. Classical Quant. Gravity 34, 044001. https://doi.org/10.1088/1361-6382/aa51f4 (2017) arXiv:1607.08697.

7. Kawamura, S. et al. Current status of space gravitational wave antenna DECIGO and B-DECIGO. Progress Theor. Exp. Phys. 2021, 05A105. https://doi.org/10.1093/ptep/ptab019 (2021).

8. Amaro-Seoane, P. et al. Laser Interferometer Space Antenna. arXiv e-prints. https://doi.org/10.48550/arXiv.1702.00786 (2017).
9. Bailes, M. et al. Gravitational-wave physics and astronomy in the 2020s and 2030s. *Nat. Rev. Phys.* **3**, 344–366. https://doi.org/10.1038/s42254-021-00303-8 (2021).

10. Braginsky, V. B., Kardashev, N. S., Polnarev, A. G. & Novikov, I. D. Propagation of electromagnetic radiation in a random field of gravitational waves and space radio interferometry. *Nuovo Cimento B Serie 105*, 1141–1158. https://doi.org/10.1007/BF02827323 (1990).

11. Gwinn, C. R., Eubanks, T. M., Pyne, T., Birkinshaw, M. & Matsakis, D. N. Quasar proper motions and low-frequency gravitational waves. *Astrophys. J.* **485**, 87–91. https://doi.org/10.1086/304424 (1997) arXiv:astro-ph/9610086.

12. Kaiser, N. & Jaffe, A. Bending of light by gravity waves. *Astrophys. J.* **484**, 545. https://doi.org/10.1086/304357 (1997).

13. Gwinn, C. R., Eubanks, T. M., Pyne, T., Birkinshaw, M. & Matsakis, D. N. Quasar proper motions and low-frequency gravitational waves. *Astrophys. J.* **485**, 87–91. https://doi.org/10.1086/304424 (1997) arXiv:astro-ph/9610086.

14. Book, L. G. & Flanagan, E. E. Astrometric effects of a stochastic gravitational wave background. *Phys. Rev. D* **83**, 024024. https://doi.org/10.1103/PhysRevD.83.024024 (2011).

15. Mihaylov, D. P., Moore, C. J., Gair, J. R., Lasenby, A. & Gilmore, G. Astrometric effects of gravitational wave backgrounds with non-Einsteinian polarizations. *Phys. Rev. D* **97**, 124058. https://doi.org/10.1103/PhysRevD.97.124058 (2018) arXiv:1804.00660.

16. Damour, T. & Esposito-Farèse, G. Gravitational-wave versus binary-pulsar tests of strong-field gravity. *Phys. Rev. D* **58**, 042001. https://doi.org/10.1103/PhysRevD.58.042001 (1998).

17. Crosta, M., Lattanzi, M. G. & Spagna, A. Prospects for astrometric measurements of gravitational waves from stellar sources with GAIA. *Baltic Astron.* **8**, 239–251. https://doi.org/10.1515/astro-1999-0209 (1999).

18. Schutz, B. F. Astrometric and timing effects of gravitational waves. in (Klioner, S. A., Seidelmann, P. K. & Soffel, M. H. Eds.) *Relativity in Fundamental Astronomy: Dynamics, Reference Frames, and Data Analysis*, vol. 261, 234–239. https://doi.org/10.1017/S1743921309990457 (2010).

19. Moore, C. J., Mihaylov, D. P., Lasenby, A. & Gilmore, G. Astrometric search method for individually resolvable gravitational wave sources with gaia. *Phys. Rev. Lett.* **119**, 261102. https://doi.org/10.1103/PhysRevLett.119.261102 (2017).

20. Klioner, S. A. Gaia-like astrometry and gravitational waves. *Classical Quant. Gravity* **35**, 045005. https://doi.org/10.1088/1361-6382/aa9f57 (2018).

21. Mignard, F. & Klioner, S. Analysis of astrometric catalogues with vector spherical harmonics. *Astron. Astrophys.* **547**, A59. https://doi.org/10.1051/0004-6361/201219927 (2012) arXiv:1207.0025.

22. Crosta, M., Geralico, A., Lattanzi, M. G. & Vecchiato, A. General relativistic observable for gravitational astrometry in the context of the gaia mission and beyond. *Phys. Rev. D* **96**, 104030. https://doi.org/10.1103/PhysRevD.96.104030 (2017).

23. Lattanzi, M. G. Astrometric cosmology. *Memorie SAIt* **83**, 1033 (2012).

24. Gai, M. et al. The astra project: A doorway to future astrometry. in (Geyl, R. & Navarro, R. Eds.) *Advances in Optical and Mechanical Technologies for Telescopes and Instrumentation IV*. https://doi.org/10.1117/12.2562671 (SPIE, 2020).

25. Gai, M. et al. Progress on the astrometric gravitation probe design. in (Lystrup, M., Batalha, N., Tong, E. C., Siegler, N. & Perrin, M. D. Eds.) *Space Telescopes and Instrumentation 2020: Optical, Infrared, and Millimeter Wave*. https://doi.org/10.1117/12.2562644 (SPIE, 2020).

26. Poisson, E. & Will, C. M. *Frontmatter*, i–iv (Cambridge University Press, 2014).

27. Thorne, K. S. Gravitational radiation. in (Hawking, S. W. & Israel, W. Eds.) *Three Hundred Years of Gravitation*, 330 (Cambridge University Press, 1989).

28. Colpi, M. & Sesana, A. *Gravitational Wave Sources in the Era of Multi-Band Gravitational Wave Astronomy*, 43–140 (WORLD SCIENTIFIC, 2017).

29. Bernacca, P. L. et al. Hubble space telescope astrometric observations of pre-main sequence stars from the HIPPARCOS program. *Astron. Astrophys.* **278**, L47–L50 (1993).

30. Bernacca, P. L., Lattanzi, M. G., Porro, I., Neuhaeuser, R. & Bucciarelli, B. Hubble space telescope astrometric observations of pre-main sequence stars from the HIPPARCOS Program. II. *Astron. Astrophys.* **299**, 933 (1995).

31. Steffen, A. et al. A dynamical mass and distance for a pms star: The binary ntt 045251+3016. *Proc. Int. Astron. Union* **200**, 19P (2000).

32. Bucciarelli, B., Holfeltz, S. T., Lattanzi, M. G., Taff, L. G. & Vener-Saavedra, P. C. The photometric properties of the hst astrometer fine guidance sensor. *Publ. Astron. Soc. Pacific* **106**, 417. https://doi.org/10.1086/133394 (1994).

33. Bosco, A., Bacchetta, A., Saponara, M. & Criado, G. S. *Euclid pointing performance: operations for the Fine Guidance Sensor reference star catalogue*, AIAA 2018–2415 (American Institute of Aeronautics and Astronautics, 2018). https://doi.org/10.2514/6.2018-2415.

34. Riva, A. et al. LEGOLAS: localizing evidence of gravitational waves by observations of light source astrometric signature. in (Jr., J. M. O., Clampin, M. C. & MacEwen, H. A. Eds.) *Space Telescopes and Instrumentation 2010: Optical, Infrared, and Millimeter Wave*, vol. 7731, 77311T. https://doi.org/10.1117/12.857161. International Society for Optics and Photonics (SPIE, 2010).

35. Lindegren, L. et al. Gaia Early Data Release 3. The astrometric solution. *Astron. Astrophys.* **649**, 2. https://doi.org/10.1051/0004-6361/202039709 (2021) arXiv:2012.03380.

36. Crosta, M. Astrometry in the 21st century. From Hipparchus to Einstein. *Nuovo Cimento Rivista Serie 42*, 443–510. https://doi.org/10.1393/ncr/i2019-10164-2 (2019).

37. Riva, A. et al. Rafter: Ring astrometric field telescope for exo-planets and relativity. in (Lystrup, M., Batalha, N., Tong, E. C., Siegler, N. & Perrin, M. D. Eds.) *Space Telescopes and Instrumentation 2020: Optical, Infrared, and Millimeter Wave*, https://doi.org/10.1117/12.2576806 (SPIE, 2020).

38. Sun, L.-Z. et al. A new telescope with three fields of view to measure the orientation parameters of the Moon and terrestrial planets. *Res. Astron. Astrophys.* **21**, 040. https://doi.org/10.1088/1674-4527/21/2/40 (2021).

39. Verbiest, J. P. W. et al. The international pulsar timing array: First data release. *Monthly Notices R. Astron. Soc.* **458**, 1267–1288. https://doi.org/10.1093/mnras/stw347 (2016) arXiv:1602.03640.

40. Wilhelm, M. J. C., Korol, V., Rossi, E. M. & D'Onghia, E. The milky way's bar structural properties from gravitational waves. *Monthly Notices R. Astron. Soc.* **500**, 4958–4971. https://doi.org/10.1093/mnras/staa3457 (2020).

41. Arca Sedda, M. et al. The missing link in gravitational-wave astronomy: Discoveries waiting in the decihertz range. *Classical Quant. Gravity* **37**, 215011. https://doi.org/10.1088/1361-6382/abb5c1 (2020) arXiv:1908.11375.

42. Aggarwal, N. et al. Challenges and opportunities of gravitational-wave searches at mhz to ghz frequencies. *Living Rev. Relativity*. https://doi.org/10.1007/s41114-021-00032-5 (2021).

43. Burke-Spolaor, S. et al. The astrophysics of nanohertz gravitational waves. *Astron. Astrophys. Rev.* https://doi.org/10.1007/s00159-019-0115-7 (2019).

## Acknowledgements

MTC and MGL acknowledge INAF for the Theory Grant "Gravitational astrometry tools for peering gravitational waves in the fabric of spacetime", INAF RF 2023.
**Author contributions**

M.C. conceived the theoretical observation equations. M.G.L. conceived the measurement procedure. M.C., M.G.L. designed the concept of the astrometric GW antenna and wrote the manuscript. C.L.P., M.G. contributed to the manuscript. All authors reviewed the manuscript and contributed to discussions.

**Competing interests**

The authors declare no competing interests.

**Additional information**

Correspondence and requests for materials should be addressed to M.C.

Reprints and permissions information is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_5.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2024