
FERMILAB-PUB-22-670-T
Freezing In Vector Dark Matter Through Magnetic Dipole Interactions
Gordan Krnjaic,1, 2, 3, ∗Duncan Rocha,4, †and Anastasia Sokolenko1, 3, ‡
1Theoretical Physics Department, Fermi National Accelerator Laboratory, Batavia, IL 60510
2Department of Astronomy and Astrophysics, University of Chicago, Chicago, IL 60637
3Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637
4Department of Physics, University of Chicago, Chicago, IL 60637
(Dated: October 14, 2022)
We study a simple model of vector dark matter that couples to Standard Model particles via
magnetic dipole interactions. In this scenario, the cosmological abundance arises through the freeze-
in mechanism and depends on the dipole coupling, the vector mass, and the reheat temperature.
To ensure cosmological metastability, the vector must be lighter than the fermions to which it
couples, but rare decays can still produce observable 3γfinal states; two-body decays can also
occur at one-loop with additional weak suppression, but are subdominant if the vector couples
mainly to light fermions. For sufficiently heavy vectors, induced kinetic mixing with the photon
can also yield additional two body decays to lighter fermions and predict indirect detection signals
through final state radiation. We explore the implications of couplings to various flavors of visible
particles and emphasize leptophilic dipoles involving electrons, muons, and taus, which offer the
most promising indirect detection signatures through 3γ,e+e−γ, and µ+µ−γdecay channels. We
also present constraints from current and past telescopes, and sensitivity projections for future
missions including e-ASTROGAM and AMEGO.
I. INTRODUCTION
While the evidence for the existence dark matter (DM)
is overwhelming, its microscopic properties remain elu-
sive (see [1] for a historical review). Since there are few
clues about its non-gravitational interactions, it is cur-
rently not known how DM was produced in the early uni-
verse or when that production took place. Thus, there is
great motivation to identify and test all predictive mech-
anisms for this key epoch in the history of the universe.
Cosmological “freeze-in” is among the simplest and
most predictive DM production mechanisms [2,3]. In
this scenario, DM is initially not present at reheating
when the hot radiation bath is first established. Rather,
its density builds up through ultra-feeble interactions
with Standard Model (SM) particles and production
halts when this process becomes Boltzmann suppressed.
Since these reaction rates are sub-Hubble, the DM never
equilibrates with the SM, so unlike freeze-out, there is no
need to deplete the large thermal entropy with additional
DM annihilation when equilibrium is lost. Freeze-in pro-
duction ends when the temperature of the universe cools
below either the mass of either the DM or the SM species
to which it couples, whichever is greater.
It is well known that dark photons A0can be produced
via freeze-in through a kinetic mixing interaction with
the SM photon [4]. Since A0are also unstable and decay
through this same interaction, only the mA0<2memass
∗krnjaicg@fnal.gov
†drocha@uchicago.edu
‡sokolenko@kicp.uchicago.edu
range can be cosmologically metastable to provide a DM
candidate. In this range, the DM decays via A0→3γre-
actions and predicts a late time X-ray flux uniquely spec-
ified by the A0mass, once the kinetic mixing parameter
is fixed to obtain the observed DM abundance. However,
this tight relationship between abundance and flux has
been used to sharply constrain this simple model with ob-
servations of X-ray lines, extragalactic background light,
and direct detection via absorption [4,5]. Collectively,
these probes have eliminated nearly all viable parame-
ter space for vector DM produced through freeze-in via
kinetic mixing interactions.
In this paper we generalize dark photon freeze-in to
allow for the possibility that its main interaction with
visible matter is a magnetic dipole coupling to charged
fermions, instead of kinetic mixing with the photon.
Since the dipole operator has mass dimension 5, the
freeze-in abundance is sensitive to the reheat tempera-
ture. Thus, unlike kinetic mixing, there is a parametric
separation between the production rate at early times
and the decay rate at late times; the former depends on
the reheat temperature and the latter does not, so it is
possible to achieve the observed DM abundance with a
much feebler coupling to SM particles and, thereby, open
up viable parameter space for direct and indirect detec-
tion.
This paper is organized as follows: in Sec. II we de-
scribe the model, in Sec. III we calculate the A0abun-
dance via freeze in, in Sec. IV we present structure for-
mation limits, in Sec. Vwe explore the indirect detection
constraints and future projections for this model, and in
Sec. VI we offer some concluding remarks.
arXiv:2210.06487v1 [hep-ph] 12 Oct 2022