Neutrino Mass, Leptogenesis, and Dark Matter from The Dark Sector with $U(1)_{D}$

I suggest a new extension of the SM by introducing a dark sector which has several new particles and a local $U(1)_{D}$ symmetry. The dark particles bring about the new and interesting physics beyond the SM. The model can generate the tiny neutrino mass by a hybrid see-saw mechanism, achieve the leptogenesis at the TeV scale, and account for the cold dark matter. All of the three things collectively arise from the dark sector. In particular, it is very feasible to test the model predictions and probe the dark sector in near future experiments.


I. Introduction
The standard model (SM) of the particle physics has successfully accounted for all kinds of the physics at or below the electroweak scale, refer to the reviews in Particle Data Group [1], but it can not explain the three important issues: the tiny neutrino mass [2], the matter-antimatter asymmetry [3], and the cold dark matter (CDM) [4]. Many theories have been suggested to solve these problems. The tiny neutrino mass can be generated by the seesaw mechanism [5] or origin from the loop-diagram radiative generation [6]. The baryon asymmetry can be achieved by the thermal leptogenesis [7] or the electroweak baryogenesis [8]. The CDM candidates are possibly the sterile neutrino [9], the lightest supersymmetric particle [10], the axion [11], and so on. These theories ordinarily focus on only one of the three problems. In recent years, some inspired ideas attempt to find some connections among the neutrino mass, the baryon asymmetry, and the CDM, for example, the lepton number violation at the super-high scale can lead to the neutrino mass and the leptogenesis [12], the neutrino mass and the leptogenesis can also be implemented by the super-heavy scalar triplet [13], the asymmetric CDM is related to the baryon asymmetry [14], and some sophisticated models unifying them into a frame [15]. Although many progresses on these fields have been made all the time, an universal and convincing theory is not established as yet.
What is exactly a realistic theory beyond the SM? The universe harmony and the nature unification are a common belief of mankind. It is very possible that there is a common origin of the tiny neutrino mass, the matter-antimatter asymmetry and the CDM, which relates the three things to each other. Therefore, a new theory beyond the SM should be capable of accounting for the three things collectively, on the other hand, this theory should keep such principles as the simplicity, the fewer number of parameters, and being feasible and promising to be tested in future experiments. If one theory is excessive complexity and unable to be tested, it is unbelievable and infeasible. Based on these considerations, I suggest a new extension of the SM. The model introduces a dark sector beyond the SM, which contains a few of new particles and has a local gauge symmetry of U(1) D , in particular, the dark sector provides a common origin of the above three things so that the model can completely account for them. Finally, it is very feasible to test the model and probe the dark sector by the TeV-scale colliders, the underground detectors, and the cosmic ray search.
The remainder of this paper is organized as follows. I outline the model in Section II. Section III and Section IV are respectively discussions of the leptogenesis and the dark matter. Section V is numerical results and discussions of the model test. Section VI is devoted to conclusions.

