Search

We apply deep learning techniques to the late-time turbulent regime in a post-inflationary model where a real scalar inflaton field and the standard model Higgs doublet interact with renormalizable couplings between them. After inflation, the inflaton decays into the Higgs through a trilinear coupling and the Higgs field subsequently thermalizes with gauge bosons via its $SU(2)\times U(1)$ gauge interaction. Depending on the strength of the trilinear interaction and the Higgs self-coupling, the effective mass squared of Higgs can become negative, leading to the tachyonic production of Higgs particles. These produced Higgs particles would then share their energy with gauge bosons, potentially indicating thermalization. Since the model entails different non-perturbative effects, it is necessary to resort to numerical and semi-classical techniques. However, simulations require significant costs in terms of time and computational resources depending on the model used. Particularly, when $SU(2)$ gauge interactions are introduced, this becomes evident as the gauge field redistributes particle energies through rescattering processes, leading to an abundance of UV modes that disrupt simulation stability. This necessitates very small lattice spacings, resulting in exceedingly long simulation runtimes. Furthermore, the late-time behavior of preheating dynamics exhibits a universal form by wave kinetic theory. Therefore, we analyze patterns in the flow of particle numbers and predict future behavior using CNN-LSTM (Convolutional Neural Network combined with Long Short-Term Memory) time series analysis. In this way, we can reduce our dependence on simulations by orders of magnitude in terms of time and computational resources., Comment: 13 pages, 5 figures; references added

We consider the effects of fragmentation on the post-inflationary epoch of reheating. In simple single field models of inflation, an inflaton condensate undergoes an oscillatory phase once inflationary expansion ends. The equation of state of the condensate depends on the shape of the scalar potential, $V(\phi)$, about its minimum. Assuming $V(\phi) \sim \phi^k$, the equation of state parameter is given by $w = P_\phi/\rho_\phi = (k-2)/(k+2)$. The evolution of condensate and the reheating process depend on $k$. For $k \ge 4$, inflaton self-interactions may lead to the fragmentation of the condensate and alter the reheating process. Indeed, these self-interactions lead to the production of a massless gas of inflaton particles as $w$ relaxes to 1/3. If reheating occurs before fragmentation, the effects of fragmentation are harmless. We find, however, that the effects of fragmentation depend sensitively to the specific reheating process. Reheating through the decays to fermions is largely excluded since perturbative couplings would imply that fragmentation occurs before reheating and in fact could prevent reheating from completion. Reheating through the decays to boson is relatively unaffected by fragmentation and reheating through scatterings results in a lower reheating temperature., Comment: 31 pages, 10 figures; v2 matches version published in JCAP

We study inflationary models based on a non-minimal coupling of a singlet scalar to gravity, focussing on the preheating dynamics and the unitarity issues in this regime. If the scalar does not have significant couplings to other fields, particle production after inflation is far less efficient than that in Higgs inflation. As a result, unitarity violation at large non-minimal couplings requires a different treatment. We find that collective effects in inflaton scattering processes during preheating make an important impact on the unitarity constraint. Within effective field theory, the consequent upper bound on the non-minimal coupling is of order a few hundreds., Comment: 16 pages, 7 figures

We study postinflationary scalar dark matter production via its non-minimal coupling to gravity. During the inflaton oscillation epoch, dark matter is produced resonantly for a sufficiently large non-minimal coupling $\xi\gtrsim 5$. We find that backreaction on the curvature and rescattering effects typically become important for the values of $\xi$ above $30$, which invalidate simple estimates of the production efficiency. At large couplings, the dark matter yield becomes almost independent of $\xi$, signifying approximate quasi-equilibrium in the inflaton-dark matter system. Although the analysis gets complicated by the presence of apparent negative energy in the Jordan frame, this behaviour can be regularized by introducing mild dark matter self-interaction. Using lattice simulations, we delineate parameter space leading to the correct dark matter relic abundance., Comment: 18 pages, 7 figures; v2: matches version published in JCAP (added minor clarifications, fixed typos, added references)

We discuss a minimal extension of the Standard Model (SM) where a single real scalar field serves as both inflaton and dark matter. The corresponding Lagrangian contains the renormalizable interactions of the inflaton field. Quantum effects generally induce a non-minimal coupling to gravity which facilitates inflation consistent with the PLANCK constraints. A large fraction of the inflaton quanta produced after inflation must be converted into the SM radiation reheating the Universe and the rest remains dark matter today. We consider thermal and non-thermal production of inflaton dark matter. In the non-thermal case, we take into account collective effects with the help of lattice simulations. Combining analytic and numerical results with the unitarity consideration, we find that the inflaton dark matter model is viable only in the thermal case where the inflaton mass is near half the Higgs mass., Comment: 6 pages, 2 figures, to be published in the Proceedings of DISCRETE 2020-2021: 7th Symposium on Prospects in the Physics of Discrete Symmetries, 29th November - 3rd December 2021, Bergen, Norway

