A Supersymmetric Theory of Baryogenesis and Sterile Sneutrino Dark Matter from B Mesons

Low-scale baryogenesis and dark matter generation can occur via the production of neutral B mesons at MeV temperatures in the early Universe, which undergo CP-violating oscillations and subsequently decay into a dark sector. In this work, we discuss the consequences of realizing this mechanism in a supersymmetric model with an unbroken U(1)R symmetry which is identified with baryon number. B mesons decay into a dark sector through a baryon number conserving operator mediated by TeV scale squarks and a GeV scale Dirac bino. The dark sector particles can be identified with sterile neutrinos and their superpartners in a type-I seesaw framework for neutrino masses. The sterile sneutrinos are sufficiently long lived and constitute the dark matter. The produced matter-antimatter asymmetry is directly related to observables measurable at B factories and hadron colliders, the most relevant of which are the semileptonic-leptonic asymmetries in neutral B meson systems and the inclusive branching fraction of B mesons into hadrons and missing energy. We discuss model independent constraints on these experimental observables before quoting predictions made in the supersymmetric context. Constraints from astrophysics, neutrino physics and flavor observables are studied, as are potential LHC signals with a focus on novel long lived particle searches which are directly linked to properties of the dark sector. a rX iv :1 90 7. 10 61 2v 1 [ he pph ] 2 4 Ju l 2 01 9


Introduction
The nature and genesis of dark matter (DM) and the dynamical origin of the matter-antimatter asymmetry of the Universe are two outstanding mysteries that the Standard Model of Particle Physics (SM) cannot explain. Meanwhile, the SM also suffers from aesthetic or fine-tuning problems, especially the gauge hierarchy problem. The hope of solving these puzzles has driven theorists and experimentalists alike to search for new physics (NP).
Many particle physics models have been proposed to explain the origin of the DM relic abundance, measured to be Ω DM h 2 = 0.1200 ± 0.0012 [1], i.e. roughly 26% of the critical energy density of the Universe [1,2]. However, searches for DM at colliders and at designated direct detection experiments together with studies of the possible indirect effects of DM in astrophysical observations have yet to shed light on its nature, thereby severely constraining many scenarios.
The origin of the primordial matter-antimatter asymmetry Y B ≡ (n B −nB)/s = (8.718 ± 0.004)×10 −11 , inferred from measurements of the Cosmic Microwave Background (CMB) [1,2] and Big Bang Nucleosynthesis (BBN) [3,4], requires a dynamical explanation: baryogenesis. Many mechanisms of baryogenesis (satisfying the three Sakharov conditions [5]: C and CP Violation (CPV), baryon number violation, and departure from thermal equilibrium) have been proposed, but remain experimentally challenging to verify due to the inaccessibly high scales and very massive particles involved.
A new mechanism for low-scale baryogenesis and DM production was introduced in [6] 1 . Contrary to the standard lore that baryogenesis is difficult to test experimentally, the mechanism of [6] predicts distinct signals at B factories and hadron colliders. In particular, a large positive enhancement of the charge/semileptonic-leptonic asymmetry in neutral B meson oscillation systems is required, as are new decay modes of charged and neutral B mesons into baryons and missing energy together with exotic decays of b-flavored baryons.
Weak scale supersymmetry (SUSY) offers a theoretically well motivated solution to the gauge hierarchy problem with a plethora of implications for LHC searches. However, the lack of discoveries at the LHC has pushed the Minimal Supersymmetric Standard Model (MSSM) to ever more tuned corners of parameter space [13], prompting theorists to introduce new variations of weak scale SUSY where experimental constraints may be evaded or relaxed. For instance, supersymmetric models with Dirac gauginos [14,15] and an exact U (1) R symmetry (R-SUSY) [16,17] are phenomenologically appealing [18]. Exact U (1) R symmetries -under which particles within the same supersymmetric multiplet carry different charges-are motivated by top down constructions. By forbidding Majorana gaugino masses and supersymmetric a-terms, the majority of collider constraints are circumvented in models of R-SUSY. Therefore, such extensions of the MSSM are intriguing in their own right, but are arguably even more motivated if the deviation from minimality accommodates a solution to other outstanding problems of the SM, such as DM and baryogenesis.
In the present work, we realize the mechanism of baryogenesis and DM production introduced in [6] in an R-SUSY setup where an exact U (1) R symmetry is identified with U (1) B baryon number 2 . In this baryogenesis scenario, neutral B mesons are produced at MeV temperatures in the early Universe, undergo CP-violating oscillations and subsequently decay to a Dirac gaugino (which is charged under the R-symmetry and therefore baryon number) through a four-fermion interaction obtained by integrating out heavy squarks. Successful baryogenesis requires at least one squark to have a mass of O(1 TeV), and critically a Dirac gaugino with mass of a few GeV (typically smaller than 4 GeV in order for B meson decay to be kinematically allowed). Current constraints do not allow the existence of so light gluinos or charginos, but the bounds do not apply to the color neutral bino [13].
The collider phenomenology of various R-SUSY setups has been greatly studied [21][22][23][24][25][26][27], along with constraints from low energy flavor observables and electric dipole moments [28]. However, previous work on R-SUSY models focused on scenarios with heavy Dirac gauginos (for instance [18] considers flavor constraints on a model where all the gauginos share the same mass scale of order 500 GeV − TeV). Realizing the mechanism of [6] motivates us to study the phenomenology of a complementary region of R-SUSY parameter space -one where multiple gaugino mass scales are present. This allows the bino to be light enough for baryogenesis purposes while the other gauge superpartners are heavier and avoid collider constraints. This scenario can be realized, for instance, if different SUSY breaking scales are assumed [29].
In the present work we further extend the R-SUSY setup to include a sterile neutrino and sneutrino with GeV scale masses, and discuss their role in providing a mass for SM neutrinos via a type-I seesaw [30][31][32][33]. In this scenario, the bino decays predominantly into a sterile neutrino-sneutrino pair. The sterile neutrino decays are generically long lived and may be searched for at SHiP [34], FASER [35], CODEXb [36], MATHUSLA [37] and in the ATLAS muon tracker [38], while the sneutrino is stable on cosmological time scales and constitutes the DM. This opens up a variety of astrophysical constraints and signals which are discussed and assessed.
While the existence of an R-symmetry is theoretically well motivated, it is generically expected that this symmetry is not realized exactly, as supergravity considerations lead to a gravitino whose mass is not U (1) R symmetric. Supergravity breaking of the U (1) R symmetry introduces, for instance, Majorana sparticle masses generated from the conformal anomaly along with aterms and small soft squark masses, resulting in additional constraints (see [39] for a detailed discussion). We quantify the degree to which the R-symmetry can be broken while still being able to accommodate the baryogenesis and DM production mechanism from B mesons. This paper is organized as follows: in Sec. 2 we review the key ingredients and results of [6], i.e. the generic model and parameter space needed to achieve baryogenesis and DM from B mesons. In Sec. 3 we present new model independent constraints on the experimental observables within our framework. Next, Sec. 4 introduces the model of R-SUSY with Dirac gauginos and right-handed neutrino multiplets, in addition to discussing the degree to which we require the R-symmetry to be exact. In Sec. 5 we present our results for flavor observables within the SUSY model -namely the contributions to the semileptonicleptonic asymmetries in neutral B meson oscillations and various other flavor constraints. Next, in Sec. 6 we discuss what region of parameter space is able to generate the observed matter-antimatter asymmetry and the DM relic abundance. This is followed by a discussion of additional constraints on the model parameters arising from neutrino physics, DM stability, supersymmetry breaking and neutron stars. In Sec. 7 various signals at colliders and B factories that can probe this model are discussed. Finally, we conclude in Sec. 8 by summarizing the allowed parameter space of the model that realizes baryogenesis. Some further details of the computations are given in the appendices.