II. Model
The model introduces several new particles and a local gauge symmetry U(1) D beyond the SM sector, which are in the dark sector. In addition, it obeys the global symmetry U(1) B−L , i.e., the difference between the baryon number and the lepton one is conserved.
are all singlets under the SM group. The particles in the dark sector have all non-vanishing D numbers, while all of the SM particles have no D numbers. Note that N R is filled into the SM sector but N L belongs to the dark sector, N R and N L will combine to form a heavy neutral Dirac fermion after the model symmetry breakings. χ R is a neutral Majorana fermion, it will become the CDM. The dark neutral scalars φ 1 and φ 2 are applied to implement the spontaneous breakings of U(1) D and U(1) B−L , while the dark doublet scalar Φ plays a key role in the generations of the neutrino mass and the matter-antimatter asymmetry. Finally, the model has also a hidden Z 2 symmetry, it is defined by the following discrete transform, and all the gauge fields remain unchanged, where f L,R denote the left-handed and righthanded fermions in Tab. 1. The last line of Tab. 1 lists the Z 2 parity of each field. Note that N L and χ R have the same gauge quantum numbers but they have opposite Z 2 parities. By virtue of the assignment of Tab. 1, it is easily verified that all of the chiral anomalies are completely cancelled in the model, namely the model is anomaly-free. Under the above symmetries, the invariant Lagrangian of the model is composed of the three parts of the gauge kinetic energy terms, the Yukawa couplings and the scalar potentials. The gauge kinetic energy terms are where f and S respectively denote all kinds of the fermions and scalars in Tab. 1. g 0 and Z ′ µ are the gauge coupling coefficient and gauge field which are associated with the local U(1) D symmetry. τ i is the Paul matrices and the other notations are self-explanatory.
The Yukawa couplings are where C is the charge conjugation matrix and CN L T = N c R , χ T R C = χ c L . Note that the Z 2 symmetry of Eq. (1) forbids the explicit mass term N L Mχ R even though it satisfies all the gauge symmetries. This will guarantee the stability of χ R since it can not mix with the other fermions. These coupling parameters in Eq. (3), Y e , Y 1 , etc., are all 3 × 3 complex matrices in the flavour space, the leading elements of them take such moderate values as In addition, we can however choose such flavour basis in which Y e , Y N , Y χ are simultaneously diagonal matrices, namely the mass eigenstate basis (see the following Eq. (14)), thus Y 1 and Y 2 certainly contain some irremovable complex phases, they eventually become CP -violating sources in the lepton sector in comparison with one in the quark sector. Eq. (3) will give rise to all kinds of the fermion masses after the scalar fields developing their non-vanishing vacuum expectation values. In particular, the spontaneous breakings of U(1) D and U(1) B−L will bring about significant results, the Y 1 and Y 2 terms can jointly lead to the tiny neutrino mass and the successful leptogenesis, while the Y χ term will generate the CDM and yield correct annihilation cross-section.
Here we assume that there is no CP violation in the scalar sector, in addition, in Eq. (4) the self-interaction of each scalar is stronger but the interactions among them are feeble, so in Eq. (5) those self-coupling parameters are far larger than those interactive coupling parameters. µ 2 φ 1 < 0 and |µ 2 φ 1 | ≫ [µ 2 φ 2 , µ 2 H , µ 2 Φ ] will lead that φ 1 first develops a nonvanishing vacuum expectation value at high scale, and then the other scalars are induced to develop non-vanishing vacuum expectation values at low scale via their coupling to φ 1 , eventually the symmetry breakings proceed along the chain of the following Eq. (9). µ 2 Φ > 0 means that µ Φ is actually the original mass of the dark doublet scalar Φ. In short, the conditions of Eq. (5) are natural and reasonable, they can sufficiently guarantee the vacuum stability and the spontaneous breakings of the model symmetries in the proper order.
The vacua of the spontaneous breakings are along the directions of the neutral component of each scalar field. We can rigorously solve the V S minimum for Eq. (4), and then derive the vacuum configurations as follows, where the four vacuum expectation values are determined by the tadpole equations as follows, A physics solution of the vacuum expectation values which fulfil the conditions of Eq (5) is a hierarchy such as In fact, v H = 246 GeV has been fixed by the electroweak physics, v 2 will be determined by the CDM physics, v 1 and v Φ will jointly be determined by the tiny neutrino mass and the successful leptogenesis.
According to the assignment of Tab. 1 and the relations of Eq. (8), the model symmetries are spontaneously broken step by step through the breaking chain as follows, At the first step, φ 1 ∼ 10 6 GeV breaks the local U(1) (B−L)+D and the discrete Z 2 but the global U(1) (B−L)−D is kept as a residual symmetry. The φ 0 1 component becomes a massive real scalar boson around the v 1 scale, while the pseudo-scalar Goldstone boson G ′0 is eaten by the massless Z ′ µ , the latter becomes a massive gauge boson through the Higgs mechanism. In addition, the Y N term in Eq. (3) generates a heavy Dirac fermion mass around the v 1 scale by a combination of the neutral N L and N R . At the second step, φ 2 ∼ 10 2 GeV violates the global U(1) (B−L)−D , φ 0 2 and G 0 in Eq. (6) respectively become a massive real scalar and a massless Goldstone boson, in addition, χ R obtains a Majorana mass around the v 2 scale, it will become the CDM because of its characteristics (which will be discussed at IV Section). Note that the Y ′ N term in Eq. (3) also generates a Majorana mass of N L around the v 2 scale, but it is far smaller than the Dirac mass of N, so this can not change the nature of N as a Dirac fermion. At the third step, the electroweak breaking is accomplished by H ∼ 10 2 GeV, the SM particles obtain their masses around the electroweak scale. It should be stressed that the B − L − D violation occurs at the time close to the electroweak breaking due to v 2 ∼ v H . Lastly, the above three breakings can also induce the neutral component of the dark doublet Φ developing a smaller Φ ∼ 0.5 MeV via the λ 7 feeble coupling term in Eq (4), of course, v Φ is too small to make an effect on the nature of the heavy Φ whose mass is ∼ µ Φ .
After the above symmetry breakings are completed, these components of G ′0 , G ′′0 , H − in Eq. (6) have been transformed into the gauge sector to generate M Z ′ , M Z , M W through the Higgs mechanism, therefore the scalar sector now includes four CP -even neutral , and the squared mass of the charged Φ ∓ , are together given as follows, where is obviously an approximately diagonal matrix because those non-diagonal elements are far smaller than those diagonal elements, of course, this arises from only feeble coupling among the scalar fields, so we can safely neglect the mixing among the CP -even bosons, for example, the mixing angle between φ 0 2 and h 0 is ∼ on account of Eq. (5) and Eq. (8), therefore the neutral and charged components of the dark doublet Φ actually keep a degenerating mass in despite of their tiny splits.
In the gauge sector, the squared mass matrix of the three neutral gauge bosons (B µ , W 3 µ , Z ′ µ ) and the squared mass of the charged gauge bosons W ± µ are given by where the photon field A µ is massless because the determinant of M 2 N GB is vanishing. The mixing angle between Z ′ µ and B µ is ∼ , obviously, they are nearly zero, so we can leave them out. The mixing angle between B µ and W 3 µ is cos θ w = , which is namely the weak mixing angle of the SM. In addition, to the SM values, which is very difficult to be detected. In the fermion sector, the mass matrix of the neutral fermions is written as where M ν , M N , M χ are given by Eq. (14) below, and ν ′ , N ′ are new fermion fields after the flavour rotation. Because of v 2 Y ′ N ≪ v 1 Y N , N L and N R actually combine into a heavy Dirac fermion. In addition, χ R has no mixing with the other neutral fermions due to the Z 2 symmetry in Eq. (1), this leads that it eventually becomes the CDM. At the low energy, the heavy Dirac fermion N has decoupled and it can be integrated out from the two terms of Y 1 and Y 2 in Eq. (3), thus we can obtain an effective coupling of the left-handed doublet lepton According to the assignment of Tab. 1, this effective coupling explicitly violates both one unit of B − L number and one unit of D number, namely violates two unit of B − L + D number but conserves B − L − D number. In fact, Eq. (13) is exactly a common origin of both the tiny neutrino mass and the leptogenesis.
We can now give all of the model particle masses as follows, M h 0 is namely the Higgs boson mass of the SM, which has been measured as M h 0 ≈ 125 GeV [1]. M Φ is about several TeVs, which is close to its original mass µ Φ . The tiny Majorana mass of the SM neutrino is generated by Eq. (12) This is easily fulfilled by choosing some suitable values of the coupling parameters in Eq. (14). However, the mass relations of Eq. (15) will successfully lead to the leptogenesis and the CDM. Finally, it is worth emphasizing that there are not any super-high scale physics in the model.

