CP in the dark

We build a model containing two scalar doublets and a scalar singlet with a specific discrete symmetry. After spontaneous symmetry breaking, the model has Standard Model-like phenomenology, as well as a hidden scalar sector which provides a viable dark matter candidate. We show that CP violation in the scalar sector occurs exclusively in the hidden sector, and consider possible experimental signatures of this CP violation. In particular, we study contribution to anomalous gauge couplings from the hidden scalars.


JHEP11(2018)091
sources of CP violation exist in the scalar sector of the model, but only in the "dark" sector. This hidden CP violation will mean that the SM-like scalar, h, behaves almost exactly like the SM Higgs boson, and in particular (unless contributions from a high number of loops are considered) h has couplings to gauge bosons and fermions which are exactly those of a scalar. This is all the more remarkable since the CP violation of the proposed model is explicit. The extra particle content of the model, as advertised, is simpler than the model of [31], consisting of two Higgs doublets (both of hypercharge Y = 1) and a real singlet (Y = 0). This is sometimes known as the Next-to-2HDM (N2HDM), and was the subject of a thorough study in [36]. The N2HDM version considered in this paper uses a different discrete symmetry than the symmetries considered in [36], designed, as will be shown, to produce both dark matter and dark CP violation. The paper is organised as follows: in section 2 we will introduce the model, explaining in detail its construction and symmetries, as well as the details of spontaneous symmetry breaking that occurs when one of the fields acquires a vacuum expectation value (vev). In section 3 we will present the results of a parameter space scan of the model, where all existing constraints -both theoretical and experimental (from colliders and dark matter searches) -are taken into account; deviations from the SM behaviour of h in the diphoton channel, stemming from the existence of a charged scalar, will be discussed, as will the contributions of the model to dark matter observables; in section 4 we will discuss how CP violation arises in the dark sector, and how it might have a measurable impact in future colliders. Finally, we conclude in section 5.

The scalar potential and possible vacua
For our purposes, the N2HDM considered is very similar to that discussed in ref. [36], in that the fermionic and gauge sectors are identical to the SM and the scalar sector is extended to include an extra doublet and also a singlet scalar field -thus the model boasts two scalar doublets, Φ 1 and Φ 2 , and a real singlet Φ S . As in the 2HDM, we will require that the Lagrangian be invariant under a sign flip of some scalar fields, so that the number of free parameters of the model is reduced and no tree-level flavour-changing neutral currents (FCNC) occur [37,38]. The difference between the current work and that of [36] consists in the discrete symmetry applied to the Lagrangian -here, we will consider a single Z 2 symmetry of the form With these requirements, the most general scalar potential invariant under SU(2) × U(1) is given by

