We study the Schr\"odinger-Poisson (SP) method in the context of cosmological large-scale structure formation in an expanding background. In the limit $\hbar \to 0$, the SP technique can be viewed as an effective method to sample the phase space distribution of cold dark matter that remains valid on non-linear scales. We present results for the 2D and 3D matter correlation function and power spectrum at length scales corresponding to the baryon acoustic oscillation (BAO) peak. We discuss systematic effects of the SP method applied to cold dark matter and explore how they depend on the simulation parameters. In particular, we identify a combination of simulation parameters that controls the scale-independent loss of power observed at low redshifts, and discuss the scale relevant to this effect.

Dark matter evolution during the process of cosmological structure formation can be described in terms of a one-particle irreducible effective action at a characteristic scale $k_m$ and a loop expansion below this scale, based on the effective propagators and vertices. We calculate the form of the effective vertices and compute the bispectrum of density perturbations within a one-loop approximation. We find that the effective vertices play a subdominant role as compared to the effective viscosity and sound velocity that modify the (inverse) propagators. For the bispectrum we reproduce the results of standard perturbation theory in the range where it is applicable, and find a slightly improved agreement with N-body simulations at larger wavenumbers.

We consider dark matter (DM) with very weak couplings to the standard model (SM), such that its self-annihilation cross section is much smaller than the canonical one, $\langle\sigma v\rangle_{\chi\chi} \ll 10^{-26}\mathrm{cm}^3/\mathrm{s}$. In this case DM self-annihilation is negligible for the dynamics of freeze-out and DM dilution is solely driven by efficient annihilation of heavier accompanying dark sector particles provided that DM maintains chemical equilibrium with the dark sector. This chemical equilibrium is established by conversion processes which require much smaller couplings to be efficient than annihilation. The chemical decoupling of DM from the SM can either be initiated by the freeze-out of annihilation, resembling a co-annihilation scenario, or of conversion processes, leading to the scenario of conversion-driven freeze-out. We focus on the latter and discuss its distinct phenomenology.

Non-thermalized dark matter is a cosmologically valid alternative to the paradigm of weakly interacting massive particles. For dark matter belonging to a $Z_2$-odd sector that contains in addition a thermalized mediator particle, dark matter production proceeds in general via both the freeze-in and superWIMP mechanism. We highlight their interplay and emphasize the connection to long-lived particles at colliders. For the explicit example of a colored t-channel mediator model we map out the entire accessible parameter space, cornered by bounds from the LHC, big bang nucleosynthesis and Lyman-alpha forest observations, respectively. We discuss prospects for the HL- and HE-LHC.

A minimal model of Cold Dark Matter (CDM) is a very massive particle with only gravitational interactions, also called Planckian Interacting Dark Matter (PIDM). Here we consider an extension of the PIDM framework by an unbroken $U(1)$ gauge symmetry under which the PIDM is charged, but remains only gravitationally coupled to the Standard Model (SM). Contrary to "hidden charged dark matter", the charged PIDM never reaches thermal equilibrium with the SM. The dark sector is populated by freeze-in via gravitational interactions at reheating. If the dark fine-structure constant $\alpha_D$ is larger than about $10^{-3}$, the dark sector thermalizes within itself, and the PIDM abundance is further modified by freeze-out in the dark sector. Interestingly, this largely reduces the dependence of the final abundance on the reheating temperature, as compared to an uncharged PIDM. Thermalization within the dark sector is driven by inelastic radiative processes, and affected by the Landau-Pomeranchuk-Migdal (LPM) effect. The observed CDM abundance can be obtained over a wide mass range from the weak to the GUT scale, and for phenomenologically interesting couplings $\alpha_D\sim 10^{-2}$. Due to the different thermal history, the charged PIDM can be discriminated from "hidden charged dark matter" by more precise measurements of the effective number of neutrino species $N_{\rm eff}$.