We consider dark matter production during the inflaton oscillation epoch. It is conceivable that renormalizable interactions between dark matter and inflaton may be negligible. In this case, the leading role is played by higher dimensional operators generated by gravity and thus suppressed by the Planck scale. We focus on dim-6 operators and study the corresponding particle production in perturbative and non-perturbative regimes. We find that the dark matter production rate is dominated by non-derivative operators involving higher powers of the inflaton field. Even if they appear with small Wilson coefficients, such operators can readily account for the correct dark matter abundance., Comment: 16 pages, 2 figures

We study the possibility that, after inflation, the inflaton reaches thermal equilibrium with the Standard Model thermal bath and eventually freezes-out in the non-relativistic regime. When the inflaton decay is the sole source of (non-thermal) dark matter, its relic density is automatically suppressed. We delineate parameter space leading to the correct dark matter abundance. The model allows for a significant Higgs-inflaton coupling which may lead to invisible Higgs decay into inflaton pairs at the LHC., Comment: 14 pages, 3 figures; v2: matches version published in PRD (extended discussions, fixed typos, updated references)

We study scalar dark matter production and reheating via renormalizable inflaton couplings, which include both quartic and trilinear interactions. These processes often depend crucially on collective effects such as resonances, backreaction and rescattering of the produced particles. To take them into account, we perform lattice simulations and map out parameter space producing the correct (non-thermal) dark matter density. We find that the inflaton-dark matter system can reach a quasi-equilibrium state during preheating already at very small couplings, in which case the dark matter abundance becomes independent of the inflaton-dark matter coupling and is described by a universal formula. Dark matter is readily overproduced and even tiny values of the direct inflaton couplings can be sufficient to get the right composition of the Universe, which reaffirms their importance in cosmology., Comment: 21 pages, 7 figures; v2: version accepted for publication in JCAP; v3: fixed a typo in eq. 3.45, added a title for ref. 39

We examine an intriguing possibility that a single field is responsible for both inflation and dark matter, focussing on the minimal set-up where inflation is driven by a scalar coupling to curvature. We study in detail the reheating process in this framework, which amounts mainly to particle production in a quartic potential, and distinguish thermal and non-thermal dark matter options. In the non-thermal case, the reheating is impeded by backreaction and rescattering, making this possibility unrealistic. On the other hand, thermalized dark matter is viable, yet the unitarity bound forces the inflaton mass into a narrow window close to half the Higgs mass., Comment: 16 pages, 4 figures

We show that the weakly interacting massive particle (WIMP) scenario of proton-philic spin-dependent inelastic dark matter can still provide a viable explanation of the observed DAMA effect in compliance with the constraints from other experiments. We also show that, although the COSINE-100 collaboration has recently tested the DAMA effect using the same target material, for the time being the comparison between DAMA and COSINE-100 still depends on the particle-physics model., Comment: 4 pages, 1 figure, contribution to the proceedings of the 16th TAUP conference, Sept. 9-13 2019, Toyama, Japan

Assuming a standard Maxwellian for the WIMP velocity distribution, we obtain the bounds from null WIMP search results of 59.5 days of COSINE-100 data on the DAMA/LIBRA-phase2 modulation effect within the context of the non-relativistic effective theory of WIMP-nucleus scattering. Here, we systematically assume that one of the effective operators allowed by Galilean invariance dominates in the effective Hamiltonian of a spin-1/2 dark matter (DM) particle. We find that, although DAMA/LIBRA and COSINE-100 use the same sodium-iodide target, the comparison of the two results still depends on the particle-physics model. This is mainly due to two reasons: i) the WIMP signal spectral shape; ii) the expected modulation fractions, when the upper bound on the time-averaged rate in COSINE-100 is converted into a constraint on the annual modulation component in DAMA/LIBRA. We find that the latter effect is the dominant one. For several effective operators the expected modulation fractions are larger than in the standard spin-independent or spin-dependent interaction cases. As a consequence, compatibility between the modulation effect observed in DAMA/LIBRA and the null result from COSINE-100 is still possible for several non-relativistic operators. At low WIMP masses such relatively high values of the modulation fractions arise because COSINE-100 is mainly sensitive to WIMP-sodium scattering events, due to the higher threshold compared to DAMA/LIBRA. A next COSINE analysis is expected to have a full sensitivity for the 5$\sigma$ region of DAMA/LIBRA., Comment: 22 pages, 6 figures, 5 tables. Updated to published version

