A Simple Testable Model of Baryon Number Violation: Baryogenesis, Dark Matter, Neutron-Antineutron Oscillation and Collider Signals

We study a simple TeV-scale model of baryon number violation which explains the observed proximity of the dark matter and baryon abundances. The model has constraints arising from both low and high-energy processes, and in particular, predicts a sizable rate for the neutron-antineutron ($n-\bar{n}$) oscillation at low energy and the monojet signal at the LHC. We find an interesting complementarity among the constraints arising from the observed baryon asymmetry, ratio of dark matter and baryon abundances, $n-\bar{n}$ oscillation lifetime and the LHC monojet signal. There are regions in the parameter space where the $n-\bar{n}$ oscillation lifetime is found to be more constraining than the LHC constraints, which illustrates the importance of the next-generation $n-\bar{n}$ oscillation experiments.


I. INTRODUCTION
Baryon number (B), an accidental global symmetry of the Standard Model (SM) Lagrangian, must be broken to dynamically generate the observed baryon asymmetry of the universe [1]. At the level of effective theories for B violation, the ∆B = 1 operators are typically of dimension-6 (or dimension-5 in supersymmetric models), so the amplitude scales as Λ −2 (or Λ −1 ), where Λ is the associated new physics scale. The stringent limits on the classic ∆B = 1 process of proton decay [2,3] imply that Λ 10 15 GeV, or the underlying new physics must be at the Grand Unified Theory (GUT) scale or above [4]. On the other hand, the ∆B = 2 operators are of dimension-9, so the amplitude scales as Λ −5 and the associated new physics scale can be much lower, even in the few TeV range [5], while satisfying the experimental constraints on ∆B = 2 processes like di-nucleon decay [2,6] and neutron-antineutron (n−n) oscillation [7][8][9]. This makes it plausible to envisage a TeV-scale model of ∆B = 2 that can be tested in laboratory experiments, while simultaneously solving the cosmological puzzle of baryogenesis, and possibly as a bonus, the apparent coincidence of baryon and dark matter (DM) abundances.
Recently, a minimal TeV-scale extension of the SM has been proposed to this effect that explains the proximity of baryon and DM abundances [10,11]. This model introduces B-violating interactions via a set of O(TeV) color-triplet scalars X α and a singlet Majorana fermion ψ that are coupled only to right-handed (RH) quarks. Baryogenesis occurs via out-of-equilibrium decays of the X α 's. The new fermion ψ becomes stable, hence a viable DM candidate, when it has approximately the same mass as the proton. Obtaining the correct DM relic abundance requires a non-thermal mechanism as the DM annihilation rate turns out to be much smaller than the thermal value of σ ann v = 3 × 10 −26 cm 3 s −1 . The late decay of a modulus field can be the common source of non-thermal DM and baryogenesis [12]. In this model, the B-violating interactions do not generate any new contribution to the electric dipole moment of neutron or electron even at two-loop level.
Since the stability of DM in this model is tied to the that of proton, no additional symmetry is needed to ensure that the DM candidate is stable. As a result, the model naturally predicts monojet/monotop signals at the LHC and the parameter space of this model is currently being probed by the LHC data [13,14]. In addition, B violation at the TeV scale results in potentially observable low-energy effects, most notably double proton decay and n −n oscillation.
In this work, we perform a detailed study of this simple, predictive model to carve out the allowed parameter space that can be probed in the near future. We show that successful production of DM and baryon asymmetry via modulus decay sets interesting bounds on the modulus mass range and the DM coupling to quarks. We also show that the parameter space of the model can be investigated at the LHC via monojet, dijet and a pair produced dijet resonance searches. The monojet signal in the allowed parameter space is mostly produced by a b quark fusion with d and s quarks where the b quark partly arises from the gluon splitting which allows a b jet to be associated with the monojet signal. There are regions in the parameter space where we find the n −n oscillation lifetime to be most constraining. Thus, the LHC and n −n oscillation provide a complementary search strategy to investigate this simple model of baryogenesis and DM at high-and low-energy frontiers, respectively.
The plan of the paper is as follows: In section II, we describe the model. Section III explains the common origin of DM and baryon abundance in this model and the ensuing constraints on the parameter space. In Section IV, we discuss the model prediction for the n −n oscilla-tion lifetime. Section V analyzes the LHC constraints from monojet, dijet and paired dijet searches, and their complementarity with the low-energy and cosmological constraints. Our conclusions are given in Section VI.