Planckian Interacting Dark Matter (PIDM) is a minimal scenario of dark matter assuming only gravitational interactions with the standard model and with only one free parameter, the PIDM mass. PIDM can be successfully produced by gravitational scattering in the thermal plasma of the Standard Model sector after inflation in the PIDM mass range from TeV up to the GUT scale, if the reheating temperature is sufficiently high. The minimal assumption of a GUT scale PIDM mass can be tested in the future by measurements of the primordial tensor-to-scalar ratio. While large primordial tensor modes would be in tension with the QCD axion as dark matter in a large mass range, it would favour the PIDM as a minimal alternative to WIMPs. Here we generalise the previously studied scalar PIDM scenario to the case of fermion, vector and tensor PIDM scenarios, and show that the phenomenology is nearly identical, independent of the spin of the PIDM. We also consider the specific realisation of the PIDM as the Kaluza Klein excitation of the graviton in orbifold compactifications of string theory, as well as in models of monodromy inflation and in Higgs inflation. Finally we discuss the possibility of indirect detection of PIDM through non-perturbative decay.

The Lyman-$\alpha$ forest is a valuable probe of dark matter models featuring a scale-dependent suppression of the power spectrum as compared to $\Lambda$CDM. In this work, we present a new estimator of the Lyman-$\alpha$ flux power spectrum that does not rely on hydrodynamical simulations. Our framework is characterized by nuisance parameters that encapsulate the complex physics of the intergalactic medium and sensitivity to highly non-linear small-scale modes. After validating the approach based on high-resolution hydrodynamical simulations for $\Lambda$CDM, we derive conservative constraints on interacting dark matter models from BOSS Lyman-$\alpha$ data on large scales, k<0.02(km/s)^(-1), with the relevant nuisance parameters left free in the model fit. The estimator yields lower bounds on the mass of cannibal dark matter, where freeze-out occurs through 3-to-2 annihilation, in the MeV range. Furthermore, we find that models of dark matter interacting with dark radiation, which have been argued to address the $H_0$ and $\sigma_8$ tensions, are compatible with BOSS Lyman-$\alpha$ data.