JHEP11(2018)091
where, with the exception of A, all parameters in the potential are real. As for the Yukawa sector, we consider all fermion fields neutral under this symmetry. As such, only the doublet Φ 1 couples to fermions, and the Yukawa Lagrangian is therefore where we have only written the terms corresponding to the third generation of fermions, with the Yukawa terms for the remaining generations taking an analogous form. The lefthanded doublets for quarks and leptons are denoted by Q L and L L , respectively; t R , b R and τ R are the right-handed top, bottom and τ fields; andΦ 1 is the charge conjugate of the doublet Φ 1 . Notice that since the two doublets have the same quantum numbers and are not physical (only the mass eigenstates of the model will be physical particles), the potential must be invariant under basis changes on the doublets. This is a well-known property of 2HDMs, which the N2HDM inherits: any unitary transformation of these fields, Φ i = U ij Φ j with a 2 × 2 unitary matrix U , is an equally valid description of the theory. Though the theory is invariant under such transformations, its parameters are not and undergo transformations dependent on U . A few observations are immediately in order: • Since only Φ 1 has Yukawa interactions it must have a vev to give mass to all charged fermions. 1 In fact the Yukawa sector of this model is identical to the one of the SM, and a CKM matrix arises there, as in the SM.
• The fact that all fermions couple to a single doublet, Φ 1 , automatically ensures that no scalar-mediated tree-level FCNC occur, as in the 2HDM with a Z 2 symmetry [37,38].
• The Z 2 symmetry considered eliminates many possible terms in the potential, but does not force all of the remaining ones to be real -in particular, both the quartic coupling λ 5 and the cubic one, A, can be a priori complex. However, using the basis freedom to redefine doublets, we can absorb one of those complex phases into, for instance, Φ 2 . We choose, without loss of generality, to render λ 5 real.
A complex phase on A renders the model explicitly CP violating. Considering, for instance, the CP transformation of the scalar fields given by we see that such a CP transformation, to be a symmetry of the potential, would require all of its parameters to be real. Notice that the CP transformation of the singlet trivially does not involve complex conjugation as Φ S is real. In fact, this is a well-known CP property of singlet fields [39]. One point of caution is in order: the complex phase of A is not invariant under the specific CP transformation of eq. (2.4), but by itself that does not prove that the model is explicitly CP violating. In fact, one could consider some form of generalized CP (GCP) transformation involving, other than complex conjugation of the fields, also doublet 1 And neutrinos as well, if one wishes to consider Dirac mass terms for them. The model can only be said to be explicitly CP  violating if there does not exist any CP transformation under which it is invariant. So, conceivably, though the model breaks the CP symmetry defined by the transformation of eq. (2.4), it might be invariant under some other one. The point is moot, however: as we will see ahead, the vacuum of the model we will be considering is invariant under the CP transformation of eq. (2.4) (and the Z 2 symmetry of eq. (2.1)), but the theory has CP violation. Thus the CP symmetry was broken to begin with, and hence the model is explicitly CP violating.
Let us consider now the possibility of spontaneous symmetry breaking in which only the Φ 1 doublet acquires a neutral non-zero vev: Given the structure of the potential in eq. (2.2), the minimisation conditions imply that this is a possible solution, with all scalar components of the doublets (except the real, neutral one of Φ 1 ) and the singlet equal to zero, provided that the following condition is obeyed: Since all fermion and gauge boson masses are therefore generated by Φ 1 , it is mandatory that v = 246 GeV as usual. At this vacuum, then, it is convenient to rewrite the doublets in terms of their component fields as where h is the SM-like Higgs boson, with interaction vertices with fermions and gauge bosons identical to those expected in the SM (the diphoton decay of h, however, will differ from its SM counterpart). The mass of the h field is found to be given by and since m h = 125 GeV, this fixes the value of one of the quartic couplings, λ 1 0.258. The neutral and charged Goldstone bosons G 0 and G + , respectively, are found to be massless as expected, and the squared mass of the charged scalar H + is given by Finally, the two neutral components of the doublet Φ 2 , ρ and η, mix with the singlet component Φ s ≡ s yielding a 3 × 3 mass matrix, withλ 345 = λ 3 + λ 4 − λ 5 and λ 345 = λ 3 + λ 4 + λ 5 . There are therefore three neutral scalars other than h, which we call h 1 , h 2 and h 3 , in growing order of their masses. This mass matrix can then be diagonalized by an orthogonal rotation matrix R, such that

JHEP11(2018)091
and the connection between the original fields and the mass eigenstates is given by The rotation matrix R can be parametrized in terms of three angles, α 1 , α 2 and α 3 , such that where for convenience we use the notation c i = cos α i , s j = sin α j . Without loss of generality, we may take the angles α i in the interval [−π/2 , π/2]. In the following we discuss several phenomenological properties of this model. The vacuum preserves the Z 2 symmetry. As a result, the physical eigenstates emerging from Φ 2 and Φ S , i.e. the charged scalar H ± and the neutral ones h 1 , h 2 and h 3 , carry a quantum number -a "dark charge" equal to −1 -which is preserved in all interactions, to all orders of perturbation theory. In the following we refer to these four eigenstates as "dark particles". On the other hand, the SM-like particles (h, the gauge bosons and all fermions) have "dark charge" equal to 1. The preservation of this quantum number means that dark particles must always be produced in pairs while in their decays they must always produce at least one dark particle. Therefore, the lightest of these dark particles -which we will choose in our parameter scans to be the lightest neutral state, h 1 -is stable. Thus, the model provides one dark matter candidate.
The model indeed shares many features with the Inert version of the 2HDM, wherein all particles from the "dark doublet" are charged under a discrete symmetry, the lightest of which is stable. The main difference with the current model is the mixing that occurs between the two neutral components of the doublet and the singlet due to the cubic coupling A, which can be appreciated from the mass matrix of eq. (2.9). In what concerns the charged scalar, though, most of the phenomenology of this model is equal to the Inert 2HDM.