II. THE MODEL
We start with the SM gauge group and add renormalizable terms with B violation. Gauge invariance then requires the introduction of new colored fields. A minimal setup includes two iso-singlet color-triplet scalars X α (where α = 1, 2) with hypercharge +4/3, which allows (potentially) B-violating interaction terms X * α d c d c in the Lagrangian (in the 2-component Weyl fermion notation), where d (with hypercharge −2/3) stands for the RH down-type quark and d c (with hypercharge +2/3) for the corresponding anti-quark. At least two X flavors are needed to produce a baryon asymmetry from the interference of tree-and loop-level decays of X governed by the X * α d c d c term. However, this term alone is not sufficient to generate the asymmetry, because the net asymmetry vanishes after summing over all flavors of d c in the final and intermediate states [15]. We therefore add other renormalizable terms X α ψu c involving a singlet Majorana fermion ψ (with hypercharge 0), which plays the role of DM in this model and also has an important role in n −n oscillation, as discussed below. The interaction Lagrangian is given by [11,16,17] where i, j are the quark flavor indices (the color indices are omitted for simplicity). Note that due to color antisymmetry, only λ αij with i = j are non-zero. Also note that without the presence of the Majorana mass term of ψ in Eq. (1), one could have simply assigned B = −2/3 to X α and B = +1 to ψ, so that the terms X * d c d c and Xψu c (and obviously X * X) do not violate baryon number. The Majorana mass termψ c ψ breaks B by two units in this model. For m ψ m Xα , which we assume to be the case for our subsequent discussion, one can integrate out X α to write an effective 4-fermion interaction ψu c i d c j d c k . This will induce the decays ψ → p+e − +ν e and ψ →p+e + +ν e , as long as m ψ > m p + m e , where m p and m e are the masses of proton and electron, respectively. On the other hand, the same 4-fermion interaction will induce the (rapid) proton decay p → ψ + e + + ν e if m p > m ψ + m e , which is unacceptable. Therefore, the only viable scenario to make ψ absolutely stable is when m p − m e ≤ m ψ ≤ m p + m e [11] (see also Ref. [18] for some related discussion). The remarkable point is that the stability of ψ within this mass window is due to that of the proton, with no need to introduce any additional symmetry. This leads to an important property of the DM in this model, i.e. m ψ m p .
We would like to note that the spin-independent direct detection cross-section of ψ is small, σ SI ≤ 10 −16 − 10 −15 pb for m X ∼ O(TeV) [11]. This is due to the coupling of ψ to a particular chirality of up-type quarks, which suppresses the spin-independent scattering cross section as m −8 X . On the other hand, the spin-dependent crosssection is only suppressed by m −4 X and could be larger: σ SD ≤ 10 −6 − 10 −5 pb for m X ∼ O(TeV), but still much below the current experimental limits [19].
Similarly, the indirect detection signals from DM annihilation in this model will be very hard to disentangle from the overwhelming astrophysical background due to the low DM mass. Thus, as we will show below, the LHC experiments at the high-energy frontier and the n −n oscillation experiments at the low-energy frontier provide the only effective ways of testing this DM model.