III. Leptogenesis
The model can account for the baryon asymmetry through the leptogenesis. Below the v 1 scale but above the v 2 scale, B − L + D is violated by φ 1 , but B − L − D is conserved, moreover, is anomaly-free. The dark doublet Φ at the TeV scale has two decay modes on the basis of the model couplings and the (15) relations, (i) the two-body decay Φ → H + φ 2 via the λ 7 term in Eq. (4), (ii) the three-body decay Φ → l α + l β + H via Eq. (13). Fig. 1 shows the tree and loop diagrams of Φ → l α + l β + H. Note that the Figure 1: The tree and loop diagrams of the decay Φ → l α + l β + H violating B − L + D but conserving B − L − D, this decay has the CP asymmetry and is out-of-equilibrium, which eventually leads to the matter-antimatter asymmetry.
three-body decay is mainly mediated via N 1 as shown Fig. 1, the diagrams in which the decay is mediated via N 2,3 can be neglected because they are greatly suppressed due to In term of the quantum number assignment in Tab. 1, explicitly, this process of Fig. 1 Fig. 1 has a CP asymmetry and is out-of-equilibrium, but the decay of Φ → H + φ 2 has no these features. However, the two-body decay rate is much larger than the three-body one, therefore the total decay width of Φ is approximately equal to the decay width of Φ → H + φ 2 .
Because the two coupling matrices of Y 1 and Y 2 contain the CP -violating sources, the decay rate of Φ → l α + l β + H is different from one of its CP -conjugate process Φ → l α + l β + H through the interference between the tree diagram and the loop one. The CP asymmetry of the two decay rates is defined and calculated as follows, A careful calculation shows that the imaginary part of the loop integration factor of the (b) diagram is derived from the three-point function Im[(C 0 +C 12 )(M 2 lα , s 12 12 , where s 12 = (p H + p l β ) 2 , but the (c) diagram has in fact no contribution to ε because the imaginary part of its three-point function is vanishing. According to the discussions in Section II, there are Y 1 ∼ Y 2 ∼ 0.01 and λ 7 ∼ 10 −7 , then we can roughly estimate ε ∼ 10 −8 by use of Eq. (8) and Eq. (15), it is exactly a reasonable and suitable value for the leptogenesis.
A further calculation shows that the decay rate Γ ± in Eq. (16) is smaller than the Hubble expansion rate of the universe, namely where M P l = 1.22 × 10 19 GeV and g * is the effective number of relativistic degrees of freedom. Therefore the decay process of Fig. 1 is actually out-of-equilibrium. At the scale of T = M Φ , the relativistic states include all of the SM particles, and φ 2 and χ i in the dark sector, so one can figure out g * = 114 in Eq. (17). By this time, we have completely demonstrated that the decay process of Fig. 1 is able to satisfy Sakharov's three conditions [16]. The above discussions are now integrated together, as a result, the decay of Fig. 1 certainly generates the following asymmetries [17], where s is the entropy density and κ is a dilution factor. The dilution is mostly from the inverse decay. For the three-body inverse decay, the dilution effect is very weak if the departure from thermal equilibrium is severe, so we can take κ ≈ 1 in Eq. (18). Note that the dilution effect from N 1 → Φ + l α is almost nothing because N 1 has early decoupled and its number density is exponentially suppressed by M Φ ∼ 20 compared to the Φ one. At the high energy, the SM sector and the dark sector are connected each other via the heavy particles Z ′ µ , N i , φ 0 1 , Φ mediating. As the universe temperature falls below M Φ , all of the heavy particles have completely decayed and decoupled, thus the connection between the SM sector and the dark sector is gradually discontinued, eventually, at the low energy the SM sector and the dark sector are isolated from each other, the surviving φ 1 and χ i in the dark sector constitute a dark world. As a consequence, the Y B−L asymmetry is totally deposited in the SM sector because φ 1 and χ i in the dark sector have vanishing B − L numbers (see Tab. 1), while the Y D asymmetry is totally deposited in the dark sector because all of the SM particles have no D numbers (see Tab. 1). In other words, at the low energy, the SM sector has the B − L asymmetry without the D asymmetry, while the dark sector has the D asymmetry without the B − L asymmetry. However, the total Y B−L−D asymmetry in the two sectors is kept to be zero.
When the universe temperature drops to the energy scale v 2 ∼ v H ∼ 250 GeV, the global B − L − D symmetry is broken because φ 2 developing φ 2 violates two unit of D number. The original complex φ 2 is now decomposed into the two neutral real φ 0 2 , G 0 which are namely their own antiparticles, and also χ R becomes a neutral Majorana fermion, namely χ = χ c . Consequently, the Y D asymmetry between dark particles and dark antiparticles is totally erased, or rather it is automatically vanishing. Although the total B − L − D conservation in the two sectors is violated, the SM sector always conserves the B −L number, therefore the Y B−L asymmetry in the SM sector is unaffected and unchanged. At this temperature of v 2 ∼ v H ∼ 250 GeV, obviously, the electroweak sphaleron process can fully put into effect in the SM sector [18], thus a part of the Y B−L asymmetry is converted into the baryon asymmetry. This is given by the following relation, where c s = 28 79 is the sphaleron conversion coefficient. 7.04 is a ratio of the entropy density to the photon number density. 6.2 × 10 −10 is the current value of the baryon asymmetry [19]. Note that χ i , φ 0 2 , G 0 in the dark sector do not at all participate in the sphaleron process since they are all singlets under the SM group and isolated from the SM sector. When the universe temperature falls below T ∼ 100 GeV, the sphaleron process is closed and the baryon asymmetry is kept up to the present day. Finally, it should be stressed that the leptogenesis is elegantly achieved at the TeV scale in the model.

IV. Dark Matter
When the global symmetry of U(1) B−L−D is broken by φ 2 at the scale v 2 ∼ v H , the dark particles χ i , φ 0 2 , G 0 are at first in thermal equilibrium in the dark world. At a later time φ 0 2 can completely decay into a pair of χ i or G 0 . However χ i is a stable particle because it is protected by the two factors. i) The model gauge symmetries prevent it from coupling to the other particles except φ 2 (see Eq. (3)), so it can not decay. ii) The Z 2 symmetry in Eq. (1) forbids the explicit mass term N L Mχ R , so it can not mix with the other fermions. Therefore χ i is a stable WIMP. Among the three generation of χ 1,2,3 , a pair of heavier χ 2,3 mainly annihilate into a pair of the lightest χ 1 via the G 0 mediator, as shown (a) in Fig. 2. At the last step a pair of χ 1 annihilate into a pair of G 0 by the two modes of (b) and (c) in Fig. 2. A careful analysis shows that the annihilation cross-section of χ 2,3 in (a) diagram is much larger than one of χ 1 in (b) and (c) diagrams, therefore, the χ 2,3 decoupling is much earlier than the χ 1 one, and the relic abundance of χ 2,3 is much smaller than one of χ 1 . Consequently, χ 1 becomes the principal particle of the CDM, while χ 2,3 only bears a tiny part of the CDM budget, in addition, G 0 becomes a dark background radiation. In short, χ 1 is a desirable candidate of the CDM because its natures and relic abundance are very well consistent with ones of the CDM. χ 1 becomes a non-relativistic particle when the temperature falls below M χ 1 . It has two annihilation channels. (i) χ 1 + χ 1 → G 0 + G 0 via the χ 1 t-channel mediation, as shown (b) in Fig. 2, note that the χ 1 u-channel mediation is not shown in Fig. 2 but it is included in the following calculation. (ii) χ 1 + χ 1 → G 0 + G 0 via the φ 0 2 s-channel mediation, as shown (c) in Fig. 2. The total annihilation rate of χ 1 is calculated as follows, where a and b are determined by the σv equation, v is a relative velocity of two annihilating particles. Eq. (20) clearly shows that the annihilation cross-section essentially arises from a p-wave contribution. In view of Eq. (8) and Eq. (15), the thermal average on the annihilation cross-section is actually σv ∼ 10 −9 GeV −2 , which is exactly a weak interaction cross-section. This naturally reproduces the so-called "WIMP Miracle" [20].
As the universe temperature decreasing, the annihilation rate of χ 1 becomes smaller than the Hubble expansion rate of the universe, then the annihilation is out-of equilibrium and χ 1 is decoupling. The freeze-out temperature is determined by After χ 1 is frozen out, its numbers in the comoving volume has no change any more. The current relic abundance of χ 1 can be calculated by the following equation [20], 0.12 is the current abundance of the CDM [21]. Obviously, both M χ 1 and v 2 are jointly in charge of the final results of Eq. (21) and Eq. (22). Provided M χ 1 ∼ 60 GeV and v 2 ∼ 250 GeV, then the solution of Eq. (21) is x ∼ 1 25 , namely the freeze-out temperature is T f ∼ 2.3 GeV. At this temperature the relativistic particles include photon, gluon, ν 0 , e − , µ − , u, d, s and G 0 , so we can figure out g * (T f ) = 62.75. Finally, we can correctly reproduce Ω χ 1 h 2 ∼ 0.12 by Eq. (20) and Eq. (22).
Obviously, the G 0 decoupling is exactly at the same temperature as the CDM χ 1 decoupling. Since G 0 is massless and a relativistic decoupling, nowadays it should become a dark background radiation, which is analogous to the CMB photon in the visible world. Because the G 0 decoupling is much earlier than the neutrino decoupling and the photon one, its effective temperature is lower than the neutrino effective temperature and the CMB photon temperature. As a result, the current abundance of G 0 , Ω G 0 , is smaller than the neutrino abundance Ω ν ≈ 1.7 × 10 −3 and the CMB photon abundance Ω γ ≈ 5 × 10 −5 , refer to the review of cosmological parameters in [1]. However, G 0 is in the dark sector and does not interact with the SM matters, so we can not detect it through the ordinary methods.
Two CDM χ 1 can interact each other through the long range exchange of G 0 , but its effective potential is a repulsive force since χ 1 is its own antiparticle, therefore there are not any bound states for the CDM χ 1 , they can only happen elastic scattering. When its reaction rate is smaller than the universe expansion rate, this elastic scattering will be frozen out and closed. The frozen-out temperature is determined by where v ≈  (20) are a weak interaction cross-section. Therefore, the elastic scattering between the CDM χ 1 is actually frozen out at the same time when they are decoupling. Thereafter the CDM χ 1 are completely free particles except the gravitational influence.
In conclusion, the model can simply account for the CDM, in particular, naturally explain the "WIMP Miracle".

