Arxiv:2011.06115V3 [Hep-Ph] 11 Feb 2021 Niﬁcant Challenges

Making the Universe at 20 MeV

Gilly Elor1, ∗ and Robert McGehee2, 3, 4, † 1Department of Physics, University of Washington, Seattle, WA 98195, U.S.A. 2Leinweber Center for Theoretical Physics, Department of Physics, University of Michigan, Ann Arbor, MI 48109, USA 3Berkeley Center for Theoretical Physics, University of California, Berkeley, CA 94720, USA 4Theory Group, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA We present a testable mechanism of low-scale baryogenesis and dark matter production in which neither baryon nor lepton number are violated. Charged D mesons are produced out-of-equilibrium at tens of MeV temperatures. The D mesons quickly undergo CP-violating decays to charged pions, which then decay into dark-sector leptons without violating lepton number. To transfer this lepton asymmetry to the baryon asymmetry, the dark leptons scatter on additional dark-sector states charged under lepton and baryon number. Amusingly, this transfer proceeds without electroweak sphalerons, which are no longer active at such low scales. We present two example models which can achieve this transfer while remaining consistent with current limits. The required amount of CP violation in charged D meson decays, while currently allowed, will be probed by colliders. Additionally, the relevant decays of charged pions to dark-sector leptons have been constrained by the PIENU and PSI experiments and will be further explored in upcoming experiments.

I. INTRODUCTION teries of the Standard Model of particle physics (SM) and are discoverable in the near-future. The standard model of inflationary cosmology pre- While the mechanism of baryogenesis is necessary to dicts a Universe born with equal parts matter and anti- explain the origin of the complex visible structures we matter, necessitating a dynamical mechanism to gener- observe today, such structures only constitute roughly ate an asymmetry which seeds the complex structures 5% of the energy budget of the Universe. The SM does observed today. The required primordial baryon asym- not explain the nature and origin of dark matter (DM), metry of the Universe (BAU) is inferred to be the gravitationally inferred component of matter which makes up roughly 26% of the energy of the Universe obs −11 YB ≡ (nB − nB¯ )/s = (8.718 ± 0.004) × 10 , (1) [1, 2]. Experimental searches for DM at colliders and direct detection experiments, together with studies of the from measurements of the Cosmic Microwave Back- possible indirect effects of DM in astrophysical observa- ground (CMB) [1, 2] and light element abundances af- tions, have yet to shed light on its nature. ter Big Bang Nucleosynthesis (BBN) [3, 4]. Discovering Many particle physics models have been proposed to baryogenesis, the mechanism responsible for generating explain the nature and origin of DM. However, with this asymmetry, is therefore critical to understanding our the simplest scenarios becoming ever more constrained, very existence. richer dark or hidden sectors containing multiple parti- A mechanism of baryogenesis must satisfy the three cles with new interactions and symmetries become more Sakharov conditions [5]; C and CP Violation (CPV), interesting.2 Such dark sectors open up a host of new re- baryon number violation, and departure from thermal constructable cosmological histories [21–23], which may equilibrium. Many mechanisms of baryogenesis have be tested by colliders [24–26], direct detection and neu- been proposed, including the perennial favorites: elec- trino experiments [27–32], and indirect searches [33–35]. troweak baryogenesis [6–16] and leptogenesis [17]. But, Moreover, an interesting subset of those models also ex- concrete realizations of these mechanisms encounter sig- plain the BAU. For instance, in many models of Asym- arXiv:2011.06115v3 [hep-ph] 11 Feb 2021 nificant challenges. Electroweak baryogenesis models of- metric Dark Matter [36–39], DM carries a conserved ten predict electric dipole moments of electrons, neu- charge whose asymmetry is tied to the BAU in a uni- trons, and atoms which are ruled out by experiments [18]. fied framework that explains both asymmetries (e.g., [40], On the other hand, leptogenesis models typically occur and references therein). at high scales and involve very massive particles, thereby In this work, we explore a novel scenario where a dark- making experimental confirmation unlikely.1 Therefore, sector state is charged under lepton number. Assuming exploring novel baryogenesis mechanisms is well motivated, especially if they address other outstanding mys-

2 While a rich dark sector may seem less compelling due to its complexity, it is a well motivated scenario from a top-down perspec- ∗ [email protected] tive. Furthermore, the SM displays signiﬁcant richness, present † [email protected] authors included, despite its meager contribution to the energy 1 See [19] for an interesting proposal, as well as [20] and references budget of the Universe. It would not be too surprising to discover therein for testable seesaw scenarios which achieve baryogenesis. a dark sector with similar complexity. 2

FIG. 1. Summary of the mechanism by which a lepton asymmetry is produced from late-time production of charged D± mesons. Here we consider CP-violating decays of the D± mesons into final states involving an odd number of charged pions. The charged pions decay into dark- and visible-sector leptons without violating lepton number, producing equal and opposite visible- and dark-sector asymmetries. late-time production at temperatures of order 20 MeV, constructing models of low-scale baryogenesis can be a mesons which undergo CP-violating decays may then challenge and there are only a few working examples (see subsequently have their decay products quickly undergo e.g [45, 46]). Furthermore, recent proposals for solutions lepton-number-conserving decays into dark leptons. In to the gauge hierarchy problem such as N naturalness [47] this way, an equal and opposite lepton asymmetry is gen- and cosmological relaxation [48] require the BAU to be erated between the visible and dark sectors. In particu- generated at a low scale. lar, CP violation in charged D± meson decays followed by If one holds out hope that the requisite CP violation prompt decays of charged pions to light, MeV-GeV scale for baryogenesis exists in the SM, one is also inevitably (dark) leptons may be used to generate such an asymme- led to consider mechanisms at such low scales. It is often try. Intriguingly, this asymmetry is directly linked to SM claimed that there is not enough CP violation within the observables, making this mechanism testable at current SM alone to provide for the baryon asymmetry, regard- and upcoming experiments (see Fig. 1 for a summary). less of the baryogenesis mechanism. However, there are While a late-time production of a lepton asymmetry potentially abundant and untapped sources of CP viola- may be interesting in its own right, to explain the BAU, tion in QCD resonances: meson oscillations [49–51] and the lepton asymmetry must generate a baryon asymme- meson decays, as in this work. Thus, there’s a relatively try. We achieve this by minimally extending the dark sec- unexplored swath of theory space in which the SM alone tor to include low-scale, dark scattering processes which provides the necessary CP violation via mesons, allow- produce an equal and opposite baryon asymmetry in the ing for different realizations of Mesogenesis. The focus dark and visible sectors using the initial lepton asym- on D mesons in this work is a natural sequel to the B metry.3 The SM baryon asymmetry is Frozen-In via meson work started in [49–51] and predicts a unique set these dark-sector scatterings. In summary, we present of experimental observables. here a novel, testable, mechanism of low-scale baryogen- This paper is organized as follows. First in Sec. II, we esis and DM production utilizing SM D± meson decays introduce the mechanism. Next in Sec. III, we present the at late times, effectively making the Universe as we know details by which baryogenesis is achieved; we solve a set it at 20 MeV. In contrast with previous mechanisms such of Boltzmann equations for the lepton and baryon asym- as high-scale leptogenesis, this does not involve lepton- metry and demonstrate that the BAU can be achieved in or baryon-number violation and does not require Elec- light of known limits on the CP violation and branching troweak sphalerons. fractions of D± mesons. We also discuss the way in which One of the most remarkable features of this model is the correct DM relic abundance can be achieved. Next the ability to achieve baryogenesis, as well as the pro- in Sec. IV, we present two models and demonstrate that duction of DM, at such low temperatures. Reasonable they can accommodate a sizeable dark-sector scattering assumptions may lead one to conclude that a baryogen- to produce the BAU. We conclude with a discussion of esis mechanism, regardless of the source of CP violation, possible extensions, additional variations of Mesogenesis, must set the asymmetry by T & 38 MeV [44]. Thus, and other future directions in Sec. V. App. A contains a detailed derivation of the Boltzmann Equations. In App. B, we tabulate the relevant D± decay modes and the current limits on their branching fractions and CP 3 For other models which transfer an asymmetry from the dark asymmetries. sector to the SM to realize baryogenesis, see e.g. [41–43]. 3