Baryogenesis and Dark Matter from B Mesons
Here we review the mechanism of low scale baryogenesis and DM production from the oscillations and subsequent decay of neutral B mesons, which was proposed in [6]. An inflaton-like 3 field Φ with mass m Φ decays out of thermal equilibrium in the early Universe, at temperatures of order T R ∼ 1 − 100 MeV. The main decay products are assumed to be quarks and antiquarks that hadronize producing neutral B 0 q=s,d = |b q mesons and antimesonsB 0 , which then undergo CP-violating oscillations before decaying partially into a dark sector 4 .
A minimal set of four new particles is required to generate the necessary interactions between the visible and dark sectors (see [6,39] for details). First consider an electrically charged (Q e = −1/3), baryon number carrying (B = −2/3), color triplet scalar which couples to SM quarks -foreshadowing, we denote this particle byq R . Such a particle couples to the SM quarks through the fully SU (3) c antisymmetric termq * Rū b c , and various other flavor combinations. Additionally, we introduce a Dirac fermion ψ carrying baryon number −1, so that couplings of the formq Rψ s c are allowed. To avoid collider constraints, the new colored scalar is assumed to have a mass around the TeV scale, allowing us to integrate it out arriving at an effective Hamiltonian valid a low energies: (2.1) To simplify the notation, we have chosen a particular combination of light quark flavors to showcase the mechanism, and have parametrized the effective coupling by λ eff . It will however be true that any other combination of couplings that yields an effective four-fermion interaction involving b and ψ and two other light quarks 5 is equally suitable. For a sufficiently light ψ, with mass smaller than about 4 GeV, the effective operator in Eq (2.1) allows theb quark within the B meson to decay; B 0 q → ψ + Baryon + X, where X parametrizes mesons or other additional SM particles. Since ψ carries baryon number, this decay is baryon number conserving. However, an asymmetry n ψ −n ψ can be generated due to the CP-violating nature of the neutral meson oscillations. If ψ is allowed to decay back into visible sector particles, the generated asymmetry will be erased. However, the asymmetry may be preserved if ψ dominantly decays into stable DM particles. This is made possible by the existence of a dark scalar baryon φ with baryon number −1 and a dark Majorana fermion ξ. If a discrete Z 2 symmetry is assumed under which the dark sector particles transform as ψ → ψ, φ → −φ and ξ → −ξ, we can write a Yukawa operator mediating the decay of ψ into the dark sector, where an overall negative baryon number is stored. A positive baryon asymmetry is thus produced in the visible sector without globally breaking baryon number. Note 3 The details of an inflation model are beyond the scope of this work. This role can also be played by a modulus or any other late-time decaying particle. 4 Additional mesons may also be produced, but are uninteresting for the current setup as only B mesons can give rise to baryogenesis in this way due to kinematic considerations. 5 More precisely, any up-plus down-type pair of first or second generation quarks, in order to preserve hypercharge and for the B meson to be able to decay through the effective operator. We refer to Table III in [6] for an exhaustive list of options. that the scalar φ carries baryon number and as such its mass is constrained to be greater than ∼ 1.2 GeV [40] 6 . In this paper we consider a scenario in which the scalar constitutes the DM (the fermion will turn out to be long lived but unstable on cosmological time scales). The combination of kinematics, proton decay, dark matter stability, and neutron star stability leads us to consider a parameter space where where m n is the neutron mass and m B is the B meson mass. The asymmetry and DM abundance can be determined by solving a system of coupled Boltzmann equations [6,41]. Doing this, the comoving baryon asymmetry is found to scale as (2.4) Here, Br(B 0 → B + X) is the branching fraction of the B 0 meson to a baryon B and X, which stands for the dark sector particles plus potentially other mesons or SM particles. The "reheat temperature" T R is taken to be the temperature at which the Φ decays become significant, that is, when 3H(T R ) = Γ Φ , where Γ Φ is the decay width of Φ. Furthermore, a q is the leptonic charge asymmetry (which in this setup is effectively equivalent to the semileptonic-leptonic asymmetry a q sl ), an experimental observable which parametrizes CP violation in the B 0 s, d systems. An important point is that, as in neutrino physics, neutral B meson oscillations can only occur in a coherent system. Additional interactions with the mesons can act to decohere the oscillations [42,43] by "measuring" the system, thereby suppressing the CP violation and diminishing the potential to generate an asymmetry. In the early Universe, decoherence is caused by the scattering of e ± off B 0 mesons due to their charge asymmetry. For CP violation to effectively happen, B mesons must oscillate at a rate similar to or faster than the e ± B 0 q → e ± B 0 q scattering in the plasma. The factor of α q T R , ∆m Bq in Eq. (2.4) is a measure of decoherence at the temperature T R for each of the q = s, d systems.
The kinematically allowed and currently experimentally unconstrained parameter space that accommodates the measured baryon asymmetry and DM relic abundance was mapped out in [6]. In order to reproduce the observed value of Y B ∼ 8.7 × 10 −11 , the experimental observables in Eq. (2.4) must, after scanning over the (m Φ , Γ Φ ) parameter space, lie within the ranges we conclude that the amount of CP violation predicted by the SM is typically not enough for successful 6 The authors of [40] showed that the observation of neutron stars with masses greater than 2M constrains the mass of a dark particle carrying baryon number to be greater than the neutron chemical potential inside the star, which is about 1.2 GeV. Otherwise, if a light dark baryon did exist, a process such as n + n → DM + DM could occur within the star. As a consequence, neutrons would be replaced by DM particles and the neutron Fermi pressure would diminish. This would destabilize the neutron star, leading to gravitational collapse. baryogenesis. However, given the current experimental constraints [4] a s, exp sl = (−0.6 ± 2.8) × 10 −3 and a d, exp we see that there is room for new physics to enhance the CP violation and accommodate the required values in Eq. (2.5). We will explore this in more detail in Sec. 3. The other relevant quantity is the branching ratio of B mesons into a baryon and dark sector states which would contribute to missing energy in an experiment. To the best of our knowledge, to the date an inclusive search for such a decay has not been performed. Given the experimental constraints on similar processes like B → Kνν, which are in the O(10 −5 ) level, we expect that existing and future data from BaBar, Belle (-II) and LHCb may be in position to fully test the parameter range that allows for successful baryogenesis. The symmetric component of the DM is generally overproduced in this set-up. Therefore, additional interactions allowing the DM to annihilate are required in order to deplete its abundance. The results of [6] show that an annihilation cross section is needed in order to reproduce the observed DM density Ω DM h 2 ∼ 0.12.
One of the main appeals of this mechanism is its susceptibility to be tested in terrestrial experiments such as flavor facilities and accelerators. Indeed, experimental constraints on a sl can directly constrain the mechanism. We explore this model-independent bounds in Sec. 3. Next, we construct a realization of this mechanism, i.e. (2.1) and (2.2), in a supersymmetric framework. Having a full model will allow us to consistently calculate the new physics contributions to a sl and the branching fraction Br(B 0 → B + X), together with any other related flavor observables that the new particles and couplings may modify. In doing so, we will be able to identify the phenomenologically viable region of parameter space where baryogenesis can successfully be achieved, and develop the best approach to experimentally test the mechanism in high and low energy experiments.

Model Independent Constraints
We now conduct a model independent study to further quantify the required new physics (NP) contributions to realize baryogenesis and DM from B mesons. This will set the stage for successfully embedding this mechanism in a full model. As seen in (2.4), the baryon asymmetry depends critically on the semileptonic asymmetres a q sl and the branching fraction Br B 0 → B + X of neutral B mesons to baryons and dark sector states. As discussed above, the SM prediction for the semileptonic asymmetries is not enough to generate the observed baryon asymmetry required for baryogenesis, and in this section we study the ability for general NP contributions to enhance this asymmetry, taking into account all the relevant experimental constraints. The amount of NP introduced is constrained by measurements of different B 0 q oscillation parameters, in turn limiting the size of the enhancement of the semileptonic-leptonic asymmetry which is consistent with experimental data. As a consequence of these bounds, we obtain a lower limit on the required branching fraction in order to generate the measured baryon asymmetry of the Universe. In this section we remain agnostic about the origin of the NP contributions and perform the study in a data driven and model independent way.

The Semileptonic Asymmetries
Experimental measurements of the oscillation parameters in neutral B meson mixing constrain the size of NP effects that introduce new CP violation into these systems. We now quantify such constraints. Recall that the evolution of the neutral B 0 q −B 0 q system is described by where the Hamiltonian 2×2 matrix H q osc = M q −iΓ q contains off diagonal terms which result in CP violation in the neutral meson mixing. Writing the off-diagonal elements as the mass and decay rate differences between the two propagation eigenstates are given by ∆M q 2 |M q 12 | and ∆Γ q 2 |Γ q 12 | cos φ q 12 , where we have defined To get to Eq. (3.3) we have assumed that |Γ q 12 | |M q 12 |, which is well established experimentally. Note that although both M q 12 and Γ q 12 can be complex, only the phase difference between the two is physical. This phase is the main quantity controlling the CP violation in the system. In the SM, where the off-diagonal elements of the oscillation Hamiltonian are generated by box diagrams, nonzero φ q M and φ q Γ arise from the phase in the CKM matrix [45][46][47][48].
The NP contributions that are relevant for baryogenesis only modify M q and leave Γ q unchanged from its SM value. In particular, there is no extra source of direct CP violation, which measures differences in the tree-level decay of B 0 → f with respect to the CP conjugateB 0 →f . In addition, no new common final state to which both B 0 andB 0 can decay is introduced. Therefore, the CP-violating observables that may be modified by NP contributions are: • CP violation in mixing, which is described by the mixing phase φ q 12 . CP violation in mixing can be measured via the semileptonic CP asymmetries, also known as CP asymmetries in flavor-specific decays. The asymmetry is defined as where the final state is such that the decays B 0 q → f andB 0 q →f are forbidden.
• Mixing-induced CP violation, also known as CP violation in interference, which arises from interference between mixing and decay. Here, we consider a final state f to which in principle both B 0 and B 0 can decay, and define This form of CP violation is controlled by the phase φ q , which is defined as In order to assess the extent to which NP can impact these observables, it is customary to define [44] M q 12 = M q,SM

12
· |∆ q | e iφ q ∆ , Γ q 12 = Γ q,SM 12 . (3.8) We reiterate that we are only allowing modifications in M q 12 , but not Γ q 12 , which is fixed to its SM value. The complex parameter ∆ q is a measure of the effects of NP. This parametrization has the advantage that each of ∆ q and φ q ∆ are individually related to the observables ∆M q and φ q . It also allows to write a simple expression for the semileptonic asymmetries in the B 0 s,d systems. For each q = s, d we have It is however more convenient for us to parametrize M q,NP 12 in a way that makes the relation with the parameters of the NP model more transparent. For that purpose, we define Note that ∆ = 1 + α. The NP effects are thus normalized with respect the SM values as The sign of the contribution to a q sl depends on the relative alignment of the NP phase with respect to the SM one (recall that a positive contribution is required for successful baryogenesis, see Sec. 2). In the SM, the CP violation can be traced back to the single physical complex phase in the CKM matrix. Using the Wolfenstein parametrization, we obtain the SM predictions for the phases of the mass splittings, The experimental measurements of the oscillation period of neutral B mesons (which constraints ∆M q ) and the CP violation in interference (which constrains φ q ) can be turned into exclusion limits on the size of NP contributions. We use the most up-to-date averages by the HFLAV group [49] to derive our constraints.
A key ingredient for the analysis are the SM predictions for the values of the relevant quantities ∆M q,SM , ∆Γ q,SM , φ q,SM 12 , φ q,SM and a q,SM sl , which we obtain from [44]. The results of our analysis are shown in Fig. 1 for the B d and B s systems. The main conclusion is that the combination of measurements pushes the absolute value of any NP contributions to be roughly one order of magnitude smaller than the SM values. This statement is slightly dependent on the phase of the NP contributions. Additionally, for small enough values of M q,NP

12
, the phase φ q,NP M is essentially unconstrained. In Fig. 1, we have also overlaid contours showing the value of the a q sl arising from the NP parameters. As expected, phases in the (0, π) range enhance a q sl with respect to its SM value, whereas phases   [41]. In order to account for the uncertainty in the calculation coming from the lack of precise knowledge of the scattering rate Γ of e ± off B mesons in the primordial plasma, we conservatively allow for up to a factor of 4 uncertainty in Γ. This uncertainty mostly stems to the lack of knowledge about the charge radius of the neutral B mesons. Right: Branching ratio of B mesons to visible baryons and dark antibaryons necessary to reproduce the baryon asymmetry of the Universe, also as a function of the reheating temperature. We use the largest experimentally allowed values for the semileptonic asymmetries and give the results for 3 different values of the mass of Φ. We use different shades of the colors to highlight uncertainties coming from the scattering rate as explained above and, additionally, from potential variations in the fragmentation ratio of b quarks to B 0 s and B 0 d mesons [4] (we use a benchmark ratio of 1 : 4 but to be conservative allow it to be as large as 1 : 2). In this range of parameters, branching ratios in the 0.001 − 0.1 range can lead to successful baryogenesis.
in the (π, 2π) range give a negative contribution to the asymmetry. It is important to note that in the B 0 d system the asymmetry is bound to be negative, while in the B 0 s system it can be positive. More precisely, we obtain the following 95% cl ranges for the semileptonic asymmetries 7 : This is a much tighter region than the one quoted in Eq. (2.7), which is an average of only direct measurements of the asymmetries. Of course, the one we have derived here is reliant on the extra assumptions mentioned above regarding the NP contributions.

