CP violation in multi-Higgs supersymmetric models

We consider supersymmetric extensions of the standard model with two pairs of Higgs doublets. We study the possibility of spontaneous CP violation in these scenarios and present a model where the origin of CP violation is soft, with all the complex phases in the Lagrangian derived from complex masses and vacuum expectation values (VEVs) of the Higgs fields. The main ingredient of the model is an approximate global symmetry, which determines the order of magnitude of Yukawa couplings and scalar VEVs. We assume that the terms violating this symmetry are suppressed by powers of the small parameter ε_PQ = O(m_b/m_t). The tree-level flavor-changing interactions are small due to a combination of this global symmetry and a flavor symmetry, but they can be the dominant source of CP violation. All CP-violating effects occur at order ε_PQ² as the result of exchange of almost decoupled extra Higgs bosons and/or through the usual mechanisms with an almost real CKM matrix. On dimensional grounds, the model gives ε_K ≈ ε_PQ² and predicts for the neutron electric dipole moment (and possibly also for ε_K¹) a suppression of order ε_PQ² with respect to the values obtained in standard and minimal supersymmetric scenarios. The predicted CP asymmetries in B decays are generically too small to be seen in the near future. The mass of the lightest neutral scalar, the strong CP problem, and possible contributions to the Z decay into b quarks (the R_b puzzle) are also briefly addressed in the framework of this model.