V. Numerical Results and Discussions
We now show some concrete numerical results of the model. All of the SM parameters have been fixed by the current experimental data [1]. Some new parameters in the model can be determined by the current data of the tiny neutrino mass, the baryon asymmetry, and the CDM abundance. For the sake of simplicity, we only choose a set of typical values in the parameter space such as where I use M Φ as an independent input parameter instead of  13 ] are taken as the two precise values in order to fit Ω χ 1 h 2 and η B respectively, while the rest of the parameters are roughly fixed to their order of magnitudes. Now put Eq. (24) into the foregoing equations, we can correctly reproduce the desired results, they are in agreement with the current experimental data very well [1]. Here we only give the upper bound of neutrino mass which is assumed as m ν 3 . The full experimental data of the neutrino masses and mixing angles can completely be fitted by choosing suitable texture of the matrix In addition, we can calculate out Γtree H ≈ 0.07 by Eqs. (16)- (17), this clearly demonstrates that the decay of Fig. 1 is indeed severely out-ofequilibrium. Fig. 3 shows the three curves of M φ 0 2 versus M χ 1 for the three cases of v 2 = 300 GeV, v 2 = 250 GeV, and v 2 = 200 GeV, each point of the curves can correctly fit Ω χ 1 h 2 ≈ 0.12. Note that the other parameters are not involved in this fitting. Evidently, a reasonable and moderate value of v 2 is in the range of 200 GeV v 2 300 GeV. If v 2 is too high, the χ 1 annihilation cross-section will be too small, then the χ 1 relic abundance will be overclose. Conversely, if v 2 is too low, the χ 1 annihilation cross-section will be too large, then the χ 1 relic abundance will be deficiency. Therefore v 2 should be around the electroweak scale v H . In addition, the curves clearly indicate that the φ 0 2 mass is probably in the area of 100 GeV M φ 0 2 250 GeV, and the CDM χ 1 mass is probably in the range of 40 GeV M χ 1 90 GeV. In short, the future experimental search for χ 1 and φ 0 2 should focus on this parameter space of Fig. 3.
In the end, we simply discuss the test of the model. Some new particles can be produced at the TeV-scale colliders. The relevant processes are as follows, At the present LHC [22], we have a chance to search Φ and Φ through two gamma photon fusion if the collider energy can reach their masses, but this detection is very difficult because its cross-section is too small. A better way to produce Φ and Φ is at the e − + e + or p + p colliders via the gamma photon s-channel mediation as long as the center-of-mass energy is enough high, for instance, the future colliders such as CEPC and ILC have some potentials to achieve this goal [23]. Only if Φ and Φ are produced, firstly, we can directly test the leptogenesis mechanism of the model by the decay asymmetry of Φ → l α + l β + H and Φ → l α + l β + H. Secondly, this can indirectly shed light on the neutrino mass origin, Figure 4: The indirect detection of the CDM χ 1 through a search for the high-energy gamma photon and Goldstone boson in the cosmic rays, and its mass can accurately be measured by the energy relation E γ = E G 0 = M χ 1 .
namely the model see-saw mechanism. Thirdly, we can probe the dark particles φ 0 2 and G 0 by Φ → H + φ 0 2 and Φ → H + G 0 . Lastly, φ 0 2 can decay into a pair of the CDM χ 1 or G 0 , by which we can measure the χ 1 mass and find the Goldstone boson. All kinds of the final state signals are very clear in the decay chain of Φ and Φ. These search are possibly a intriguing direction in the future collider experiments.
Of course, the model can also be tested by some non-collider experiments. An indirect detection for the CDM χ 1 is a search for the high-energy gamma photon and Goldstone boson in the cosmic rays [24], they arise from an annihilation process of the CDM χ 1 in the dark galactic halo as shown Fig. 4. Since both gamma photon and Goldstone boson are massless, their energy are fixed as E γ = E G 0 = M χ 1 due to conversation of energy and momentum. If the gamma ray whose energy is 40 − 90 GeV is found, this is not only a definite signal of the process of Fig. 4, but also it can tell us the accurate mass of the CDM χ 1 . A direct detection is very difficult by means of scattering off nuclei at the underground detectors such as DAMA, XENON, etc. [25], but it is not impossible. The diagram of the scattering process is similar to Fig. 4 but the photon line becomes an internal one and it is connected with a proton as an external line. In short, it will be very large challenges to actualize the above-mentioned experiments, this needs the researchers make a great deal of efforts, however, its scientific significance is beyond all doubt. We will give an in-depth discussion on the model test in another paper.