II. THE MECHANISM We consider decays of charged pions into dark and SM leptons that proceed through an effective operator of the We now introduce the mechanism of baryogenesis and form DM from D± mesons. First, we assume the late decay of 1 h ¯ µ ih¯ i an inflaton-like scalar field Φ into quarks and anti-quarks O = dΓ u `dΓµ` + h.c. , (3) Λ2 when the temperature of the Universe was roughly tens of MeV. In particular, we assume the decay occurs at where ` is a SM charged lepton and Γµ represents all pos- temperatures in the range TBBN . TR . TQCD, so that sible distinct Lorentz tensors. The UV model from which the produced quarks hadronize but do not spoil the pre- the operator in Eq. (3) arises depends on the Lorentz dictions of BBN. TR is the “reheat temperature” corre- structure. For instance, a scalar operator could arise from sponding to the time at which Φ decays. Such a field a charged scalar mediator similar to [55], while a vector Φ may arise naturally out of multi-field inflation models, operator could arise from a new vector of a left-right sym- or may be identified as a flavon in the context of flavor metric model e.g. [28]. Depending on the UV model, to theories. While such models are interesting to consider, be consistent with current constraints, the scale Λ could for generality, we remain agnostic about the nature of Φ be anywhere from hundreds of GeV to a few TeV. and simply consider its mass and decay width (mΦ, ΓΦ), The result of the fast decays, as well as relevant branching ratios introduced later. + + + The produced quarks and anti-quarks hadronize into π → `d + ` , m`d < mπ − m` , (4) an equal number of mesons and anti-mesons. By ad- justing the mass and decay width of Φ, we consider sce- along with the conjugate decays, is the generation of a narios in which D± mesons (with mass of 1.87 GeV) lepton asymmetry in the dark sector are produced out of equilibrium. Thus, the tempera- n` − n¯ ture must be low enough so that D± decay before an- d `d Y`d ≡ , (5) nihilating with other species. The D-meson lifetime is s τ = 1.5 × 109 MeV−1 [4], while the typical cross section D which is equal and opposite to a lepton asymmetry cre- for hadrons is determined by the pion mass σ ' m−2 ∼ π ated in the visible sector. Throughout this work, we use O(10 mb). Following the argument in [49], we find an up- the common co-moving yield variables Y defined as the per bound on the reheat temperature such that the D± ratio of the number density to the entropy density in the mesons decay before annihilating: SM bath. In the absence of any other lepton-charged, dark-sector states, Y = Y dark, the total lepton asym- 3.5 MeV . TR . 20 MeV . (2) `d L metry in the dark sector. But, in later sections, we intro- The lower bound of 3.5 MeV comes from the requirement duce additional dark-sector leptons in order to generate that the asymmetry generation completes before SM neu- dark 4 the baryon asymmetry, resulting in Y`d ≤ YL . Re- trino decoupling and we restrict our reheat temperatures gardless, since we never introduce lepton-violating inter- to this range [52–54]. actions, the following is always true: The D± mesons then undergo CP-violating decays into an odd number of charged pions. Since these decays oc- dark SM YL = −YL . (6) cur out of equilibrium, an asymmetry in charged pions is temporarily generated. These charged pions themselves In this way, lepton asymmetries are generated in both the quickly decay into a lighter, dark-sector Dirac fermion dark and visible sectors while conserving the total lepton `d which carries visible sector lepton number (L = +1). number of the Universe.5 Since annihilations of pions are subdominant to their de- The generated lepton asymmetry is directly related to cays for the range of temperatures in Eq. (2), these fast SM observables, pion decays are able to happen before any appreciable washout of the temporary pion asymmetry. It may be dark `d X f f YL ∝ Brπ ACP BrD+ , (7) possible to reheat at slightly higher temperatures than f those in Eq. (2) at the cost of washing out some pion asymmetry as well as depleting some initial D meson `d + + where Brπ ≡ Br (π → `d + ` ), the sum is over final abundance. However, calculating a more accurate re- states f which contain an odd number of π±, and Af heat temperature upper bound would require solving the CP Boltzmann equations for both charged D’s and pions while accounting for non-negligible scattering with the rest of the SM bath. This complexity is beyond the scope 4 of this work and we instead take the conservative upper In much of the parameter space that results in the measured baryon asymmetry, the dark lepton asymmetry is much greater bound in Eq. (2). By introducing this new, dark-sector dark and Y`d ≈ YL even after baryogenesis completes. decay channel for pions, an asymmetry can start to form 5 This mechanism does not require lepton number violation. But between the dark and visible sectors. Without it, the the presence of lepton violation, for instance in neutrino masses, generated charged pion asymmetries would wash out. will not spoil this mechanism. 4 is the CP violation observable for a given decay mode, asymmetry. Critically, this dark scattering can occur defined by through an operator which conserves the total baryon and lepton number of the Universe; a dark-sector lepton + − ¯ f Γ(D → f) − Γ(D → f) asymmetry is partially transferred to equal and opposite A = . (8) CP Γ(D+ → f) + Γ(D− → f¯) dark- and visible-sector baryon asymmetries. Schemati- cally, we consider scatterings of the form Brf ≡ Br (D+ → f) is the branching fraction of the D+ D+ `¯ + χ → χ + B , (11) decay (the relevant decay modes and the current limits d 1 2 on their branching fractions and CPV are summarized in where B is a SM baryon, and χ1 and χ2 are the gauge- `d Table. III). The current limits on Brπ may be extracted singlet, dark-sector states which may be fermions or from the limits on sterile neutrinos. Recasting and im- scalars depending on the exact dark-sector model. For posing current limits for charged pion decays into elec- possible baryon and lepton number charge assignments, trons [56, 57], we find that the allowed branching ratio see Table I. Note that the mass of a dark-sector state is not large enough to generate the requisite asymme- charged under baryon number must be greater than 1.2 try when m`d > 1 MeV. However, the branching ratio is GeV [61], but dark leptons may be considerably lighter. unconstrained for sub-MeV `d masses so that this decay Additional kinematic and stability requirements will be mode can generate the entire asymmetry. Recasting the model dependent, and we leave these details for Sec. IV. most current bound from PIENU [57, 58] for final-state Depending on the details of the dark-sector charge as- muons yields signment and the UV model, either χ1 or χ2 (or both) ± ± −6 −5 may constitute (part of) DM. A Z2 discrete symmetry Br(π → µ + MET) . 10 − 10 , will generically need to be imposed to stabilize the DM

for 15.7 MeV < m`d < 33.8 MeV , (9) and evade washing out the produced asymmetry. In Sec. IV, we describe the cosmological assumptions and which is just at the threshold of producing enough asym- possible models of the dark sector that allow for a large- metry. For lighter `d masses, constraints can be recast enough cross section to transfer the asymmetry consis- from PSI [57, 59] tent with current bounds as well as produce the measured

± ± −3 DM relic abundance. Br(π → µ + MET) . 10 ,

for 5 MeV < m`d < 15 MeV . (10) III. THE DETAILS Note that for ∼1-5 MeV, the bound on the branching −2 fraction can be as weak as 10 . Given the `d mass de- Having given a broad-brush overview of the important pendence, these bounds do not constrain the entire pa- ingredients of this mechanism in the previous section, we rameter space of interest to us; as with decays to final- move on to calculate the relevant matter contents in de- state electrons, sub-MeV `d masses lead to completely tail. We consider the generation of the (dark-sector) lep- unconstrained branching ratios. ton asymmetry, (visible-sector) baryon asymmetry, and Improved measurements of these decays will be the DM in turn. focus of upcoming searches at future experiments and as such will be able to further probe this mechanism [60]. In what follows, we will demonstrate that a large lepton A. Generating a Lepton Asymmetry asymmetry may be generated which is consistent with current experimental bounds and may be probed in the In this section, we demonstrate that a dark lepton future. 6 asymmetry equal to (or much greater than) the measured Baryogenesis is achieved by transferring the dark lep- baryon asymmetry may be generated via the processes ton asymmetry into a SM baryon asymmetry using ad- outlined in Fig. 1, postponing a discussion of how it may ditional dark-sector states and dynamics which can be be transferred to a SM baryon asymmetry to Sec. III B. rich and possibly reconstructable. In particular, we con- In order to numerically solve for the generated lepton sider `d interactions with additional dark-sector states asymmetry, we consider the coupled Boltzmann equa- (χ1 and χ2) that carry lepton- and baryon-number which tions which track the production and CP-violating de- can transfer the dark lepton asymmetry into a SM baryon cays of D± mesons into π±, which then subsequently decay into dark leptons and anti-leptons. For simplicity, we compute the generated lepton asymmetry for the range of reheat temperatures in Eq. (2) so that annihilations 6 For simplicity, we refer to this as an asymmetry transfer, since ± ± ¯ of D and π mesons can be ignored. The reheat tem- the asymmetry in `d − `d is being partially translated into an asymmetry in SM baryons and dark-sector particles. Note that perature is defined by 4H (TR) = ΓΦ, so that Eq. (2) the total lepton asymmetry in the dark (and SM) sectors does corresponds to an inflaton decay width in the range −22 −21 not change as a result, and so this is not a “transfer,” strictly ΓΦ ∈ 1 × 10 GeV, 3 × 10 GeV . Additionally, as speaking. the inflaton must be heavy enough to produce D±, its 5 mass must be in the range mΦ ∈ [5 GeV , 100 GeV]. Φ late decays to radiation so that the evolution of the Φ number density and the radiation density are governed by the interplay of the following Boltzmann equations dn Φ + 3Hn = − Γ n , (12) dt Φ Φ Φ dρ rad + 4Hρ = + Γ m n , (13) dt rad Φ Φ Φ where the Hubble parameter is given by