The Branching Fraction Br(B 0 → B + X)
Taking the values in Eq. (3.13) as reference for the maximum positive asymmetry that can be accommodated by NP contributions, we may use (2.4) to set a limit on the possible branching fraction of neutral B mesons to baryons and dark sector states Br B 0 → B + X that is required to reproduce the observed baryon asymmetry of the Universe. The baryon asymmetry is solved for by evolving a set of coupled Boltzmann equations [6,39,41] governing the evolution of the number densities of the relevant particles Φ, B,B and the DM. The solution may be found numerically, and factorizes into the following form: (3.14) In addition to the branching fraction and the semileptonic asymmetry, the asymmetry depends on m Φ and Γ Φ . The Φ width may be related to the late era reheating temperature in the usual way Γ Φ = 3H(T R ).
Following the prescription of [41], the value of the asymmetry is simply inversely proportional to m Φ and depends on the decay width through a numerically determined polynomial denoted by α q (T R ) for the B 0 q=s,d systems, whose dependence on T R is depicted in Fig. 2. Roughly speaking, α q acts to parametrize the decoherence effects of the system.
In analogy with neutrino oscillations, additional interactions with the B mesons can "measure" the state of the meson and decohere the oscillating system through the Zeno effect [51]. Such decoherence will spoil the oscillations and therefore prevent a baryon asymmetry from being generated. For the setup at hand, the only possible source of dechorence is the e ± scattering off the neutral B meson in the early Universe plasma. To a good approximation, decoherence becomes important when the scattering rate exceeds the oscillation lengths of the B 0 s or B 0 d system, which are measured to be ∆M s = 1.17 × 10 −11 GeV and ∆M d = 3.34 × 10 −13 GeV. The temperature-dependent scattering rate is roughly Γ e ± B→e ± B ∼ 10 −11 GeV (T /20 MeV) 5 up to uncertainties in the charge radius of B 0 s,d mesons. Therefore, decoherence becomes more relevant at high temperatures. More specifically, the oscillations decohere for T (15,30) MeV for the B d and B s systems, respectively. Decoherence effects may be properly treated in the Boltzmann equations through the prescription of [43]. Results are displayed on the left panel of Fig. 2, for which we assume that reheating occurs between ∼ 4 − 100MeV, that is, between BBN and hadronization. To avoid excessive decoherence, baryogenesis favors reheating temperatures in the range T R ∼ 15 − 50 MeV.
The behavior of the numerical results for α q (T R ) can be understood as follows: at high reheat temperatures (T R 50MeV), the B mesons formed from the decay of the Φ do not oscillate coherently and asymmetry production is therefore highly suppressed. As we consider lower values of T R , the scattering rate for e ± B → e ± B drops such that the B 0 s mesons are produced and oscillate coherently thereby producing a baryon asymmetry. At even lower temperatures, B 0 d −B 0 d coherent oscillations will also contribute to the asymmetry. The change of behavior of α s,d (T R ) at very low T R is the result of dilution by a greater entropy density.
Recall that for baryogenesis to proceed through B meson oscillations and subsequent decay to the dark sector, a large positive net semileptonic-leptonic asymmetry is favored. As we saw above, the NP contributions to the neutral B meson oscillating systems can enhance the semileptonic-leptonic asymmetry in B 0 s decays to at most a s sl = 8.80 × 10 −4 while not being in conflict with measurements in B physics. Meanwhile, a d sl < 0 in all regions of parameter space, so that maximizing the produced baryon asymmetry favors the smallest possible value of this negative contribution, namely a d sl = −4 × 10 −4 . At the same time, Fig. 2 shows that the maximal values of α s occurs for a reheat temperature of about 30 MeV. Taking everything into account, we conclude that in order to reproduce the observed value of Y = (8.718 ± 0.004) × 10 −11 , a sizable inclusive branching ratio of B mesons to baryons is needed. The required value for this branching ratio is displayed on the right panel of Fig. 2 for extremal values of a s,q sl , therefore representing a lower limit on the size of the branching fraction. The relatively large branching ratio is a robust prediction of the mechanism and provides a unique opportunity for high-energy experiments like LHCb [52] and Belle-II [53] to test this baryogenesis scenario.

A Supersymmetric Realization
The model independent analysis presented above allows to place constraints on the observables that play a direct role in the baryogenesis mechanism considered here. However, it is expected that other related observables can be used to indirectly test the scenario once the required minimal NP content is embedded in a consistent framework. With this aim, we construct a complete realization of the model in a supersymmetric context. The key ingredients of the theory which allow for baryogenesis to be realized are an unbroken U (1) R symmetry [16,17] and the presence of Dirac gauginos [14,15]. Supersymmetric models with Dirac gauginos and U (1) R symmetries are phenomenologically appealing [18]. By forbidding Majorana gaugino masses and supersymmetric a-terms, the majority of collider constraints that push the MSSM into ever more fine tuned regions of parameter space are avoided. Their collider phenomenology has been greatly studied [21][22][23][24][25][26], along with flavor considerations [28]. Among other applications, these models can accommodate diverse baryogenesis scenarios [19,20], sometimes at the cost of breaking the R-symmetry [21]. In this regard, realizing the mechanism of [6] is a novel application of such ingredients and, as will be discussed, motivates a rather unstudied region of parameter space: a light (GeV mass scale) Dirac bino.
Let us briefly summarize the main ideas behind our supersymmetric theory for baryogenesis and DM from B mesons. R-symmetries are known to arise naturally in supersymmetric constructions [14][15][16][17][18]28]. Supercharges and their conjugates transform in opposite ways under an R-symmetry, allowing for particles within the same supermultiplet to carry different R-charges. Although U (1) R is usually taken to be broken at low energies, an interesting possibility is for it to remain exact 8 and be identified with baryon number. Following this avenue leads to quarks and squarks having different baryon number assignments, as shown in Tables 1 and 2. In this way, a right-handed down-type squarkd R has the required charge assignments to play the role of the heavy colored scalar in Eq. (2.1). We may also identify an adequate candidate for the neutral GeV scale Dirac baryon ψ in the form of the Dirac bino. After integrating out the heavy squark, the effective four-fermion operator in Eq. (2.1) arises from a combination of R-parity violating and gauge interactions. Furthermore, the Dirac partner of the bino is by definition a singlet under the SM gauge group, allowing it to naturally couple to a dark sector. The dark sector can be minimally accommodated by introducing a generic new singlet chiral multiplet. More interestingly, this new singlet can be identified with a right-handed neutrino multiplet. In this way, the bino mediates the decay of B mesons into dark sector states. The complete set of superfield charge assignments is summarized in Table 1, and the corresponding particle fields and their properties are given in Table 2. In what follows, we will further elaborate upon the ingredients of this model needed to realize baryogenesis.

An Exact U (1) R Symmetry
The model is similar to the one studied in [21,22]. An exact R-symmetry requires a superpotential with R-charge of 2 (superspace derivatives D α have R-charge −1). Given the U (1) R charges of superfields (see Table 1), only the following terms respect the R-symmetry: The first four terms of Eq. (4.1) constitute the usual MSSM superpotential, with the addition of the "Rparity violating" (RPV) term U c i D c j D c k (the subscripts denote generation indices) which is now allowed. The rest describe the Higgs sector; the R u,d multiplets are added to generate µ terms which are forbidden in models with R-symmetry. Electroweak symmetry breaking proceeds as usual when the scalars of H u and H d acquire vacuum expectation values, while the VEVs of R u.d remain zero. Given the U (1) R charge assignments of the matter content (see Table 1), U c and D c have an R-charge of 2/3. If we explicitly write out the decomposition of the multiplets, e.g.
we see that the d-type right-handed squarks have R-charge −2/3 while the quarks 9 have R-charge 1/3 (recall that θ has R-charge 1). Since we identify baryon number with R-charge,d R has the appropriate charge assignments to be identified with the heavy colored scalar that mediates the effective interaction of Eq. (2.1). At the quark/squark level, the U c D c D c term in the superpotential produces the interactions As in Sec. 2, we have chosen a particular combination of quark flavor to keep the notation simple. In general, the role ofd R can be played by any first or second generation up-or down-type squark, which means that any of the couplings λ ij3 with i, j = 1, 2 can be relevant for the baryogenesis mechanism.

Dirac Gauginos
In Wess-Zumino gauge, the field strength of a vector superfield can be written in the usual way as where fields are expressed in terms of the   they are embedded. We work in terms of Weyl spinors so that the Dirac bino consists of two components ψ = (B, λ † s ). Here, Φ is an "inflaton-like" field that decays out of equilibrium producing b andb quarks. The heavy colored squark qR is integrated out to mediate quark-bino four fermion interactions (note that various flavors of squarks may generate this interaction, as will be discussed further below). The Dirac bino mediates decays of hadrons into the dark sector. The dark sector is embedded into a right-handed neutrino supermultiplet where the sneutrino is stable and represents our DM candidate, while the sterile neutrino generally decays on time scales of interest for collider searches.
in a similar way 10 .
With the addition of the new superfields, Dirac gaugino mass terms can be generated [18,55]. Consider the following superpotential terms: are generated (the second term of (4.6) vanishes on-shell when D 1 = 0). We denote the 4-component Dirac bino and its corresponding Dirac mass as in consonance with Eq. (2.1). Similar expressions are obtained for the other gauginos, which acquire masses m ψ2,3 = √ 2c 2,3 D/Λ 2,3 . In order to reproduce the spectrum required for baryogenesis, m ψ needs to be a few four-fermion operator GeV. At the same time, collider [13] and flavor [56,57] constraints generically push charginos and gluinos to lie above the TeV scale. If no other source of mass for gauginos is present, this implies that there must exist a hierarchy (c 1 /Λ 1 )/(c 2,3 /Λ 2,3 ) 10 −3 . Such hierarchy can be explained either by invoking some tuning of the c i coefficients or by allowing for at least two separate sources of SUSY breaking [29] so that Λ 1 Λ 2,3 .

Squark Masses and Couplings
Soft masses for squarks (and other superpartners) can be generated through an F -term SUSY breaking spurion X = θF [28] at a scale Λ F , which produces squark masses mq ∼ F * F/Λ F . Given the experimental constraints [21,26,28] and depending on the scale of gluino and chargino masses, squarks as light as a few 100 GeV can be feasible 11 . As discussed above, the role of the colored scalar participating in the baryogenesis mechanism is played by a right-handed squark, which can be up-or down-type and belong to any of the first two generations. The squark mass is assumed to be in the TeV range so that it can be safely integrated out at the GeV scale where the baryogenesis dynamics happen. It is straightforward to see how theq Rψ s term in Eq. (2.1) is generated in the supersymmetric setup. Interacting chiral matter theories have gauge interactions of the form where Q q R and g = e/ cos θ W denote the charge assignment and hypercharge gauge coupling (as customary, e is the electromagnetic charge quantum and θ W is the weak mixing angle). Whens R is integrated out, the effective four-fermion interactionB u s b is recovered. Of course, other flavor combinations for the operator are possible: they are summarized in Table 3 along with the resulting coefficient of the effective interaction. It is worth noting that the gauge coupling being fixed, the only unknown parameter that sets the strength of the effective interaction is λ ij3 .