III. COMMON ORIGIN OF DARK MATTER AND BARYON ABUNDANCE
For m ψ m p , the only kinematically available DM annihilation channel is ψψ → u c u c mediated by t-channel X α . The corresponding annihilation rate is given by For m Xα ∼ O(TeV), even |λ α1 | ∼ 1 results in a rate that is much smaller than the thermal annihilation rate of 3 × 10 −26 cm 3 s −1 for obtaining the correct relic density. This implies that thermal freeze-out leads to overproduction of DM in this model, in accordance with the Lee-Weinberg bound [20]. This requires a non-thermal production mechanism to explain the correct DM abundance. One possible scenario is the late decay of a scalar field φ that reheats the universe to a sub-GeV temperature. A concrete example of this scalar field φ is the moduli, which typically arise in string compactification and supersymmetric models [12]. They are displaced from the minimum of their potential in the early universe and start oscillating when the Hubble expansion rate is H ∼ m φ (where m φ is mass of the modulus). Moduli are long lived due to their gravitationally suppressed couplings to other fields and dominate the energy density of the universe before decaying. Moduli decay rate is given by with c φ ∼ 0.01 − 1 in typical string compactification scenarios, such as KKLT-type [21], and M Pl = 2.4 × 10 18 GeV being the reduced Planck mass. Moduli decay occurs when Γ φ ∼ H 0.33g 1/2 * T 2 /M Pl . Thus, using Eq. (3), the reheat temperature of the universe after modulus decay is given by where g * denotes the number of relativistic degrees of freedom at T R . Preserving the success of Big Bang Nucleosynthesis (BBN) requires that T R 3 MeV. On the other hand, non-thermal DM production in our model sets an absolute upper bound of T R < 1 GeV. The corresponding range of the modulus mass from Eq. (4) is The number density of DM particles (normalized by the entropy density s) produced by φ decay is given by where Y φ ≡ n φ /s is the yield from moduli decay and Br φ→ψ ≡ n ψ /n φ denotes the number of ψ quanta produced per φ quanta (either directly or as a secondary from decay of other particles, like X α , that are produced from φ decay). Successful non-thermal production of DM via moduli decay requires that this number density match the observed DM abundance in the universe For Here the lower bound corresponds to the case when φ democratically decays to all of the particles in the model (i.e., SM particles plus X α and ψ), and the upper bound corresponds to the extreme case when φ dominantly decays to ψ (directly or indirectly). We then see from Eq. (5) that moduli decay can yield the correct DM abundance in non-thermal fashion in our model.
Regarding baryogenesis, the interference of the treeand one-loop level diagrams in out-of-equilibrium decay of X α → ψu c i , d c i d c j can generate a baryon asymmetry, provided that the couplings λ αi and/or λ αij have CPviolating phases and there are at least two X α flavors (i.e. α = 1, 2). Note that the X α decays can directly generate the baryon asymmetry, independent of the electroweak sphalerons, and therefore, can be realized at temperatures either above or below the electroweak scale. In the latter case, this provides a concrete example of the postsphaleron baryogenesis mechanism [22][23][24]. In the nonthermal scenario discussed above, the number density of X particles depends on the moduli decay rate; therefore, the final baryon asymmetry also depends on the number of X quanta produced per decay of φ, parameterized by Br φ→Xα , apart from the size of the CP asymmetry α in X α decays. We find that This should be compared with the observed baryon asymmetry η obs B 10 −10 [25]. Taking Br φ→Xα between 10 −2 (for democratic modulus decay to all degrees of freedom in the model) and 1 (for modulus predominantly decaying to X α ), and the allowed range of Y φ ∼ 10 −8 − 10 −7 [see below Eq. (6)], we find from Eq. (7) that a relatively large CP asymmetry ∼ 10 −3 − 10 −1 is needed to generate the observed baryon asymmetry.
Such large values of CP asymmetry can only be realized by a resonant enhancement mechanism, similar in spirit to the resonant leptogenesis scenario [26,27], where the self-energy graphs dominate the CP -asymmetry for quasi-degenerate X α 's, i.e. when the mass difference is of the same order as their decay width: where Γ X is the average decay width of X 1 and X 2 . In our model, see (1), the tree-level decay width of X α is given by The flavored CP asymmetry is then given by [16,28] with α, β = 1, 2 and α = β. In the exact resonant limit (8), it can be seen from Eq. (10) that the regulator (second term on the RHS) goes as m X /Γ X , where we have defined m X as the average mass of X 1,2 . From Eq. (9), we know that Γ X ∝ m X , so the CP asymmetry in the resonance limit becomes independent of m X . This in turn implies that the baryon asymmetry given by Eq. (7) also becomes independent of m X , as long as m φ m X , so that the moduli decay is not kinematically suppressed. In what follows, we work in the resonance limit (8) to calculate the CP asymmetry for given values of the couplings.
For simplicity, we make a few reasonable assumptions on the flavor structure of the new couplings λ αi and λ αij , taking into account the experimental constraints. In particular, the exchange of X α in combination with the Majorana mass of ψ leads to ∆B = 2 and ∆S = 2 process of double proton decay pp → K + K + in this model [16]. For m ψ ∼ O(GeV) and m X ∼ O(TeV), current experimental limits on di-nucleon decay from Super-Kamiokande [2,6] imply that |λ α1 λ α12 | ≤ 10 −6 . To satisfy this bound, we assume that |λ α12 | is very small (and drop it from our subsequent discussion), while |λ α1 | can be O(1).
Another ∆B = ∆S = 2 process is pK − →pK + , which is related to the double proton decay pp → K + K + via crossing symmetry. This is relevant for reactions involving cosmic rays and its cross-section goes like |λ α1 λ α12 | 4 m 2 ψ /s 2 at center-of-mass energies √ s 1 TeV.
The limit from double proton decay then implies a crosssection below 10 −64 cm 2 , which is in agreement with the bounds from cosmic rays experiments [29]. The model also predicts the ∆B = 2 process of n −n oscillation for which the relevant couplings are λ α1 and λ α13 . The ensuing constraints will be discussed in Section IV. The couplings λ α13 and λ α23 can also lead to ∆S = 2 processes of K 0 s − K 0 s and B 0 s − B 0 s mixing. However, color conservation does not allow any tree-level contribution to the mixing [12,30], and hence, experimental bounds on these processes are easily satisfied in all of the allowed parameter space shown below. A combination of λ αi and λ α13 , λ α23 also leads to novel monojet/monotop signals (depending on the relative size of |λ α1 |, |λ α2 | and |λ α3 |), as well as dijet events, at the LHC [34,35]. The corresponding constraints from the latest LHC data will be discussed in Section V.
To derive the cosmological constraints on the parameter space of the model, we choose a common value for all |λ αi |, denoted by |λ|. We also take |λ 1ij | = |λ 2ij | ≡ |λ ij |, with λ 12 being negligible as mentioned before. Modifications due to non-universality of these couplings will be straightforward to follow. Note that the pair of couplings (λ 1i , λ 2i ) or (λ 1ij , λ 2ij ) must have a relative phase between them, otherwise the CP asymmetry in Eq. (10) vanishes. It is worth mentioning that successful baryogenesis can occur despite |λ α12 | being very small in order to satisfy the di-nucleon decay and related bounds. In fact, we can seen from Eq. (10) that a sizable CP asymmetry can be generated as long as at least one of the |λ α13 | and |λ α23 | is large.
From Eqs. (5) and (7), we find the ratio of DM to baryon density (each normalized by the critical energy density of the universe): We note that there are two contributions to Br φ→ψ in Eq. (11): (i) due to the direct production of ψ quanta from modulus decay (since the modulus has coupling to all fields), and (ii) due to the production of X quanta from φ decay, followed by X α decay to ψ. Since m ψ m X , decay of X α produces ψ quanta with the following branching fraction: Thus, the numerator of Eq. (11) can be written as with Br Xα→ψ given by Eq. (12), which is the same for α = 1, 2 due to the simplified flavor structure considered above (the same is true for α ). From Eqs. (11) and (13), we find the following important inequality for the ratio of DM to baryon density in this model: Using Eqs. (10) and (12), the requirement that Ω DM /Ω B ≈ 5 [25] sets an absolute upper limit on the ratio of couplings |λ/λ | ≤ 1/ √ 2, where |λ | ≡ |λ 13 | 2 + |λ | 2 23 . This is shown in Fig. 1 by the intersection of the RHS of Eq. (14) (blue curve) and the observed value of Ω DM /Ω B (horizontal line). We emphasize that the upper limit on |λ/λ | is independent of the choice of m X , m φ . In Fig. 1, the green shaded region is allowed, because for values of Ω DM /Ω B [as given by the RHS of Eq. (14)] less than the observed value, we could always arrange the matching contribution from the direct modulus decay [cf. Eq. (13)].
It is clear from Eqs. (7) and (10) that obtaining the observed baryon asymmetry also imposes a lower bound on the same ratio |λ/λ |. This is because the CP asymmetry goes to zero in the limit of |λ/λ | → 0. This is shown in Fig. 2 where various constraints are plotted in the Br φ→Xα − |λ/λ | plane. Here we superimpose the Ω DM /Ω B constraint from Fig. 1 (green shaded region) with the baryogenesis constraint (blue shaded region). The upper blue curve corresponds to the minimum allowed value of Y φ , whereas the lower blue curve corresponds to the maximum allowed value of Y φ in the model. The region below the blue curves is excluded because it gives rise to a value of η B that is lower than observed value, even with the maximum possible value of CP asymmetry (i.e., with resonant enhancement). In the region above the blue curves, one can always obtain the observed value of η B by departing from the resonant condition to lower the CP asymmetry to the desired level. Two important results emerge from Fig. 2: (i) Since Br φ→Xα cannot be larger than 1, we obtain an absolute lower limit on |λ/λ | 0.02, independent of m X and m φ . The lower bound becomes more stringent for smaller values of Br φ→Xα . (ii) The requirement of successful baryogenesis also imposes an absolute lower bound on the branching ratio of modulus decay to X quanta: Br φ→Xα 2.7 × 10 −3 . This is independent of the details of any underlying string compactification model.