2 8π H = 2 (ρrad + mΦnΦ) . (14) 3MPl We assume that Φ was in equilibrium at some high temperature with the bath and as such has a number density ∝ T 3. While it may be possible to achieve this mechanism in an inflationary model where Φ is identified as the inflaton, this assumption of high-temperature equilibrium simplifies this analysis at the cost of presuming other scalars responsible for inflation. Since the focus in this section is on the lepton asymme- ¯ FIG. 2. The shaded green region corresponds to values of try, we assume a minimal dark sector with only `d and `d and do not include any additional dark-sector states or charged D meson and pion observable needed to produce a lepton asymmetry that is equal to or greater than the ob- interactions, deferring this discussion to Sec. III B. Since obs ± served baryon asymmetry YB . The contours correspond to the formation and subsequent decay of the D meson solving the Boltzmann equations Eq. (15), as summarized in and the following decay of the π± meson occurs quickly Eq. (16), with values of (mΦ,TR) = (5 GeV, 20 MeV) which (before any scattering effects can significantly change the correspond to maximizing the produced asymmetry. The dot- abundance of these mesons), the generated dark-sector ted gray line and the gray shaded region represents the current lepton asymmetry can be written simply as (for a de- limits on the pion branching fraction [56] and the limits on the tailed derivation, see App. A) sum of the D meson CP asymmetry and branching fraction, respectively. d n` − n¯ + 3H n` − n¯ = (15) dt d `d d `d

D `d X f f f sizable lepton asymmetry may be generated. We find 2 ΓΦ nΦBrπ Nπ aCP BrD+ , f f f dark `d P f Y Br N a Br + T 10 GeV L ' π f π CP D R . (16) f + obs −3 −5 where Nπ is the number of π minus the num- YB 10 3 × 10 20 MeV mΦ ber of π− in each channel labeled by f. Note that only decay modes with an odd number of If all the lepton asymmetry can be instantaneously con- charged pions contribute, as expected. Here we de- verted into a baryon asymmetry, then the SM baryon SM dark D asymmetry will be YB = YL . In practice, the fine ΓΦ ≡ ΓΦBr(Φ → c)Br(c → D) (where we account for the possibility that Φ can also populate dark-sector dark-sector dynamics need not transfer the asymme- f f f f try completely. Therefore, Eq. (16) represents a lower states). Also, aCP ≡ ACP /(1 + ACP ) ≈ ACP for most f bound on the observables such that baryogenesis can be decay channels since ACP is a small number. The sum dark obs achieved. In Fig. 2, we show contours of YL /YB is over the exclusive rates to each of the final states f for a range of values of the experimental observables listed in Table III. In this way, an asymmetry in ` is d Br`d and P N f af Brf . Also shown for reference generated, as defined in Eq. (5), that is equal and op- π f π CP D+ posite to an asymmetry generated in the visible-sector is the PSI constraint from Eq. (10) which holds when 5 MeV . m`d . 15 MeV. leptons. This asymmetry is interestingly related to ob- ± servable CP-asymmetries and branching fractions in SM Summing over the relevant D decay modes in Ta- mesons systems. Critically, note again that the total lep- ble III, we find ton number of the Universe is actually conserved, as we X f f +0.0031 N f a Br = −9.3 × 10−4 , (17) have not introduced any lepton-number-violating inter- π CP D+ −0.0039 actions. f We numerically integrate the above set of Boltzmann where the central value corresponds to taking the central P f f f `d f f equations and float the values of f Nπ aCP BrD+ , Brπ , values of both ACP and BrD+ for each decay channel. TR, and mΦ to discover the parameter space in which a The lower bound corresponds to the “lowest-reasonable” 6 value for the sum and is calculated in the following way. f To make the sum as negative as possible, we take all ACP values 1σ below their mean. For channels with values of f ACP which are still positive, we assume their correspond- f ing BrD+ is 1σ below the mean. For channels which in- f stead (now) have negative ACP , we assume their corre- f sponding BrD+ is 1σ above the mean. The upper bound in Eq. (17) is calculated in an analgous way. The measured central value is shown in dashed black in Fig. 2, while the solid gray region corresponds to the (absolute value) of the most negative possible sum. Comparing Eq. (17) to Eq. (16), it is clear that it is possible to generate a dark-sector lepton asymmetry that is orders of magnitude larger than the measured baryon asymmetry. f f Future, more precise measurements of ACP and BrD+ for the various pion decay channels in Table III will shift the gray, ruled-out region to the left. Such improvements are expected to be made by experiments such as LHCb. f While ACP are expected to be small in the SM, quanti- fying them is plagued with the usual technical challenges of the charm sector. If better SM predictions result in f ACP which are smaller than we require for this mecha- f nism, new physics contributions could also enhance ACP while keeping them within current experimental bounds. The PIENU experiment has accessed the majority of its data, and as such, an improvement in the sensitiv- FIG. 3. Numerical solutions for the yields of relevant `d ity of Brπ is unlikely. However, the relevant mass range species for a benchmark point which produces the observed could be extended as uncertainties are improved which baryon asymmetry. We take TR = 10 MeV, mΦ = 6 GeV, −15 −2 `d previously made certain areas of phase space difficult to Br (Φ → χ1χ¯1) = 0.1, hσvi = 1 × 10 GeV , and Brπ = probe. Additionally, next-generation experiments which −3 P f f f 10 , and for f Nπ aCP BrD+ we take the maximum value would improve the limit on the branching fraction are be- in Eq. (17) (although note that saturating this bound was ing proposed [60]. Note that for a given UV model gener- not required to get the right asymmetry). There are three `d distinct phases which are delineated by vertical dashed, gray ating Eq. (3), the branching ratio Brπ can be computed and will depend on the scale of the higher-dimensional lines, and are highlighted particularly in the top panel. See operator Λ. This in turn will be constrained by collider the text for details. and astrophysical searches in a model-dependent way. As a sanity check, we have performed this computation for the charged-scalar mediator scenario e.g. [55], and find evolves according to that direct constraints on the scale Λ in this model do dnχ1 not exclude any of the parameter space of Fig. 2. + 3Hnχ = ΓΦnΦBr (Φ → χ1χ¯1) (18) dt 1 − hσvi n¯ n , `d χ1 B. Generating a Baryon Asymmetry where hσvi is the thermally averaged cross section7 of the baryon-transfer process in Eq. (11). Detailed derivations We now complete baryogenesis by elucidating the de- of all of the Boltzmann equations in this section may be tails by which equal and opposite baryon asymmetries in found in App. A. the dark and visible sectors are frozen-in. We remain ag- Recall that the evolution of the asymmetry in `d in nostic about the dark-sector model which generates the Eq. (15) simply tracked the production of a lepton asym- scattering process in Eq. (11), deferring a detailed dis- metry. We modify this equation to include the relevant cussion to Sec. IV. Instead, we compute how large the cross section must be for the process in Eq. (11) to efficiently transfer the dark lepton asymmetry to the measured baryon asymmetry of the Universe. 7 While it is technically correct that this is a thermally averaged For simplicity, we take Φ to also decay to χ1 and there- cross section, the phase space distribution functions will not be fore require mχ1 < mΦ, though another scalar could in- the usual thermal Maxwell-Boltzmann distributions. Rather, stead be responsible for this late-time, out-of-equilibrium they are determined by the kinematics of the relevant decays and Hubble expansion. χ1 production. The number density of χ1 therefore 7 scattering term and obtain the evolution equation for the C. Generating the Dark Matter asymmetry in `d: d Since baryon number is never violated, the measured n` − n¯ + 3H n` − n¯ = (19) dt d `d d `d SM baryon asymmetry is always balanced by an equal D `d X f f f and opposite baryon asymmetry in the dark sector. This 2Γ nΦBr N a Br + − hσvi nχ n` − n¯ . Φ π π CP D 1 d `d dark baryon asymmetry, therefore, is always an asym- f metric component of DM, and a substantial fraction at For simplicity, we take nχ1 ∼ nχ¯1 here as both are ini- that due to the lower bound on baryon-charged masses tially produced in equal amounts from Φ decays. The of 1.2 GeV. Further, there are equal and opposite lep- Boltzmann equation for the evolution of the SM baryon ton asymmetries in the dark and visible sectors. If the asymmetry is then given by state(s) which comprise the dark-sector baryons due not d also account for this dark-sector lepton asymmetry, then (n − n ) + 3H (n − n ) = (20) dt B B B B these additional dark leptons must make up a different, asymmetric subcomponent of DM. Clearly, the details − hσvi nχ n` − n¯ . 1 d `d depend on baryon- and lepton- number assignments to Next, we turn our attention to numerically solving the set states χ1 and χ2, which we defer to the next section. of coupled Boltzmann equations for the baryon asymme- Here, we just make simple qualitative remarks about try of the Universe. how generating the correct DM abundance is relatively From Eq. (20), it is clear that the produced lepton straightforward (as compared to generating the baryon asymmetry will be efficiently transferred when the scat- asymmetry). tering rate nχ1 hσvi dominates over the Hubble expan- Perhaps the simplest scenario is to assume that the sion. Since we want the transfer to happen quickly, we dark baryon-charged state comprises almost the entirety require the cross section satisfy of DM, making it a well motivated case of completely n hσvi Y obs asymmetric DM. In this case, the lightest dark lep- χ1 B . (21) & dark ton is appreciably lighter than the dark baryon so that H(T ) T =TR Y L it makes up a negligible subcomponent of DM. Thus, The number density of χ1 is found by integrating Eq. (18) mDM ∼ 5 GeV. 3 and is roughly given by nχ1 ∼ Br (Φ → χ1χ¯1) T . As an alternative, we may also consider the case where We numerically integrate the Boltzmann equations the dark baryon state is lighter so that the other dark- Eqs. (12), (18), (19), (20) to solve for lepton and baryon sector particles must comprise the remaining relic abun- asymmetries, floating the model parameters. We find dance of DM. These could be new additional states, or that the dark scattering cross section is required to be just the dark-sector states already present to provide for greater than the baryon-asymmetry transfer. The details here become Y obs less relevant to the baryogenesis mechanism considered hσvi 10−16 GeV−2 B since there are generic, dark-sector freeze-out possibili- & Y dark L ties with viable parameter space. 10 GeV 20 MeV 10−1 × . (22) Both of the above scenarios make one important as- mΦ TR Br(Φ → χ1χ¯1) sumption which we require to be true generically. There In Fig. 3, we plot the solution of the Boltzmann equa- must exist a portal between the dark and visible sections for a benchmark point that achieves baryogenesis. tors that becomes efficient at late times (before SM neutrino decoupling) to allow the symmetric component(s) We plot the yields corresponding to Φ abundance, `d asymmetry, and SM baryon asymmetry. There are three of dark-sector states to sufficiently annihilate away, pre- clearly distinct regions in the plot, particularly empha- venting any overabundance of DM or non-negligible contribution to the relativistic degrees of freedom at BBN. sized by the top panel which zooms in on Y`d . First, a dark lepton asymmetry is produced as the inflaton begins The kinetic mixing portal involving a massive dark pho- ton fits our needs here and is a commonly used portal to decay. Soon after, `d-χ1 scatterings begin to dominate to transfer entropy out of dark sectors at late times (see and the `d asymmetry converts to a baryon asymmetry. As inflaton decay completes, the baryon asymmetry e.g. [40, 62, 63] for various such usages). One additional concern is that protons and anti- freezes in and the `d asymmetry is fixed. In the following section, we will present possible UV protons may not be able to efficiently annihilate below 8 models that can accommodate a cross section of the size 20 MeV , so that a large, symmetric baryon component in Eq. (22) while remaining consistent with present con- may freeze-in. One solution is to introduce other dark- straints. Note that the branching fraction of the inflaton sector processes which can efficiently deplete the symmetric component of ` and `¯ . As a proof of principle, into the dark sector Br (Φ → χ1χ¯1) depends on the spe- d d cific inflationary model, but can be sizable. Furthermore, hσvi nB and hσvi nχ2 are easily both less than Hubble as these particles are not sourced by Φ. Therefore, any possible washout effects are negligible. 8 We thank Seyda Ipek for pointing this out. 8

