Probing Dark Matter Clumps, Strings and Domain Walls with Gravitational Wave Detectors

Gravitational wave astronomy has recently emerged as a new way to study our Universe. In this work, we survey the potential of gravitational wave interferometers to detect macroscopic astrophysical objects comprising the dark matter. Starting from the well-known case of clumps we expand to cosmic strings and domain walls. We also consider the sensitivity to measure the dark matter power spectrum on small scales. Our analysis is based on the fact that these objects, when traversing the vicinity of the detector, will exert a gravitational pull on each node of the interferometer, in turn leading to a differential acceleration and corresponding Doppler signal, that can be measured. As a prototypical example of a gravitational wave interferometer, we consider signals induced at LISA. We further extrapolate our results to gravitational wave experiments sensitive in other frequency bands, including ground-based interferometers, such as LIGO, and pulsar timing arrays, e.g. ones based on the Square Kilometer Array. Assuming moderate sensitivity improvements beyond the current designs, clumps, strings and domain walls may be within reach of these experiments.


I. INTRODUCTION
LIGO's measurement of gravitational waves emitted from black hole binary mergers [1,2] has opened novel ways to study the properties of astrophysical objects via their gravitational interactions. After decades-lasting developments and technical improvements on interferometrical methods, an unprecedented sensitivity in the measurement of gravitational interactions has been achieved. Future ground-and space-based interferometers, such as the Cosmic Explorer [3], the Einstein Telescope [4] and LISA [5,6], are expected to further improve on this sensitivity and to complement it in different frequency ranges. In addition, this may also be supplemented by instruments such as AEDGE [7], AION [8], BBO [9] or DE-CIGO [10]. To achieve sensitivity to lower gravitational wave frequencies and smaller strains, longer interferometer arms and less background noise are needed. In the former, pulsar timing arrays, such as NANOGrav [11], PPTA [12] or EPTA [13], currently provide additional tests of gravitational physics. In the future, pulsar timing arrays might even use the Square Kilometer Array (SKA) [14] to further improve on sensitivity.
In this work, we want to explore the potential of gravitational wave interferometers to probe concentrated structures of dark matter via their gravitational interaction with the apparatus. This provides an additional science goal for these instruments, but also requires different analysis strategies. While previous studies [15][16][17][18][19][20][21][22][23][24][25][26][27][28][29] have focused mainly on clumps and primordial black holes, we consider also topological defects, such as cosmic strings and domain walls. The opportunity to search for the latter two has already been noted in [24,30] 1 . Here, 1 In [30] the main focus is on axion streams arising from the de-we aim for a more detailed analysis of the experimental sensitivity and also consider the possibility that the topological defects feature a more general equation of state, e.g. due to a non-trivial behavior of the string/domain wall network. We also take a brief look beyond localized structures and consider the possibility to measure the local dark matter density fluctuation power spectrum (see also [20]). The sensitivity to the latter is, however, somewhat limited.
Our analysis is based on the simplest effect (originally discussed in [16]), namely, that any massive object or, more precisely, any localized energy density will locally perturb the gravitational field in the vicinity of the interferometer, thereby exerting a different gravitational pull on each node of the detector. Due to the differential gravitational acceleration, this perturbation will, in turn, lead to a measurable Doppler shift in the apparatus. For instance, a similar analysis of such acceleration burst signals at space-based interferometers has been discussed for nearby asteroids and similar objects [31][32][33]. In addition to this, a signal could also arise due to the (changing) gravitational potential when a structure exists within the line of site connecting the different nodes of the gravitational wave detector [17,20,23,26,28]. It would be interesting to study this Shapiro effect [34] also for cosmic strings and domain walls. We leave this for future work.
In principle, this type of dark matter search is very general as it solely relies on measuring purely gravitational interactions 2 . Of course, the downside is that it requires struction of axion miniclusters, which are, however, rather similar to strings. 2 Signals of dark matter structures in gravitational wave interferometers can also arise from non-gravitational interactions, see, e.g., [24,[35][36][37][38][39][40][41].
the dark matter to be strongly concentrated in highly localized structures which, in turn, re-introduces a model dependence. While our study is purely phenomenological and does not rely on the origin of the structures in question, let us nevertheless mention a few possibilities allowing for such scenarios. Localized clumps in the dark matter structure can appear in a variety of models and situations. Perhaps the most obvious scenario for clumps are primordial black holes [42][43][44][45][46] (in the context of gravitational wave astronomy, see also more recently, e.g., [47][48][49][50][51][52] and [53] for a recent review). That said, already standard cold dark matter may feature clumpy structures (see, e.g., [20,22]), although these may likely be beyond the reach of near future gravitational wave detectors [20]. In addition, localized clumps can also be produced if dark matter features strong self-interactions [54][55][56][57]. Another important potential source are initial conditions featuring large inhomogeneities, a prominent example of which are axion miniclusters [58][59][60][61] (see also [62][63][64][65][66] for some more recent work). More recently, such large fluctuations have also been discussed in the context of inflationary production mechanisms [67,68] or as a consequence of a fragmentation of homogeneous fields due to their (self-)interactions [69,70]. Furthermore, macroscopic clumps could arise as solitonic objects such as, e.g., Q-balls [71][72][73][74].
In addition to localized clumps, dark matter structures might also come in form of topological defects, such as cosmic strings or domain walls. These could have formed in the early Universe [75][76][77]. In a cosmological context, topological defects, in particular dynamical networks of cosmic strings or domain walls (see, e.g., [78][79][80][81][82][83][84]), contribute to the total energy budget of the Universe. However, while their equation of state is usually negative, there is still uncertainty in the behaviour of a network of these objects (see, e.g., [85,86]). Therefore they may account for dark energy but also for dark matter [82][83][84]. This has led to a variety of ideas to investigate their experimental signatures as dark matter candidates [29,35,38,41,[87][88][89][90][91][92][93]. We will follow this more phenomenological spirit and be agnostic about the equation of state of these networks by treating it as a free parameter. This allows us to investigate their experimental imprint in gravitational wave detectors independent of the dynamics of the network.
Our discussion is structured as follows. Section II first reviews the analysis of acceleration burst signals at the LISA gravitational wave interferometer. As an initial test case, we then apply these techniques to obtain the signal power spectrum associated to localized clumps of dark matter passing by the detector. In Section III, we extend our analysis to the case of cosmic strings and domain walls. Next, Section IV gives an overview of how the same technique might also be used to measure stochastic fluctuations of the local dark matter density with LISA. In Section V, we extrapolate our results to other gravitational wave experiments, i.e. ground-based interfer-ometers and pulsar timing arrays, that are sensitive to different frequency bands. As particular examples, we examine LIGO and a future PTA using the SKA. Finally, we summarize our results and conclude in Section VI.

II. LOCALIZED CLUMPS OF DARK MATTER
Localized dark matter clumps can cause a signal in a gravitational wave interferometer as shown by a number of previous studies [15][16][17][18][19][20][21][22][23][24][25][26][27][28][29]. In this section, we will review how such a signal is generated. We focus on the effect where the clump exerts a stronger gravitational acceleration on one of the interferometer nodes and derive the associated signal power spectrum for the LISA detector. This discussion will then serve as a basis and test bed for our investigation of strings and domain walls in the next section. For a general introduction to the physics and measurement techniques of gravitational wave interferometers see, e.g., [94,95].
Although the general strategy in principle applies to any gravitational wave detector, in the present work, we will consider LISA as a prototypical example to obtain the experimental signature of gravitational perturbations caused by the above macroscopic astrophysical objects in a gravitational wave interferometer. We will later (rather naively) use the same technique to estimate the sensitivity in other gravitational wave experiments. Here, loosely speaking, we focus on an experimental setup with three distinct nodes arranged in an equilateral triangle (see Appendix A for some more details). Along the three interferometer arms with a length of about 2.5 million kilometers, the satellites exchange laser beams. Consequently, a differential acceleration due to a gravitational perturbation caused by macroscopic objects in the vicinity of the interferometer will then lead to a measurable Doppler signal.
To begin with, let us first consider the case of a single dark matter clump passing by one of the LISA satellites. This is similar to the case of the detection of asteroids treated in [31], which we will follow closely. As a simple coordinate frame, we choose the dark matter clump to be in a straight uniform motion with velocity V parallel to the y-axis. The trajectory of the clump is, in addition, confined to the xy-plane. The closest distance between the clump and the satellite located at the origin, i.e. the impact parameter, is denoted by D. For a schematic illustration of this reference frame we refer the reader to Appendix A 1. The differential acceleration of the interferometer node is then given by [31], where G denotes the gravitational constant and M is the mass of the dark matter clump. A schematic form of this acceleration burst was first considered in [16] with regard to the detection of primordial black holes with spacebased interferometers. As [31], in our analysis we will instead consider the velocity shift, v(t), i.e. the integrated gravitational acceleration. However, we will use the more appropriate signal response functions defined in [96,97] that take into account a time retardation of the signals from different nodes (see also discussion below). In principle, the gravitational field of the dark matter clump passing by the detector will exert a gravitational pull on each interferometer node separately. Due to the differential velocity shifts relative to each other, the laser beams exchanged between the satellites will be influenced by Doppler shifts. That is, we expect a time-dependent response of the detector to these velocity perturbations. This signal is parametrized by a so-called response function X(t). The signal power spectrum, that we are interested in, is given by the absolute square of the Fourier transform of the detector response 3 , (II.2) In the case of LISA, this response function typically is a linear combination of the velocity perturbations of all three interferometer nodes. Its exact form depends on which undesirable noise sources are tried to be removed from the signal spectrum. Therefore, the detector response function is not unique. For concreteness, throughout this work, we will use the so-called Michelson response function X(t) for the readout of a signal at a single detector node [31,96]. For simplicity, we will not present its exact form here. Instead, we give a detailed definition in (A.1) of Appendix A. Nevertheless, let us point out its main features. Naively, the components of the response function are given by projections of the velocity perturbations of the nodes onto the interferometer arms, Here, the n i denote the unit vectors pointing between two nodes, labelled by the opposite side of the triangle, v i is the velocity perturbation of the i-th node induced by the gravitational pull and c is the speed of light. The integer coefficients a ij take into account retardation effects along the different signal paths. For more details on this notation see Appendix A and also [96,97] for alternative response functions. In fact, in some situations the detector response function above can be simplified (see, e.g., [31]). For instance, 3 Technically, the power spectrum of a given signal X(t) is defined by the Fourier transform of the so-called auto-correlation function, P (ω) = F [(X X) (τ )] = F R dt X * (t)X(t + τ ) . The latter is, in fact, equivalent to |X(ω)| 2 , i.e. to the definition given in the main text. Throughout this work, we will denote the Fourier transform of a function f (t) byf (ω).
if the velocity perturbations of two nodes are negligible for all practical purposes, e.g., v 2 ≈ v 3 ≈ 0, it reduces to where v 1 (t) is the velocity perturbation of a single node and we have used that n 1 + n 2 + n 3 = 0. The latter condition is true, if all interferometer arms are of equal length. Carefully note that this form of X(t) is only an approximation of the detector response. Eq. (II.4) is the dominant contribution to the exact response function, if the dark matter clump approaches the detector node very closely, or, in other words, if the impact parameter is smaller than the arm length of the interferometer, D L. This is the so-called close-approach limit (see, e.g., [16,31]). In contrast, if the impact parameter is much larger than the arm length (sometimes called the tidal limit, see [16,31]), D L, the differential gravitational pull on the nodes can be qualitatively different from the pull on a single node. To see this, imagine a simple situation where two nodes are aligned on an axis perpendicular to the trajectory of the dark matter clump, with their relative distance to each other being much smaller than their distance to the clump's trajectory. In this scenario, the impact parameters of both satellites will differ by the arm length, i.e. they read D and D +L, respectively. According to (II.3), the detector response will be proportional to the differential gravitational acceleration between the two nodes, with impact parameters and D and D + L, Hence, comparing both regimes, we observe the behavior That is, the detector response in the tidal regime, D L, falls off much faster with distance than in the closeapproach limit, D L. In this regime, we therefore expect the close-approach approximation of the detector response to break down and one would need to take the exact response function into account. Nevertheless, we will argue a posteriori that, for the purpose of this section, we only have to consider the close-approach case, D L, for events with a detectable signal-to-noise ratio. Therefore, (II.4) gives a reasonable approximation to the detector response. We note however, that when the sensitivity of the experiment becomes better, i.e. the noise is reduced, events with larger impact parameters will become detectable and the close-approach approximation might have to be reconsidered.
In practice, by means of (II.4) we can compute the response of the interferometer to an arbitrary velocity perturbation. For concreteness, we have considered the velocity perturbation associated to a gravitational pull by a dark matter clump passing by a single LISA satellite in (II.1) in a specific coordinate frame. However, in principle, the dark matter clump can approach the satellite from any direction. In order to account for this, we can equivalently choose an arbitrary orientation of the LISA experiment. That is, we can parametrize the unit vector n 1 in (II.4) by n 1 = (sin ϑ cos ϕ, sin ϑ sin ϕ, cos ϑ). The angles ϑ and ϕ essentially implement the arbitrary orientation of the detector plane relative to the dark matter clump. For a detailed discussion of this, see Appendix A 1. As we are interested in the signal power spectrum (II.2), it is then natural to consider the Fourier transform 4 of the response function,X(ω). In our scenario, the latter is given by (cf. [31], adding suitable time delays which give the L dependent factor) where K i is the i-th modified Bessel function of the second kind. In principle, we could now obtain the signal power spectrum P (ω) by squaring this expression. Obviously, the spectrum would then depend on the orientation of the LISA detector relative to the trajectory of the dark matter clump through the angles ϑ and ϕ. As a simple approximation it is reasonable to assume that the dark matter clumps generically move in a random direction. That means, in principle, any value of ϑ and ϕ is equally likely. To account for this, we can assume a uniform distribution for both, of which we then take the average (see Appendix A 1 for details). Before proceeding we note, however, that a uniform average is only a rough approximation. It assumes that the local dark matter distribution is isotropic. While this may be true in an isotropic reference frame where the detector is at rest, the Sun, together with the experiment, is moving through the dark matter halo at a constant velocity, thereby imposing a preferred direction on the system. Therefore, strictly speaking, the distribution of the angles parametrizing the relative orientation between the dark matter and the detector plane should not be uniform. This, however, is somewhat ameliorated by the motion of LISA itself (see also [31]), which changes direction along its orbit around the sun (cf., e.g., [98]). While certainly not all directions are attained with equal probability, it nevertheless amounts to at least a partial averaging. For a more detailed discussion of this see A 4. Finally, with this caveat in mind, we arrive at the angular-averaged signal power spectrum associated to the 4 In this work, we use the symmetric convention of the Fourier transform,f (ω) = (2π) −1/2 R dt f (t) exp(iωt).
gravitational pull of a dark matter clump on an interferometer node [31], For a few examples of the typical shape of P (ω) we refer the reader to Fig. 1. Note that there, for reasons that will become clear momentarily, we show the signal power spectral density, 4ωP (ω). Obviously, the signal power spectrum is only useful for an experimental test, if the desired signal can be distinguished from the background noise that the experiment is subject to. In general, at a gravitational wave interferometer the background noise is characterized by a noise power spectrum, commonly defined by ñ(f )ñ * (f ) = 1 2 δ(f − f )S n (f ) (for a comprehensive review see, e.g., [99]). For future experiments such as LISA this is not yet completely settled. In principle, different estimates can lead to quite different results for the detection rate. For concreteness, we will use a recent estimate [100] as our benchmark,  Note that, here, we have added an additional factor of 4(ωL/c) 4 /(1 + (ωL/c) 2 ) 2 compared to [100]. This essentially acts as an estimate of a transfer function to convert the original strain spectrum to the equivalent of our signal spectrumX 5 . Furthermore, we have divided by 2π in order to match our (symmetric) convention of the Fourier transform to the convention typically used in signal processing. We further remark that other choices of a noise power spectrum might also be reasonable, for example, the ones presented in [98,101].
Finally, as a measure for distinguishing a signal from the background noise, we can define the signal-to-noise ratio by comparing the signal-to the noise power spectrum over the range of all frequencies (see, e.g., [94]), (II.12) Here, one factor of 2 arises from the definition of the noise power spectral density given above. A second factor of 2 reflects the fact that we are considering single-sided power spectra only [31,94,99]. The form on the very right hand side of the equation is particularly useful for a quick estimation of the signal-to-noise ratio from the usual logarithmic plots of the sensitivity, as it is obtained from an integral of the logarithm over a dimensionless ratio between the signal, 4ωP (ω), and the noise, S n (ω), power spectral density. Before continuing, let us remark that the signal power spectrum P given in (II.8) remains constant with decreasing frequency due to the constant velocity of the satellite in the asymptotic future. Nevertheless, this does not lead to an infinite signal-to-noise ratio when normalizing to our choice of the noise power spectrum, S n . Instead, it remains finite, even when integrated over all frequencies. As was already noted in [31], alternative noise power spectra, however, might not share this feature and hence require the introduction of an experimentspecific lower frequency cut-off. For LISA imposing a cutoff ω c ∼ 10 −4 Hz due to experimental limitations provides a relatively conservative estimate of the lower end of the frequency band that the detector is sensitive to (see, e.g., [31]). However, in this paper we use the closeapproach approximation of dark matter clump encounters with the interferometer. This is valid for impact parameters smaller than the size of the experiment, D L. Therefore, we will use an even more conservative cut-off, that is essentially determined by the characteristic time of flight of a dark matter clump through the detector volume, ω c ∼ 2πV /L ∼ 10 −3 Hz.
In Fig. 1 we show the angular-averaged signal power spectral density, 4ωP (ω) with P (ω) given in (II.8), associated to dark matter clumps traversing the detector ating a strain signal of a gravitational wave in the channel X(t), yielding a transfer function proportional to 4 sin 4 (ωL/c). In this case our factor corresponds to the envelope of this conversion factor. Note that this approximation is correct up to a possible constant factor of order O(1 − 10). volume and compare it to the experimental sensitivity of LISA. As an example, we consider events with an impact parameter of D = 50, 000 km. The coloured lines correspond to the signal, while the black solid line shows the noise power spectrum S n (ω) given in (II.9), which essentially determines the overall experimental sensitivity of the detector. As indicated earlier, the signal-tonoise ratio is the logarithimic integral over the ratio of the two plotted quantities and we can roughly read this figure in such way, that we have a chance to distinguish a signal from the detector noise, whenever the coloured signal exceeds the black noise spectrum for a sufficient range of frequencies. Also note that the signal power spectral density shown here, is uniformly averaged over the angles parametrizing the relative orientation between the interferometer and the trajectory of the dark matter clump (cf. Appendix A 1). As explained above, this approximates the typical size of the signal. The signal shape and strength of individual events will be different. For instance, in extreme cases, we expect a different signal from a clump with normal incidence to the detector plane as from one that traverses the detector volume almost parallel to it. Furthermore, as we have already pointed out before, the uniform average can only be seen as an approximation. This is because there is preferred direction given by the Sun, together with the detector, moving through the dark matter halo at a constant velocity (see Appendix A 4) as well as the rotation of the detector itself not covering all angles equally. In general, the signal power spectrum receives large contributions from low frequencies, while it quickly drops in the high frequency regime. As pointed out earlier, the constant value at low frequencies is due the constant velocity component of the LISA satellite in the asymptotic future. Consequently, the sensitivity may benefit from background noise that is reduced in the low frequency tail of the spectrum. This, however, will also require going beyond the close approach approximation. Keeping this caveat in mind, below we investigate the benefits of improving at low frequencies by considering two different low-frequency cutoffs.
Having reviewed all necessary aspects of distinguishing a possible signal induced by a localized clump of dark matter traversing the detector volume of LISA from background noise, let us now quantify the discovery potential for these clumps inside the dark matter halo of our Galaxy. The mass of the dark matter clumps, M DM , determines the characteristic distance between them, d ∼ (M DM /ρ DM ) 1/3 . Therefore, it controls the average rate of encounters of a clump with one of the satellites at or below a given impact parameter D, where Φ DM is the effective dark matter flux at velocity v DM , given by Φ DM ∼ ρ DM v DM /M DM . That is, naively, η is the rate at which, on average, a dark matter clump of mass M DM and velocity v DM passes through a surface of radius D. Therefore, applied to our scenario, we can use it to estimate the rate at which we expect a dark matter clump to induce a signal in the interferometer. Note, however, that using the average signal strength in the calculation of the required impact parameter and thereforė η is only an approximation. In general, clumps passing by at a suitable angle may already give rise to a signal at somewhat larger distances than indicated by the average signal strength and, similarly, closer encounters are needed for other angles. However, we expect that this simplistic treatment nevertheless captures the effect reasonably well. Furthermore, we assume the velocity distribution of the dark matter clumps inside the halo to be a (simplified) superposition of the Sun moving through the Galaxy at v ≈ 220 km/s [102] and the dark matter velocity having a uniformly random direction 6 while its magnitude is Maxwell-Boltzmann distributed with a root mean square of v rms = 3/2v ≈ 270 km/s (see, e.g., [103]). For a discussion of possible caveats of this approximation and how it is implemented in practice, see Appendix A 4. In addition, we fix the dark matter energy density to ρ DM ≈ 0.39 GeV/cm 3 [104]. Note that these values can come with relative uncertainties of up to 25 %, which could alter our results by a similar amount. We As long as this is smaller than the size of LISA we expect the close-approach approximation to be valid. We use the same baseline parameters for the dark matter density and velocity as in Fig. 1.
remark, however, that the approximations we employ as well as the differences in the noise power spectra probably cause larger uncertainties.
In Fig. 2 we show the average gravitational interaction rate as a parametric function of the signal-to-noise ratio as well as the mass of the dark matter clumps. As pointed out before, here, we only take impact parameters into account that are smaller than the arm length of LISA, D L, i.e. we consider the close-approach regime. This is indicated by the black dashed line in the bottom panel (with the appropriate scale being on the right hand side), which illustrates typical impact parameters at which the signal-to-noise ratio exceeds one, SNR 1, for a given mass of the dark matter clump. We can see that this is always smaller than the typical size of LISA and therefore the close-approach approximation of the detector response function (II.4) is valid in the mass regime we show here. This approximation has to be reconsidered, if the sensitivity of the experiment is improved. In this case, dark matter clumps at larger impact parameters can still be detected. The impact parameter shown in Fig. 2 also gives an upper bound on the size of the dark matter clumps for which our treatment of the clumps as point-like is valid (most of the signal rate originates from the largest detectable impact parameters).
In general, we observe that with increasing mass of the dark matter clumps, the signal-to-noise ratio is enhanced. However, at the same time, the event rate of encounters is reduced due to the reduced effective dark matter flux, Φ DM ∼ M −1 DM . That is, at higher masses there is a balance between the increased signal-to-noise ratio and the reduced interaction rate. In particular, we find an optimal detection potential that passes a minimal detection threshold, SNR 1, for dark matter clumps of mass M DM ≈ 10 10 kg which could be observed at LISA approximately every 200 years on average. However, we note that this threshold of a signal-to-noise ratio, SNR 1, is rather optimistic. In practice, as the signal itself has to be distinguished from other sources, the detection threshold may be significantly higher, of the order of SNR 10 (e.g. this is on the lower end of the typical signals considered in the Mock LISA Data Challenges [105,106]). As shown in Fig. 2, this will reduce the average detection rate by an order of magnitude while shifting the optimal detection potential to dark matter clumps of higher masses.
Unfortunately, the sensitivity is not yet on a desirable level for a near future discovery potential. Let us note, however, that the latter value has to be understood as a relatively conservative estimate, which could also be significantly higher. As already mentioned, the signal-tonoise ratio depends on the noise power spectrum which features a comparably large uncertainty. For the relatively broad signal spectra associated to the dark matter clumps, in particular the low frequency tails of the spectrum may be relevant. Therefore, the signal-to-noise ratio may vary significantly, if the overall noise spectrum was smaller, to some degree also at low frequencies. We try to illustrate this through the bars in the top panel of Fig. 2. For these, we push the close-approach approximation closer to the boundary of its validity by introducing a low-frequency cutoff of ω c ≈ 5 × 10 −4 Hz as compared to the one imposed by the close-approach approximation, ω c ≈ 10 −3 Hz. Let us stress, however, that this really has to be understood only as a very first estimate of potential improvements. Finally, we note that, in light of the very rough nature of our estimates, the results presented Fig. 2 are reasonably in line with but perhaps somewhat more pessimistic than earlier works, in particular [24], with regards to the detection of dark matter clumps LISA.
While the obtained rate is still rather low, it nevertheless brings us much closer to a desired level, so that we can hope that further improvements both in the detector as well as in the analysis might allow for a detection in a reasonable time frame. Indeed, one such improvement could be in the analysis. For example, already the authors of [15] noted that the dark matter clump inter-action is inelastic and therefore differs from the elastic interaction of a gravitational wave exploited in current detection strategies. Further improvements are expected from a more detailed analysis of the time structure of potential signals as discussed in [28].

III. TOPOLOGICAL DEFECTS
In this section, we now want to go beyond dark matter clumps and investigate the detection of structures such as cosmic strings and domain walls 7 . As already mentioned in the introduction, these topological defects might have been produced in the early Universe [75][76][77], but their use as dark matter requires a non-trivial behavior, e.g. interacting networks, of such objects. We do not address how such a network is formed or how it can be made to satisfy the constraints imposed by the properties of dark matter (or dark energy) but simply assume its presence with a given density. This is in the spirit of and follows the completely phenomenological approach to study their detection, also pursued in [29,35,38,41,[87][88][89][90][91][92][93] in the context of various different detection techniques. We briefly note that our discussion also applies to localized structures of "ordinary" dark matter with a string-like or domain wall-like geometry, an example of which are the string-like axion streams that were already investigated with regard to LISA in [30].
The gravitational properties of topological objects can significantly differ from those of non-relativistic matter. Cosmologically, this is due to a different equation of state, relating the pressure and the energy density of the cosmic fluid, w = p/ρ. This then results in a typical scaling behavior ρ ∼ a −3(1+w) of the average energy density on cosmological scales.
For our signals we are interested also in the behavior in the vicinity of individual objects. Very naively, this can be obtained from Poisson's equation for the Newtonian gravitational potential. For a fluid this is sourced by the combination of pressure and energy (ρ + 3p) (see, e.g., [109]), and the gravitational potential satisfies i.e. it explicitly depends on the equation of state of the source. However, it is not obvious that the gravitational field of topological objects such as strings and domain walls behaves locally as suggested by their global equation of state. That said, in good approximation this is nevertheless true for static strings as well as domain walls, as shown by the results of [110,111] and discussed below for each case in the respective subsection. To model the local field of more complicated networks of topological structures we therefore simply assume that it can be approximated by Eq. (III.1) with the equation of state parameter given by its cosmological value.
The equation of state as well as their dimension for cosmic vacuum strings and domain walls leads to significantly different gravitational fields sourced by these objects compared to non-relativistic clumps of matter. In the following subsections, we want to discuss their gravitational properties and the related imprints they might leave at LISA.

A. Cosmic strings
In order to obtain the gravitational potential of a cosmic string, we can solve (III.1) in a cylindrically symmetric space for an energy density that is distributed along an infinite string, e.g. ρ = µδ(x)δ(y). Here, µ is the tension of the string, i.e. the energy stored per unit length. This yields a gravitational potential that grows logarithmically with distance, φ ∼ log (r 0 /r), and sources the gravitational field where r denotes the radial distance to the string. In general, we can directly use (III.2) to obtain the signal strength at LISA. However, before we proceed, let us mention a noteworthy special case. If we consider a single static vacuum string, the equation of state is given by w = −1/3, and its energy density dilutes as ρ S ∼ a −2 [112]. At the same time, this implies that its gravitational field vanishes, indicating that a static cosmic string does not couple to matter [110,113]. Nevertheless, this picture can change, if the vacuum string starts to be dynamical. For instance, a string moving at a velocity β = v/c has a modified equation of state, w = 2/3β 2 − 1/3 [112]. Going beyond such scenario, the situation can deviate from this even more, if multiple strings are considered. In particular, if the strings interact with each other to form a network, the corresponding equation of state can drastically change. Therefore it is possible that they can contribute to the dark matter or dark energy component of the Universe (see, e.g., [82]).
The naive reasoning from above suggests, that we would not expect any gravitational pull on the interferometer by a static vacuum string at all. However, it was shown that, due to the globally non-trivial conical spacetime geometry sourced by the string, it will still attract massive objects around it [114]. In fact, the gravitational acceleration of a mass m in the vicinity of the string is given by [114] where κ ≈ 1/32 for small Gµ/c 2 . In this case, the gravitational field falls of as g ∼ r −2 in contrast to the r −1 asymptotics of the general configuration (III.2). Naively, this can be understood as follows. An infinite, straight and static cosmic string sources a conical spacetime geometry that can lead to double copies, i.e. mirror images, of nearby objects [110,113]. In this sense, (III.3) can be understood as the gravitational field sourced by the effective mass, (Gµ/c 2 )m. Hence, we can interpret the gravitational field of the string as similar to the usual attractive force between two massive objects. Cosmological considerations, such as the isotropy of the cosmic microwave background, provide bounds on the dimensionless string tension, which are typically of the order Gµ/c 2 10 −6 (see, e.g., [115,116]). Therefore, we expect a tiny acceleration of a test mass in the gravitational field of a static cosmic string.
In the following, we will distinguish between signals due to the gravitational pull on the interferometer generated by a static vacuum string and by an interacting network of cosmic strings that accounts for dark matter or dark energy. The former will be characterized by the gravitational field (III.3), while the latter is given by (III.2) with a general equation of state.

Signal power spectrum
Given the gravitational field of a cosmic string, we can proceed analogously to the case of dark matter clumps and determine the gravitational acceleration that each node of the interferometer experiences in the vicinity of a string. For concreteness, let us choose a coordinate frame, in which the (infinite) string is parallel to the z-axis and the LISA satellite located at the origin is, initially at t 0 = 0, at a minimum distance D to the string. We furthermore assume that the string is uniformly moving at velocity V . Due to the additional internal orientation of the string as compared to spherical clumps, in the frame where the string motion is confined to the yz-plane, its velocity can have a component in the y-as well as the z-direction, i.e. V y = V sin θ and V z = V cos θ, respectively. In other words, the string has an additional inclination angle when approaching the interferometer node. For a schematic illustration see the bottom panel of Fig. 12 in Appendix A. In this reference frame, in the gravitational field of a static cosmic string (III.3), a test mass is subject to an acceleration of while in the field of an interacting string network (III.2) it reads Note that in both expressions only the velocity component V y = V sin θ appears. This is because only the radial distance to the string, in this reference frame determined by the x-and y-coordinate, enters the gravitational field. Along the lines of our discussion in Section II, in order to obtain the frequency spectrum of the detector response, we now need to determine the Fourier transform of the velocity perturbations associated to the accelerations of the detector nodes. As we will again argue a posteriori, we consider the regime where the strings traverse the detector volume in the close vicinity of a single node, i.e. the close-approach limit, D L. Therefore, the detector response can be approximated by (II.4), which in the frequency domain readsX(ω) In addition to the orientation of the string and its motion relative to the node of the interferometer, the arbitrary orientation of the detector plane has yet to be implemented. In the close-approach approximation, we can take this into account by parametrizing n 1 accordingly, n 1 = (sin ϑ cos ϕ, sin ϑ sin ϕ, cos ϑ). For a detailed discussion of this see Appendix A 2.
By squaring the Fourier transform of the response function, |X(ω)| 2 , we obtain the signal power spectrum induced by a cosmic string in the vicinity of the interferometer. At this point, the latter still depends on the relative orientation between the string and the detector plane, parametrized by ϑ and ϕ. As already noted for the case of clumps, we assume that any orientation occurs equally likely (a discussion of the validity of this simplistic approximation and additional details are given in Appendix A 4). Therefore, we take the uniform average over both (cf. Appendix A 2). Finally, the angular-averaged signal induced by the gravitational pull of a static cosmic string is given by where G is the Meijer G-function. For an interacting network of cosmic strings with an arbitrary equation of state, it reads (III.7) Here, K 1 denotes the first modified Bessel function of the second kind.
In the panels of Fig. 3, we illustrate an example of the angular-averaged signal power spectral densities induced by a static cosmic string (top) and an interacting string network with w = 0 (bottom) and compare it to the experimental sensitivity of LISA. In particular, we show different values of the string tension and consider, as an example, events with a fixed impact parameter of D = 100, 000 km. Moreover, we have assumed Here, we have chosen an impact parameter of D = 100, 000 km. Therefore, as the strings are taken to account for dark matter, these signal events are very rare and occur every 10,000,000 (blue) to 1,000 (green) years on average. The black lines correspond to the sensitivity of LISA and the vertical dashed lines to two different values of the used low-frequency cutoff. In both panels we have assumed a Maxwell-Boltzmann distribution for the velocity of the strings with root mean square vrms ≈ 270 km/s in combination with the LISA experiment moving through the galaxy at v ≈ 220 km/s. a (simplified) superposition of a Maxwell-Boltzmann distribution for the velocity of the strings with root mean square v rms ≈ 270 km/s in combination with the LISA experiment moving through the galaxy at v ≈ 220 km/s. We present the practical implementation of this in Appendix A 4.
Similar to the dark matter clumps, we find that the signal receives large contributions from low frequencies, while it quickly drops in the high frequency regime. The signal induced by a static cosmic string is orders of magnitude smaller compared to the signal caused by interacting cosmic strings with w = 0 as appropriate for the dark matter component of the Universe. This is due to the fact that the gravitational perturbation of the former is sourced by the tiny mirror mass, (Gµ/c 2 )m. In fact, the suppression already indicates that only an interacting network of cosmic strings is within experimental reach of LISA. Average gravitational interaction rate of cosmic strings with the interferometer as a function of the signalto-noise ratio. We distinguish between a network of static cosmic strings (top) and an interacting network of dynamical strings (bottom), that constitutes the dark matter or dark energy component of the universe. In the top panel, the colors denote typical values of the dimensionless string tension, Gµ/c 2 , while we have fixed the energy density to the dark matter one, ρS ≈ 0.39 GeV/cm 3 . In the bottom panel, we have chosen a string tension of Gµ/c 2 ≈ 10 −23 and fixed the energy density and equation of state for dark matter as (ρDM ≈ 0.39 GeV/cm 3 , w = 0) and for dark energy as (ρDE ≈ 3.2 keV/cm 3 , w = −1). In both panels we have assumed a Maxwell-Boltzmann distribution for the velocity of the strings with root mean square vrms ≈ 270 km/s in combination with the LISA experiment moving through the galaxy at v ≈ 220 km/s. Given the angular-averaged signal power spectra, let us now estimate the discovery potential for cosmic strings in the vicinity of the gravitational wave interferometer. Assuming that the energy stored in the cosmic string network is determined by a single scale d, that is the characteristic distance between the strings, the energy density is given by ρ S ∼ µ/d 2 . In addition, this scale will also determine the average interaction rate of a string with the detector. Since we have averaged over the angles describing the orientation of the string, for simplicity, we can assume that the strings are all aligned and moving into the same direction perpendicular to their orientation. Projecting into the plane orthogonal to the strings then effectively reduces the problem of estimating the flux of cosmic strings to computing the flux of pointparticles through a given surface element. Here, however, we have projected out one dimension and are thus dealing with the flux through a line element. Accordingly, the rate of cosmic strings approaching the detector at velocity V and impact parameter D can be estimated to beη Note that, in line with the angular-averaged signal power spectral density shown in Fig. 4,η is an estimate of the average rate at which a cosmic string of velocity V is passing by the detector at a distance D. Therefore, it is only an approximation, as it does not take the shape of the signal into account, which, for instance, can crucially depend on the relative orientation between the string and the interferometer (this is similar to the discussion of dark matter clumps in Section II). Nevertheless, we do not expect any significant changes of our conclusions, if this complication was taken into account more carefully.
The panels of Fig. 4 show the average gravitational interaction rate of a cosmic string network of static (top) and interacting strings for w = 0 and w = −1 (bottom; for the yellow curve we also choose the energy density to be that of dark energy, ρ DE ≈ 3.2 keV/cm 3 [117], which will use throughout this work) with the LISA interferometer as a function of the signal-to-noise ratio. The bars illustrate the alternative low-frequency cutoff for the noise power spectrum of the detector (cf. Section II for details). For static cosmic strings, we find that, due to their tiny gravitational field, the rate of sufficiently strong interactions with the detector is basically negligible, even for the most optimistic values of the string tension. Naively, one could of course consider situations where one tries to probe the gravitational field with more massive test objects, thereby increasing the source of the field. This, however, is completely determined by the experimental setup, such as LISA in our case, and not a free parameter. Consequently, the tiny gravitational field renders the signal of a static cosmic string by its gravitational pull on the interferometer unobservable.
The more interesting case is an interacting network with an equation of state w = 0, suitable for being the dominant component of the dark matter. Due to the greatly enhanced gravitational field, the situation is significantly improved, as we demonstrate in the bottom panel of Fig. 4. For instance, we find that a network with strings of tension Gµ/c 2 ≈ 10 −23 could induce a signal in the interferometer with a signal-to-noise ratio of about 10 every 10,000 years on average. Taking a more optimistic low-frequency cutoff this could even be improved, as indicated by the bars in the figure. Therefore, an overall gain in experimental sensitivity could allow to gain a sizeable factor in the signal event rate, giving hope that already moderate improvements will allow for a detection. As a consequence, similar to the case of dark matter clumps, our analysis strategy would benefit from improvements in the low frequency regime. Note that, in order to obtain a decent interaction rate with a signal-to-noise ratio still greater than one, we have, to some extent, chosen the string tension close to an optimal value, Gµ/c 2 ≈ 10 −23 . We will show the overall experimental sensitivity of LISA for different string tensions in Fig. 10 of Section VI.
We conclude that the overall discovery potential of a cosmic string network with the LISA interferometer depends on its dynamics, i.e. its equation of state. Our estimates of the signal-to-noise ratio indicate that a network of static vacuum strings appears not to be observable, while an interacting network with an equation of state appropriate for dark matter may be closer to the experimental reach of LISA. Furthermore, our comparison of different low-frequency cutoffs suggests that the overall discovery potential may be increased by improvements in the low frequency regime.

B. Domain walls
In order to obtain the gravitational field sourced by a domain wall, we can solve Poisson's equation (III.1) for an energy density confined to an infinite plane, e.g., ρ = σδ(x). Here, σ denotes the surface tension of the domain wall and we have neglected a possibly finite thickness. In this background, the gravitational potential grows linearly with the distance, φ ∼ r. Therefore, the gravitational field of a domain wall with an arbitrary orientation is constant, Here, n is the unit vector normal to the plane parametrizing the domain wall and the sign ensures that the field always points towards (or, as we will see momentarily, even away from) the wall 8 . That is, the gravitational field of a domain wall is pointing in the direction normal to it and, in particular, is independent of the distance to the wall. Similar to cosmic strings, a network of domain walls can have different equations of state, depending on its dynamics. For instance, the equation of state of a static domain wall is given by w = −2/3, such that its energy density dilutes as ρ DW ∼ a −1 , while a domain wall moving at a velocity β = v/c obeys an equation of state of w = β 2 − 2/3 [112]. Intriguingly, according to (III.9), this implies that the gravitational field of a static domain wall is repulsive rather than attractive [110,111]. Nevertheless, complicated dynamics of an interacting network of domain walls may lead to a very different equation of state. Therefore, it can as well serve as an exotic candidate for dark matter or dark energy [82][83][84] with a corresponding equation of state.
In the following, we want to derive the signal due to the gravitational pull (or push, in the case of a repulsive potential) of a domain wall travelling through the LISA interferometer.

Signal power spectrum
The gravitational field of a domain wall is independent of the distance to it. In contrast to the detection of dark matter clumps or cosmic strings, this means that all three nodes of the interferometer will experience the same acceleration due to the presence of a domain wall. Therefore, in order to obtain a differential acceleration and hence a signal in the interferometer, the domain wall has to traverse the space between the different satellites. That is, it has to completely separate one spacecraft from the other two, thereby accelerating them into opposite directions.
The velocity perturbation that each node picks up due to the constant acceleration in the gravitational field of the domain wall reads v(t) = ∓2πG(1 + 3w)σtΘ (t) n , (III.10) where n denotes the normal vector of the domain wall and the different signs follow the conventions of (III.9). Here, for simplicity, we assume that the wall starts traversing the interferometer at a time t 0 = 0. Hence, the Θ-function implements the fact that, due to the equal acceleration of the three nodes, we do not expect a signal, if the domain wall does not separate the individual nodes from each other. Similarly, we will assume that the signal ceases to exist when the domain wall has traversed the detector volume completely. In other words, for simplicity, we consider a signal induced by a domain wall, that starts travelling through the detector by passing the first node, then passing the second, which thereafter gets accelerated in the opposite direction, and finally traverses the third node, after which it ceases to exist. That is, strictly speaking, we view the signal as caused by the acceleration burst instead of the individual velocity perturbations. When considering the latter, obviously, the signal will not cease to exist after the domain wall has passed the last satellite, as there is still a differential velocity shift between all three nodes. In principle, this will, in addition, lead to a persistent deformation of the LISA triangle. In general, the situation is even more complicated, when considering a complex network of interacting domain walls which traverse the detector volume at different times and directions. A correct and thorough treatment would require an involved numerical simulation of this scenario. Here, we will not consider this layer of complexity. Due to the constant gravitational field, the signal that a domain wall will induce in the interferometer when traversing the detector volume involves each of its three nodes. Therefore, the detector response has to be parametrized by the exact response function given In contrast to the case of clumps and strings, the lower one is chosen at a smaller frequency as our calculation for the domain walls has no close-approach limitation in the present case. Here, we assume that the domain wall is moving at a velocity of v/c = 10 −3 . If the domain wall network is to cover the dark matter component of the Universe, the signal events shown here are quite infrequent. On average, we expect them every 1.5 (blue), 150 (yellow) and, in the extreme case, 15,000 (green) years.
in (A.1), where the v i are now given by (III.10). The terms of the detector response function can be evaluated explicitly, e.g., in a reference frame where each node is located on a coordinate axis, r i = L/ √ 2e i , and we parametrize the domain wall by the unit vector n = (sin θ cos φ, sin θ sin φ, cos θ). We can then proceed by considering the Fourier transform of the response function, i.e. of the velocity perturbations, and uniformly averaging over the angles (θ, φ) in order to obtain the angular-averaged signal power spectrum. For a detailed discussion of the geometrical aspects of this, see Appendix A 3. We also note that, similar to our discussion of localized dark matter clumps, a uniform average is only an approximation. We discuss its validity and more details in Appendix A 4.
In Fig. 5 we illustrate the angular-averaged signal power spectral density due to the differential gravitational acceleration by a domain wall traversing the interferometer and compare it to the experimental sensitivity of LISA. In particular, we consider domain walls of different surface tensions and fix their typical velocity to v/c = 10 −3 . We find that, compared to the case of dark matter clumps or cosmic strings (cf. Fig. 1 and Fig. 4), the signal caused by a domain wall is significantly enhanced in the high frequency regime of LISA's characteristic frequency range. One factor contributing to this is that we assumed an infinitely thin domain wall and point like nodes of the experiment. In practice, both have a finite thickness that should lead to a faster drop off at large frequencies. However, this does not play any significant role in our determination of the sensitivity.
Furthermore, we can estimate the average rate of gravitational interactions with the LISA detector that we expect for a given network of domain walls. The energy density that is stored in a network is essentially determined by the surface tension of the walls σ and the characteristic distance d between them, ρ DW ∼ σ/d. In order to estimate the rate at which we expect them to approach the detector, we can consider the simplified situation where the domain walls of the network are all parallel to each other 9 . The overall rate is then simply given by the inverse time interval between two consecutive walls travelling through the detector volume. That is, if the domain walls move at a certain velocity v, the rate at which we expect them to approach the detector can be roughly estimated bẏ Note that, in line with the angular-averaged signal power spectral density shown in Fig. 5,η is merely an estimate of the average rate at which a domain wall is passing through the detector and is therefore only an approximation. In particular, it does not take the shape of the signal caused by a specific domain wall into account, which, for instance, can crucially depend on the relative orientation between the domain wall and the detector (this is similar to the our discussion of dark matter clumps in Section II). As an extreme example, a domain wall that is exactly parallel to the interferometer plane would not leave any experimental imprint in the detector. While this does not make any significant difference for the results we show in this section, we present results of using a somewhat improved approach in Fig. 10 of Section VI. There we account for this additional angular dependence. In particular, a thorough treatment requires to weight the signal events according to their relative orientation with respect to the detector. For instance, as a simple schematic example, this can be written (see also Appendix A 3) (III.12) where s is the required signal-to-noise ratio. That is, naively, only signal events are taken into account for configurations which lead to a detectable signal at the interferometer. Nevertheless, this does not qualitatively alter the results presented in this section. Fig. 6 shows the average gravitational interaction rate of a domain wall network as a function of the signal-tonoise ratio. In the top panel, we consider the case of a static domain wall network, while in the bottom panel, we consider an interacting network of domain walls that constitute the dark matter or dark energy component of the Universe. Similar to Fig. 2, the bars illustrate an alternative choice of the low-frequency cutoff for the noise power spectrum of the detector given in (II.9). However, as our estimate of the signal induced by a domain wall does not rely on the close-approach approximation, the cutoff can be shifted towards even lower frequencies.
For both setups, we observe that with increasing signalto-noise ratio the event rate reduces according to a linear power law. This is because, here, we consider the angular-averaged signal-to-noise ratio, SNR ∝ σ, as well as interaction rate,η ∝ σ −1 , which are solely determined by the domain wall surface tension σ. That is, each data point shown in the figure corresponds to a specific value of the latter. In general, as expected, the overall discovery potential for static and interacting domain walls is comparable, as long as the energy density is fixed to that of dark matter. Intriguingly, it is also somewhat better than in the case of dark matter clumps or cosmic strings. For instance, we expect a signal due to a domain wall traversing the interferometer with a signal-to-noise ratio of about ten every 10 to 100 years on average. Taking into account potential improvements due to a more optimistic noise spectrum, this rate might also be higher. In this case the same signal may even be expected almost every 1 to 10 years on average. Therefore, in the future, such a domain wall network could certainly be within experimental reach of LISA.
Let us close this discussion with a few words of caution. The estimates presented in this section are subject to some simplifications we have made in our derivation. Our results are based on the assumption that the domain wall can be parametrized by an infinite plane of vanishing thickness. This assumption, however, might not always be fully justified in particle physics models that admit domain wall solutions. Moreover, and perhaps more importantly, we have only considered the situation of a domain wall inducing a signal that ceases to exist once the wall has traversed the entire detector volume. In principle, one would also have to include the remnant velocities of the nodes relative to each other once the domain wall has passed, which most likely would modify the signal power spectrum at low frequencies. In addition, this would also change the triangular detector geometry persistently. Nevertheless, our relatively conservative estimates for gravitational interactions of domain wall networks with the LISA interferometer provide hope for a future discovery potential.

IV. STOCHASTIC FLUCTUATIONS OF THE DARK MATTER DENSITY
In addition to strongly localized energy densities, such as dark matter clumps or topological defects travelling through the Universe, the gravitational potential in the vicinity of the interferometer can also be perturbed by stochastic fluctuations of the dark matter density, as already discussed in [20] for the example of pulsar timing arrays. To linear order, these gravitational perturbations satisfy where δφ denotes the perturbation of the gravitational potential, while δ = δρ/ρ is the relative fluctuation of the mean densityρ. In principle, (time dependent) perturbations of the gravitational potential are precisely what gravitational wave experiments try to measure. In this section, we aim to estimate the sensitivity of LISA to these perturbations caused by fluctuations of the dark matter density. For simplicity, we consider a simple density fluctuation, oscillating in space, that travels through the interferometer. That is, after decomposing δ in Fourier space, the Here, we have fixed the normalization of the dark matter power spectrum to P = 1. We furthermore assume a velocity of v ≈ 220 km/s at which the detector is moving through the density wave, while for the dark matter density we use the same value as in Fig. 1.
gravitational response to a Fourier mode with wave vector k is given byδ where k = |k|.
We can now determine the detector response to a specific Fourier mode of the density fluctuation by considering the contribution to the gravitational perturbation of each mode separately, δg(x) = d 3 k δg k (x), where the perturbation of the gravitational field associated to a mode with wave vector k reads Naively, this corresponds to a situation, where we select a plane wave of specific wave length and "freeze it", as there is no dynamical wave evolution. Instead, we obtain a detector response to this density wave, if the interferometer moves through the perturbation of the gravitational field. To some extent, this setup is similar to a network of domain walls travelling through the interferometer, discussed in Section III B. Assuming a detector node is moving through the density perturbation of wave vector k at a constant velocity v, i.e. x(t) = vt + x 0 , the frequency spectrum of the associated gravitational acceleration is given bỹ (IV.4) Note that the right most term denotes the δ-distribution, not to be confused with the density fluctuation. This represents that, if the detector moves through a plane wave at a constant velocity, there is only one frequency contributing, which depends on the angle between the wave vector of the density perturbation and the motion of the detector.
The above gravitational acceleration will, in turn, lead to a velocity perturbation of each node. Similar to the discussion of domain walls in Section III B, we can plug each velocity perturbation into the detector response function, X(t), to obtain the response of the interferometer to a given Fourier mode of a density perturbation, X k (ω). Furthermore, in order to determine the detector response to a superposition of density fluctuations, we can then sum over all Fourier modes. The absolute square of this finally yields the signal power spectrum associated to the linear superposition of density perturbations, which schematically reads where, similar to the previous sections, n i denote the unit vectors pointing between two nodes, labelled by the opposite side of the triangle, and r i is the initial position of each node. The coefficients c ij (ω) encode the linear combination of the velocity perturbations of each interferometer node in the detector response (A.1). Therefore, they also include phases of the form exp(iωL/c) due to retardation effects. That is, naively, the signal response of the detector to a superposition of density fluctuations with different wavelengths is weighted according to the dark matter density power spectrum. For the latter we use the conventional definition (see, e.g., [118]) As a simple example of the density power spectrum, we consider a broken power law 10 , (IV.7) 10 In particle physics models of very light dark matter, such a spectrum and the wavelengths of interest to us may be achieved quite naturally, see, e.g., [67,68,119].
where P is a normalization constant and n is the spectral index. Obviously, P(k) exhibits a peak P at a characteristic scale k , which we treat as a free parameter. Indeed, for the purpose of this work, we assume that it obtains its maximum within the characteristic frequency band of LISA, ω . The corresponding wavelength is then of the order Again we average over angles (cf. Appendices A 3 and A 4). We illustrate the angular-averaged signal power spectral density associated to a dark matter density fluctuation traversing the interferometer in Fig. 7 and compare it to the overall experimental sensitivity of LISA. The signal power spectrum shown is normalized to unity, i.e. P = 1. We also show different choices of the spectral index n as well as different peak positions k of the dark matter power spectrum. For simplicity, we have fixed the velocity at which the detector is moving through the density wave to v ≈ 220 km/s. For density fluctuations of order 1 we find that, even for optimistic choices of the peaks, the overall signal power spectrum is suppressed by many orders of magnitude compared to the noise. As the signal scales linearly with the power spectrum of the fluctuations, P , we can see that enormous fluctuations are needed for a signal to be detectable. Indeed, this result is, to some extent, expected from our discussion of localized dark matter clumps. Naively, clumps correspond to a fluctuation of significant overdensity compared to the background. However, as we have shown in Section II, these are scarcely within experimental reach of the LISA detector. In summary, we therefore do not expect any detectable experimental signature of stochastic fluctuations of the dark matter density at LISA, unless the overdensities are very large (at which point they may actually resemble more the localized structures we have already discussed).

V. EXTRAPOLATION TO OTHER GRAVITATIONAL WAVE EXPERIMENTS
The general strategy we have presented in the previous sections is not limited to the LISA interferometer, but in principle applies to any experimental apparatus that is sensitive to gravitational perturbations. In this section, we want to extrapolate our results to other gravitational wave experiments. In particular, we aim to obtain the detectable rate corresponding to the gravitational pulls of dark matter clumps and topological defects in groundbased interferometers, for example LIGO [120,121], and pulsar timing arrays, for instance utilizing SKA [14]. Due to their different characteristic sizes compared to LISA, these detectors are sensitive to gravitational interactions in other frequency regimes and therefore, ultimately, to other energy densities of dark matter objects.

A. LIGO
With its pioneering measurements, LIGO has strongly advanced the field of gravitational wave astronomy [1,2]. For the purpose of this work, we will treat it as a classical Michelson interferometer, which aims to measure phase shifts between laser beams sent across two orthogonal arms. Since the characteristic length of both arms is L ≈ 4 km, it is sensitive to gravitational perturbations at frequencies of approximately 1 Hz to 1000 Hz [120,121].
In general, while the basic idea of our survey for LISA also applies to LIGO, there are a few differences due to the detector geometry. Most importantly, LIGO is a ground-based interferometer. Naively, its detector nodes are freely hanging mirrors suspended from a static laboratory frame, in stark contrast to the freely moving satellites of LISA. Therefore, while the detector response at LISA was determined by relative velocity shifts between the satellites, at LIGO we can only expect a signal due to sufficiently short gravitational acceleration bursts on the mirrors. Consequently, as a simplified detector response to a gravitational perturbation, we consider the differential acceleration between the four mirrors [94], where x and y denote the arms of the interferometer. In addition, being a ground-based interferometer, LIGO is subject to different sources of background noise, such as seismic motion, for example. To take the different background noise into account, we model the noise power spectrum according to the sensitivity curve of the design sensitivity of advanced LIGO [121]. Since we parametrize the detector response X(t) in terms of differential accelerations, we then convert the strain-to an acceleration noise by multiplying with (ω 2 L/2) 2 (see, e.g., [95]). With these definitions, we can repeat the analysis of the previous sections. We summarize our results in Fig. 8, where we show the average gravitational interaction rate of localized dark matter clumps as well as string and domain wall networks with LIGO as a function of the signal-to-noise ratio. Note that here, due to the small characteristic size of the detector as compared to LISA, we cannot use the close-approach approximation of the detector response, D L, in every region of the parameter space we are interested in. For simplicity, and in order to also take the regime D L into account, we extrapolate our results by multiplying the gravitational acceleration by a factor L/D in this limit (see (II.6) for a qualitative explanation).
In general, for localized dark matter clumps, we observe a behaviour qualitatively similar to the LISA analysis shown in Fig. 2. As one may expect, we find that LIGO with its much smaller size is sensitive to clumps with comparatively low masses. We find that dark matter clumps of masses 10 3 kg, could be observed every 2,000 years on average. More optimistically, given that these estimates crucially depend on the noise spectrum as well Here the color indicates the equation of state of the cosmic string network. We furthermore assume strings of tension Gµ/c 2 ≈ 10 −28 . (c) As in (b) but for a domain wall network travelling at a fixed velocity, v/c = 10 −3 (as noted in Fig. 6, the combination of interaction rate, density and velocity fixes the domain wall tension). All remaining parameters for dark matter and dark energy, as well as for the velocity distribution of clumps and strings, are chosen as in Fig. 4. as the frequency cut-offs used to determine the signalto-noise ratio (in addition to the other approximations we employed), these values could also be a factor of 10 or even higher. Bearing this in mind, our results are in rough agreement with [24]. An overview of the expected signal event rate for different masses of the clumps is shown in Fig. 10.
We obtain similar results for networks of topological defects. As an example, we show a cosmic string network with tension Gµ/c 2 ≈ 10 −28 , for which we find that LIGO could observe a corresponding signal with a signalto-noise ratio above one approximately every 20,000 years on average. Note that these values depend on the string tension and velocity distribution we assume. The signals induced by a domain wall network traversing the interferometer are similar to the ones we obtained at LISA. Quantitatively, we expect a possible signal with a signalto-noise ratio above one every 1,000 to 10,000 years on average.
In summary, LIGO can measure sufficiently short gravitational acceleration bursts caused by localized clumps of dark matter or networks of cosmic strings and domain walls. Due to LIGO's smaller size and correspondingly different frequency range compared to LISA, it is sensitive to objects of typically smaller masses as well as string and domain wall tensions. Indeed, with the same analysis strategy, further experimental improvements are required to lift the discovery potential of LIGO for these dark matter structures to a level that allows for a detection in a reasonable time frame.

B. Pulsar timing arrays
Another type of experiment, aiming towards a measurement of a gravitational waves, are pulsar timing arrays (PTAs). These exploit the fact that the time of arrival of radiation bursts from pulsars across the Universe can be precisely predicted. Gravitational waves that traverse the space between a certain set of pulsars and the observer on Earth will perturb the correlations between the different times of arrival, thereby generating a signal in the corresponding detector. Due to the long lines of sight between the pulsars and the telescopes, PTAs are able to probe gravitational wave spectra in very low frequency regimes, typically of the order 10 −9 Hz to 10 −6 Hz (see, e.g., [11-14, 122, 123]).
In this section we want to explore the discovery potential of PTAs with respect to the astrophysical structures discussed in the previous sections, i.e. dark matter clumps (note the previous works [17,18,20,21,23,[26][27][28]) as well as networks of cosmic strings and domain walls. For simplicity, we will do so by viewing a PTA as a gravitational wave interferometer of the same type as common ground-or space-based interferometers at very large scales, i.e. as an interferometer with arm lengths of several thousand light-years. As a prototypical example, we will consider a replica of LISA, that is an interferometer made up by three nodes at equal distance to each other, but with an arm length of L ≈ 1000 ly. However, the sensitivity of PTAs is usually based on the observation of multiple pulsars. To use this, the signal must exist in all Earth-pulsar combinations. For our search strategy this requires that the dark matter structure is affecting the velocity of the Earth. Compared to LISA, where all nodes can be treated equally, this reduces the detectable rate by a factor of 3. Although this is a crude simplification of the experimental techniques used for gravitational wave measurements with PTAs, we still expect reasonable estimates of prospective signal strengths for the purpose of this work. In particular, in terms of average signal event rates, we expect to benefit from the largely increased detector volume.
As we are essentially considering a LISA experiment at large scales, our analysis strategy is similar to the one presented in Sections II and III. The only major difference is the sensitivity of the experiments in different frequency regimes. We implement this by modelling the noise power spectrum according to the typical sensitivity curve of a PTA. As a prototype, we use the noise power spectrum of a PTA utilizing the future Square Kilometer Array (SKA) [14] as shown in [99], and extrapolate it to high frequencies.
Here, similar to (II.9), we again introduced a transfer function to adapt the raw strain spectrum to our signal spectrum. Furthermore, we added an additional factor of (2π) −1 in order to match the Fourier conventions. Due to the limited observation time of SKA, we introduce a cut-off for frequencies below f c ≈ 1.5 × 10 −9 Hz in the computation of the corresponding signal-to-noise ratio. A summary of our results is given in Fig. 9, where we show the average gravitational interaction rate of localized dark matter clumps as well as string and domain wall networks as a function of the signal-to-noise ratio at SKA. Note that here, similar to the analysis of LISA, we consider the close-approach regime of the localized dark matter clumps and cosmic strings. That is, we only consider events with impact parameters smaller than the size of the detector, D L. Because of the large size of the experiment, this condition is basically fulfilled in any region of the parameter space we are interested in.
For localized clumps of dark matter, we find that we benefit from the increased detector volume, allowing for very massive clumps while still retaining decent interaction rates. As expected, the qualitative behaviour we observe is similar to the detection of dark matter clumps with LISA shown in Fig. 2. Due to the increased detector volume, we expect a discovery potential for signals induced dark matter clumps with masses of M DM ≈ 10 −9 M at a signal-to-noise ratio of SNR ∼ 1 every 300 years on average. Apparently, this is a slightly more pessimistic estimate compared to what was found in earlier works. For example, in [18] the original estimate of a detection rate for primordial black hole dark matter with PTAs is somewhat higher. Nevertheless, bearing in mind the crude approximations we have applied in our derivation, our results are in rough agreement. We find a similar discovery potential for networks of cosmic strings as well as domain walls, constituting the dark matter or dark energy budget of the Universe, respectively. In particular, a cosmic string network featuring a larger string tension appears to benefit from the experimental setup at large scales. Here, we show an example of strings with a tension of Gµ/c 2 ≈ 10 −18 . For this network we expect a signal with a signal-to-noise ratio of about 1-10 every 1000 to 10,000 years on average. For a domain wall network with the same equation of state, we find an expected signal with a signal-to-noise ratio of 10 5 every few million years on average. We clearly do not need such a high signal-to-noise ratio. However, to preserve the validity of our analysis we require to have only a single domain wall traversing the detector volume at the same time. This, in turn, implies a lower bound on the surface tension of the domain wall network. As each data point shown is uniquely determined by a specific surface tension, this forbids us to go beyond the shown signal-to-noise ratios and corresponding interaction rates. Nevertheless, keeping these caveats in mind, we can extrapolate the given data. This raises the hope that we can have events every 10 to 100 years with a signal-to-noise ratio above one, SNR 1.
In summary, we find that, in a scenario where the dark matter properties optimally fit the experimental sensitivity, dark matter clumps, cosmic strings and domain walls, may be observed at SKA with a slightly higher rate than LIGO while being comparable to LISA. All categories benefit from the increased detector volume, such that PTAs are sensitive to more massive objects than smaller gravitational wave interferometers. We will give an overview of this feature in Section VI, which then also makes explicit the complementarity that exists between these different types of experiments. However, carefully note that the values we present in this section are all based on the assumption that a measurement of gravitational perturbations with a PTA is the same as with the LISA detector, but at large scales. The discovery potential of PTAs for strings and domain walls using a more accurate treatment of its experimental techniques certainly merits further investigation.

VI. CONCLUSIONS
In this paper, we have studied the potential of gravitational wave interferometers to measure gravitational perturbations caused by the presence of macroscopic dark matter objects in the vicinity of the detector. Any localized energy density passing by sufficiently close to the interferometer will exert a gravitational pull on each of its nodes, and hence cause a differential acceleration. This acceleration leads to Doppler shift signal in the interferometer that can be measured.
The objects we have considered in the present work include localized clumps of dark matter (see also [15][16][17][18][19][20][21][22][23][24][25][26][27][28][29]), topological defects, such as cosmic strings and domain walls, as well as stochastic fluctuations (cf. [20]) of the dark matter density. As our baseline example, we have examined the LISA experiment for which we have given the signal power spectrum associated to the presence to each of these sources in the vicinity of the detector. Based on this we then looked 11 at LIGO and a pulsar timing array using SKA that are sensitive to complementary frequency ranges.
Our results are summarized in Fig. 10. For each experiment, we show contours of the average gravitational 11 Note that our analysis strategy for LIGO and PTAs corresponds essentially to an extrapolation of LISA to different scales. Therefore, these estimates need to be taken with a bit of caution.
interaction rate of events with a signal-to-noise ratio greater than one, SNR 1, as a function of the energy density of each dark matter structure. Due to their different characteristic sizes, LIGO, LISA and SKA complement each other very well in the sense that they are sensitive to different masses and tensions of dark matter clumps and topological defects, respectively.
In the most sensitive regimes of each experiment, we find a prospective signal with a signal-to-noise ratio SNR 1, every 10,000 up to 100 years on average. Note that the most striking signature is a domain wall traversing LISA which might even be expected almost annually. However, we also note that a signal-to-noise ratio SNR 1 is in fact, at best, a minimum requirement for a signal to be detected. As the latter still needs to be distinguished from various other sources, a significantly higher signal-to-noise threshold is probably more realistic. As an example, we illustrate a detection threshold of SNR 10 by the dash-dotted lines in Fig. 10. While this is not yet on a desirable level for a near future discovery potential, let us remark that the signal estimates we show here are based on the close-approach approximation, imposing a relatively high low-frequency cutoff. Indeed, crucially, the signal-to-noise ratio can vary significantly, if a different low-frequency cutoff is taken into account (cf. Section II for a detailed discussion). This indicates that improvements of the experimental sensitivity as well as the theoretical analysis, notably in the low-frequency region, may lead to a sizable enhancement of the detection rate. Keeping this in mind, localized clumps of dark matter as well as cosmic strings and domain walls may still be within experimental reach of LIGO, LISA and PTAs. In contrast, as clumps of dark matter clumps naively constitute a critical overdensity of dark matter, stochastic fluctuations of the latter most likely cannot be measured above background noise. Overall, this clearly requires more sensitive future gravitational wave experiments.
In addition to the acceleration burst signals we have studied in this work, a signal could also arise due to the Shapiro effect [34], i.e. from the changing gravitational potential due to a dark matter structure within the line of sight connecting different nodes of a gravitational wave detector (see [17,20,23,26,28]). It would be interesting to investigate the corresponding signal associated to the presence of cosmic strings or domain walls, which we leave for future work.
In summary, not only localized clumps of dark matter but also cosmic strings or domain walls are close to the experimental reach of gravitational wave interferometers. Current and future gravitational wave experiments, such as LIGO, LISA and PTAs, are sensitive to gravitational perturbations due to the presence of these objects in the vicinity of the detector. These experiments are complementary to each other, as the different characteristic sizes and time-scales of the detectors make them sensitive to different parameter regions of the gravita- where we used a uniform average, here, the determination of the total rate takes into account the relative orientation between the domain walls and the experiment (see (III.12)). In all panels, a dashed line illustrates a naive extrapolation of our analysis. In (a) and (b) we extrapolate outside the close-approach regime by rescaling the differential acceleration with an additional tidal factor L/D (see (II.6)). In (c) we extrapolate the signal to a regime, where, on average, there is more than one domain wall traversing the detector volume at the same time.
tional sources. Already moderate improvements in the detector noise and analysis may yield interesting discovery potential to intriguingly exotic dark matter objects such as cosmic strings and domain walls. The shape and strength of a signal at LISA (or any other gravitational wave interferometer) induced by a gravitational perturbation due to a massive object in the vicinity of the detector, of course, depends on the distance, velocity and orientation of the latter with respect to the experiment. In this appendix, we aim to define the relevant geometrical quantities that enter the derivation of the signal power spectrum associated to these events.
In general, LISA is a space-based gravitational wave interferometer involving three distinct nodes, which are arranged in an equilateral triangle of side length L ≈ 2.5 × 10 6 km. Following the notation of [97], we schematically illustrate the general experimental setup of LISA in Fig. 11. Each pair of nodes exchanges laser beams, such that, in principle, there are six functions, U 1,2,3 and V 1,2,3 , that encode Doppler shifts due to gravitational perturbations by massive objects in the vicinity of the experiment. The desired detector response to an acceleration burst is then given by a suitable linear combination of these functions. Here, we use the so-called Michelson response function for the readout of a signal at a single detector node [96], where we have assumed that the arms of the interferometer are of equal length L. The components of the response function are given by projections of the velocity perturbations onto the interferometer arms, and cyclic permutations thereof. Here, the n i denote the unit vectors pointing between two nodes, labelled by the opposite side of the triangle, v i is the velocity perturbation of the i-th node induced by the gravitational pull and c is the speed of light. We note that other response function are also possible, see, e.g., [96,97]. From the response function (A.1) it is clear that any gravitational pull exerted by a massive object in the vicinity of the interferometer nodes has to be projected into the detector plane. For example, in the extreme case, where only one satellite is accelerated perpendicular to the detector plane, e.g. n 2 · v 1 = n 3 · v 1 = 0, the object will not leave any signature in the interferometer 12 . Therefore, the detector response will depend on how an object traverses the detector volume, i.e. on its orientation relative to the detector plane.
In the following, we want to define the relevant geometrical quantities describing this relative orientation of different macroscopic astrophysical objects we aim to probe with LISA. For simplicity, we will only consider the closeapproach limit, where the object passes by an interferometer node with an impact parameter smaller than the characteristic size of the detector, D L. As explained in the main text, in this case the gravitational perturbation of two of the three interferometer nodes can be 12 Strictly speaking, this is not true, because the two other nodes will also experience a gravitational pull which is not perpendicular to the detector plane. Nevertheless, at large distances between the source and these nodes this effect is negligible. 11. Schematic overview of the LISA interferometer.
Here, ni are unit vectors connecting pairs of satellites. The Ui and Vi encode the possible Doppler shifts of the laser beams that are exchanged between the detector nodes.
neglected, such that the detector response function can be approximated by adding suitable time delays to [31] where v i (t) is the velocity perturbation of the i-th node, that dominates compared to the two others. In other words, i is the interferometer node with smallest impact parameter with respect to the object traversing the detector volume. Hence, we do not sum over the indices in (A.4).

Spherical clumps
Massive spherical objects, such as localized clumps of dark matter, are in a sense the most symmetric configuration when they pass by one of the LISA nodes. That is, in the close-approach limit, their lack of an internal orientation allows to parametrize their motion relative to the detector by a velocity vector v and an impact parameter D, i.e. the closest distance in an encounter between the massive clump and a detector node. The former is completely determined by its magnitude v and an arbitrary direction given in terms of two angles, i.e. v = v (sin θ cos φ, sin θ sin φ, cos θ). Clearly, not all of these four parameters will enter the detector response. In fact, since the gravitational force between the spherical clump and the LISA satellite is only determined by their relative distance, we can consider the projection into the plane spanned by the satellite and the trajectory of the clump. This effectively removes two degrees of freedom, such that we are left with the relative velocity v and the impact parameter D. As a particular example discussed in the main text, we can choose the spherical clump to be in a straight uniform motion with velocity v parallel to the y-axis at an initial distance D to the satellite. The clump is furthermore confined to the xy-plane. We illustrate this scenario in the top panel of Fig. 12.
The above considerations determine the gravitational pull exerted by a massive clump on a single LISA satel-lite. As pointed out earlier, in order to determine the corresponding detector response, this gravitational acceleration burst has to be projected into the detector plane. The latter can be parametrized by an arbitrary unit normal vector, n T = (sin ϑ cos ϕ, sin ϑ sin ϕ, cos ϑ). In our analysis, the orientation of the detector plane is implemented in the detector response function, which in the close-approach regime, D L, is given by (A.4). That is, we can take this orientation into account by parametrizing the "dominant" unit vector of the LISA triangle accordingly, n i = (sin ϑ cos ϕ, sin ϑ sin ϕ, cos ϑ). Note that, strictly speaking, we are slightly abusing notation here. Obviously, the angles parametrizing the normal vector of the detector plane and the unit vector connecting two nodes of the triangle are not the same. Nevertheless, since we will average over these angles later, we denote them by the same symbol to avoid an overload of notation.
In summary, in the close-approach limit, there are four geometrical degrees of freedom in total that enter the detector response to a localized massive clump travelling through the interferometer. In particular, the clump's velocity v, the impact parameter D as well as the orientation of the detector plane (ϑ, ϕ) completely determine the signal at LISA. That is, the detector response is a function of these parameters, X(t) = X(t, v, D, ϑ, ϕ).
Finally, since we assume a locally isotropic situation, i.e. the clumps can approach the interferometer from any direction equally likely, we uniformly average over the orientation of the detector plane in order to obtain the signal power spectrum from the detector response, Overall, the signal power spectrum then still depends on the velocity of the clump as well as the impact parameter of the encounter. Let us close this discussion with a few words of caution. Strictly speaking, the uniform average we have employed above, is not fully justified. This is because the configuration we consider is not strictly isotropic. Instead, there is a preferred direction in the system, given by the Sun, together with the detector, moving through the Universe. In this sense, a uniform average is only an approximation. We present a more detailed discussion of this in Appendix A 4. by line elements, thereby having an additional orientation themselves. Obviously, when determining the detector response to a gravitational acceleration caused by a cosmic string, this orientation has to be taken into account.
In general, an infinite string can be parametrized by a straight line, γ(s) = x 0 + sn γ , where n γ denotes an arbitrary unit vector. In addition, the string can move in a direction with a certain velocity v relative to the nodes of the interferometer. Again, similar to the case of spherical clumps, not all of these parameters will enter the detector response to a gravitational perturbation in the close-approach regime. In fact, the gravitational field of the string only depends on the radial distance to the source, such that we can consider the projection into the plane perpendicular to the string. As in the main text, we can choose a coordinate frame, in which the infinite string is parallel to the z-axis and the satellite located at the origin is, initially at t 0 = 0, at a minimum distance D to the string. We can then assume that the string is uniformly moving in a random direction in the yz-plane with velocity v, i.e. v y = v sin θ and v z = v cos θ, respectively. Indeed, this reflects the fact that the string has an additional internal orientation as compared to spherical objects such as clumps. This situation is depicted in the top panel of Fig. 12. Therefore, the gravitational acceleration burst induced by a cosmic string on a single LISA node depends on three parameters, namely the relative velocity v, the impact parameter D as well as the direction of motion relative to the string orientation, parametrized by θ. Note that, equivalently, we could also choose a reference frame where the satellite is moving uniformly and the satellite is at rest.
Finally, as pointed out in the previous section, the overall gravitational acceleration has to be projected into the detector plane, parametrized by the angles ϑ and ϕ. Therefore, in summary, LISA's detector response to a gravitational pull by cosmic string (in the close-approach regime) is a function of five geometrical parameters in total, X(t) = X(t, v, D, θ, ϑ, ϕ). In an isotropic Universe, the signal power spectrum is finally given by a uniform average over all arbitrary orientations involved, In total, the signal power spectrum depends on the velocity of the string as well as the impact parameter of the encounter. However, we note that, similar to the case of dark matter clumps, a uniform average over all possible directions might not be fully justified, see Appendix A 4.

Domain walls
When determining the detector response of LISA to a gravitational potential sourced by an energy density localized on an infinite plane, i.e. a domain wall, additional degrees of freedom compared to a spherical clump or an infinite string have to be taken into account.
Geometrically, the plane parametrizing a domain wall can be described by the algebraic equation n · (r − r 0 ) = 0, where r 0 is an arbitrary point in the plane and n denotes the unit vector normal to it. Nevertheless, as we will see momentarily, the exact signal shape caused by a domain wall involves fewer geometrical parameters than, e.g., spherical objects or cosmic strings. This is due to the fact that its gravitational field only induces a signal at LISA, if the domain wall is located in between the detector nodes, thereby separating them from each other 13 . Therefore, we only have to consider a situation, where the triangle spanned by the LISA satellites intersects an infinite plane. If the domain wall, or equivalently LISA, is 13 As discussed in the main text, this is because LISA only measures differential accelerations between the satellites. However, the gravitational field sourced by a domain wall does not depend on the distance, but is constant everywhere. Therefore, signals are only generated for configurations where the satellites are accelerated into opposite directions.
moving at a certain velocity, this line of intersection will move, too, until it has completely passed the detector plane. We illustrate this in the bottom panel of Fig. 12.
The only geometrical quantities that enter the detector response function X(t) in this scenario are, in fact, the orientation of the domain wall with respect to the triangle spanned by the LISA satellites as well as its relative velocity. Without loss of generality, the former can be completely parametrized by, e.g., the unit vector normal to the plane, n = (sin θ cos φ, sin θ sin φ, cos θ), while we can assume the latter to point into the normal direction, v = vn, (cf. Fig. 12). Accordingly, the detector response will be a function of these parameters only, X(t) = X(t, v, θ, φ). Finally, similar to the previous sections, in a locally isotropic dark matter distribution, the signal power spectrum associated to the gravitational perturbation by a domain wall traversing the detector volume is given by the uniform average over all possible orientations relative to the detector, That means, the overall signal power spectrum will only depend on the velocity of the domain wall relative to the LISA detector. However, we also note that, similar to the case of dark matter clumps, a uniform average over all possible directions might not be fully justified, as we will discuss in the following subsection.

Velocity distribution of dark matter
In the previous subsections we have illustrated how our estimate of the signal power spectrum accounts for the relative orientation between the detector plane of LISA and the source of the gravitational perturbation. In particular, we have assumed that the dark matter can approach the interferometer from any direction equally likely and hence uniformly averaged over the solid angle which parametrizes the latter (see, e.g., (A.5)). Naively, this partly follows from the Maxwell-Boltzmann distributed velocities of the dark matter structures. However, strictly speaking, this is not fully justified for the following reason.
In a naive approximation, it is reasonable to assume that the dark matter inside the halo surrounding our Galaxy behaves like an ideal gas of non-interacting particles and therefore roughly follows a Maxwell-Boltzmann distribution (see, e.g., [124]), for some normalization v 0 , which, from a microscopic perspective, is determined by the dark matter mass and the temperature of the gas. The Maxwell-Boltzmann distribution (A.8) of the dark matter inside the halo of our Galaxy yields an isotropic uniform distribution for the direction in which the dark matter structures are moving. That is, dark matter can approach the experiment from any direction equally likely. In the previous sections, this feature is taken into account by a uniform average of the angles parametrizing the relative orientation between the detector plane of LISA and the trajectory of the dark matter structure (see, e.g., (A.5)). Obviously, this is true in an isotropic reference frame where an observer is at rest inside the dark matter halo of the Galaxy. However, in practice, the Sun, together with the detector, is moving through the halo at a constant velocity of v ≈ 220 km/s [102], thereby imposing a preferred direction on the system. That means that not every direction in an encounter occurs equally likely, such that the average over these directions should not be uniform. Instead, to correctly account for this one would need to weight the velocity in each direction according to the normal distribution (A.8) with the appropriate velocity shift by v , for example in At first glance the situation looks even worse, because, in addition, LISA is moving on a complicated orbit around the Sun (see, e.g., Fig. 4 in [98]). However, this composite motion of the detector might turn out to be a blessing in disguise [31], as it does not impose a single preferred direction but (at least naively) periodically changes the latter. Hence, in order to account for the relative orientation between the experiment and the dark matter trajectory, a uniform average average over the orientation might indeed be closer to the experimental scenario than singling out only one preferred direction [31]. In practice, as an approximation, we therefore take a uniform average over the solid angle accounting for the direction (see, e.g., (A.5)) and weight the detector response according to the probability distribution Here, we try to approximate the motion of the detector through the dark matter halo with v ≈ 220 km/s [102] and finally normalize to the dark matter rms velocity of the latter, v 0 = v rms / √ 3 with v rms ≈ 270 km/s (see, e.g., [103]). Let us remark that, for the purpose of this work, we do not expect any large quantitative changes if a more accurate estimation of the dark matter velocity distribution with respect to the detector was performed.