VI. Conclusions
In summary, I suggest a new extension of the SM by introducing the dark sector with the local U(1) D symmetry. The particles in the dark sector have all non-vanishing D numbers, while all of the SM particles have no D numbers. The model also conserves the global B − L symmetry and the hidden discrete Z 2 one. The three symmetries are together broken by φ 1 at the scale of thousands of TeVs, but the global B − L − D is kept as a residual symmetry. This breaking gives rise to heavy neutral gauge boson M Z ′ and neutral Dirac fermion M N in the dark sector. When the universe temperature is close to the electroweak scale, the global B − L − D is violated by φ 2 , this generates the CDM χ 1 mass and leads to the "WIMP Miracle". The dark doublet scalar Φ with several TeVs mass can decay into two left-handed doublet leptons and one doublet Higgs of the SM, this process can elegantly achieve the leptogenesis at the TeV scale. The tiny neutrino mass is generated by the hybrid see-saw mechanism, it is suppressed by both the heavy M N and the small Φ induced from the feeble scalar coupling. In brief, the model with fewer parameters is a simple and natural extension of the SM, it can collectively account for the tiny neutrino mass, the matter-antimatter asymmetry, and the CDM. In particular, the model gives some interesting predictions, for example, the leptogenesis at the TeV scale, the CDM χ 1 with dozens GeVs mass, the dark neutral scalar boson φ 0 2 with 100 − 250 GeV mass, the dark background radiation of Goldstone bosons with a tiny abundance, all of them are possibly probed by the TeV collider experiments, the underground detectors, and the cosmic ray search. In short, new physics of the dark sector beyond the SM sector are very attractive and worth researching in depth.