Guided by the non-relativistic effective field theory of interactions between Weakly Interacting Massive Particles (WIMPs) of spin 1/2 and nuclei we study direct detection exclusion plots for an example of non-standard spin-dependent interaction and compare it to the standard one. We analyze an extensive list of 15 existing experiments including the effects of momentum dependence and isospin violation. In our analysis, we fixed the dark matter velocity distribution to a Maxwellian., Comment: Proceeding for ICHEP 2018, International Conference on High Energy Physics, 4-11 July 2018, Seoul, Korea

Since the DAMA/LIBRA collaboration released updated results from their search for the annual modulation signal expected from Dark Matter (DM) scattering in their NaI detectors, we have fitted the updated DAMA result for the modulation amplitudes for a Weakly Interacting Massive Particle (WIMP) signal, parameterizing the interaction with nuclei in terms of the most general effective Lagrangian for a WIMP particle spin 1/2, systematically assuming dominance of one of the 14 possible interaction terms and a standard Maxwellian for the WIMP velocity distribution. We find that most of the couplings of the non-relativistic effective Hamiltonian can provide a better fit compared to the standard Spin Independent interaction case, and with a reduced fine-tuning of the three parameters (WIMP mass, WIMP-nucleon effective cross-section and ratio between the WIMP-neutron and the WIMP-proton couplings). In addition, effective models for which the cross section depends explicitly on the WIMP incoming velocity can provide a better fit of the DAMA data at large values of $m_{\chi}$ compared to the standard velocity-independent cross-section due to a different phase of the modulation amplitudes. All the best fit solutions are in tension with exclusion plots of both XENON1T and PICO60., Comment: 4 pages, 2 figures, Contribution to the Proceedings of the 39th International Conference on High Energy Physics, ICHEP 2018, Seoul, Korea

We show that the Weakly Interacting Massive Particle scenario of proton-philic spin-dependent inelastic Dark Matter (pSIDM) can still provide a viable explanation of the observed DAMA modulation amplitude in compliance with the constraints from other experiments after the release of the DAMA/LIBRA-phase2 data and including the recent bound from COSINE-100, that uses the same $NaI$ target of DAMA. The pSIDM scenario provided a viable explanation of DAMA/LIBRA--phase1 both for a Maxwellian WIMP velocity distribution and in a halo-independent approach. At variance with DAMA/LIBRA-phase1, for which the modulation amplitudes showed an isolated maximum at low energy, the DAMA/LIBRA-phase2 spectrum is compatible to a monotonically decreasing one. Moreover, due to its lower threshold, it is sensitive to WIMP-iodine interactions at low WIMP masses. Due to the combination of these two effects pSIDM can now explain the yearly modulation observed by DAMA/LIBRA only when the WIMP velocity distribution departs from a standard Maxwellian. In this case the WIMP mass $m_{\chi}$ and mass splitting $\delta$ fall in the approximate ranges 7 GeV $\lesssim m_{\chi}\lesssim$ 17 GeV and 18 keV$\lesssim\delta\lesssim$29 keV. The recent COSINE-100 bound is naturally evaded in the pSDIM scenario due to its large expected modulation fractions., Comment: 8 pages, 4 figures. COSINE-100 bound added. Discussion extended. Updated to published version

Assuming for Weakly Interacting Massive Particles (WIMPs) a Maxwellian velocity distribution in the Galaxy we provide an assessment of the sensitivity of existing Dark Matter (DM) direct detection (DD) experiments to operators up to dimension 7 of the relativistic effective field theory describing dark matter interactions with quarks and gluons . In particular we focus on a systematic approach, including an extensive set of experiments and large number of couplings, both exceeding for completeness similar analyses in the literature. The relativistic effective theory requires to fix one coupling for each quark flavor, so in principle for each different combination the bounds should be recalculated starting from direct detection experimental data. To address this problem we propose an approximate model-independent procedure that allows to directly calculate the bounds for any combination of couplings in terms of model-independent limits on the Wilson coefficients of the non-relativistic theory expressed in terms of the WIMP mass and of the neutron-to-proton coupling ratio $c^n/c^p$. We test the result of the approximate procedure against that of a full calculation, and discuss its possible pitfalls and limitations. We also provide a simple interpolating interface in Python that allows to apply our method quantitatively., Comment: 23 pages, 11 figures, 2 tables. Updated to published version. The companion code can be downloaded from https://github.com/NRDD-constraints/NRDD

