Wilson coefficient of the effective Hamiltonian as defined in Eq. (\[NPH\]) is \[o2coeff\] C\^[NP]{}(m\_t) = . where $\hat{m}_S^2 = \hat{m}_1^2 + \left(\hat{\lambda}_1 + \hat{\lambda}_2\right) v^2/4$. This leads to the mass bound \[upperlimitmass\] \_S\^2 && (98 [GeV]{})\^4 |’\_D|\^2 \_3 + (174 [GeV]{})\^4 \^2\_3. in terms of the parameters defined in the potential. The mass spectrum of the new doublet is given by m\_[S\^]{}\^2 &=& \_S\^2 - \_2 ,\
m\_[S\_R\^0]{}\^2 &=& \_S\^2 + \_3 ,\
m\_[S\_I\^0]{}\^2 &=& \_S\^2 - \_3 , We show the masses of the neutral scalars for this model in Fig.7.

This model has eight new scalar doublets. Nevertheless, precision electroweak constraints are satisfied (when the Higgs is fixed to be $m_h = 96^{+29}_{-24} \, {\rm GeV}$ for $m_S \gtrsim 100 \, {\rm GeV}$) when $|m_I - m_\pm| < 50 \, {\rm GeV}$ [@Burgess:2009wm]. Conversely, custodial $\rm SU(2)$ violation in such a light scalar doublet leading to a positive contribution to $\Delta T$ can raise the allowed mass of the Higgs in EWPD [@Peskin:2001rw; @Burgess:2009wm].

In this model the coupling constant analogous to $\eta_U$ is naturally of order $(m_b/v)^2 \sim 10^{-3}$ and the $\bar{B} \rightarrow X_S \, \gamma$ and neutron EDM effects of the model are suppressed as phenomenologically 