Field L B Field L B Field Spin L B Z2 Mass

χ1 1 0 χ1 1 1 φc 0 0 −2/3 +1 & 1 TeV χ2 0 -1 χ2 0 0

`d 1/2 1 0 +1 O(10 − 140 MeV) χ1 0 1 χ1 0 0

χ2 1 0 χ2 -1 -1 ψB 1/2 0 −1 +1 & 1.2 GeV

TABLE I. Possible baryon and lepton charge assignments for TABLE II. Dark-sector states which interact directly with the dark-sector states χ1 and χ2. Any baron-charged state must SM. ψB is a dark-sector baryon introduced in this section to be heavier than 1.2 GeV. (Left) The two states dark χ2 and generate interactions between the dark sector and SM baryons χ1 involved in the scattering Eq.(11) may be charged under through Eq. (24). SM baryon or lepton number. In this case, DM is multi- component with contributions from χ2 and χ1.(Right) One of the states involved in the scattering Eq.(11) may be a lep- χ2. Depending on the charge assignment, either χ1, χ2 tobaryon while the other uncharged under SM lepton and or both will constitute (at least part of) DM. baryon number. In this case, the leptobaryon can be the sin- gle component of DM. Regardless of the charge assignment of χ1 and χ2, a coupling between a dark-sector baryon and SM fields must be generated. To do so, we simply invoke the model one can imagine the extreme limit where such processes from [49] and introduce the following interactions which are allowed by all the symmetries:10 deplete all `d and leave only the tiny necessary asymmetric amount of `¯ . These would then only freeze-in d L ⊃ −y φ∗u¯ dc − y φ ψ¯ dc + h.c. . (23) SM baryons and not anti-baryons, avoiding the problem uidj c i j ψdk c B k altogether. Here, ψB is a dark-sector Dirac fermion carrying baryon number B = −1, and φc is a colored scalar mediator with baryon number B = −2/3. Integrating out the heavy φc IV. THE MODELS mediator leads to the following effective baryon-number conserving four fermion interaction Thus far, we have remained agnostic about the na- y c ¯c ture of the dark-sector fields participating in the baryon- Leff = u¯ djd ψB , (24) M 2 i k generating process in Eq. (11). The baryon asymmetry φc has been computed independently of the details of the where we have defined y ≡ yuidj yψB dk . Note that the col- dark-sector model and we have found that the dark lep- ored mediator mass and couplings are constrained to be ton asymmetry can be efficiently transferred to a SM √ Mφc y & 1 TeV to be consistent with collider bounds baryon asymmetry provided the dark-sector cross sec- (for details, see [51] and references therein). At low tion is sizable - as given by Eq. (22). We now turn our scales, this generates an effective mass mixing between attention to the details of the dark sector and possible SM baryons and the dark-sector baryon. Note that since viable toy models. ψB may couple to protons and neutrons through the op- There are several minimal variations of the dark-sector erator in Eq. (24), the stability of baryonic matter must field content which will suffice. The states χ1 and χ2 be ensured kinematically by mψB > 1.2 GeV. The field may be identified with a dark baryon and dark lepton as content is given in Table. II. summarized on the left side of Table I. Alternatively, one Additional dark-sector states are necessary to transfer state can be neutral and the other a dark leptobaryon9, the asymmetry; the two states χ1 and χ2 as well as an- charged under both SM lepton and baryon number as other mediator. These states are odd under a discrete summarized on the right side of Table I. Note that a Z2 (while `d and ψB, which interact directly with the SM, dark-sector state carrying baryon number must have a must be even), thereby evading washout and ensuring the mass greater than 1.2 GeV in order to be consistent with stability of the dark-sector lepton and baryon asymme- the observation of old neutron stars [61], while a dark tries. In this way, the dark lepton ` may scatter lepton may be significantly lighter. For concreteness, we d ¯ ¯ enumerate two dark-sector models corresponding to two `d + χ1 → χ2 + ψB , (25) lepton- and baryon-number charge assignments of χ1 and