IV. NEUTRON-ANTINEUTRON OSCILLATION
After integrating out X α , the B-violating terms in Eq. (1) induce an effective B-violating dimension-6 operator of the form [36] 1 Due to Majorana nature of ψ, the ψu c d c d c term can give rise to n −n oscillation at the tree-level. However, because the couplings λ αij are color antisymmetric, the two down-type quarks in (15) must involve different families, which suppresses the tree-level contribution to n −n oscillation because of its dependence on the strange content of the neutron [12]. The leading ∆B = 2 operator arises from the conversion of two strange or bottom quarks to two down quarks. This will need a ∆s = 2 or ∆b = 2 effective interaction. Due to the constraints on the ∆s = 2 operator from pp → KK lifetime [2,6], the dominant contribution comes from the ∆b = 2 operator of dimension-9, which can be parameterized as m ψ m 6 This gives rise to n −n oscillation at one-loop level, with the amplitude given by [16] G nn We translate this into the oscillation lifetime τ nn ∼ (Λ 6 QCD G nn ) −1 and use the current experimental lower limit on τ nn ≥ 3×10 8 sec [7-9] to derive n−n constraints on the model parameter space. For concreteness, we fix the ratio |λ/λ | at the maximum and minimum allowed values of 1/ √ 2 and 0.02, respectively (cf. Fig. 2) to derive the constraints in the m X − λ 13 parameter space, assuming that λ 23 = 0. 1 This already excludes a sizable portion of the parameter space, as shown by the red shaded regions in Fig. 3. The future sensitivity of τ nn ≥ 5 × 10 10 sec, as predicted by next generation experiments [5], can probe even larger portion of the parameter space.