The Higgs Sector and Neutralino Mixing
The number of neutral fermions in R-symmetric models is larger than in the MSSM due to the presence of the new SU (2) doublets R u,d and the singlet and triplet fields S and T . Here we follow [26] to investigate the mixing between all these states after electroweak symmetry breaking. The relevant terms in the Higgs sector, allowed by gauge and R-Symmetries, are given by Recall that R u,d are added to generate the otherwise R-symmetry forbidden µ-terms. Electroweak symmetry breaking proceeds as usual when the scalars of H u and H d acquire vacuum expectation values v u and v d (the VEVs of R u.d remain zero). The terms in (4.9), together with the Dirac gaugino masses and the gauge potential, generate the Dirac mass mixing matrix showing that one of the states with mass µ d decouples. It is clear that if λB u = λW u = 0, there is no mixing between the bino, the neutral weakino and the higgsinos. Even if λB ,W u are nonzero, in the limit where m ψ v mW , µ u , the lightest state is mostly bino with mass m ψ . This is the case of interest to us, as we are dealing with a GeV bino but other gauginos and higgsinos are at or above the TeV scale. In this limit, the mixings of the bino are approximately given by Although not relevant for the baryogenesis dynamics, these mixings can play a role in constraining the model using flavor observables such as µ → eγ, which are discussed in Sec. 5.2.

A Dark Chiral Multiplet
As was highlighted in Sec. 2, the crucial step to generate a baryon asymmetry in the visible sector without incurring into global baryon number violation is to allow ψ to decay into a dark sector where the corresponding antibaryon number is stored. In this way, a DM abundance is generated in parallel to the baryon asymmetry.
Thus far, the model we have introduced cannot accommodate DM. However, being uncharged under the SM gauge group, the singlet S offers a natural portal to a dark sector. Given that the fermion in S is the Dirac partner of the bino, such a coupling allows ψ to decay to the dark sector. The minimal extension to take advantage of this portal includes a new chiral multiplet φ where the B = −1 dark scalar φ and the B = 0 Majorana spinor ξ can both be embedded, The only superpotential terms allowed by gauge and U (1) R symmetries are precisely (4.14) The former generates the coupling λ s φ * ξ while the latter gives masses to the dark states. Note that the invariance of (4.14) under φ ↔ −φ acts to stabilize the DM (corresponding to the Z 2 symmetry introduced in [6]). This most economic and generic extension contains all the ingredients to produce the baryon asymmetry and an asymmetric DM component. However, given that the dark states φ and ξ and the mediator ψ are rather light, a symmetric relic abundance tends to be overproduced [6] in this minimal setup. In order to partially deplete the symmetric abundance, new interactions have to be added to allow for the excess DM to annihilate 12 . In our case, there exists an attractive solution to this problem: to identify the dark sector states with a right-handed neutrino supermultiplet.

Sterile Sneutrino Dark Matter
In order to accommodate neutrino masses in the U (1) R symmetric model, the minimum requirement is the presence of a right-handed sterile neutrino multiplet N R with R-charge +1, Identifying U (1) R with U (1) B leads to the right-handed sneutrinoν R carrying baryon number −1, while the right-handed neutrino remains uncharged under baryon number. At the same time, N c R has to carry −1 lepton number. The following operators, allowed by all the symmetries (but including a ∆L = 2 Majorana neutrino mass term), can be added to the superpotential: To keep the discussion general, we allow for three flavors of sterile neutrinos (I = e, ν, τ ; J = 1, 2, 3). Note that the first and third terms are equivalent to the ones in (4.14), so the identification of N c R as the DM multiplet is direct. The first term of (4.16) may be expanded in component fields as This operator generates the three point interaction (2.2) that mediates the decay of the Dirac bino (through its right-handed partner λ s ) into the dark sector states ν R andν R . The second term in (4.16) generates Dirac neutrino masses. The SU (2) L multiplets can be expanded out as where I = e, µ, τ . Upon electroweak symmetry breaking, h 0 u = v u produces a Dirac mass for the neutrino, with m D = 2y N v u = 2y N v sin β. For simplicity, we have assumed a diagonal flavor structure and suppressed generation indices. This simple setup reproduces the type-I seesaw mechanism [30][31][32][33], whereby the light neutrinos obtain masses m ν m 2 D /M M , as long as m D M M . For the baryogenesis mechanism to be viable, at least one of the right-handed neutrinos needs to have a mass in the GeV range so that ψ can decay into it. In order to obtain a neutrino mass of m ν 0.1 eV, a small Yukawa y N 10 −8 is required, as will be further discussed in Sec. 6.3.
Aside from generating Dirac neutrino masses, the second term in (4.16) also provides the interactions necessary to avoid the overproduction of a symmetric DM component. The crucial difference to the generic case discussed in the previous section is the fact that the DM particles, being identified with sterile neutrino and sneutrino, now carry lepton number. The annihilation of right-handed sneutrinos into left-handed neutrinos can occur via t-channel exchange of a bino or a higgsino. We will study this in more detail in Sec. 6.2.

U (1) R Symmetry Breaking from Supergravity
The model presented above relies on an exact global U (1) R symmetry which is identified with baryon number. While R-symmetries are motivated by top down constructions, we do generically expect them to be broken at low scales. Supersymmetry breaking effects coming from supergravity are generically always present (however see [59] for a possible alternative), and these generate a tiny Majorana mass for the gauginos through the conformal anomaly (as well as soft squark masses and a-terms). The order parameter is the gravitino mass m 3/2 so that gaugino Majorana masses Additionally, soft squark masses (m 2 ) i j = 1 2 |m 3/2 | 2 d dt γ i j and non-zero a-terms a ijk = −m 3/2 β y ijk will be generated. Clearly, the introduction of (even tiny) a-terms, soft squark masses and Majorana gaugino masses from AMSB can spoil the phenomenological advantages of R-SUSY. Therefore, it is interesting to quantify the degree to which the R-symmetry must be exact, i.e. the size of the gravitino mass and couplings λ ijk of Eq. (4.1) that still accommodate a phenomenologically viable model. This will be further discussed in Sec. 6.5, where we will see that the most stringent constraints come from the AMBS Majorana gaugino masses that lead to dinucleon decay.

Flavor Phenomenology
After having outlined a complete model with a field content that realizes the mechanism of [6], we may begin to preform calculations within this framework in order to identify the phenomenologically viable parameter space. As we will see, generating the observed baryon asymmetry and the DM relic abundance of the Universe requires a precise range of model parameters and observables, which are summarized in Table 5. A key prediction of this scenario is the modification of some of the oscillation parameters of the neutral B meson system with respect to the values expected in the SM. Because of this, sufficiently precise measurements of the B meson system can directly probe the baryogenesis dynamics. Any NP inducing such modifications is also expected to have an effect on other ∆F = 2 observables, like oscillations in the B, K, and D neutral meson systems; and on ∆F = 1 processes such as µ → eγ and b → sγ. Making use of the supersymmetric model, these effects can be quantified within a consistent framework.
It is common lore that SUSY models with an R-Symmetry are not too severely constrained by flavor observables (see e.g. [28]). However, due to the presence of the O(GeV) bino and the λ ijk as the only RPV couplings, our mechanism highlights a different slice of parameter space from what has been previously considered. It is therefore important to re-evaluate the constraints and apply them to the relevant parameter space, as well to assess whether the model prediction for the semileptonic-leptonic asymmetry in the B s and B d systems lies within the range that can accommodate baryogenesis.
The RPV couplings λ ijk U c i D c j D c k are an obvious source of flavor violation, and may also produce CP violation through complex phases. A comprehensive study of the phenomenological implications of RPV couplings was made in [62]. The main particularity of our setup is that the most constraining processes, i.e. proton decay and nucleon-antinucleon oscillations, are forbidden as long as the R−symmetry remains exact (we discuss the implications of R−symmetry breaking in Sec. 6.5). As a consequence, the flavor violating observables described below give the strongest constraints for all the RPV couplings.
Soft masses for SM fermion superpartners can also be sources of flavor and CP violation. A particularity of the unbroken U (1) R is that left-right squark and slepton mixing is forbidden by the R-symmetry [16]. The soft scalar mass matrices can be parametrized as usual by dimensionless ratios δ ij LL ≡ (m 2 q ) ij /|m 2 q | of the off-diagonal elements to the diagonal ones m 2 q , where i, j are quark or lepton flavor indices (analogously for δ ij RR ). As we will see, CP violation can arise if these ratios have complex phases.

∆F = 2 Meson Oscillations
In addition to the flavor violation already present in the SM, our setup allows for new processes contributing to neutral meson-antimeson oscillations, where flavor is violated by two units. For baryogenesis, a NP contribution in the B 0 -B 0 system is desirable. At the same time, experimental constraints, which are strongest in the K 0 -K 0 system, must be satisfied. Studies in general supersymmetric extensions of the SM [56,57], in RPV contexts [62,63] and in U (1) R models with Dirac gauginos [28] have been performed.
Here we extend and adapt those results to our scenario. The physics of neutral meson oscillations, as well as the relevant observables, have been described in Section 3.1 in a model independent way. We now link such observables to the NP model at hand. Firstly, the modification of the mass splitting between meson and antimeson can be computed as Here, Ci encode the model-specific partonic processes, which are discussed below. Q i , that is, Q i |B 0 and similarly for K 0 and B 0 s . The operators are numbered following [57]. The ones relevant for us are where α and β are color indices. The corresponding matrix elements are 2) In our calculations, we use the values obtained in [64] for the bag parameters B i .
The decay rate difference ∆Γ is not modified in our setup, as there is no new common state to which both meson and antimeson can decay. Its value can therefore be set to the SM prediction as was done in Sec. 3. As a consequence, the previous expression Eq. (5.1) also allows to calculate CP-violating observables, which depend on the complex phase of (∼) Ci with respect to the SM counterpart, as made explicit in Eq. (3.11). We are particularly interested in contributions to the semileptonic asymmetries in the B 0 s,d systems, as these are the key observables that enter the baryogenesis mechanism.
At the partonic level, the meson oscillations arise from box diagrams. In the SM they are mediated by W bosons, but in our supersymmetric scenario new contributions arise, which can be classified into two groups depending on the origin of the flavor violation.

