
2
The coupling S6= 0 gives rise to a force on a SM
object moving at speed v:
F=εSeemQS−∇V0−∂tV+v×(∇×V),(2)
where QSis the U(1)Scharge of the SM object.
If the Vµgauge boson comprises all of the local, non-
relativistic DM (speed vdm ∼10−3), we write Vµ≡
(V0,V) exp [−iωdmt+ikdm ·x+iα] within a single co-
herence time/length,5with αan arbitrary phase, ωdm ≈
mdmp1 + v2
dm ≈mdm, and |kdm| ≈ mdmvdm. Note that
we have identified mV≡mdm. Now, ∂µVµ= 0 implies
that |V0| ∼ vdm|V|. Moreover, we have6|∇×V| ∼
mdmvdm|V|and |∇V0| ∼ mdmvdm|V0| ∼ mdmv2
dm|V|.
Additionally, |∂tV| ∼ mdm|V|. Since, by assumption,
T00
V=ρdm, it follows that7|V| ∼ √2ρdm/mdm.
Considering objects gravitationally bound to the Solar
System, we also have |v|.vdm, owing to the motion
of the Solar System relative to the galactic rest frame.
Therefore, the force can be written as
F≈iεSeemQSmdmVe−imdm(t−vdm·x)+iα
+O(v2
dm, vdm|v|)×εSeemQSmdm|V|.(3)
If the relevant SM object has mass M, and we drop the
subleading corrections,8this force causes an oscillatory
acceleration of that object [19,21,24,25]:
a≈2εSp2παemρdm
QS
Me−imdm(t−vdm·x)+iφ ˆ
V,(4)
where
ˆ
Vgives the polarization state of the DM (
ˆ
V·ˆ
V∗=
1), and we have absorbed a phase into φ.
Consider now a GW detector that operates by mea-
suring fluctuations in the proper distance between two
or more inertial test masses that define the endpoints of
one or more detector baselines. The DM-induced accel-
eration Eq. (4) causes the TMs to oscillate, leading to
a modulation of the proper length of the detector base-
line(s) that manifests as an oscillatory strain component.
5The coherence time is Tcoh ∼2π/(mdmv2
dm) and the coherence
length is λcoh ∼2π/(mdmvdm) [i.e., the de Broglie wavelength].
6We clarify that these are order of magnitude estimates for the
largest that the expressions on the lhs can be; additional geo-
metrical factors can suppress these even further.
7Technically this is only true as an average statement over many
coherence times of the DM field, as the field amplitude executes
O(1) stochastic fluctuations from one coherence time to the next
(see, e.g., Refs. [36–38]); nevertheless, it is a good figure of merit.
8The vectorial orientation of the force relative to a GW detec-
tor baseline is the relevant quantity to consider for GW detector
effects. Because the subleading corrections can be oriented dif-
ferently than the leading term, it is possible that the leading
term gives no effect, but that the subleading terms do. However,
because the subleading corrections are suppressed by at least
∼v2
dm, this can only occur in highly tuned orientations. Real
GW detector baselines evolve in orientation over time relative
to inertial space, which will always spoil the tuning required for
such cancellation to be maintained; see also Appendix B. The
subleading terms can thus always be dropped.
There are two contributions to this strain: (1) a term
arising from the finite light-travel time for a null signal
propagating between the TMs, present even if the acceler-
ation had no spatial dependence at all (i.e., vdm = 0) [18–
20,22,24,25], and (2) a term arising from the spatial
gradient of the acceleration [19,21–25].
Averaged over time, DM field orientations with respect
to the baseline, and DM momentum directions, and ig-
noring further corrections at O(v2
dm), the mean square of
the strain signal h(t)≡∆t/(2L), where ∆tis the change
to the round-trip light travel time between the TMs, can
be written as [24]:9
hh2i=hh2
1i+hh2
2i; (5)
hh2
1i ≡ 4cgeom
1H2×sin41
2mdmL; (6)
hh2
2i ≡ 1
3cgeom
2H2×(mdmLvdm)2; (7)
H2≡8π
3ε2
Sαem
ρdm
m4
dmL2QS
M2
,(8)
where cgeom
jare O(1) geometrical factors that depend on
the baseline orientation, and Lis the unperturbed base-
line length. The baseline is assumed to be approximately
fixed at least on the ∼2Lround-trip light travel time be-
tween the TMs, but it can vary secularly both in length
and orientation on longer timescales. The expression
for hh2
2iat Eq. (7) is correct in the limit mdmLvdm 1;
this is satisfied in all cases that we consider. We dis-
cuss the appropriateness of averaging over DPDM field
orientations and momentum directions at length in Ap-
pendix B, as well as a relevant correction we apply when
we find this to be incompletely justified.
Note that hh2
2i/hh2
1i ∝ (mdmLvdm)2cosec41
2mdmL.
Because mdmLvdm 1, the h2signal is only dominant
when either (1) the baseline is much shorter than the
DM Compton wavelength mdmL.vdm (never the case
for interesting mass ranges for the detectors we consider),
or (2) in certain narrow mass ranges where the round-trip
light travel time is such that the h1signal is nearly zero:
|mdmL−2πn| ∼ 2√mdmLvdm for n= 1,2, . . .. The latter
occurs only at frequencies above the peak sensitivity for
every detector we consider.
The geometrical factors for a single-baseline detector
are given by
cgeom
1=cgeom
2=1
2,[single baseline] (9)
9For detectors that are constructed with Fabry-P´erot (FP) cav-
ities in the baseline arms, these expressions are referred to the
input in the sense that the FP transfer function has been omitted.
Published GW characteristic-strain curves for detectors with FP
cavities are similarly referred to the GW input. Since the FP
cavity transfer function is the same for a GW and the DPDM
case that we consider in this work, no correction is required for
this [24]; see also Ref. [39] for further discussion of FP cavities
in GW detectors.