We show that there is a constraint on the parameter space of two Higgs doublet models that comes from the existence of the stable vortex-domain wall systems. The constraint is quite universal in the sense that it depends on only two combinations of Lagrangian parameters and does not depend on how fermions couple to two Higgs fields. Numerical solutions of field configurations of domain wall-vortex system are obtained, which provide a basis for further quantitative study of cosmology which involve such topological objects.
Contrary to the standard model that does not admit topologically nontrivial solitons, two Higgs doublet models admit topologically stable vortex strings and domain walls. We numerically confirm the existence of a topological $Z$-string confining fractional $Z$-flux inside. We show that topological strings at $sintheta_W = 0$ limit reduce to non-Abelian strings which possess non-Abelian moduli $S^2$ associated with spontaneous breakdown of the $SU(2)$ custodial symmetry. We numerically solve the equations of motion for various parameter choices. It is found that a gauging $U(1)_Y$ always lowers the tension of the $Z$-string while it keeps that of the $W$-string. On the other hand, a deformation of the Higgs potential is either raising or lowering the tensions of the $Z$-string and $W$-string. We numerically obtain an effective potential for the non-Abelian moduli $S^2$ for various parameter deformations under the restriction $tanbeta=1$. It is the first time to show that there exists a certain parameter region where the topological $W$-string can be the most stable topological excitation, contrary to conventional wisdom of electroweak theories. We also obtain numerical solutions of composites of the string and domain walls in a certain condition.
The Two Higgs Doublet Model predicts the emergence of 3 distinct domain wall solutions arising from the breaking of 3 accidental global symmetries, $Z_2$, CP1 and CP2, at the electroweak scale for specific choices of the model parameters. We present numerical kink solutions to the field equations in all three cases along with dynamical simulations of the models in (2+1) and (3+1) dimensions. For each kink solution we define an associated topological current. In all three cases simulations produce a network of domain walls which deviates from power law scaling in Minkowski and FRW simulations. This deviation is attributed to a winding of the electroweak group parameters around the domain walls in our simulations. We observe a local violation of the neutral vacuum condition on the domain walls in our simulations. This violation is attributed to relative electroweak transformations across the domain walls which is a general feature emerging from random initial conditions.
The Two Higgs Doublet Model (2HDM) with spontaneously broken $Z_2$ symmetry predicts a production of domain walls at the electroweak scale. We derive cosmological constraints on model parameters for both Type-I and Type-II 2HDMs from the requirement that domain walls do not dominate the Universe by the present day. For Type-I 2HDMs, we deduce the lower bound on the key parameter $tanbeta > 10^5$ for a wide range of Higgs-boson masses $sim$ 100 GeV or greater close to the Standard Model alignment limit. In addition, we perform numerical simulations of the 2HDM with an approximate as well as an exact $Z_2$ symmetry but biased initial conditions. In both cases, we find that domain wall networks are unstable and, hence, do not survive at late times. The domain walls experience an exponential suppression of scaling in these models which can help ameliorate the stringent constraints found in the case of an exact discrete symmetry. For a 2HDM with softly-broken $Z_2$ symmetry, we relate the size of this exponential suppression to the soft-breaking bilinear parameter $m_{12}$ allowing limits to be placed on this parameter of order $mu$eV, such that domain wall domination can be avoided. In particular, for Type-II 2HDMs, we obtain a corresponding lower limit on the CP-odd phase $theta$ generated by QCD instantons, $theta stackrel{>}{{}_sim} 10^{-11}/(sinbeta cosbeta)$, which is in some tension with the upper limit of $theta stackrel{<}{{}_sim} 10^{-11}$--$10^{-10}$, as derived from the non-observation of a non-zero neutron electric dipole moment. For a $Z_2$-symmetric 2HDM with biased initial conditions, we are able to relate the size of the exponential suppression to a biasing parameter $varepsilon$ so as to avoid domain wall domination.
We show that a topological Nambu monopole exists as a regular solution for a large range of parameters in two Higgs doublet models, contrary to the standard model admitting only non-topological Nambu monopoles. We analyze a Higgs potential with a global $U(1)$ symmetry and a discrete symmetry $mathbb{Z}_2$. The monopole is attached by two topological $Z$ strings ($Z$ flux tubes) from both sides. Despite of a trivial second homotopy group, the discrete symmetry $mathbb{Z}_2$ together with a non-trivial first homotopy group for $Z$ strings topologically ensures the topological stability. After analytically constructing an asymptotic form of such a configuration, we explicitly construct a solution of the equation of motion based on a 3D numerical simulation, in which magnetic fluxes spherically emanating from the monopole at large distances are deformed in the vicinity of the monopole. Since the monopoles are expected to be abundant in the present universe, they might be observed in the current monopole searches.
We study Two-Higgs-Doublet Models (2HDM) where Abelian symmetries have been introduced, leading to a drastic reduction in the number of free parameters in the 2HDM. Our analysis is inspired in BGL models, where, as the result of a symmetry of the Lagrangian, there are tree-level scalar mediated Flavour-Changing-Neutral-Currents, with the flavour structure depending only on the CKM matrix. A systematic analysis is done on the various possible schemes, which are classified in different classes, depending on the way the extra symmetries constrain the matrices of couplings defining the flavour structure of the scalar mediated neutral currents. All the resulting flavour textures of the Yukawa couplings are stable under renormalisation since they result from symmetries imposed at the Lagrangian level. We also present a brief phenomenological analysis of the most salient features of each class of symmetry constrained 2HDM.