  We study the phenomenology of electric dipole moments (EDMs) induced in
various scalar leptoquark models. We consider generic leptoquark couplings to
quarks and leptons and match to Standard Model effective field theory. After
evolving the resulting operators to low energies, we connect to EDM experiments
by using up-to-date hadronic, nuclear, and atomic matrix elements. We show that
current experimental limits set strong constraints on the possible CP-violating
phases in leptoquark models. Depending on the quarks and leptons involved in
the interaction, the existing searches for EDMs of leptons, nucleons, atoms,
and molecules all play a role in constraining the CP-violating couplings. We
discuss the impact of hadronic and nuclear uncertainties as well as the
sensitivities that can be achieved with future EDM experiments. Finally, we
study the impact of EDM constraints on a specific leptoquark model that can
explain the recent $B$-physics anomalies.