We present a comprehensive analysis of top-philic Majorana dark matter that interacts via a colored t-channel mediator. Despite the simplicity of the model -- introducing three parameters only -- it provides an extremely rich phenomenology allowing us to accommodate the relic density for a large range of coupling strengths spanning over six orders of magnitude. This model features all `exceptional' mechanisms for dark matter freeze-out, including the recently discovered conversion-driven freeze-out mode, with interesting signatures of long-lived colored particles at colliders. We constrain the cosmologically allowed parameter space with current experimental limits from direct, indirect and collider searches, with special emphasis on light dark matter below the top mass. In particular, we explore the interplay between limits from Xenon1T, Fermi-LAT and AMS-02 as well as limits from stop, monojet and Higgs invisible decay searches at the LHC. We find that several blind spots for light dark matter evade current constraints. The region in parameter space where the relic density is set by the mechanism of conversion-driven freeze-out can be conclusively tested by R-hadron searches at the LHC with 300\,fb$^{-1}$.

We perform a quantitative comparison between N-body simulations and the Schr\"odinger-Poisson system in 1+1 dimensions. In particular, we study halo formation with different initial conditions. We observe the convergence of various observables in the Planck constant h-bar and also test virialization. We discuss the generation of higher order cumulants of the particle distribution function which demonstrates that the Schr\"odinger-Poisson equations should not be perceived as a generalization of the dust model with quantum pressure but rather as one way of sampling the phase space of the Vlasov-Poisson system -- just as N-body simulations. Finally, we quantitatively recover the scaling behavior of the halo density profile from N-body simulations.

The Standard Model could be self-consistent up to the Planck scale according to the present measurements of the Higgs mass and top quark Yukawa coupling. It is therefore possible that new physics is only coupled to the Standard Model through Planck suppressed higher dimensional operators. In this case the WIMP miracle is a mirage, and instead minimality as dictated by Occam's razor would indicate that dark matter is related to the Planck scale, where quantum gravity is anyway expected to manifest itself. Assuming within this framework that dark matter is a Planckian Interacting Massive Particle, we show that the most natural mass larger than $0.01\,\textrm{M}_p$ is already ruled out by the absence of tensor modes in the CMB. This also indicates that we expect tensor modes in the CMB to be observed soon for this type of minimal dark matter model. Finally, we touch upon the KK graviton mode as a possible realization of this scenario within UV complete models, as well as further potential signatures and peculiar properties of this type of dark matter candidate. This paradigm therefore leads to a subtle connection between quantum gravity, the physics of primordial inflation, and the nature of dark matter.

Chemical equilibrium is a commonly made assumption in the freeze-out calculation of coannihilating dark matter. We explore the possible failure of this assumption and find a new conversion-driven freeze-out mechanism. Considering a representative simplified model inspired by supersymmetry with a neutralino- and sbottom-like particle we find regions in parameter space with very small couplings accommodating the measured relic density. In this region freeze-out takes place out of chemical equilibrium and dark matter self-annihilation is thoroughly inefficient. The relic density is governed primarily by the size of the conversion terms in the Boltzmann equations. Due to the small dark matter coupling the parameter region is immune to direct detection but predicts an interesting signature of disappearing tracks or displaced vertices at the LHC. Unlike freeze-in or superWIMP scenarios, conversion-driven freeze-out is not sensitive to the initial conditions at the end of reheating.

Treating dark matter at large scales as an effectively viscous fluid provides an improved framework for the calculation of the density and velocity power spectra compared to the standard assumption of an ideal pressureless fluid. We discuss how this framework can be made concrete through an appropriate coarse-graining procedure. We also review results that demonstrate that it improves the convergence of cosmological perturbation theory.

The measured (central) values of the Higgs and top quark masses indicate that the Standard Model (SM) effective potential develops an instability at high field values. The scale of this instability, determined as the Higgs field value at which the potential drops below the electroweak minimum, is about $10^{11}$ GeV. However, such a scale is unphysical as it is not gauge-invariant and suffers from a gauge-fixing uncertainty of up to two orders of magnitude. Subjecting our system, the SM, to several probes of the instability (adding higher order operators to the potential; letting the vacuum decay through critical bubbles; heating up the system to very high temperature; inflating it) and asking in each case physical questions, we are able to provide several gauge-invariant scales related with the Higgs potential instability.

We present a new analytic approach to describe large scale structure formation in the mildly non-linear regime. The central object of the method is the time-dependent probability distribution function generating correlators of the cosmological observables at a given moment of time. Expanding the distribution function around the Gaussian weight we formulate a perturbative technique to calculate non-linear corrections to cosmological correlators, similar to the diagrammatic expansion in a three-dimensional Euclidean quantum field theory, with time playing the role of an external parameter. For the physically relevant case of cold dark matter in an Einstein--de Sitter universe, the time evolution of the distribution function can be found exactly and is encapsulated by a time-dependent coupling constant controlling the perturbative expansion. We show that all building blocks of the expansion are free from spurious infrared enhanced contributions that plague the standard cosmological perturbation theory. This paves the way towards the systematic resummation of infrared effects in large scale structure formation. We also argue that the approach proposed here provides a natural framework to account for the influence of short-scale dynamics on larger scales along the lines of effective field theory.

We use time-sliced perturbation theory (TSPT) to give an accurate description of the infrared non-linear effects affecting the baryonic acoustic oscillations (BAO) present in the distribution of matter at very large scales. In TSPT this can be done via a systematic resummation that has a simple diagrammatic representation and does not involve uncontrollable approximations. We discuss the power counting rules and derive explicit expressions for the resummed matter power spectrum up to next-to leading order and the bispectrum at the leading order. The two-point correlation function agrees well with N-body data at BAO scales. The systematic approach also allows to reliably assess the shift of the baryon acoustic peak due to non-linear effects.

Following an approach of Matarrese and Pietroni, we derive the functional renormalization group (RG) flow of the effective action of cosmological large-scale structures. Perturbative solutions of this RG flow equation are shown to be consistent with standard cosmological perturbation theory. Non-perturbative approximate solutions can be obtained by truncating the a priori infinite set of possible effective actions to a finite subspace. Using for the truncated effective action a form dictated by dissipative fluid dynamics, we derive RG flow equations for the scale dependence of the effective viscosity and sound velocity of non-interacting dark matter, and we solve them numerically. Physically, the effective viscosity and sound velocity account for the interactions of long-wavelength fluctuations with the spectrum of smaller-scale perturbations. We find that the RG flow exhibits an attractor behaviour in the IR that significantly reduces the dependence of the effective viscosity and sound velocity on the input values at the UV scale. This allows for a self-contained computation of matter and velocity power spectra for which the sensitivity to UV modes is under control.

We derive a non-perturbative equation for the large scale structure power spectrum of long-wavelength modes. Thereby, we use an operator product expansion together with relations between the three-point function and power spectrum in the soft limit. The resulting equation encodes the coupling to ultraviolet (UV) modes in two time-dependent coefficients, which may be obtained from response functions to (anisotropic) parameters, such as spatial curvature, in a modified cosmology. We argue that both depend weakly on fluctuations deep in the UV. As a byproduct, this implies that the renormalized leading order coefficient(s) in the effective field theory (EFT) of large scale structures receive most of their contribution from modes close to the non-linear scale. Consequently, the UV dependence found in explicit computations within standard perturbation theory stems mostly from counter-term(s). We confront a simplified version of our non-perturbative equation against existent numerical simulations, and find good agreement within the expected uncertainties. Our approach can in principle be used to precisely infer the relevance of the leading order EFT coefficient(s) using small volume simulations in an `anisotropic separate universe' framework. Our results suggest that the importance of these coefficient(s) is a $\sim 10 \%$ effect, and plausibly smaller.