10 Such models can arise in, for instance, supersymmertic theories 9 We thank Ann Nelson for this suggestion. [50]. 9 √ with charges of χ1 and χ2 chosen such that this process cross section. The color mediator mass Mφ y has been ¯ c conserves baryon, and lepton number. ψB subsequently set to saturate the collider bound of order 1 TeV. Note mixes into a SM baryon through Eq. (24). Which fields that Eq. (28) holds when mΦ & mχ2 + mB and mΦ & make up the DM depend upon further details of the dark- mχ1 ; as the inflaton populates both `d and χ1 in this sector model which we now explore. setup, and the energy available in the scattering will be of order mΦ. The DM will consist of the asymmetric parts of both χ1 Model 1: DM as Scalar Baryons and Leptons and χ2. We will need additional dark-sector interactions to annihilate away any symmetric part of the DM, as it In this model, we take χ1 to have (L = 1,B = 0) while will generically be overproduced, and to obtain the cor- χ2 has (L = 0,B = −1), corresponding to assignments in rect relic abundance. Such a set-up is simple to achieve the top left two row of Table I. We take both χ1,2 to be and there exists a host of dark-sector production mecha- scalars and introduce a Dirac fermion mediator ξ, which nisms that can deplete the asymmetry (for instance, see can be MeV scale. We make all three odd under a Z2. discussion in [49]). The DM consists of both (asymmetric components of) χ2 and χ1, stabilized under the discrete Z2 symmetry, as to maintain the equal and opposite asymmetries in the dark Model 2: DM as Fermionic Leptobaryons and visible sectors. As discussed above, χ1 is produced through Φ decays In this second example model, we take χ to have ± 1 and `d, generated from the π decay through Eq. (3), (L = 0,B = 0) while χ2 has (L = −1,B = −1), corre- then scatters off χ1 producing χ2 and the dark-sector sponding to charge assignments in the bottom right two fermion ψB as in Eq. (25). This scattering is mediated rows of Table I. We take both χ1,2 to be fermions and by the dark Dirac fermion ξ and is generated through the introduce a scalar mediator ΦL with L = 1 which may following baryon and lepton number conserving Yukawa be light. As in the first toy model, we take all three odd interactions allowed by the 2 symmetry: Z under a Z2. The χ2 are leptobaryons in this model and ¯ ¯ could be, for instance, a neutrino multiplet in a super- L ⊃ yb ψBξχ2 + yl `dξχ1 + h.c. . (26) symmetric model with an exact R-symmetry identified ξ (which can be relatively light) mediates the s-channel with baryon number, similar to [50] (e.g. a right handed ¯ process in Eq. (25). An intermediate ψB is produced sterile neutrio multiplet had a νR with B = 0, L = 0 which then quickly mixes into a SM baryon through andν ˜R with B = −1, L = 1). As per Model 1, we Eq. (24). This scattering transfers the dark lepton asym- have a dark fermionic baryon field ψB coupling to the metry to the SM baryon asymmetry. SM by the same UV construction as Eq. (24). The fol- Stability of baryonic matter is ensured kinematically lowing ∆B = 0 = ∆L Lagrangian is allowed by all of the by mχ2 > 1.2 GeV, as with any dark-sector state charged symmetries under baryon number. Meanwhile χ can have a sub-GeV 1 ¯ ¯ mass, unless otherwise restricted by kinematics. In gen- L ⊃ yb ψBχ2ΦL + yl `dχ1ΦL + h.c. . (29) eral the Yukawa coupling could induce decays of ψ into B This generates the scattering Eq. (25) mediated by the the dark sector. Since ψB transforms into SM baryons via the operator in Eq. (24), we require dark scalar lepton ΦL. As with the first model, the scattering cross section in this setup is easily large enough to

mχ2 + mξ > mψB > mB , (27) accommodate baryogenesis. In this model, DM is always (partially) comprised of which also ensures the stability of SM baryons. asymmetric χ2. If m` < mχ , it will also have an asym- We have computed the thermally averaged cross sec- d 2 metric component of `d. As in the previous model, we do tion corresponding to Eq. (11) for the s-channel scatter- not illustrate explicitly how to obtain the remainder of ing in this model and have confirmed that it can easily the relic abundance, nor do we detail how to remove sym- be sizable enough to satisfy Eq. (22). We leave a thor- metric components sufficiently. The possibilities here are ough exploration of the corresponding parameter space quite generic to dark-sector DM production and are thus to future work (and a more detailed UV embedding) and decoupled from the baryogenesis mechanism at hand. simply present, as a proof of principle, the following result: −15 −2 2 V. DISCUSSION hσvi ' 10 GeV (yl yb) (28) 10 MeV 20 GeV 10 GeV × . m m m We have introduced a novel, low-scale mechanism for `d χ1 χ2 generating the BAU from the late-time production of ± ± Here for simplicity we have taken mψB ∼ 5 GeV. We have mesons. D mesons decay in a CP-violating way to π , also fixed the dark-sector mediator mass to be mξ = 10 which in turn decay into a dark-sector state charged un- MeV— a heavier mediator will result in a slightly smaller der lepton number. The processes are out-of-equilibrium, 10 occur at tens of MeV, and generate equal and oppo- Eq. (3), thereby generating a sizable lepton asymmetry of site dark- and visible-sector lepton asymmetries. Addi- 0 s,d s,d order YL ∝ Br B → Mesons + X × ASL, where ASL is tional dark-sector states charged under SM baryon num- 0 the semi-leptonic asymmetry in the Bs,d systems. Simi- ber scatter with the dark leptons to transfer this dark larly, one could use the CP violation in neutral D0 meson lepton asymmetry into the observed SM baryon asym- oscillations. Although measurements by LHCb [66] have metry. Since we never explicitly violate lepton or baryon shown there is little evidence of sizable CP violation in number, these dark states also (partially) comprise DM. D0−D¯ 0 mixing, there exist some possible caveats [67, 68]. The matter-antimatter asymmetry is related to exper- We leave explaining the BAU through new iterations of imental observables. The measurement of the CP viola- Mesogenesis to future work. ± f tion in charged D decays, ACP , will be improved upon Generating a lepton asymmetry at low scales can be at for instance LHCb. The branching fraction of charged interesting in its own right, e.g. to resolve astrophysi- pions into charged leptons and missing energy, which has cal anomalies [69]. However, applications usually require been constrained by the PIENU experiment, will be fur- a larger lepton asymmetry than can be achieved with ther probed at future experiments. D± decays alone. To make progress in such a direc- We have presented two simple, dark-sector models tion, it may be interesting to consider the cumulative which can efficiently transfer the dark lepton asymmetry effect of multiple CP-violating meson systems (as dis- and achieve baryogenesis. However, we have remained cussed above). We leave such an avenue to future work. agnostic about many of the details of the dark sector It is a well motivated, albeit difficult, exercise for the- and have not UV completed these models. Such comple- orists to make the Universe at 20 MeV. That the finger- tions are the subject of future work and will likely open prints of such an extreme creation may hide within the up additional model-dependent, complimentary probes. SM itself compels us to discover whether Nature did it Another option is to introduce a complex, dark-sector too. gauge group which transfers the asymmetry through a dark Sphaleron process. We leave the details of this in- triguing possibility to future work. ACKNOWLEDGMENTS The study of this mechanism within the context of spe- cific flavor and inflationary models is left to future work. We thank Doug Bryman, Tim Cohen, Jeff Dror, We have focused here on explaining the BAU through Miguel Escudero, David McKeen, Seyda Ipek, Ann Nel- D± decays to π±. Depending on the flavor structure, one son, Jan Stube, Xabier Cid Vidal, Mark Williams for could also consider D± decays to Kaons which then decay useful conversations. We thank Doug Bryman and En- to dark-sector leptons. Current limits allow for a sizable rique Fernandez Martinez for pointing out an error in branching fraction of charged Kaons to muons given by the recast PIENU bound. We thank Miguel Escudero Br(K± → µ± + X) < 10−3 − 5 × 10−6 [64, 65]. Thus, for useful conversations regarding the numerics and re- we could repeat the calculations of this work using the casting sterile neutrino bounds. We thank Tim Cohen currently allowed values for CP-violating decay channels and Seth Koren for comments on the draft. GE is thank- of D± which involve an odd number of Kaons. A benefit ful to Ann Nelson for always encouraging the pursuit of using Kaons is the possibility of reheating at temper- of creative ideas. GE is supported by the U.S. Depart- atures above 20 MeV since washout-inducing Kaon an- ment of Energy, under grant number DE-SC0011637. nihilations stop at higher temperatures than their pion RM was supported in part by NSF grant PHY-1915314 counterparts. and the U.S. DoE Contract DE-AC02-05CH11231, and Φ can also decay into neutral mesons, such as B0, is currently supported by the U.S. DoE under grant de- which undergo CP-violating oscillations (this was lever- sc0007859. GE thanks the Berkeley Center for Theoreti- aged in the mechanism of [49] to explain the BAU). The cal Physics and Lawrence Berkeley National Laboratory produced B0 could decay into a dark lepton through for their hospitality during the completion of this work.