Parameter scan, the diphoton signal and dark matter observables
With the model specified, we can set about exploring its available parameter space, taking into account all of the existing theoretical and experimental constraints. We performed a vast scan over the parameter space of the model (100.000 different combinations of the parameters of the potential of eq. (2.2)), requiring that: • The correct electroweak symmetry breaking occurs, and the correct value for the mass of the observed Higgs boson is obtained; as already explained, this is achieved by requiring that v = 246 GeV in eq. (2.6) while at the same time the parameters of the model are such that eqs. (2.5) and (2.7) are satisfied.

JHEP11(2018)091
• By construction, all tree-level interactions and vertices of the Higgs particle h are identical to those of the SM Higgs boson. As a consequence, all LHC production cross sections for h are identical to the values expected in the SM. Additionally, all decay widths of h, apart from the diphoton case to be treated explicitly below, are identical to their SM values up to electroweak corrections. This statement holds as we require the h 1 mass to be larger than roughly 70 GeV, to eliminate the possibility of the decay h → h 1 h 1 (when this decay channel is open it tends to affect the branching ratios of h, making it difficult to have h be SM-like).
• The quartic couplings of the potential cannot be arbitrary. In particular, they must be such that the theoretical requirements of boundedness from below (BFB) -that the potential always tends to +∞ along any direction where the scalar fields are taking arbitrarily large values -and perturbative unitarity -that the model remains both perturbative and unitary, in all 2 → 2 scalar scattering processes -are satisfied. The model considered in the current paper differs from the N2HDM discussed in ref. [36] only via the cubic coefficient A, so the tree-level BFB and perturbative unitarity constraints described there (in sections 3.1 and 3.2) are exactly the ones we should use here.
• The constraints on the scalar sector arising from the Peskin-Takeuchi electroweak precision parameters S, T and U [40][41][42] are required to be satisfied in the model. Not much of the parameter space is eliminated due to this constraint, but it is still considered in full.
• Since the charged scalar H ± does not couple to fermions, all B-physics bounds usually constraining its interactions are automatically satisfied. The direct LEP bound of m H ± > 90 GeV assumed a 100 % branching ratio of H ± to fermions, so that this constraint also needs not be considered here.
• Since all scalars apart from h do not couple to fermions, no electric dipole moment constraints need be considered, this despite the fact that CP violation occurs in the model.
With these restrictions, the scan over the parameters of the model was such that: • The masses of the neutral dark scalars h 1 and h 2 and the charged one, H ± , were chosen to vary between 70 and 1000 GeV. The last neutral mass, that of h 3 , is obtained from the remaining parameters of the model as explained in [36].
• The mixing angles of the neutral mass matrix, eq. (2.12), were chosen at random in the interval −π/2 and π/2.
All other parameters of the model can be obtained from these using the expressions for the masses of the scalars and the definition of the matrix R. The scan ranges for the quartic couplings are chosen larger than the maximally allowed ranges after imposing unitarity and BFB. Therefore, all of the possible values for these parameters are sampled. We have used the implementation of the model, and all of its theoretical constraints, in ScannerS [47].
N2HDECAY [48], a code based on HDECAY [49,50], was used for the calculation of scalar branching ratios and total widths, as in [36].
As we already explained, the tree-level interactions of h are identical to the ones of a SM Higgs boson of identical mass. The presence of the charged scalar H ± , however, changes the diphoton decay width of h, since a new loop, along with those of the W gauge boson and charged fermions, contributes to that width. This is identical to what occurs in the Inert model, and we may use the formulae of, for instance, ref. [22]. Thus we find that the diphoton decay amplitude in our model is given by (3.1) where the sum runs over all fermions (of electric charge Q f and number of colour degrees of freedom N c,f ) and A 0 , A 1/2 and A 1 are the well-known form factors for spin 0, 1/2 and 1 particles (see for instance refs. [51,52]). The charged Higgs contribution to the diphoton amplitude in eq. (3.1) changes this decay width, and therefore the total decay width, hence all branching ratios, of h with respect to the SM expectation. However, the diphoton decay width being so small compared to the main decay channels for h (bb, ZZ and W W ), the overall changes of the total h width are minimal. In fact, numerical checks for our allowed parameter points have shown that the branching ratios of h to bb, ττ , ZZ and W W change by less than 0.05% compared to the corresponding SM quantities -therefore, all current LHC constraints for the observed signal rates of h in those channels are satisfied at the 1σ level.
As for the branching ratio into two photons, it can and does change by larger amounts, as can be appreciated from figure 1. In that figure we plot the ratio of the branching ratio of h into two photons to its SM value as a function of the charged Higgs mass. Comparing these results to the recent measurements of the h → γγ signal rates 2 µ γγ from ref. [53], we see that our model can accommodate values well within the 2σ interval. The lower bound visible in figure 1 emerges from the present experimental lower limit from [53] at 2σ. The experimental upper limit, however, is larger than the maximum value of ∼ 1.2 possible in our model. The latter results from the combination of BFB and unitarity bounds which constrain the allowed values of the coupling λ 3 . The lowest allowed value for λ 3 , which JHEP11(2018)091 Thus we see that the model under study in this paper is perfectly capable of reproducing the current LHC data on the Higgs boson. Specific predictions for the diphoton signal rate are also possible in this model -values of µ γγ larger or smaller than unity are easily accommodated, though they are constrained to the interval 0.917 µ γγ 1.206. As the parameter scan was made taking into account all data from dark matter searches, we are comfortable that all phenomenology in that sector is satisfied by the dark particles.
Let us now study how the model behaves in terms of dark matter variables. Several experimental results put constraints on the mass of the dark matter (DM) candidate, and on its couplings to SM particles. The most stringent bound comes from the measurement of the cosmological DM relic abundance from the latest results from the Planck Collaboration [45], (Ωh 2 ) obs DM = 0.120 ± 0.001. The DM relic abundance for our model was calculated with MicrOMEGAs [44]. In our scan we accepted all points that do not exceed the value measured by Planck by more than 1σ. This way, we consider not only the points that are in agreement with the DM relic abundance experimental values but also the points that are underabundant and would need further dark matter candidates to saturate the measured experimental value.
Another important constraint comes from direct detection experiments , in which the elastic scattering of DM off nuclear targets induces nucleon recoils that are being measured by several experimental groups. Using the expression for the spin-independent DM-nucleon cross section given by MicrOMEGAs, we impose the most restrictive upper bound on this cross section, which is the one from XENON1T [46,54]. In the left panel of figure 2 we use the parameter scan previously described to compute dark matter observables. We show the points that passed all experimental and theoretical constraints in the relic abundance versus dark matter mass plane. We present in pink the points that saturate the relic abundance, that is the points that are in the interval between +1σ and −5σ around the central value, and in violet the points for which the relic abundance is below the measured value. It is clear that there are points in the chosen dark matter mass range that saturate the relic density. In the right panel we present the spin-independent nucleon dark matter scattering cross section as a function of the dark matter mass. The upper bound (the grey line) represents the latest XENON1T [46,54] results. The pink points in the right plot show that even if the direct bound improves by a few orders of magnitude there will still be points for the entire mass range where the relic density is saturated.

