Can CP be conserved in the two-Higgs-doublet model?

We study the conditions under which the CP violation in the quark mixing matrix can leak into the scalar potential of the real two-Higgs-doublet model (2HDM) via divergent radiative corrections, thereby spoiling the renormalizability of the model. We show that any contributing diagram must involve 12 Yukawa-coupling insertions and a factor of the $U(1)_{PQ}$-breaking scalar potential parameter $λ_5$, thereby requiring at least six loops; this also implies that the 2HDM with only softly-broken $U(1)_{PQ}$ is safe from divergent leaks of CP violation to all orders. We demonstrate that additional symmetries of the six-loop diagrams in the type I and II 2HDMs guarantee that all of the divergent CP-violating contributions cancel at this order. We also show that these symmetries are violated at seven loops and enumerate the classes of diagrams that can contribute to CP-violating divergences, providing evidence that the real 2HDM is theoretically inconsistent starting at the seven-loop level.