APPENDICES

Appendix A: Boltzmann Equations for the Lepton and Baryon Asymmetries

Here we present a detailed derivation of the Boltzmann equations presented in Sec. III. It is useful to solve the set of Boltzmann equations in terms of temperature T rather than time. In order to do so, one can write [70–72]

4 30 dT −4Hg?,ST + π2 ΓΦmΦnΦ = 3 , (A1) dt T (4g? + T dg?/dT ) which follows from energy conservation and is valid for energies above the neutrino decoupling temperature T > 3 MeV, as neutrinos and other light degrees of freedom are still coupled to the plasma. The number of relativistic 11 species g? (g?,S) contributing to the energy (entropy) density is given in [73]. The number density Boltzmann equations for charged D mesons are

dnD+ D+ dnD− D + 3Hn + = Γ n − Γ + n + + 3Hn − = Γ n − Γ + n − , (A2) dt D Φ Φ D D dt D Φ Φ D D D+ + where ΓΦ is the rate of decay of Φ to a ﬁnal state with 1 D (we assume no ﬁnal states with multiple) and ΓD+ is the total D± decay rate. We assume that all decays above are much faster than SM annihilations, which is true for D+ D− D TR < 20 MeV. Note that ΓΦ = ΓΦ ≡ ΓΦ since there is no CP violation in Φ decays. The number density Boltzmann equations for the charged pions are:

dnπ+ + + + + + 3Hn + = n + Γ + Br D → π + other + 2Γ + Br D → 2π + other dt π D D D − + + ΓD− Br D → π + other nD− − Γπ+ nπ+ , (A3)

dnπ− − − − − + 3Hn − = n − Γ − Br D → π + other + 2Γ − Br D → 2π + other dt π D D D + − + ΓD+ Br D → π + other nD+ − Γπ+ nπ− , (A4) ± where Γπ+ is the total π decay rate. All SM annihilation terms that would appear are negligible relative to the decay terms thanks to the low reheating temperature. The number of different decay terms are due to grouping by the number of final state charged pions, where we are only interested in decays with up to two pions of the same charge (see Tab. III). We’re also neglecting Φ decays into lighter quarks which could hadronize into pions. For intuition and simplicity, we first solve for the lepton asymmetry generated in the dark sector in the case that the additional dark-sector interactions that give rise to Eq. (11) are absent. Then, the number density Boltzmann equations for the dark leptons are

dn dn¯ `d + + `d − − ¯ + 3Hn` = Γ + Br π → ` `d n + , + 3Hn¯ = Γ + Br π → ` `d n − . (A5) dt d π π dt `d π π + + We assume any possibly present dark-sector annihilations are slow, `d is stable, and back scatters of ` `d → π are slow. ¯ To ﬁnd the generated lepton asymmetry in the dark sector, we simply take the diﬀerence of the above `d/`d Boltzmann equations above:

d + + n` − n¯ + 3H n` − n¯ = Γ + Br π → ` `d (n + − n − ) . (A6) dt d `d d `d π π π To simplify things analytically, we assume the Φ’s, produced D’s, and the subsequently-produced π’s decay quickly, which is approximately true as long as H ΓΦ ΓD Γπ. The right hand side of the above Boltzmann equation becomes + + D + + Γπ+ Br π → ` `d (nπ+ − nπ− ) = ΓΦ nΦBr π → ` `d (A7)  X + + − − X + + − − ×  Br D → 1π + f − Br D → π + f¯ + 2 Br D → 2π + f − Br D → 2π + f¯ f f  X + − − + − Br D → π + f − Br D → π + f¯  . f In the above, the sums over final states f do not include any additional charged pions. Rather, the number of charged pions we’re considering in a given sum is explicitly highlighted in the decay channel branching ratio. After some partial cancellations between channels with both signs of charged pions, we find that each channel is multiplied by its + − f net ‘+’ charge in π s vs π s, which we define as Nπ , and as expected: d D `d X f + − ¯ n` − n¯ + 3H n` − n¯ = Γ nΦBr N Br D → f − Br D → f . (A8) dt d `d d `d Φ π π f

f To simplify this further, we rewrite the branching ratio diﬀerences in terms of the observable ACP (as relevant to the D± decays): f + − ¯ + 2ACP f f Br D → f − Br D → f = Br D → f f ≡ 2BrD+ aCP . 1 + ACP 12

Thus, the generated dark-sector lepton asymmetry simpliﬁes to d D `d X f f f n` − n¯ + 3H n` − n¯ = 2Γ nΦBr N a Br + . (A9) dt d `d d `d Φ π π CP D f

Next, we modify the story by allowing the interactions in Eq. (11) to cause a net dark baryon asymmetry (and dark therefore, a net SM baryon asymmetry) to form. With this introduction, notice that Y`d 6= YL since χ1 or χ2 will have lepton number. The χ1’s andχ ¯1’s are populated by Φ decays and their Boltzmann equations are

dnχ1 + 3Hnχ = ΓΦnΦBr (Φ → χ1χ¯1) − hσvi n¯ nχ , dt 1 `d 1 dn χ¯1 + 3Hn = Γ n Br (Φ → χ χ¯ ) − hσvi n n , (A10) dt χ¯1 Φ Φ 1 1 `d χ¯1

χ2, B, and their conjugates initially have negligible abundances and we assume their abundances are always less than the abundances of `d, χ1, and their conjugates while these scattering processes are active. We have thus neglected the reverse scattering terms that would contribute to the above. The major modiﬁcation to the previous story (prior to the inclusion of the process in Eq. (11)) is the addition of ¯ scattering terms in the Boltzmann equations for `d and `d. Eq. (A5) becomes

dn`d + + + 3Hn = Γ + Br π → ` ` n + − hσvi n n , dt `d π d π `d χ¯1 dn¯ `d − − ¯ + 3Hn¯ = Γ + Br π → ` `d n − − hσvi n¯ nχ . (A11) dt `d π π `d 1 Taking the diﬀerence, we ﬁnd that the generated dark-sector lepton asymmetry may now be used to source equal and opposite dark-sector and SM baryon asymmetries: d D `d X f f f n` − n¯ + 3H n` − n¯ = 2Γ nΦBr N a Br + − hσvi nχ n` − n¯ . (A12) dt d `d d `d Φ π π CP D 1 d `d f

Note that nχ¯1 ≈ nχ1 for all times. Finally, the number density Boltzmann equation for the SM baryon asymmetry is simply