Cosmological perturbations of sufficiently long wavelength admit a fluid dynamic description. We consider modes with wavevectors below a scale $k_m$ for which the dynamics is only mildly non-linear. The leading effect of modes above that scale can be accounted for by effective non-equilibrium viscosity and pressure terms. For mildly non-linear scales, these mainly arise from momentum transport within the ideal and cold but inhomogeneous fluid, while momentum transport due to more microscopic degrees of freedom is suppressed. As a consequence, concrete expressions with no free parameters, except the matching scale $k_m$, can be derived from matching evolution equations to standard cosmological perturbation theory. Two-loop calculations of the matter power spectrum in the viscous theory lead to excellent agreement with $N$-body simulations up to scales $k=0.2 \, h/$Mpc. The convergence properties in the ultraviolet are better than for standard perturbation theory and the results are robust with respect to variations of the matching scale.

We discuss the renormalization of the initial value problem in Nonequilibrium Quantum Field Theory within a simple, yet instructive, example and show how to obtain a renormalized time evolution for the two-point functions of a scalar field and its conjugate momentum at all times. The scheme we propose is applicable to systems that are initially far from equilibrium and compatible with non-secular approximation schemes which capture thermalization. It is based on Kadanoff-Baym equations for non-Gaussian initial states, complemented by usual vacuum counterterms. We explicitly demonstrate how various cutoff-dependent effects peculiar to nonequilibrium systems, including time-dependent divergences or initial-time singularities, are avoided by taking an initial non-Gaussian three-point vacuum correlation into account.

We compute non-linear corrections to the matter power spectrum taking the time- and scale-dependent free-streaming length of neutrinos into account. We adopt a hybrid scheme that matches the full Boltzmann hierarchy to an effective two-fluid description at an intermediate redshift. The non-linearities in the neutrino component are taken into account by using an extension of the time-flow framework. We point out that this remedies a spurious behaviour that occurs when neglecting non-linear terms for neutrinos. This behaviour is related to how efficiently short modes decouple from long modes and can be traced back to the violation of momentum conservation if neutrinos are treated linearly. Furthermore, we compare our results at next to leading order to various other methods and quantify the accuracy of the fluid description. Due to the correct decoupling behaviour of short modes, the two-fluid scheme is a suitable starting point to compute higher orders in perturbations or for resummation methods.

We analyze the dark matter power spectrum at three-loop order in standard perturbation theory of large scale structure. We observe that at late times the loop expansion does not converge even for large scales (small momenta) well within the linear regime, but exhibits properties compatible with an asymptotic series. We propose a technique to restore the convergence in the limit of small momentum, and use it to obtain a perturbative expansion with improved convergence for momenta in the range where baryonic acoustic oscillations are present. Our results are compared with data from N-body simulations at different redshifts, and we find good agreement within this range.

Recently reported tentative evidence for a gamma-ray line in the Fermi-LAT data is of great potential interest for identifying the nature of dark matter. We compare the implications for decaying and annihilating dark matter taking the constraints from continuum gamma-rays, antiproton flux and morphology of the excess into account. We find that higgsino and wino dark matter are excluded, also for nonthermal production. Generically, the continuum gamma-ray flux severely constrains annihilating dark matter. Consistency of decaying dark matter with the spatial distribution of the Fermi-LAT excess would require an enhancement of the dark matter density near the Galactic center.