Contributions from RPV couplings
The corresponding Feynman diagrams are shown in Fig. 7 in Appendix A. The first two contain four RPV vertices with λ couplings, while the lower two diagrams involve mixed contributions with two RPV vertices and two gauge/Higgs vertices with the corresponding CKM factors. These diagrams were calculated in [63] for the (s)top couplings λ 3jk , here we generalize the computation to include also contributions from first two generation (s)quarks. The explicit analytic expressions are given in (A.1). Flavor changing squark mass mixing is neglected for this computation.
As can be seen in Eq. (A.1), one-loop diagrams involving RPV interactions generate the operatorsQ 1 , Q 4 and Q 5 . The former is associated with the box diagram involving four λ couplings and therefore scales asC 1 ∼ λ 4 /m 2 q , where λ 4 symbolizes the appropriate product of four potentially different λ ( ) ijk couplings (see Eq. (A.1) for the explicit combinations that arise). For simplicity, we assume universal squark masses and compute the loop integrals in the limit where mq is much larger than any other mass 13 . The other two operators are associated with diagrams with only two RPV couplings, while the other two vertices correspond to weak interactions where flavor violation is due to CKM insertions. Therefore, they both scale as C 4,5 ∼ λ 2 V 2 /m 2 q , where V 2 represents the appropriate CKM factors (see Eq. (A.1) for explicit expressions). These diagrams also carry a factor of m ui m uj /m 2 W , which introduces an extra suppression unless only top quarks run in the loop (i.e. i = j = 3). Note that the RPV couplings generically source CP violation in the neutral meson mixing because the λ couplings can be complex. 13 We do not take into account renormalization group running effects for the operators (∼) Q i from the squark scale down to the hadronic scale. We estimate this approximation to be good within a factor of ∼ 2 for the coefficients (∼) Ci and therefore do not expect a large effects on the limits presented in Table 4    s systems constrain different combinations of couplings. For each product, the first limit assumes that the NP contribution is aligned with the SM one, while the second one assumes maximal CP-violating phase. A line indicates that the derived limit is above the unitarity bound on the couplings. The limits are computed assuming a universal squark mass mq = 1 TeV, and scale with mq (diagonal ones) or m 2 q (off-diagonal ones).
Using the experimental measurements of the neutral meson oscillation parameters, (∆M, φ d,s ) for B 0 d,s and (∆M, K ) for K 0 , together with the SM predictions for such observables 14 , we can obtain constraints on the products of RPV couplings. The resulting limits are compiled in Table 4. For all of them, we assume a value of mq = 1 TeV. The limits on the diagonal of each table are dominated by theC 1 contribution with four λ factors, which means that they scale with mq. However, the off-diagonal entries in the tables, which rely on the C 4,5 contributions, scale with m 2 q and are weaker due to the CKM and/or quark mass suppressions explained in the previous paragraph. The exception are the bottom right entries of the tables, which correspond to diagrams mediated by top quarks for which the contributions fromC 1 and C 4,5 are of similar size for mq = 1 TeV. Each of the entries in the table quotes two limits: the left one corresponds to the case where the phase of the NP contribution to M 12 is aligned with the SM one; the limit in that case comes from measurements of ∆M . The right one assumes that the NP phase is orthogonal to the SM one and therefore enhances CP violation, the limits in that case come from φ d,s and K . For the general case where the phases are neither aligned nor orthogonal, the limit lies between the two extremal values quoted in Table 4.
In the general case where the RPV couplings are complex and source CP violation, contributions to the semileptonic asymmetries in the B 0 d and B 0 s are expected. Indeed, as Fig. 1 shows and was discussed in Sec. 3.1, NP effects can significantly enhance a d,s sl with respect to their SM value while satisfying all other constraints. For instance, the benchmark values used for the baryogenesis calculations are obtained if any one product of complex RPV couplings saturates the corresponding limit in Table 4.

Contributions from squark mass mixings
Even if the product of RPV couplings is too small to generate sizeable flavor violation, off-diagonal terms in the squark mass matrix can also induce neutral meson oscillations. The bino-mediated diagrams are shown in Fig. 9 in App. A for the B 0 d system, the case of B 0 s and K 0 are analogous. Note that for Dirac gauginos the number of contributing diagrams is limited compared to the case of the MSSM, as the Dirac partners do not couple to quarks so that chirality flips in the internal gaugino propagators are not possible 15 . Furthermore, the U (1) R symmetry prohibits the presence of a-terms, i.e. no left-right squark mass mixing is allowed.
The authors of [28] considered the ∆F = 2 flavor constraints arising from box diagrams with gluinos, which due to the stronger coupling give the largest contribution if the masses of all gauginos are assumed to be similar. This is not necessarily the case in the regime where m ψ ∼ 1 GeV is much lighter than gluinos and winos. The analytic expressions for the contribution of a light bino to the oscillation parameters defined in Sec. 3 are given in App. A. Parametrically, the contribution of a light GeV bino to ∆M is of the same order as that of a gluino which is about 10 times heavier than the squarks.
CP violation arises when ∆M NP is a complex quantity. This happens if the off-diagonal elements of the squark mixing matrix are complex, and thus so are δ LL and δ RR . The CP-violating observables 15 For heavy gauginos that can be integrated out, the effective dimension five operators liked Rs L d † R sL are absent, and the transition happens through more suppressed dimension six operators likedL∂µs L d † L γ µ sL. For the case of the light bino this extra suppression is not present, but the chirality structure of the diagrams is in any case different than for Majorana gauginos. defined in Sec. 3 can therefore receive NP contributions. With this, the experimental constraints detailed in Sec. 3, together with the analogous ones for the kaon system (∆M K and K , see [44]) can be used to place constraints on the squark mixings. Figure 3 shows the excluded ranges of values for the diagonal and (generally complex) off-diagonal elements of the squark soft mass matrix, mq and δ. The non-trivial dependence of the constraints on the phase of δ comes from the interference with the SM contributions, i.e. the alignment between the NP effects and their SM counterparts.
The values of the semileptonic asymmetries as a function of the model parameters are shown as contour lines in Fig. 3. For squark masses in the TeV scale and O(1) absolute values of the mixings δ sb and δ bd , a broad interval of phases can lead to semileptonic asymmetries compatible with baryogenesis. At the same time, the stringent bounds in the kaon sector restrict the size of δ ds to be 10 −2 unless the phase of this mixing is close to zero (or π/2). This is compatible with the baryogenesis mechanism as δ ds doesn't play a role in B meson oscillations, but points towards some flavor structure that particularly suppresses any flavor violation which does not involve the third generation; see for instance [65,66] for constructions where such flavor structures arise.

∆F = 1 Observables
Processes in which flavor is violated by only one unit can also probe the model parameters relevant for baryogenesis. We limit the discussion to two observables, b → sγ and µ → eγ, as they offer the best experimental prospects. Although the latter is the best tested experimentally, the former is directly related to parameters relevant for the baryogenesis mechanism, making it worth studying in detail. For this purpose we mainly adapt the results of [28] to our setup.

Contributions from RPV couplings
The Feynman diagrams contributing to the flavor-violating decay b → sγ are shown in Fig. 8 in Appendix A, along with the full analytic expression for the decay rate, given in Eq. (A.3). Assuming a universal mass mq for all squarks, we obtain Comparing with the experimental measurement [67] Br(B → X s γ) = (3.55 ± 0.24 ± 0.09)×10 −4 and the SM prediction [68] Br(B → X s γ) = (3.15 ± 0.23) × 10 −4 , we see that this process does not place any meaningful constraint on the RPV couplings for squark masses around or above the TeV scale.

Contributions from squark and slepton mass mixings
Flavor violating decays can also be triggered by off-diagonal elements in the mass matrix of squarks (for b → sγ) or sleptons (for µ → eγ). Given the particularities of the R-symmetric model, i.e. the absence of L − R sfermion mixing and the inert nature of the Dirac partner of the gauginos, the bino can only mediate two kind of diagrams contributing to this process. They are depicted in Fig. 10 for b → sγ, the ones for µ → eγ are analogous. The first diagram requires an external chirality flip while the second one involves mixing between the appropriate Dirac higgsino and the bino, as discussed in Sec. 4.4. The flavor-violating decay of the b quark is most important because it directly probes the couplings relevant for baryogenesis. Evaluating Eq. (A.9) in the small tan β limit, we obtain a branching ratio This should be compared with the experimental measurement [67] Br(B → X s γ) = (3.55 ± 0.24 ± 0.09) × 10 −4 , which is compatible with the SM prediction [68] Br(B → X s γ) = (3.15 ± 0.23) × 10 −4 . The corresponding exclusion on the b − s squark mixing is shown in the top right panel of Fig. 3. The computation of the flavor-violating muon decay is analogous to the one for the b → sγ transition with the obvious substitutions, see Eq. (A.14) for the corresponding expression. Evaluating it again in the small tan β limit, we find This process is extremely suppressed in the SM, so this contribution is constrained by the stringent experimental upper limit Br(µ + → e + γ) < 4.2 × 10 −13 (90% cl) set by the MEG experiment [69]. For sleptons at the TeV scale, the mass mixing in the e − µ sector is bound to be δ eµ < 4 × 10 −3 , which again indicates that any flavor violation that only involves the first two generations should be suppressed.

Theoretical and Phenomenological Considerations
We now discuss constraints on the model parameters that arise from the requirement of achieving the measured matter-antimatter asymmetry and DM abundance. The model independent results of Sec. 3 are applied to the R-SUSY model of Sec. 4, as is the analysis of flavor observables performed in Sec. 5. With this, we find that successful baryogenesis requires squark masses in the scale of 1 − 4 TeV and at least one of the λ ij3 RPV couplings to be O(1). Both these requirements stem from the necessity for a large branching fraction of (visible) baryon number violating decays of neutral B mesons. The observed DM abundance is easily accounted for with order one couplings of ψ to the sterile neutrino multiplet. We further discuss constraints that apply to parameters that do not directly control baryogenesis. Such effects are rather a byproduct of realizing the mechanism in the R-SUSY model. At the phenomenological level, we discuss bounds coming from neutrino physics and astrophysics. On the more theoretical side, supergravity considerations on the size of R-symmetry breaking constrain the gravitino mass and some RPV couplings.