V. COLLIDER SIGNALS
There is a novel monojet (or monotop, depending on the flavor structure of the λ coupling) signal at the LHC from the on-shell production of X and its subsequent decay to ψu c through the λ-coupling: pp → X ( * ) → ψu c i , where the jet recoils against the DM particle and its transverse momentum has a Jacobian peak near one half of the resonance energy √ŝ = m X [34]. The cross section for this kind of monojet signal does not suffer from the high p T cut in typical LHC monojet searches for DM, where the jet comes from the initial state radiation. However, since we only have λ 13 and λ 23 couplings sizable in our model for the Xd c i d c j type interactions, one of the initial states in the quark fusion process must include the b-quark, which has a relatively smaller parton distribution function (PDF) inside the proton, as compared to the light quarks. This PDF suppression is overcome by considering a gluon in the initial state (with much larger PDF), which splits into bb and one of the b's then fuses with the d quark in the other proton to produce X. 2 The additional b-quark in the final state can either be part of the inclusive monojet search (which just requires one high p T jet) or could be used to tag the signal events. In addition, the colored X particles can be pair-produced via gluon-gluon fusion, each of which can subsequently decay into the DM and a u c quark. If one of the final state quarks satisfies the monojet selection criteria, this would also contribute to the total monojet signal. To sum up, our monojet production cross-section has two separate contributions, namely, a single X production from the b − d or b − s fusion, including the possibility of gluon  Fig. 1) and baryogenesis (cf. Fig. 2) constraints, respectively. The shaded regions are ruled out at 95% CL (except for the n −n constraint, which is at 90% CL).
splitting with an additional b jet in the final state, and pair production of X via gluon fusion. In the gluon fusion case, the additional b jet does not exist.
A recent dedicated analysis from CMS [13,14] puts constraints on the λ − λ parameter space using the √ s = 13 TeV monojet data with at least one jet with p T > 100 GeV, |η| < 2.5, E miss T > 250 GeV and no leptons. In this analysis, the monojet production is assumed to have occurred via X production due to d − s quark fusion and gg fusion. However, as we have noted earlier, λ 12 is small due to the double proton decay limit in our case and the monojet is produced from the b − d, b − s and gluon fusions via single and pair production of X. In any case, the CMS limit from the dedicated analysis can be directly applied to our case. We first use the total cross-section limit as a function of m X as given in Ref. [14] which is shown by the green, dot-dashed line in the m X − λ 13 plane in Fig. 3 labeled as monojet (CMS). The CMS analysis covers the mass range 750 GeV -2.5 TeV. Here we have considered two different values of |λ/λ | = 1/ √ 2 (left panel) and 0.02 (right panel), as for n −n oscillation discussed above. It is interesting to see the complementarity between the low-energy n −n constraint (red solid line) and the high-energy monojet constraint. In particular, for smaller X masses, the monojet constraint is more stringent, while for m X 3 TeV, the production cross section at the LHC is kinematically suppressed, and the n −n process serves as a better probe in this part of the parameter space. For a smaller |λ/λ |, the monojet constraint becomes significantly weaker, as the Br X→ψu c is suppressed compared to Br X→d c d c . This is illustrated in the right panel of Fig. 3.
The ATLAS monojet analysis [39] provides a constraint only on the effective theories of monojet production with upper limits on the cross-section for various missing energy cuts. Since this is not a dedicated analysis for the model we are using where the production cross-section depends on the mass scale of m X , it is not possible to use ATLAS data in our case unlike what we have done with the CMS analysis. We expect, however, that if a similar analysis is done on the model using the ATLAS data, the constraint in the m X −λ 13 plane will be similar to the CMS constraint since both CMS and AT-LAS have produced similar limits for monojet analysis on effective theories.
Apart from the monojet signal, the on-shell production of X could also give a distinct dijet resonance through the λ coupling: pp → X ( * ) → d c i d c j . We follow the recent CMS analysis [40] (see also Ref. [41] for the corresponding ATLAS search) for dijet resonance search (above 600 GeV) with √ s = 13 TeV LHC data and select the dijet events passing the cuts p T > 30 GeV, |η| < 2.5, m jj > 450 GeV and H T > 250 GeV. Requiring our dijet signal cross section times acceptance to be compatible with the observed 95% CL upper limit for the quarkquark case [40], we derive the constraint shown by the purple, dashed line in Fig. 3, labeled as dijet. It is clear that for larger |λ/λ |, the monojet constraint is better than the dijet one, while for smaller |λ/λ |, the dijet constraint is better. This is simply due to the interplay between the corresponding branching ratios of X, namely, Br X→ψu c and Br X→d c d c .
The colored X particles can also be pair-produced at a hadron collider: pp → XX * , purely through QCD interactions, independent of the λ, λ couplings. The subsequent decay of X → d c d c will then lead to a distinct paired dijet resonance. Following the recent AT-LAS analysis for pair-produced resonances in four-jet final states [42], we select events with at least four jets, each with p T > 120 GeV and |η| < 2.4 and compare the signal cross section with the observed 95% CL upper limit to obtain a lower limit on the mass of X. Note that this depends on [Br X→d c d c ] 2 , which in turn only depends on the ratio |λ/λ | in our case. Therefore, for a given |λ/λ |, the lower limit on the X mass is independent of λ , as shown by the vertical orange, dotted line in Fig. 3. This The shaded areas are ruled out at 95% CL (except for the n −n constraint, which is at 90% CL).
limit becomes stronger for a smaller |λ/λ |. This model can also give rise to monotop final states with a RH top quark. This feature can be utilized to distinguish the model from the SM single top background which is consistent with left-handed top quarks [35].
To further illustrate the complementarity between the high-energy constraints from LHC and low-energy constraints from n −n oscillation, along with the DM and baryogenesis requirements, we show in Fig. 4 the λ − λ parameter space for two benchmark values of the X mass, m X = 1 and 2 TeV. As noted above, the DM-to-baryonratio requires |λ/λ | < 1/ √ 2, whereas successful baryogenesis requires |λ/λ | > 0.02, independent of m X , as shown by the gray and brown dotted lines, respectively. The n −n constraint is shown by the red solid line, whereas the LHC constraint from monojet is shown by the green dashed line and from dijet by the purple dotdashed curves. The LHC constraints have been obtained with the recent √ s = 13 TeV data following the procedure described above. As in Fig. 3, we find that for smaller X masses, the LHC constraints are better, while for larger X mass, the n −n constraint becomes more effective in probing the parameter space of this model.