JHEP11(2018)091
Thus we see that the model under study in this paper can fit, without need for fine tuning, the existing dark matter constraints. Next we will study the rise of CP violation in the dark sector.

CP violation in the dark sector
As we explained in section 2, the model explicitly breaks the CP symmetry defined in eq. (2.4). Notice that the vacuum of the model which we are studying -wherein only Φ 1 acquires a vev -preserves that symmetry. Therefore, if there is CP violation (CPV) in the JHEP11(2018)091 interactions of the physical particles of the model, it did not arise from any spontaneous CPV, but rather the explicit CP breaking mentioned above. 3 There are several eventual experimental observables where one could conceivably observe CPV. For instance, a trivial calculation shows that all vertices of the form Zh i h j , with i = j, are possible. These vertices arise from the kinetic terms for Φ 2 where from eq. (2.6) we obtain, in terms of the neutral components of the second doublet, where g is the SU(2) L coupling constant and θ W is the Weinberg angle. With the rotation matrix between field components and neutral eigenstates defined in eq. (2.12), we easily obtain (i, j =1,2,3) Notice that no vertices of the form Zhh i are possible. This is not due to any CP properties, however, but rather to the conservation of the Z 2 quantum number. Thus observation of such decays or production mechanisms (all three possibilities for Z → h j h i , i = j, would have to be confirmed) could serve as confirmation of CPV in the model, though the non-observability of the dark scalars would mean they would only contribute to missing energy signatures. Both at the LHC and at future colliders, hints on the existence of dark matter can appear in mono-Z or mono-Higgs searches. The current model predicts cascade processes such as qq (e + e − ) → Z * → h 1 h 2 → h 1 h 1 Z and qq (e + e − ) → Z * → h 1 h 2 → h 1 h 1 h 125 , leading to mono-Z and mono-Higgs events, respectively. This type of final states occurs in many dark matter models, regardless of the CP-nature of the particles involved. Therefore, these are not good processes to probe CP-violation in the dark sector. However, though CPV occurs in the dark sector of the theory, it can have an observable impact on the phenomenology of the SM particles. A sign of CPV in the model -possibly the only type of signs of CPV which might be observable -can be gleaned from the interesting work of ref. [34] (see also ref. [  boson vertices ZZZ and ZW + W − were considered. A Lorentz structure analysis of the ZZZ vertex, for instance [55][56][57][58], reveals that there are 14 distinct structures, which can be reduced to just two form factors on the assumption of two on-shell Z bosons and massless fermions, the off-shell Z being produced by e + e − collisions. Under these simplifying assumptions, the ZZZ vertex function becomes (e being the unit electric charge) where p 1 is the 4-momentum of the off-shell Z boson, p 2 and p 3 those of the remaining (onshell) Z bosons. The dimensionless f Z 4 form factor is CP violating, but the f Z 5 coefficient preserves CP. In our model there is only one-loop diagram contributing to this form factor, shown in figure 3. As can be inferred from the diagram there are three different neutral scalars circulating in the loop -in fact, the authors of ref. [34] showed that in the 2HDM with explicit CPV (the C2HDM) the existence of at least three neutral scalars with different CP quantum numbers that mix among themselves is a necessary condition for non-zero values for f Z 4 . Notice that in the C2HDM there are three diagrams contributing to f Z 4other than the diagram shown in figure 3, the C2HDM calculation involves an additional diagram with an internal Z boson line in the loop, and another, with a neutral Goldstone boson G 0 line in the loop. In our model, however, the discrete Z 2 symmetry we imposed forbids the vertices ZZh j and ZG 0 h i (these vertices do occur in the C2HDM, being allowed by that model's symmetries), and therefore those two additional diagrams are identically zero. In [34] an expression for f Z 4 in the C2HDM was found, which can easily be adapted to our model, by only keeping the contributions corresponding to the diagram of figure 3. This results in where α is the electromagnetic coupling constant and the LoopTools [59] function C 001 is used. The f 123 factor denotes the product of the couplings from three different vertices, given in ref. [34] by

JHEP11(2018)091
where the e i,j,k (i, j, k = 1, 2, 3) factors, shown in figure 3, are related to the coupling coefficients that appear in the vertices Zh i h j (in the C2HDM they also concern the ZG 0 h i and ZZh i vertices, cf. [35]). With the conventions of the current paper, we can extract these couplings from eq. (4.2) and it is easy to show that where the simplification that led to the last line originates from the orthogonality of the R matrix. We observe that the maximum value that f 123 can assume is (1/ √ 3) 3 , corresponding to the maximum mixing of the three neutral components, ρ, η and Φ S ≡ s. This is quite different from what one expects to happen in the C2HDM, for instance -there one of the mixed neutral states is the observed 125 GeV scalar, and its properties are necessarily very SM-like, which implies that the 3 × 3 matrix R should approximately have the form of one diagonal element with value close to 1, the corresponding row and column with elements very small and a 2 × 2 matrix mixing the other eigenstates. 4 Within our model, however, the three neutral dark fields can mix as much or as little as possible.
In figure 4 we show, for a random combination of dark scalar masses (m h 1 80.5 GeV, m h 2 162.9 GeV and m h 3 256.9 GeV) the evolution of f Z 4 normalized to f 123 , 5 with p 2 1 , the 4-momentum of the off-shell Z boson. This can be compared with figure 2 of ref. [34], where we see similar (if a bit larger) magnitudes for the real and imaginary parts of f Z 4 , despite the differences in masses for the three neutral scalars in both situations (in that figure, the masses taken for h 1 and h 3 were, respectively, 125 and 400 GeV, and several values for the h 2 mass were considered). As can be inferred from figure 4, f Z 4 is at most of the order of ∼ 10 −5 . For the parameter scan described in the previous section, we obtain, for the imaginary part of f Z 4 , the values shown in figure 5. We considered two values of p 2 1 (corresponding to two possible collision energies for a future linear collider). The imaginary part of f Z 4 (which, as we will see, contributes directly to CP-violating observables such as asymmetries) is presented as a function of the overall coupling f 123 defined in eq. (4.6). We in fact present results as a function of f 123 /(1/ √ 3) 3 , to illustrate that indeed the model perfectly allows maximum mixing between the neutral, dark scalars. Figure 5 shows that the maximum values for |Im(f Z 4 )| are reached for the maximum mixing scenarios. We also highlight in red the points for which the dark neutral scalars h i have masses smaller than 200 GeV. The loop functions in the definition of f Z 4 , eq. (4.4), have a complicated dependence on masses (and external momentum p 1 ) so that an analytical demonstration is not possible, but the plots of figure 5 strongly imply that choosing all dark scalar masses small yields smaller values for |Im(f Z 4 )|. Larger masses, and larger mass splittings, seem to be required for larger |Im(f Z 4 )|. A reduction on the maximum values of |Im(f Z 4 )| (and |Re(f Z 4 )|) with increasing external momentum is observed (though that variation is not linear, as can be appreciated from figure 4). A reduction of the maximum values of |Im(f Z 4 )| (and |Re(f Z 4 )|) when the external momentum tends to infinity is also observed.

JHEP11(2018)091
The smaller values for |Im(f Z 4 )| for the red points can be understood in analogy with the 2HDM. The authors of ref. [34] argue that the occurrence of CPV in the model implies a non-zero value for the basis-invariant quantities introduced in refs. [60,61], in particular for the imaginary part of the J 2 quantity introduced therein. Since Im(J 2 ) is proportional to the product of the differences in mass squared of all neutral scalars, having all those scalars with lower masses and lower mass splittings reduces Im(J 2 ) and therefore the amount of CPV in the model. Now, in our model the CPV basis invariants will certainly be different from those of the 2HDM, but we can adapt the argument to understand the behaviour of the red points in figure 5: those red points correspond to three dark neutral scalars with masses lower than 200 GeV, and therefore their mass splittings will be small (compared to the remaining parameter space of the model). In the limiting case of three degenerate dark scalars, the mass matrix of eq. (2.9) would be proportional to the identity matrix and therefore no mixing between different CP states would occur. With this analogy, we can understand how regions of parameter space with larger mass splittings between the dark neutral scalars tend to produce larger values of |Im(f Z 4 )|. Experimental collaborations have been probing double-Z production to look for anomalous couplings such as those responsible for a ZZZ vertex [62][63][64][65][66][67][68][69][70]. The search for anomalous couplings in those works uses the effective Lagrangian for triple neutral vertices proposed in ref. [55], parametrised as where γZZ vertices were also considered. In this equation, F µν is the electromagnetic tensor, Z µν = ∂ µ Z ν − ∂ ν Z µ andZ µν = µνρσ Z ρσ /2. The f Z 4 coupling above is taken to be a constant, and as such it represents at most an approximation to the f Z 4 (p 2 1 ) of eq. (4.4). Further, the analyses of the experimental collaborations mentioned above take this coupling to be real, whereas the imaginary part of f Z 4 (p 2 1 ) is the quantity of interest in many interesting observables. With all that under consideration, latest results from LHC [70] already probe the f Z 4 coupling of eq. (4.7) to order ∼ 10 −3 , whereas the typical magnitude of f Z 4 (p 2 1 ) (both real and imaginary parts) is ∼ 10 −5 . We stress , however, that the two quantities cannot be directly compared, as they represent very different approaches to the ZZZ vertex. A thorough study of the experimental results of [70] using the full expression for the ZZZ vertex of eq. (4.3) and the full momentum (and scalar masses) dependence of the form factors is clearly necessary, but beyond the scope of the current work.
The crucial aspect to address here, and the point we wish to make with the present section, is that f Z 4 (p 2 1 ) is non-zero in the model under study in this paper. Despite the fact that the neutral scalars contributing to the form factor are all dark particles, CP violation is therefore present in the model and it can indeed be "visible" to us, having consequences in the non-dark sector. We also analysed other vertices, such as ZW + W −there CPV form factors also arise, also identified as "f Z 4 ", and for our parameter scan we computed it by once again adapting the results of ref. [34] to our model. In the C2HDM three Feynman diagrams contribute to this CP-violating form factor (see figure 17 in [34]) but in our model the Z 2 symmetry eliminates the vertices h i W + W − and h i G + W − , so only JHEP11(2018)091 Interestingly, this form factor is larger, by roughly a factor of ten, than the corresponding quantity in the ZZZ vertex (though still smaller than the corresponding C2HDM typical values). This is illustrated in figure 6, where we plot the imaginary part of f Z 4 as given by eq. (4.8) for p 2 1 = (450 GeV) 2 , having obtained non-zero values. Therefore CPV also occurs in the ZW + W − interactions in this model, though presumably it would be no easier to experimentally establish than for the ZZZ vertex. The point we wished to make does not change, however -if even a single non-zero CPV quantity is found, then CP violation occurs in the model.
As an example of a possible experimental observable to which the form factors f Z 4 for the ZZZ interactions might contribute, let us take one of the asymmetries considered in ref. [34], using the techniques of ref. [71]. Considering a future linear collider and the process e + e − → ZZ, taking cross sections for unpolarized beams σ λ,λ for the production of two Z bosons of helicities λ andλ (assuming the helicity of the Z bosons can be determined), the asymmetry A ZZ 1 is defined as with θ the angle between the electron beam and the closest Z boson with positive helicity, β = 1 − 4m 2 Z /p 2 1 denoting the velocity of the produced Z bosons and the function F 1 (β, θ) is given in appendix D of ref. [34]. Choosing the two points in our parameter scan with largest (positive) and smallest (negative) values of Im f Z 4 (p 2 1 ) for p 2 1 = (450 GeV) 2 , we obtain the two curves shown in figure 7. Clearly, the smallness (∼ 10 −5 ) of the f Z 4 form factor renders the value of this asymmetry quite small, which makes its measurement challenging. This raises the possibility that asymmetries involving the ZW + W − vertex might be easier to measure than those pertaining to the ZZZ anomalous interactions, since we have shown that f Z 4 is typically larger by a factor of ten in the former vertex compared to the latter one. To investigate this possibility, we compared A ZZ 1 , considered above, with the A W W 1 asymmetry defined in eq. (5.21) of ref. [34]. A direct comparison of the maximum values of A W W 1 and A ZZ 1 shows that for some regions of parameter space the former quantity can indeed be one order of magnitude larger than the latter one; but that is by no means a generic feature, since for other choices of model parameters both asymmetries can also be of the same order. Notice that both asymmetries show a quite different √ s dependence.

Conclusions
We presented a model whose scalar sector includes two Higgs doublets and a real singlet. A specific region of parameter space of the model yields a vacuum which preserves a discrete symmetry imposed on the model -thus a charged scalar and three neutral ones have a "dark" quantum number preserved in all interactions and have no interactions with JHEP11(2018)091 fermions. The lightest of them, chosen to be a neutral particle, is therefore stable and a good dark matter candidate. The first doublet yields the necessary Goldstone bosons and a neutral scalar which has automatically a behaviour almost indistinguishable from the SM Higgs boson. A parameter scan of the model, imposing all necessary theoretical and experimental constraints (including bounds due to relic density and dark matter searches, both direct and indirect) shows that the SM-like scalar state indeed complies with all known LHC data for the Higgs boson -some deviations may occur in the diphoton signal rate due to the extra contribution of a charged scalar to the involved decay width, but we have shown such deviations are at most roughly 20% of the expected SM result when all other constraints are satisfied, and this is still well within current experimental uncertainties.
The interesting thing about the model presented in this paper is the occurrence of explicit CP violation exclusively within the dark matter sector. A complex phase allowed in the potential forces the neutral components of the second (dark) doublet to mix with the real singlet to yield three neutral eigenstates, none of which possesses definite quantum numbers. Signals of this CP violation would not be observed in the fermion sector (which, by the way, we assume is identical to the one of the SM, and therefore has the usual CKMtype source of CP violation) nor in the interactions of the SM-like scalar -protected as it is by the unbroken Z 2 symmetry, and by the mass ranges chosen for the dark scalars, h will behave like a purely CP-even SM-like scalar, even though the CP symmetry of the model is explicitly broken in the scalar sector as well! Can the model then be said to be CP violating at all? The answer is yes, as an analysis of the contributions from the dark sector to the ZZZ vertex demonstrates. Even though the dark particles have no direct fermion interactions and could elude detection, their presence could be felt through the emergence of anomalous triple gauge boson vertices. Though we concentrated mainly on ZZZ vertices we also studied ZW + W − interactions, but our main purpose was to show CPV is indeed occurring. Direct measurements of experimental observables probing this CPV are challenging: we have considered a specific asymmetry, A ZZ 1 , built with ZZ production cross sections, but the magnitude of the CPV form factor f Z 4 yields extremely small values for that asymmetry, or indeed for other such variables we might construct. Direct measurements of ZZ production cross sections could in theory be used to constraint anomalous ZZZ vertex form factors -and indeed several experimental collaborations, from LEP, Tevatron and LHC, have tried that. But the experimentalists' approach is based on constant and real form factors, whereas model-specific expressions for f Z 4 such as those considered in our work yield quantities highly dependent on external momenta, which boast sizeable imaginary parts as well. Thus a direct comparison with current experimental analyses is not conclusive.
The other remarkable fact is the amount of "damage" the mere inclusion of a real singlet can do to the model with two doublets. As repeatedly emphasised in the text, the model we considered is very similar to the Inert 2HDM -it is indeed simply the IDM with an added real singlet and a tweaked discrete symmetry, extended to the singlet having a "dark charge" as well. But whereas CP violation -explicit or spontaneous -is entirely impossible within the scalar sector of the IDM, the presence of the extra singlet produces a completely different situation. That one obtains a model with explicit CPV is all the more JHEP11(2018)091 remarkable when one considers that the field we are adding to the IDM is a real singlet, not even a complex one. Notice that within the IDM it is even impossible to tell which of the dark neutral scalars is CP-even and which is CP-odd -all that can be said is that those two eigenstates have opposite CP quantum numbers. The addition of a real singlet completely changes the CP picture.
The occurrence of CP violation in the dark matter sector can be simply a matter of curiosity, but one should not underestimate the possibility that something novel might arise from it. If the current picture of matter to dark matter abundance is indeed true and the observed matter is only 5% of the total content of the universe, then one can speculate how CP violation occurring in the interactions of the remainder matter might have affected the cosmological evolution of the universe. We reserve such studies for a follow-up work.