Baryogenesis
As was highlighted in Sec. 3.2, successful baryogenesis requires a relatively large branching ratio of B mesons into final states containing a single baryon and dark sector particles. In the supersymmetric setup, such decays can be mediated by any of the four-fermion operators listed in Table 3. Assuming that the decay is dominated by a single one of them yields a branching ratio where ∆m = m B 0 − m B − m X is the mass splitting between B meson and final states, and X represents the dark particles together with any other mesons. Note that Eq. (6.1) is computed using a spectator approximation [39], which breaks down when the mass splitting ∆m 2 GeV. In general, given the hadron spectrum of interest here, we expect 2 GeV ∆m 3 GeV for the operators of interest. For final states involving charm quarks, ∆m may be too small to reliably use the spectator quark approximation and an analytic estimation of the branching ratio isn't feasible.
For hypercharge coupling g (m b ) ∼ 0.34 and assuming degenerate squark masses, achieving a large enough branching fraction for baryogenesis where we have used the maximum positive value of (a s sl ) max = 4.1 × 10 −4 and the minimum negative of (a d sl ) min = −9.0 × 10 −5 to maximize the asymmetry. It is also important to take into account the lower bound on squark masses of ∼ 1 TeV, which comes from collider constraints [6]. Because of this, it is clear that attaining a large enough branching fraction to accommodate the observed baryon asymmetry favors at least one of the RPV couplings being large, together with (at least one) squark mass in the TeV scale 16 , that is Here the range in the allowed squark mass parameter space takes into account enhancement to the coupling from the renormalization group running of the anomaly dimension of the 4 fermi operator 17 .

Dark Matter Relic Abundance
The mechanism of [6] tends to generically overproduce a symmetric DM component, therefore requiring additional dark sector interactions to deplete its abundance. The R-SUSY model presented here naturally provides interactions that can regulate the overproduction of DM. As discussed in Sec. 4.6, the scalar sneutrinoν R dark sector particle is stable on cosmological time scales and therefore serves as a DM candidate. At the same time and as long as mν R > m ν R , a sneutrino-antisneutrino pair can annihilate into right-handed neutrinos. After that, the ν R decays as usual into SM leptons, thus contributing to the depletion of the symmetric component of DM. The annihilation cross section forν Rν R → ν R ν R is dominated by the exchange of light binos. The thermally averaged cross section is given by With this, we can easily leverage the results of [6] to find that the observed DM relic abundance is reproduced for values of the superpotential coupling λ N of the order In what follows we take this as the benchmark value for λ N .

Generation of Neutrino Masses
As was discussed in Sec. 4.6, the SM neutrinos ν acquire a tiny mass through their Dirac mass mixing with the Majorana right-handed neutrinos ν R . This is nothing but the well-known type-I seesaw mechanism [30][31][32][33]. As usual, diagonalization of the neutrino mass matrix with terms from Eq. (4.19) and Eq. (4.18) leads to a large eigenvalue m ν R ∼ M M and a small one m ν = 4 y 2 N v 2 sin 2 β/m ν R , as long as the Yukawa y N sin β 16 The parameter space for squark masses can be opened up a little by increasing λ ij3 up to where perturbativity begins to break down (for instance [62] computes the renormalization group running in R-SUSY models and obtains λ 123 < 1.25). 17 A detailed RGE study is beyond the scope of the current work. is sufficiently small. Successful baryogenesis requires the existence of at least one right-handed neutrino with a mass m ν R ∼ GeV. At the same time, current experimental data constrain the mass of any of the SM neutrinos to be m ν 0.12 eV [34,70,71]. Fulfilling these two conditions bounds ν R 's contribution to the seesaw, constraining the sterile neutrino Yukawa coupling to be 18 This translates into an upper limit for the mixing between the active and the sterile neutrino, which is given as usual by the ratio of the Dirac to Majorana mass terms, corresponding to |U | ≡ y N v sin β/m ν R in terms of our model parameters. We find |U | The usual phenomenological constraints for type-I seesaw models with GeV-scale right-handed neutrinos apply to our scenario. The strongest limits come from a combination of beam dump, LEP and B factory searches (see [72] for a recent review). Meaningful bounds can also be obtained from displaced vertices searches at LHC [73]. In the most simplistic scenario which we discuss here, the mixing between sterile and active neutrino flavors is expected to be small, |U | 2 ∼ 10 −11 , which is beyond the reach of current experimental setups. As a particular feature, our model predicts that a relatively large number of sterile neutrinos are produced in B meson decays. As was discussed in Sec. 6.1, a new decay channel for B mesons with branching ratio of 10 −3 is required for baryogenesis. The bino produced in this process decays predominantly to a right-handed neutrino and sneutrino pair, constituting a large new production cross section for sterile neutrinos. Experiments like LHCb and other B factories and SHiP [34], where ∼ 7 · 10 13 B mesons are expected to be produced, are best posed to take advantage of this production channel, potentially allowing to place meaningful constraints in the parameters of the model. We leave a more detailed investigation of this interesting possibility for future work.

Dark Matter Stability
The stability of the right-handed sneutrino DM is a consequence of kinematic relations and the action of a stabilizing symmetry which is a combination of baryon and lepton number. This stability is however not absolute, as decays into three light (anti)quarks and a SM neutrino are possible via a higher dimensional operator. The diagram for the process is shown in Fig. 4. The first vertex corresponds to the coupling of the Dirac bino to the dark sector sterile neutrino multiplet through Eq. (4.17). The second is an effective four-fermion operator analogous to the ones shown in Table 3, but with a different flavor combination, namely The decay is thus induced by the λ 112 RPV coupling, which is not directly relevant for baryogenesis. The full processν R → ψ * + ν ( * ) R →ū +d +s + ν SM is mediated by an off-shell bino and is further suppressed by the small mixing between the active and the sterile neutrino, so the decay rate is expected to be very small. Let us nevertheless estimate it to obtain a constraint on the RPV coupling to light quarks from the requirement of DM stability.
The hadronized final state must contain a baryon and a strange quark, so the lowest mass final state in the decay contains a Λ baryon, which has a mass m Λ 1.1 GeV. As the DM mass is constrained to be in the range 1.2 − 2.7 GeV, the maximum momentum for any final state particle is around 1.1 GeV. This is too small for the simple approximation of considering 3 free light quarks as decay products to be good. However, we can leverage the fact that the available energy is low enough so that final states with fewer particles, perhaps only one hadron, should dominate the width. We therefore compute the rate for the decaỹ ν R → Λ + ν. A rough estimate can be obtained by using dimensional analysis to match uds ↔ 4πf 3 π Λ. The lifetime for the two-body decay is then where ∆m 2 = m 2 ν R −m 2 Λ and the rest of the notation is as described in previous sections. A most conservative bound can be obtained by requiring that this lifetime is larger than the age of the Universe, 4.3 × 10 17 sec. Using the benchmark parameter values indicated in Eq. (6.9), we obtain a constraint on the RPV coupling, λ 112 10 −5 . (6.10) The operator of Eq. (6.8) also leads to a (visible) baryon number non-conserving decay of the sterile neutrino, ν R →ν Rūds . As one sterile neutrino is produced for each DM antibaryon during baryogenesis, such a process could washout the generated baryon asymmetry if its rate is too large. However, the width of this decay channel is set by the same RPV coupling which is bound to be small by DM stability in Eq. Eq. (6.9)). Because of this, its rate is negligible compared to the usual baryon number preserving SM right-handed neutrino decays. Even if the bound Eq. (6.10) is satisfied and the DM is cosmologically stable, the sneutrino may have a non-negligible decay width into antiquarks and a neutrino. As discussed above, the most likely decay mode is expected to be into Λ + ν. If such process occurs inside a large volume detector such as as Super-Kamiokande [74], the pions and the annihilation of thep andn from the Λ decay could produce a detectable signal. Furthermore, observations of the antiproton flux in cosmic rays can be used to the constrain the rate of DM annihilation or decay into antimatter [75][76][77][78][79]. However, the antiprotons produced in our setup have very low energies, which makes modelling their propagation in the galactic medium a difficult task. On top of that, the background contribution to the low-energy antiproton flux from astrophysical sources is hard to model [80,81]. A careful analysis would therefore be needed to apply any such constraint to our scenario.

How Exact Must the R-Symmetry Be?
As discussed in Sec. 4.7, AMSB [60,61] effects coming from supergravity generate small Majorana gaugino Masses, soft squark masses and a-terms (trilinear scalar couplings), all of which are of order the gravitino mass m 3/2 . While other supergravity scenarios exist in which these contributions are suppressed (see e.g. [59]), it is interesting to study generic constraints on the gravitino mass and associated bounds on the RPV couplings.
The Majorana mass term in Eq. (4.21) allows for two neutrons to decay into 2 kaons (dinucleon decay) or neutron-antineutron oscillations via a Majorana mass insertion. In a nuclear environment, the antineutrons annihilate with other nucleons leading to the decay of the nucleus. This process is severely constrained by the lack of observation of such decays. Flavor structure may also severely suppress this process [39], with the most constrained operators involving only up, down and strange quarks. Such operators are generated via the combination of the λ 112 operator and a Majorana mass for the bino, or by any of the other λ ijk in combination with a Majorana bino mass and flavor violation from the weak interactions or squark mixing parameters. Among the RPV couplings, λ 112 is the most constrained as it can lead to dinucleon decay without any squark mixing or weak interaction insertion. Adapting the results of [39], where the dinucleon decay bound is translated into the Λ − Λ mixing rate δ ΛΛ , to the R-SUSY setup presented here, we find where m ψ is the Dirac mass of the bino given in Eq. (4.7), λ 112 is the RPV coupling as defined in Eq. (4.1) and M 1 is the bino Majorana mass, e.g. generated by AMSB (see Eq. (4.21)). Given that in our model the R-symmetry is identified with baryon number, the gravitino is a baryon and therefore subject to the constraints coming from neutron stars as discussed in Sec. 6.6, i.e. it must satisfy m 3/2 1.2 GeV. Taking this into account and using the vanilla AMSB scenario without any additional suppression, Eq. (6.11) implies the bound λ 112 5 × 10 −6 (6.12) for O(1 TeV) squark masses. Note that this is similar to the bound from DM stability discussed in Sec. 6.4. Additionally, the squark mass matrix receive a contribution from AMSB as discussed in Sec. 4.7. This effect should comply with the limits obtained in Sec. 4.3. Finally, AMSB generates non-zero a-terms which can also contribute to di-nucleon decay. This effect is however sub-leading to the one discussed above arising from the Majorana bino mass. All in all, in the scenario where the R-symmetry is slightly broken by supergravity effects, avoiding a too large dinucleon decay rate places strong constraints on some combinations of RPV and flavor violating couplings.