We re-examine the case of anapole dark matter as an explanation for the DAMA annual modulation in light of the DAMA/LIBRA-phase2 results and improved upper limits from other DM searches. If the WIMP velocity distribution is assumed to be a Maxwellian, anapole dark matter is unable to provide an explanation of the DAMA modulation compatible with the other searches. Nevertheless, anapole dark matter provides a better fit to the DAMA-phase2 modulation data than an isoscalar spin-independent interaction, due to its magnetic coupling with sodium targets. A halo-independent analysis shows that explaining the DAMA modulation above 2 keVee in terms of anapole dark matter is basically impossible in face of the other null results, while the DAMA/LIBRA-phase2 modulation measurements below 2 keVee are marginally allowed. We conclude that in light of current measurements, anapole dark matter does not seem to be a viable explanation for the totality of the DAMA modulation., Comment: 16 pages, 6 figures. Updated to published version

Assuming for Weakly Interacting Massive Particles (WIMPs) a Maxwellian velocity distribution in the Galaxy we explore in a systematic way the relative sensitivity of an extensive set of existing and projected Dark Matter (DM) direct detection experiments to each of the 14 couplings that parameterize the most general non-relativistic (NR) effective Hamiltonian allowed by Galilean invariance for a contact interaction driving the elastic scattering off nuclei of WIMPs of spin 1/2. We perform our analysis in terms of two free parameters: the WIMP mass $m_{\chi}$ and the ratio between the WIMP-neutron and the WIMP-proton couplings $c^n/c^p$. We include the modified signal spectral shape due to non-standard interactions when it is needed in the determination of the bound, such as in the case of background subtraction or of the application of the optimal-interval method. For each coupling, in the $m_{\chi}$-$c^n/c^p$ plane we provide contour plots of the most stringent 90 % C.L. bound on the WIMP-nucleon cross section and show the experiment providing it. We also introduce NRDD_constraints, a simple interpolating code written in Python that allows to obtain the numerical value of the bound as a function of the WIMP mass $m_{\chi}$ and of the coupling ratio $c^n/c^p$ for each NR coupling. We find that 9 experiments out of the 14 present Dark Matter searches considered in our analysis provide the most stringent bound on some of the effective couplings for a given choice of $(m_{\chi},c^n/c^p)$: this is evidence of the complementarity of different target nuclei and/or different combinations of count-rates and energy thresholds when the search of DM is extended to a wide range of possible interactions., Comment: 40 pages, 20 figures, 5 tables. XENON1T bound updated. Python code to access numerical values in tables introduced. Section added with discussion on relativistic models. Updated to published version

The DAMA/LIBRA collaboration has recently released updated results from their search for the annual modulation signal expected from Dark Matter (DM) scattering in their NaI detectors. We have fitted the updated DAMA result for the modulation amplitudes in terms of a Weakly Interacting Massive Particle (WIMP) signal, parameterizing the interaction with nuclei in terms of the most general effective Lagrangian for a WIMP particle spin up to 1/2, systematically assuming dominance of one of the 14 possible interaction terms, and assuming for the WIMP velocity distribution a standard Maxwellian. We find that most of the couplings of the non-relativistic effective Hamiltonian can provide a better fit compared to the standard Spin Independent interaction case, and with a reduced fine-tuning of the three parameters (WIMP mass, WIMP-nucleon effective cross-section and ratio between the WIMP-neutron and the WIMP-proton couplings). Moreover, effective models for which the cross section depends explicitly on the WIMP incoming velocity can provide a better fit of the DAMA data at large values of $m_{\chi}$ compared to the standard velocity-independent cross-section due to a different phase of the modulation amplitudes. All the best fit solutions are in tension with exclusion plots of both XENON1T and PICO60., Comment: 17 pages, 5 figures, 1 table. One high-energy control bin added to the fit, results of Table 1 updated accordingly, no change in the conclusions with the exception of model c5. Some references added and more details provided in the discussion. Updated to published version