Although proposed long ago, the search for internal bremsstrahlung signatures has only recently been made possible by the excellent energy resolution of ground-based and satellite-borne gamma-ray instruments. Here, we investigate thoroughly the current status of internal bremsstrahlung searches in light of the results of direct dark matter searches and in the framework of a minimal mass-degenerate scenario consisting of a Majorana dark matter particle that couples to a fermion and a scalar via a Yukawa coupling. The upper limits on the annihilation cross section set by Fermi-LAT and H.E.S.S. extend uninterrupted from tens of GeV up to tens of TeV and are rather insensitive to the mass degeneracy in the particle physics model. In contrast, direct searches are best in the moderate to low mass splitting regime, where XENON100 limits overshadow Fermi-LAT and H.E.S.S. up to TeV masses if dark matter couples to one of the light quarks. In our minimal scenario we examine carefully the prospects for GAMMA-400, CTA and XENON1T, all planned to come online in the near future, and find that: (a) CTA and XENON1T are fully complementary, with CTA most sensitive to multi-TeV masses and mass splittings around 10%, and XENON1T probing best small mass splittings up to TeV masses; and (b) current constraints from XENON100 already preclude the observation of any spectral feature with GAMMA-400 in spite of its impressive energy resolution, unless dark matter does not couple predominantly to light quarks. Finally, we point out that, unlike for direct searches, the possibility of detecting thermal relics in upcoming internal bremsstrahlung searches requires, depending on the concrete scenario, boost factors larger than 5-10.

We discuss the convergence of cosmological perturbation theory. We prove that the polynomial enhancement of the non-linear corrections expected from the effects of soft modes is absent in equal-time correlators like the power or bispectrum. We first show this at leading order by resumming the most important corrections of soft modes to an arbitrary skeleton of hard fluctuations. We derive the same result in the eikonal approximation, which also allows us to show the absence of enhancement at any order. We complement the proof by an explicit calculation of the power spectrum at two-loop order, and by further numerical checks at higher orders. Using these insights, we argue that the modification of the power spectrum from soft modes corresponds at most to logarithmic corrections. Finally, we discuss the asymptotic behavior in the large and small momentum regimes and identify the expansion parameter pertinent to non-linear corrections.

The latest XENON100 data severely constrains dark matter elastic scattering off nuclei, leading to impressive upper limits on the spin-independent cross-section. The main goal of this paper is to stress that the same data set has also an excellent \emph{spin-dependent} sensitivity, which is of utmost importance in probing dark matter models. We show in particular that the constraints set by XENON100 on the spin-dependent neutron cross-section are by far the best at present, whereas the corresponding spin-dependent proton limits lag behind other direct detection results. The effect of nuclear uncertainties on the structure functions of xenon isotopes is analysed in detail and found to lessen the robustness of the constraints, especially for spin-dependent proton couplings. Notwithstanding, the spin-dependent neutron prospects for XENON1T and DARWIN are very encouraging. We apply our constraints to well-motivated dark matter models and demonstrate that in both mass-degenerate scenarios and the minimal supersymmetric standard model the spin-dependent neutron limits can actually override the spin-independent limits. This opens the possibility of probing additional unexplored regions of the dark matter parameter space with the next generation of ton-scale direct detection experiments.

Over the last years both cosmic-ray antiproton measurements and direct dark matter searches have proved particularly effective in constraining the nature of dark matter candidates. The present work focusses on these two types of constraints in a minimal framework which features a Majorana fermion as the dark matter particle and a scalar that mediates the coupling to quarks. Considering a wide range of coupling schemes, we derive antiproton and direct detection constraints using the latest data and paying close attention to astrophysical and nuclear uncertainties. Both signals are strongly enhanced in the presence of degenerate dark matter and scalar masses, but we show that the effect is especially dramatic in direct detection. Accordingly, the latest direct detection limits take the lead over antiprotons. We find that antiproton and direct detection data set stringent lower limits on the mass splitting, reaching 19% at a 300 GeV dark matter mass for a unity coupling. Interestingly, these limits are orthogonal to ongoing collider searches at the Large Hadron Collider, making it feasible to close in on degenerate dark matter scenarios within the next years.

We present general constraints on dark matter stability in hadronic decay channels derived from measurements of cosmic-ray antiprotons.We analyze various hadronic decay modes in a model-independent manner by examining the lowest-order decays allowed by gauge and Lorentz invariance for scalar and fermionic dark matter particles and present the corresponding lower bounds on the partial decay lifetimes in those channels. We also investigate the complementarity between hadronic and gamma-ray constraints derived from searches for monochromatic lines in the sky, which can be produced at the quantum level if the dark matter decays into quark-antiquark pairs at leading order.