Dark Matter Capture in Neutron Stars
The extremely dense environment found inside neutron stars has been utilized as a laboratory for constraining DM interactions and properties. Neutron stars, the remnants of a collapsed core of a giant star, are prevented from further gravitational collapse into a black hole by Fermi pressure of neutrons and electrons. DM that interacts with neutrons and electrons can upset this balance leading to destabilization of the star. Many studies have been conducted to investigate such effects and the constraints they imply [82][83][84][85][86][87][88]. As an example of an application, in Sec. 2 we have used the results of [40] to restrict our parameter space to dark baryon masses greater than ∼ 1.2 GeV, as lighter dark baryons would lead to a process that would deplete the neutrons (and therefore the Fermi pressure) via conversion into DM within the star. We now examine constraints on our model arising from the possibility that DM is captured inside the core of a neutron star. Scattering processes of DM on electrons or neutrons within the neutron star imparts the DM's kinetic energy on the target such that the emitted DM particle may not have enough energy to escape the star. DM that accumulates within the core of a neutron star due to this loss of kinetic energy will begin to self gravitate, eventually overcoming the Fermi pressure, thereby destabilizing the star and leading to gravitational collapse [89,90].
In our framework, DM-neutron scattering proceeds through two higher dimensional four-fermion operators coupling to light quarks, and as such has a negligible rate. However, DM-electron scattering e −ν R → e −ν R occurs via a one loop diagram involving a bino, W boson and sterile neutrinos. The cross section for this process is given approximately by where on the right hand side we have fixed GeV mass dark sector particles and TeV scale squarks along with other benchmark values as discussed above. Note that since we consider bosonic DM there is no contribution to the Fermi pressure from the captured DM, and we simply recast the bound from [91] on the scattering cross section that would result in enough DM being captured to cause it to self gravitate. To avoid this, the DM-electron cross-section needs to be smaller than ∼ 10 −58 cm 2 (for 1 GeV DM mass). Therefore, black hole formation is avoided for small enough sterile neutrino oscillation parameters, (6.14) Recall from Eq. (6.6) that in a type-I seesaw framework, y N sin β should not exceed ∼ 10 −8 in order to accommodate constraints on the maximal size of SM neutrino masses. Therefore, in the simplest type-I seesaw framework DM capture in neutron stars proceeds at a low enough rate and does not induce the collapse of the star. The accumulation of DM in the core of the neutron star is even less constrained in our framework, as the captured DM carries antibaryon number and can therefore annihilate with neutrons within the star. Such annihilation process may proceed through several different channels in our model. Let us consider as an example the annihilation to a sterile neutrino and a kaon, which is kinetically favored and proceeds through t-channel exchange of a Dirac bino. We estimate the cross section for the parton level process ν R + udd → ν R + ds as where we have used our usual benchmark values. For given annihilation and capture rates C a and C c , which measure the probability of a single particle annihilating or being captured over a unit time interval, the number Nν R of DM particles captured by a neutron star evolves according to The capture rate C c depends on the scattering rate of DM with neutron star constituents. As shown above, in our case the scattering with electrons dominates over the one with neutrons. Using Eq. (6.13) and the results of [89], the capture rate is estimated to be C c ∼ 10 23 yr −1 σ eν R 10 −55 cm 2 . (6.17) Assuming that DM particles at the core of the neutron star are relativistic and that the density of surrounding neutrons is similar to nuclear density n s , the annihilation rate is simply given by C a ∼ σ inc annh n s , where σ inc annh is the inclusive cross section for DM annihilation with neutrons, which receives a large contribution from Eq. (6.15). The number of DM particles that is required for bosonic DM to form a black hole [89] is about Nν R ≈ 10 38 . To avoid the accumulation of too many DM particles inside the neutron star, we can simply require C a Nν R > C c . This condition translates into a lower bound on the annihilation cross section, which in our case reads This condition is easily satisfied by the annihilation process described above.

Signals at Colliders and B Factories
In this section we identify the main features of the model that are likely to be tested in current and upcoming terrestrial experiments, with focus on colliders and B factories. As we will see, a rich phenomenology including exotic B decays and long lived particle searches provides an opportunity to fully test this baryogenesis scenario. More precisely, the large ∼ 10 −3 branching ratio of B mesons to baryons and missing transverse energy (E miss T ) Eq. (6.1), along with the requirement of 1 − 4 TeV squark masses (Eqs. (6.2) and (6.3)), result in multiple correlated signals arising from processes occurring at different vertices within LHC detectors.

Semileptonic Asymmetries
An avenue to probe the baryogenesis scenario is through the semileptonic asymmetries. As was discussed in Sec. 3.1, an enhancement of a s sl and/or a d sl with respect to their SM values is required to reproduce the observed baryon asymmetry. This amplification of the CP violation in the B meson system can be tested as the still relatively loose experimental bounds tighten up on the SM predictions. Given stronger constraints on the semileptonic asymmetries, larger values of the branching ratio of B to a baryon and missing energy than the ones quoted in Eq. (6.2) would be required for the viability of the baryogenesis mechanism, strengthening the signals of the searches described below. Therefore, experimental efforts in this direction complement the searches proposed below in testing the dynamics of the model.

Exotic B Meson Decays at B Factories
As was discussed in Sec. 6.1, successful baryogenesis requires a new decay mode of B mesons with a relatively large branching ratio of ∼ 10 −3 . The decay products include a single baryon and dark sector states (a sneutrino DM particle and a sterile neutrino) carrying baryon number, potentially along with some mesons. This apparently baryon number violating decay can be looked for in B factories like BABAR, Belle [92] and Belle II [53]. Depending on what precise combination of coupling and light squark dominates the process (see Eq. (6.1)), different baryons and mesons are expected to appear as decay products. Generically, protons and strange or charmed baryons will be produced. A more experimentally challenging situation would arise if a large fraction of decays occurs into a neutron and dark particles, as such a virtually invisible final state would be extremely difficult to detect. However, all of the possible operators contributing to the decay contain either an up quark or a charm quark, and either a down quark or a strange quark. The only operator that would not produce either charmed or strange particles in the final state is isospin symmetric, so averaging over both the B 0 and B + decays it is valid to assume an equal number of decays containing protons and neutrons.
A dedicated search for exotic B decays with a final state containing a baryon and missing energy has not been performed to this date. However, an inclusive branching ratio Br (B → B + anything) = 6.8 ± 0.6% was reported in [93] (see also [94]), which is large compared to the known exclusive modes. A loose limit Br (B → B + X) 1 − 2 % can be derived assuming that the (also reported [95,96]) modes Br (B → pp + anything) = 2.47 ± 0.23% and Br B → Λp/Λp + anything = 2.5 ± 0.4 % do not contain any invisible fermion in the final state. A dedicated search in BABAR or Belle II is highly desirable and has the potential to probe the full parameter space, given the sensitivity 19 of similar searches like Br (B → Λpνν) < 3 · 10 −5 (90% cl) using BABAR data [97]. An additional signal of interest arises in the case when the mixing of the sterile neutrino with SM neutrinos is large enough for it to decay into SM states within the detection volume. This possibility is discussed in Sec. 7.5.

Exotic decays of b hadrons at the LHC
At the LHC, bb pairs are copiously produced via strong interactions. These b-quarks subsequently hadronize and decay. Through the effective four-fermion operators in Table 3, one of the b-quarks may undergo the visible baryon number violating decay b →ūdψ (or b →ūsψ). Recall that ψ is the Dirac bino (Eq. (4.7)), which subsequently decays into two dark sector states: a stable sneutrino DM particle and a sterile neutrino. If the sterile neutrino is sufficiently long-lived, both particles leave the detector contributing to missing transverse energy of the event. The rate for this exotic b decay is directly related to the baryogenesis dynamics (see Eq. (6.1)) and is expected to be significant, corresponding to a branching fraction of ∼ 10 −3 .
At the partonic level, the process of interest is shown in Fig. 5 (left) and can be summarized as b b →būq ψ , where q = s , b , (7.1) 19 We thank S. Robertson for his input on this matter. and ψ subsequently decays into a DM particle and a sterile neutrino, both of which are assumed to be stable at detector scales (the possibility of the sterile neutrino decaying before leaving the detector is discussed in Sec. 7.5). The corresponding signature at the LHC is an event with a single jet along with significant missing transverse energy (E miss T ) [98]. The fact that the single jet is produced from a b-quark allows for the use of b-tagging [99], which provides a further distinctive feature of the process.

LHC Searches for j + E miss
T A similar final state to the one discussed in Sec. 7.3 can be generated in pp collisions through the effective four-fermion operators given in Table 3. As an example, we consider the process u d → ψb, which involves couplings relevant for baryogenesis. The corresponding diagram is shown in Fig. 5 (right). At the partonic level, we can estimate the production cross section as where √ s is the center-of-mass energy of the process and we have neglected all fermion masses. Similarly to the exotic b-hadron decay, this process results in a distinctive final state consisting of a single b-jet plus missing transverse energy corresponding to the decay of the bino into dark sector particles. The generalization to the other RPV operators λ ijk is straightforward and would motivate more general monojet searches [98].