Proton-philic spin-dependent inelastic Dark Matter is a Weakly Inelastic Massive Particle scenario that can explain the DAMA yearly modulation effect in compliance with the constraints from other direct detection Dark Matter searches. We obtain updated ranges for its parameters both in a halo-independent approach and adopting a truncated Maxwellian for the WIMP velocity distribution constructing approximate frequentist confidence intervals from an effective chi-square including, among others, the latest experimental constraints from XENON1T, PANDAX-II, SuperCDMS and CDMSlite. In the halo-independent analysis we have implemented the dependence on the WIMP velocity distribution through a step-wise parameterization of the halo function. For the calculation of the profile likelihood we used a Markov Chain Montecarlo generator that required a large number of evaluations of the expected rates. To reduce computational time we introduced expected rates expressions suitable for a fast evaluation through tabulation and interpolation. Our frequentist analysis confirms the present viability of the pSIDM scenario as a possible explanation of the DAMA effect. In terms of the WIMP mass $m_{\chi}$ and of the mass splitting $\delta$ we find the 1-sigma ranges $\mbox{12.5 GeV} \le m_{\chi}\le \mbox{15.7 GeV}$, $\mbox{22.1 keV} \le \delta \le \mbox{26.1 keV}$ for the halo--independent analysis and $\mbox{11.4 GeV} \le m_{\chi}\le \mbox{13.6 GeV}$, $\mbox{24.4 keV} \le \delta \le \mbox{27.0 keV}$ for the Maxwellian case. We find that a full year of data taking of XENON1T should allow to start probing the pSIDM scenario., Comment: 22 pages, 8 figures

We show that compatibility between the DAMA modulation result (as well as less statistically significant excesses such as the CDMS Silicon effect and the excess claimed by CRESST) with constraints from other experiments can be achieved by extending the analysis of direct detection data beyond the standard elastic scattering of a WIMP off nuclei with a spin--dependent or a spin--independent cross section and with a velocity distribution as predicted by the Isothermal Sphere model. To do so we discuss several new approaches for the analysis of Dark Matter direct detection data, with the goal to remove or reduce its dependence on specific theoretical assumptions, and to extend its scope: the factorization approach of astrophysics uncertainties, the classification and study of WIMP-nucleon interactions within non--relativistic field theory, inelastic scattering and isovector-coupling cancellations including subdominant two-nucleon NLO effects. Typically, combining two or more of these ingredients can lead to conclusions which are very different to what usually claimed in the literature. This shows that we are only starting now to scratch the surface of the most general WIMP direct detection parameter space., Comment: 6 pages, 3 figures, to appear in the proceedings of the 14th Marcel Grossmann meeting, Roma, July 2015

Guided by non-relativistic Effective Field Theory (EFT) we classify the most general spin-dependent interactions between a fermionic Weakly Interacting Massive Particle (WIMP) and nuclei, and within this class of models we discuss the viability of an interpretation of the DAMA modulation result in terms of a signal from WIMP elastic scatterings using a halo-independent approach. We find that, although several relativistic EFT's can lead to a spin-dependent cross section, in some cases with an explicit, non-negligible dependence on the WIMP incoming velocity, three main scenarios can be singled out in the non-relativistic limit which approximately encompass them all, and that only differ by their dependence on the transferred momentum. For two of them compatibility between DAMA and other constraints is possible for a WIMP mass below 30 GeV, but only for a WIMP velocity distribution in the halo of our Galaxy which departs from a Maxwellian. This is achieved by combining a suppression of the WIMP effective coupling to neutrons (to evade constraints from xenon and germanium detectors) to an explicit quadratic or quartic dependence of the cross section on the transferred momentum (that leads to a relative enhancement of the expected rate off sodium in DAMA compared to that off fluorine in droplet detectors and bubble chambers). For larger WIMP masses the same scenarios are excluded by scatterings off iodine in COUPP., Comment: 24 pages, 4 figures

We discuss within the framework of effective four-fermion scalar interaction the phenomenology of a Weakly Interacting Massive Particle (WIMP) Dirac Dark Matter candidate which is exothermic (i.e. is metastable and interacts with nuclear targets down-scattering to a lower-mass state) and $Ge$-phobic (i.e. whose couplings to quarks violate isospin symmetry leading to a suppression of its cross section off Germanium targets). We discuss the specific example of the CDMS-II Silicon three-candidate effect showing that a region of the parameter space of the model exists where WIMP scatterings can explain the excess in compliance with other experimental constraints, while at the same time the Dark Matter particle can have a thermal relic density compatible with observation. In this scenario the metastable state $\chi$ and the lowest-mass one $\chi^{\prime}$ have approximately the same density in the present Universe and in our Galaxy, but direct detection experiments are only sensitive to the down-scatters of $\chi$ to $\chi^{\prime}$. We include a discussion of the recently calculated Next-to-Leading Order corrections to Dark Matter-nucleus scattering, showing that their impact on the phenomenology is typically small, but can become sizable in the same parameter space where the thermal relic density is compatible to observation., Comment: 19 pages, 7 figures. References added, updated to published version