B_(s,d)⁰ to ell^+ell^- Decays in Two-Higgs Doublet Models

We study the rare leptonic decays $B_{s,d}^0 \to \ell^+\ell^-$ within the general framework of the aligned two-Higgs doublet model. A complete one-loop calculation of the relevant short-distance Wilson coefficients is presented, with a detailed technical summary of the results. The phenomenological constraints imposed by present data on the model parameters are also investigated.