LHC Searches for Long Lived Particles
In the discussion of the previous section we have assumed the Dirac bino to decay into dark sector states which leave the detector as missing transverse energy, but this might not always be the case. In fact, in the model presented above the dark sector particles are identified with a sneutrino and a right-handed neutrino. The sneutrino is sufficiently stable to be a good DM candidate and therefore invisible in any high energy experiment. The sterile neutrino ν R is however unstable and expected to decay back into SM particles. Interestingly, for the range of parameters of interest for baryogenesis, the sterile neutrino can decay with a lifetime that is long on collider time scales, thereby providing exciting implications for long lived particle searches at MATHUSLA [37], FASER [35], CODEX-b [36], and the ATLAS muon tracker [38]. The right-handed neutrino can in principle decay into DM and three light quarks through an off-shell bino and an RPV coupling as shown in Fig. 6 (left). However, this process is related to the same coupling that is required to be small by DM stability (as discussed in Sec. 6.2). As a consequence, this decay channel is constrained to have a rate that is too small for detection purposes 20 . Nevertheless, the neutral and charged current decay of the sterile neutrinos through a mass mixing to SM neutrinos may in general be long lived 21 . In particular, the processes ν R → q q l and ν R → l + l − ν (see the center and left diagrams in Fig. 6) are long lived with decay lengths of order 1 − 10 6 meters for values of the Yukawa couplings that lie within the requisite range to accommodate SM neutrino masses. Therefore, depending on the exact parameters and experimental setup, the sterile neutrino might either leave the detector as missing energy or decay in a displaced vertex, motivating the use of long lived particle detection techniques.
In the framework presented here, the prompt production and subsequent long lived decay of the sterile neutrino would produce novel signals at the LHC. Let us exemplify the particular case involving b quarks, which is the interesting one for baryogenesis. As discuss in the previous sections, two different production channels are possible, each one of them producing a signature with some distinctive features.
On the one hand, the creation of a b quark and a Dirac bino from a pp collision can proceed as in Fig. 5 (right) through the ψ u d b operator in Table 3, with the cross section estimated in Eq. (7.2). The bino then promptly decays into a sterile neutrino and DM through the operator in Eq. (4.17). This is followed by the long lived decay of the sterile neutrino as discussed above, which results in lepton jet (for the CC case). The full parton level process is therefore In this case, the long lived particle is centrally produced. A potential study in, for instance, the ATLAS detector would involve a mono-jet signal (as discussed in Sec. 7.4) and an associated long lived decay leading to a lepton jet [101]. Note that our signal is novel in that it contains a single lepton within the jet. On the other hand, QCD production of a bb pair can be followed by the decay of one of the quarks through the same four-fermion operator ψ u d b, resulting in a forward production of quarks, DM and a sterile neutrino as shown in Fig. 5 (left). The prompt b decay signal can be searched for as discussed in Sec. 7.3 by triggering on the associately produced b quark. The long lived sterile neutrino can then decay through the charged or neutral current channels just discussed. At the partonic level, this corresponds to the process b →ūd ψ →ūdν R ν R →ūdν R qq l , orūdν R l + l − ν . (7.4) The kinematics of these decays would be primed for forward searches such as FASER. We leave a detailed analysis of all these interesting collider signals to future work. 20 Additional frameworks of low-scale baryogenesis through meson and baryon oscillations contain decays which are long lived on collider time scales and can produce three light quarks. For instance, in the set-up of [39], baryogenesis proceeds through the oscillation of charmed and beautiful baryons which creates an asymmetry in a Majorana fermion χ. χ can be produced at the LHC, and would be relatively long lived before decaying into light quarks χ → u d s. Note that this scenario does however not accommodate DM. 21 Constructing models which produce collider signals from long lived sterile neutrino decays has been a subject of much study [100] Parameter Description Achieve Baryogenesis [6] m ν R Majorana neutrino mass 1.2 − 2.7 GeV Sec. 2, see also [6] mν R sneutrino mass 1.2 − 2.7 GeV Sec. 2, see also [6] Br Br of B → νν + Baryon 4 × 10 −4 − 0.1 Fig. 2 Fig. 1, Sec. 3.1, Fig. 3

Summary and Outlook
In this work we have constructed a supersymmetric realization of the mechanism proposed in [6], where the baryon asymmetry and the DM abundance of the Universe are produced from neutral B meson oscillation and subsequent decay into a dark sector. We have introduced a model of R-SUSY in which an exact U (1) R symmetry is identified with U (1) B baryon number (such that particles within the same multiplet are charged differently under baryon number). Heavy squarks of mass 1 − 2 TeV are integrated out to generate baryon number conserving four-fermion operators coupling SM quarks to a light Dirac bino (about 2-4 GeV in mass) which carries baryon number. The existence of an exact R-symmetry forbids the existence of Majorana masses for gauginos, which are instead of Dirac type. This is a crucial feature of the model, as the Majorana gauginos and the absence of R-parity allow for neutron-antineutron oscillations and dinucleon decays, both of which processes are very tightly constrained experimentally. That said, it is well known that supergravity effects generically break the R-symmetry and generate Majorana gaugino masses (along with soft squark masses and a-terms). Taking this into account, we have quantified the degree to which the R-symmetry needs to be exact in the face of the experimental constraints.
The baryogenesis mechanism starts when b-quarks and antiquarks are produced from the late decay of an "inflaton-like" TeV scale particle at temperatures T R ∼ 10 − 50 MeV. The quarks then hadronize into neutral B 0 s,d mesons which oscillate and quickly decay into a baryon and a Dirac bino via the aforementioned operator. The Dirac bino further decays into dark sector particles, which are identified with a righthanded sterile neutrino and its baryon number carrying superpartner -a sterile sneutrino which is stable on cosmological time scales and and therefore constitutes the DM.
We have shown that the R-SUSY model presented here reproduces the measured baryon asymmetry and DM abundance of the Universe while being consistent with experimental constraints from colliders and flavor observables. In particular, this mechanisms motivates a study of a yet unexplored region of the R-SUSY model parameter space -that of a light Dirac bino. Indeed, this work represents the first excursion into the exploration of this slice of parameter space.
Focusing first on the flavor observables that are relevant for the baryogenesis dynamics, we have shown that the NP modifications to the semileptonic-leptonic asymmetries a s,d sl in the neutral B 0 s,d −B 0 s,d systems can be large enough to generate a sufficient baryon asymmetry while being allowed by current experimental constraints. In our model, the extra CP violation that enhances the semileptonic asymmetries can be traced back to complex phases in the RPV couplings λ ijk or in the off-diagonal elements of the squark mass matrix. Another robust prediction of the mechanism is the existence of a new decay channel of B mesons into a baryon and dark sector states, which needs to have a large branching fraction 10 −3 . In the supersymmetric setup, this is only possible if at least one of the squarks has a mass close to 1 TeV while the RPV couplings λ ij3 have a flavor structure that allows at least one of the combinations 11, 12, 21 or 22 to be O(1). The R-SUSY NP modifications to additional flavor observables (most notably K 0 −K 0 oscillations and the branching ratio of b → sγ and µ → eγ) further constrain couplings that, while not directly relevant for baryogenesis, give a hint towards the flavor structure of the model. The values of the parameters for which R-SUSY successfully realizes baryogenesis are summarized in Table. 5.
As a new feature of the supersymmetric model, the dark sector particle content may be elegantly realized as a sterile neutrino supermultiplet. The Dirac bino then decays predominantly into the sterile neutrino and the cosmologically stable sterile (scalar) sneutrino, which carries baryon number −1 and constitutes the DM. In this way, equal and opposite matter-antimatter asymmetries are generated in the visible and dark sectors without violating the global baryon number of the Universe. Dark sector interactions between the sneutrino and neutrino act to deplete the otherwise overproduced symmetric abundance of the sneutrino DM. These partially invisible decays of the bino (and, by extension, of b-flavored hadrons) can be looked for at B-factories like Babar and Belle II or at the LHC, by exploiting mono-jet searches and b-tagging techniques. In addition to that, the decays of the sterile neutrino into SM fermions may be long lived on collider time scales and can be searched for at experiments like SHiP, MATHUSLA, FASER, CODEX-b, and the ATLAS and CMS muon trackers. Additional constraints on model parameters are derived from measurements of neutrino oscillations and by exploiting astrophysical observations, most importantly the existence of stable neutron stars.
To conclude, we have presented a supersymmetric model that can achieve low-scale baryogenesis and DM production without baryon number violation. In contrast to the standard lore of high scale baryogenesis, our setup enjoys a plethora of possible experimental signatures ranging from measurements of flavor observables to long lived decays at colliders. Such exciting possibilities leave the door open to fully test this baryogenesis scenario in the not too distant future.

Appendices A Analytic Expressions for Flavor Violating Observables
Here we give the analytic expressions and some more details on the evaluation of the NP contributions to the flavor violating observables discussed in Sec. 5. We distinguish two kinds of contributions depending on whether the light bino plays a role in the process or not.

A.1 Contributions from RPV Couplings
The R−parity violating couplings λ ijk can mediate flavor changing transitions and induce CP-violating effects through their complex phases. A comprehensive study of RPV SUSY phenomenology was done in [62], to which we refer for any further details not discussed here. ∆F = 2 processes. Neutral meson oscillation parameters are modified by the existence of diagrams like the ones depicted in Fig. 7 for the case of the B 0 d system. The corresponding box diagrams were computed in [63] for the (s)top couplings, here we generalize the results including the two lighter generations. One important difference in our setup is the absence of left-right squark mass mixing, which is forbidden by the R-symmetry. Taking this into account, we find Here, V ij denotes the corresponding CKM matrix element and C rs refers to the mixing matrix of positive charginos. The loop integrals are given by .

(A.2)
In the limit where one of the masses M is much larger than the others, we have I 4 ∼ 1/M 2 and I 2 ∼ 1/(M µ) 2 , where µ is some combination of the remaining masses. Additionally, logarithmic factors with quotients of masses can appear in the expressions. Therefore, squark masses dominate the integrals and the coefficients C i are always suppressed by 1/m 2 q . ∆F = 1 processes. The RPV couplings alone can induce the flavor-violating decay b → sγ through the diagrams shown in Fig. 8. Summing all contributions, the resulting branching ratio is  where The loop function G 1 is defined in Eq. (A.12) and is evaluated at For heavy squarks, we have G 1 (x 1) 1/2.

A.2 Contributions from Squark Mass Mixing
Here we compute the contribution of a light Dirac bino to the meson oscillation parameters, following the spirit of [28]. The flavor violation is generated in this case by off diagonal elements in the squark soft mass mixing matrix. To be specific, we give the expressions for consider Comparing with [28], we see that binos give rise to fewer operators than gluinos due to the more restricted color structure of the diagrams. The loop function depends on the ratio λB = m 2 B /m 2 q and is given bỹ f 6 (λB) = 6λB(1 + λB) log λB − λ 3 B − 9λ 2 B + 9λ + 1 3(1 − λB) 5 , (A.8) which approachesf 6 → 1/3 when mB mq. To get this result we have neglected the momentum of the external quarks and applied Fierz rearrangements to obtain the desired color structures. In general, ∆M NP is a complex quantity because the off-diagonal elements of the squark mixing matrix can be complex, and thus so are δ LL and δ RR . Effects in the CP-violating observables defined in Sec. 3 are therefore expected. Note however that C 1 is always real because it only depends on the absolute value of the squark mixing parameters δ 2 = δ ·δ . Analogous formulas, with the obvious substitutions, hold for the B 0 s and K 0 systems. ∆F = 1 processes. The calculation of these processes are done following the conventions of [28]. The diagrams for b → sγ are shown in Fig. 10. Summing contributions from both diagrams, the decay width of a b quark to a s quark and a photon mediated by a light bino can be written as where B L1,R1 correspond to diagrams with an external chirality flip, while B L2,R2 come from the diagrams with internal bino-higgsino mixing, The loop functions G 1 and G 2 , which are to be evaluated at x = m 2 B /m 2 q , are given by In the case where the bino is much lighter than the squarks, we have G 1 (xB) 1/2 and G 2 (xB) 1.
The computation of the flavor-violating muon decay µ → eγ is very similar to the b → sγ transition, with the obvious substitutions. We therefore directly give the result δ LL m 2 G 1 (xB), (A.14) The loop functions as defined above should now be evaluated at x = m 2 B /m 2 .