d (nB − n ) + 3H (nB − n ) = − hσvi nχ n` − n¯ , (A13) dt B B 1 d `d where again, we have assumed that the backreaction processes are negligible due to the minuscule number densities of χ2, B, and their conjugates relative to `d, χ1 and their conjugates. 13

+ f −2 f −2 D decay mode ACP /10 BrD+ /10

0 + KS π −0.41 ± 0.09 1.562 ± 0.031

K−π+π+ −0.18 ± 0.16 9.38 ± 0.16

K−π+π+π0 −0.3 ± 0.6 ± 0.4 5.98 ± 0.08 ± 0.16∗ [74]

0 + 0 ∗ KS π π −0.1 ± 0.7 ± 0.2 6.99 ± 0.09 ± 0.25 [74]

0 + + − ∗ KS π π π 0.0 ± 1.2 ± 0.3 3.122 ± 0.046 ± 0.096 [74]

π+π0 2.4 ± 1.2 (1.247 ± 0.033) × 10−1

π+η 1.0 ± 1.5 (3.77 ± 0.09) × 10−1

π+η0(958) −0.6 ± 0.7 (4.97 ± 0.19) × 10−1

K+K−π+ 0.37 ± 0.29 (9.35 ± 0.17 ± 0.24)∗ × 10−1 [74]

φπ+ 0.01 ± 0.09 (5.70 ± 0.05 ± 0.13) × 10−1

0 + +8 7.0 −2 a a0(1450) π −19 ± 12−11 4.5−1.8 × 10

+ +4.0 −3 b φ(1680)π −9 ± 22 ± 14 4.9−1.9 × 10

π+π+π− −1.7 ± 4.2 (3.27 ± 0.18) × 10−1

a + − this only includes the subsequent decay mode in which a0(1450) → K K b this only includes the subsequent decay mode in which φ(1680) → K+K−

TABLE III. Summary of D+ decay modes which violate CP and involve an odd number of π± and therefore help to generate a dark-sector lepton asymmetry.

Appendix B: D Meson Decay Modes

+ f We summarize the relevant D decay modes, including their values of ACP and branching ratios, in Table III. All quoted values come from the latest Particle Data Group (PDG) [75], with the following exception. There are some decay modes for which PDG does not provide their own ﬁt to the branching ratio, denoted by an asterisk on the branching ratio value. For these, we use the top listed reference within PDG and cite it in our table accordingly. 14

[1] P. A. R. Ade et al. (Planck), Planck 2015 results. XIII. els, JHEP 08, 157, arXiv:1606.06719 [hep-ph]. Cosmological parameters, Astron. Astrophys. 594, A13 [21] C. Cheung, G. Elor, L. J. Hall, and P. Kumar, Origins (2016), arXiv:1502.01589 [astro-ph.CO]. of Hidden Sector Dark Matter I: Cosmology, JHEP 03, [2] N. Aghanim et al. (Planck), Planck 2018 results. VI. 042, arXiv:1010.0022 [hep-ph]. Cosmological parameters, Astron. Astrophys. 641, A6 [22] X. Chu, T. Hambye, and M. H. Tytgat, The Four Basic (2020), arXiv:1807.06209 [astro-ph.CO]. Ways of Creating Dark Matter Through a Portal, JCAP [3] R. H. Cyburt, B. D. Fields, K. A. Olive, and T.-H. Yeh, 05, 034, arXiv:1112.0493 [hep-ph]. Big Bang Nucleosynthesis: 2015, Rev. Mod. Phys. 88, [23] Y. Hochberg, E. Kuflik, T. Volansky, and J. G. Wacker, 015004 (2016), arXiv:1505.01076 [astro-ph.CO]. Mechanism for Thermal Relic Dark Matter of Strongly [4] M. Tanabashi et al. (ParticleDataGroup), Review of Par- Interacting Massive Particles, Phys. Rev. Lett. 113, ticle Physics, Phys. Rev. D98, 030001 (2018). 171301 (2014), arXiv:1402.5143 [hep-ph]. [5] A. D. Sakharov, Violation of CP Invariance, c Asym- [24] G. Elor, H. Liu, T. R. Slatyer, and Y. Soreq, Comple- metry, and Baryon Asymmetry of the Universe, mentarity for Dark Sector Bound States, Phys. Rev. D Pisma Zh. Eksp. Teor. Fiz. 5, 32 (1967), [Usp. Fiz. 98, 036015 (2018), arXiv:1801.07723 [hep-ph]. Nauk161,61(1991)]. [25] Y. Tsai, L.-T. Wang, and Y. Zhao, Dark Matter Anni- [6] V. A. Kuzmin, V. A. Rubakov, and M. E. Shaposhnikov, hilation Decay at The LHC, Phys. Rev. D 93, 035024 On the Anomalous Electroweak Baryon Number Non- (2016), arXiv:1511.07433 [hep-ph]. conservation in the Early Universe, Phys. Lett. 155B, [26] C. Cheung, G. Elor, L. J. Hall, and P. Kumar, Origins of 36 (1985). Hidden Sector Dark Matter II: Collider Physics, JHEP [7] A. G. Cohen, D. B. Kaplan, and A. E. Nelson, 03, 085, arXiv:1010.0024 [hep-ph]. WEAK SCALE BARYOGENESIS, Phys. Lett. B245, [27] J. A. Dror, G. Elor, and R. McGehee, Directly Detecting 561 (1990). Signals from Absorption of Fermionic Dark Matter, Phys. [8] A. G. Cohen, D. B. Kaplan, and A. E. Nelson, Baryogen- Rev. Lett. 124, 18 (2020), arXiv:1905.12635 [hep-ph]. esis at the weak phase transition, Nucl. Phys. B349, 727 [28] J. A. Dror, G. Elor, and R. McGehee, Absorption of (1991). Fermionic Dark Matter by Nuclear Targets, JHEP 02, [9] N. Turok and J. Zadrozny, Dynamical generation of 134, arXiv:1908.10861 [hep-ph]. baryons at the electroweak transition, Phys. Rev. Lett. [29] J. A. Dror, G. Elor, R. McGehee, and T.-T. Yu, Absorp- 65, 2331 (1990). tion of sub-MeV fermionic dark matter by electron tar- [10] N. Turok and J. Zadrozny, Electroweak baryogenesis in gets, Phys. Rev. D 103, 035001 (2021), arXiv:2011.01940 the two doublet model, Nucl. Phys. B358, 471 (1991). [hep-ph]. [11] L. D. McLerran, M. E. Shaposhnikov, N. Turok, and [30] R. Essig, M. Fernandez-Serra, J. Mardon, A. Soto, M. B. Voloshin, Why the baryon asymmetry of the uni- T. Volansky, and T.-T. Yu, Direct Detection of sub-GeV verse is approximately ∼ 10−10, Phys. Lett. B256, 477 Dark Matter with Semiconductor Targets, JHEP 05, 046, (1991). arXiv:1509.01598 [hep-ph]. [12] M. Dine, P. Huet, R. L. Singleton, Jr, and L. Susskind, [31] G. Krnjaic, Freezing In, Heating Up, and Freezing Out: Creating the baryon asymmetry at the electroweak phase Predictive Nonthermal Dark Matter and Low-Mass Di- transition, Phys. Lett. B257, 351 (1991). rect Detection, JHEP 10, 136, arXiv:1711.11038 [hep- [13] A. G. Cohen, D. B. Kaplan, and A. E. Nelson, Sponta- ph]. neous baryogenesis at the weak phase transition, Phys. [32] M. Geller and O. Telem, Self Destructing Atomic DM, Lett. B263, 86 (1991). (2020), arXiv:2001.11514 [hep-ph]. [14] A. E. Nelson, D. B. Kaplan, and A. G. Cohen, Why there [33] D. Hooper, R. K. Leane, Y.-D. Tsai, S. Wegsman, and is something rather than nothing: Matter from weak in- S. J. Witte, A systematic study of hidden sector dark teractions, Nucl. Phys. B373, 453 (1992). matter:application to the gamma-ray and antiproton ex- [15] A. G. Cohen, D. B. Kaplan, and A. E. Nelson, De- cesses, JHEP 07 (07), 163, arXiv:1912.08821 [hep-ph]. bye screening and baryogenesis during the electroweak [34] G. Elor, N. L. Rodd, T. R. Slatyer, and W. Xue, Model- phase transition, Phys. Lett. B294, 57 (1992), arXiv:hep- Independent Indirect Detection Constraints on Hidden ph/9206214 [hep-ph]. Sector Dark Matter, JCAP 06, 024, arXiv:1511.08787 [16] G. R. Farrar and M. E. Shaposhnikov, Baryon asymme- [hep-ph]. try of the universe in the standard electroweak theory, [35] P. Barnes, Z. Johnson, A. Pierce, and B. Shakya, Simple Phys. Rev. D50, 774 (1994), arXiv:hep-ph/9305275 [hep- Hidden Sector Dark Matter, Phys. Rev. D 102, 075019 ph]. (2020), arXiv:2003.13744 [hep-ph]. [17] M. Fukugita and T. Yanagida, Baryogenesis Without [36] H. Davoudiasl and R. N. Mohapatra, On Relating the Grand Unification, Phys. Lett. B174, 45 (1986). Genesis of Cosmic Baryons and Dark Matter, New J. [18] V. Andreev et al. (ACME), Improved limit on the electric Phys. 14, 095011 (2012), arXiv:1203.1247 [hep-ph]. dipole moment of the electron, Nature 562, 355 (2018). [37] K. Petraki and R. R. Volkas, Review of asymmetric [19] J. A. Dror, T. Hiramatsu, K. Kohri, H. Murayama, and dark matter, Int. J. Mod. Phys. A28, 1330028 (2013), G. White, Testing the Seesaw Mechanism and Leptoge- arXiv:1305.4939 [hep-ph]. nesis with Gravitational Waves, Phys. Rev. Lett. 124, [38] K. M. Zurek, Asymmetric Dark Matter: Theories, Sig- 041804 (2020), arXiv:1908.03227 [hep-ph]. natures, and Constraints, Phys. Rept. 537, 91 (2014), [20] P. Hernández, M. Kekic, J. López-Pavón, J. Racker, arXiv:1308.0338 [hep-ph]. and J. Salvado, Testable Baryogenesis in Seesaw Mod- 15

[39] H. Davoudiasl, D. E. Morrissey, K. Sigurdson, and [57] D. Bryman and R. Shrock, Constraints on Sterile Neutri- S. Tulin, Hylogenesis: A Uniﬁed Origin for Baryonic Vis- nos in the MeV to GeV Mass Range, Phys. Rev. D 100, ible Matter and Antibaryonic Dark Matter, Phys. Rev. 073011 (2019), arXiv:1909.11198 [hep-ph]. Lett. 105, 211304 (2010), arXiv:1008.2399 [hep-ph]. [58] R. Abela, M. Daum, G. Eaton, R. Frosch, B. Jost, P. Ket- [40] E. Hall, T. Konstandin, R. McGehee, and H. Murayama, tle, and E. Steiner, Search for an Admixture of Heavy Asymmetric Matters from a Dark First-Order Phase Neutrino in Pion Decay, Phys. Lett. B 105, 263 (1981), Transition, (2019), arXiv:1911.12342 [hep-ph]. [Erratum: Phys.Lett.B 106, 513 (1981)]. [41] J. Shelton and K. M. Zurek, Darkogenesis: A baryon [59] A. Aguilar-Arevalo et al. (PIENU), Search for heavy neu- asymmetry from the dark matter sector, Phys. Rev. D82, trinos in π → µν decay, Phys. Lett. B 798, 134980 (2019), 123512 (2010), arXiv:1008.1997 [hep-ph]. arXiv:1904.03269 [hep-ex]. [42] G. Servant and S. Tulin, Baryogenesis and Dark Mat- [60] D. Bryman, Private communication., private communi- ter through a Higgs Asymmetry, Phys. Rev. Lett. 111, cation. 151601 (2013), arXiv:1304.3464 [hep-ph]. [61] D. McKeen, A. E. Nelson, S. Reddy, and D. Zhou, Neu- [43] E. Hall, T. Konstandin, R. McGehee, H. Murayama, and tron stars exclude light dark baryons, Phys. Rev. Lett. G. Servant, Baryogenesis From a Dark First-Order Phase 121, 061802 (2018), arXiv:1802.08244 [hep-ph]. Transition, JHEP 04, 042, arXiv:1910.08068 [hep-ph]. [62] K. Harigaya, R. McGehee, H. Murayama, and K. Schutz, [44] E. W. Kolb and M. S. Turner, The Early Universe, Front. A predictive mirror twin Higgs with small Z2 breaking, Phys. 69, 1 (1990). JHEP 05, 155, arXiv:1905.08798 [hep-ph]. [45] S. Dimopoulos and L. J. Hall, Baryogenesis at the {MeV} [63] Y.-D. Tsai, R. McGehee, and H. Murayama, Resonant Era, Phys. Lett. B 196, 135 (1987). Self-Interacting Dark Matter from Dark QCD, (2020), [46] T. Cohen, D. J. Phalen, and A. Pierce, Supersymmetric arXiv:2008.08608 [hep-ph]. Baryogenesis from Exotic Quark Decays, Phys. Rev. D [64] A. V. Artamonov et al. (E949), Search for heavy neutri- + + 81, 035020 (2010), arXiv:0908.2998 [hep-ph]. nos in K → µ νH decays, Phys. Rev. D91, 052001 [47] N. Arkani-Hamed, T. Cohen, R. T. D’Agnolo, A. Hook, (2015), [Erratum: Phys. Rev.D91,no.5,059903(2015)], H. D. Kim, and D. Pinner, Solving the Hierarchy arXiv:1411.3963 [hep-ex]. Problem at Reheating with a Large Number of De- [65] E. Cortina Gil et al. (NA62), Search for heavy neutral grees of Freedom, Phys. Rev. Lett. 117, 251801 (2016), lepton production in K+ decays, Phys. Lett. B778, 137 arXiv:1607.06821 [hep-ph]. (2018), arXiv:1712.00297 [hep-ex]. [48] P. W. Graham, D. E. Kaplan, and S. Rajendran, Cosmo- [66] R. Aaij et al. (LHCb), Updated determination of D0- logical Relaxation of the Electroweak Scale, Phys. Rev. D0 mixing and CP violation parameters with D0 → Lett. 115, 221801 (2015), arXiv:1504.07551 [hep-ph]. K+π− decays, Phys. Rev. D 97, 031101 (2018), [49] G. Elor, M. Escudero, and A. Nelson, Baryogenesis and arXiv:1712.03220 [hep-ex]. Dark Matter from B Mesons, Phys. Rev. D 99, 035031 [67] Y. Nir, Lessons from BaBar and Belle measurements (2019), arXiv:1810.00880 [hep-ph]. of D0 - anti-D0 mixing parameters, JHEP 05, 102, [50] G. Alonso-Alvarez,´ G. Elor, A. E. Nelson, and H. Xiao, arXiv:hep-ph/0703235 [hep-ph]. A Supersymmetric Theory of Baryogenesis and Sterile [68] E. Golowich, J. Hewett, S. Pakvasa, and A. A. Petrov, Sneutrino Dark Matter from B Mesons, JHEP 03, 046, Implications of D0 - D¯ 0 Mixing for New Physics, Phys. arXiv:1907.10612 [hep-ph]. Rev. D76, 095009 (2007), arXiv:0705.3650 [hep-ph]. [51] G. Alonso-Alvarez,´ G. Elor, and M. Escudero, Collider [69] G. Barenboim, W. H. Kinney, and W.-I. Park, Flavor Signals of Baryogenesis and Dark Matter from B Mesons: versus mass eigenstates in neutrino asymmetries: impli- A Roadmap to Discovery, (2021), arXiv:2101.02706 [hep- cations for cosmology, Eur. Phys. J. C77, 590 (2017), ph]. arXiv:1609.03200 [astro-ph.CO]. [52] T. Hasegawa, N. Hiroshima, K. Kohri, R. S. Hansen, [70] T. Venumadhav, F.-Y. Cyr-Racine, K. N. Abazajian, and T. Tram, and S. Hannestad, MeV-scale reheating tem- C. M. Hirata, Sterile neutrino dark matter: Weak inter- perature and thermalization of oscillating neutrinos by actions in the strong coupling epoch, Phys. Rev. D 94, radiative and hadronic decays of massive particles, JCAP 043515 (2016), arXiv:1507.06655 [astro-ph.CO]. 12, 012, arXiv:1908.10189 [hep-ph]. [71] R. J. Scherrer and M. S. Turner, Primordial Nucleosyn- [53] M. Kawasaki, K. Kohri, and N. Sugiyama, Cosmological thesis with Decaying Particles. 1. Entropy Producing De- constraints on late time entropy production, Phys. Rev. cays. 2. Inert Decays, Astrophys. J. 331, 19 (1988). Lett. 82, 4168 (1999), arXiv:astro-ph/9811437. [72] S. Hannestad, What is the lowest possible reheating tem- [54] M. Kawasaki, K. Kohri, and N. Sugiyama, MeV perature?, Phys. Rev. D 70, 043506 (2004), arXiv:astro- scale reheating temperature and thermalization of neu- ph/0403291. trino background, Phys. Rev. D 62, 023506 (2000), [73] M. Laine and Y. Schroder, Quark mass thresholds in arXiv:astro-ph/0002127. QCD thermodynamics, Phys. Rev. D73, 085009 (2006), [55] A. Ibarra and S. Wild, Dirac dark matter with a arXiv:hep-ph/0603048 [hep-ph]. charged mediator: a comprehensive one-loop analysis [74] S. Dobbs et al. (CLEO), Measurement of absolute of the direct detection phenomenology, JCAP 05, 047, hadronic branching fractions of D mesons and e+ e- —> arXiv:1503.03382 [hep-ph]. D anti-D cross-sections at the psi(3770), Phys. Rev. D [56] A. Aguilar-Arevalo et al. (PIENU), Improved search for 76, 112001 (2007), arXiv:0709.3783 [hep-ex]. heavy neutrinos in the decay π → eν, Phys. Rev. D 97, [75] P. Zyla et al. (Particle Data Group), Review of Particle 072012 (2018), arXiv:1712.03275 [hep-ex]. Physics, PTEP 2020, 083C01 (2020).