We derive solutions to the Schwinger-Dyson equations on the Closed-Time-Path for a scalar field in the limit where backreaction is neglected. In Wigner space, the two-point Wightman functions have the curious property that the equilibrium component has a finite width, while the out-of equilibrium component has zero width. This feature is confirmed in a numerical simulation for scalar field theory with quartic interactions. When substituting these solutions into the collision term, we observe that an expansion including terms of all orders in gradients leads to an effective finite-width. Besides, we observe no breakdown of perturbation theory, that is sometimes associated with pinch singularities. The effective width is identical with the width of the equilibrium component. Therefore, reconciliation between the zero-width behaviour and the usual notion in kinetic theory, that the out-of-equilibrium contributions have a finite width as well, is achieved. This result may also be viewed as a generalisation of the fluctuation-dissipation relation to out-of-equilibrium systems with negligible backreaction.

If the dark matter particle is a Majorana fermion, annihilations into two fermions and one gauge boson could have, for some choices of the parameters of the model, a non-negligible cross-section. Using a toy model of leptophilic dark matter, we calculate the constraints on the annihilation cross-section into two electrons and one weak gauge boson from the PAMELA measurements of the cosmic antiproton-to-proton flux ratio. Furthermore, we calculate the maximal astrophysical boost factor allowed in the Milky Way under the assumption that the leptophilic dark matter particle is the dominant component of dark matter in our Universe. These constraints constitute very conservative estimates on the boost factor for more realistic models where the dark matter particle also couples to quarks and weak gauge bosons, such as the lightest neutralino which we also analyze for some concrete benchmark points. The limits on the astrophysical boost factors presented here could be used to evaluate the prospects to detect a gamma-ray signal from dark matter annihilations at currently operating IACTs as well as in the projected CTA.

The decay of dark matter particles which are coupled predominantly to charged leptons has been proposed as a possible origin of excess high-energy positrons and electrons observed by cosmic-ray telescopes PAMELA and Fermi LAT. Even though the dark matter itself is electrically neutral, the tree-level decay of dark matter into charged lepton pairs will generically induce radiative two-body decays of dark matter at the quantum level. Using an effective theory of leptophilic dark matter decay, we calculate the rates of radiative two-body decays for scalar and fermionic dark matter particles. Due to the absence of astrophysical sources of monochromatic gamma rays, the observation of a line in the diffuse gamma-ray spectrum would constitute a strong indication of a particle physics origin of these photons. We estimate the intensity of the gamma-ray line that may be present in the energy range of a few TeV if the dark matter decay interpretation of the leptonic cosmic-ray anomalies is correct and comment on observational prospects of present and future Imaging Cherenkov Telescopes, in particular the CTA.

We discuss a toy model where baryogenesis and cosmic acceleration are driven by a leptonic quintessence field coupled to the standard model sector via a massive mediating scalar field. It does not require the introduction of B-L-violating interactions below the inflationary scale. Instead, a B-L-asymmetry is stored in the quintessence field, which compensates for the corresponding observed baryon asymmetry.

We investigate the impact of quantum fluctuations on a light rolling quintessence field from three different sources, namely, from a coupling to the standard model and dark matter, from its self-couplings and from its coupling to gravity. We derive bounds for time-varying masses from the change of vacuum energy, finding \Delta m_e/m_e << 10^{-11} for the electron and \Delta m_p/m_p << 10^{-15} for the proton since redshift z~2, whereas the neutrino masses could change of order one. Mass-varying dark matter is also constrained. Next, the self-interactions are investigated. For inverse power law potentials, the effective potential does not become infinitely large at small field values, but saturates at a finite maximal value. We discuss implications for cosmology. Finally, we show that one-loop corrections induce non-minimal gravitational couplings involving arbitrarily high powers of the curvature scalar R, indicating that quintessence entails modified gravity effects.

We consider a baryogenesis scenario, where the difference of baryon (B) and lepton (L) number is conserved in such a way that the B-L asymmetry in the standard model sector is compensated by an asymmetry of opposite sign stored in the dark energy sector. Therefore, we introduce a toy-model in which a complex quintessence field carries a B-L asymmetry at late times. We determine the produced baryon asymmetry in the visible sector for a large range of initial conditions and find it easy to achieve a value of the observed order of magnitude. While the size of the produced baryon asymmetry depends on details of the underlying inflationary model, it turns out to be independent of the reheating temperature in many cases. We also discuss possible sources of instability like the formation of Q-Balls.