VI. CONCLUSION
We analyzed a simple TeV-scale model of B-violation leading to successful baryogenesis, which also explains the apparent coincidence of the DM and baryon energy densities as due to a common origin from moduli decay. In this model, the DM particle is a Majorana fermion whose mass is required to be in a tiny window between m p − m e and m p + m e , and its stability is linked to that of the proton with no ad-hoc discrete symmetry imposed. This naturally results in monojet/monotop signals at hadron colliders. Obtaining the observed DM-tobaryon ratio imposes an upper bound on the ratio of the couplings of the B-violating interactions involving the up-type and down-type quarks, while generating the observed baryon asymmetry imposes a lower bound on the same ratio, independent of the other model parameters. Neutron-antineutron oscillation is predicted at one-loop level, with a sizable rate that could probe a large and interesting part of the allowed parameter space. There exists a novel complementarity among the DM-to-baryon ratio, baryon asymmetry, n −n oscillation lifetime, and the LHC monojet and dijet signals in this model. The existing n −n oscillation lifetime constraint is already probing the parameter space for m X 2 TeV, whereas the LHC monojet constraints are more relevant for the lower mass range. Thus, the ongoing LHC and future n −n oscillation experiments provide a very important complementary search strategy to investigate this simple model of baryogenesis and DM.