UltraViolet freeze-in

If dark matter is thermally decoupled from the visible sector, the observed relic density can potentially be obtained via freeze-in production of dark matter. Typically in such models it is assumed that the dark matter is connected to the thermal bath through feeble renormalisable interactions. Here, rather, we consider the case in which the hidden and visible sectors are coupled only via non-renormalisable operators. This is arguably a more generic realisation of the dark matter freeze-in scenario, as it does not require the introduction of diminutive renormalisable couplings. We examine general aspects of freeze-in via non-renormalisable operators in a number of toy models and present several motivated implementations in the context of Beyond the Standard Model (BSM) physics. Specifically, we study models related to the Peccei-Quinn mechanism and Z′ portals.


Introduction
In the freeze-out paradigm of dark matter (DM), such as the much studied WIMP scenario, the DM is initially in thermal equilibrium and its abundance evolves with its equilibrium distribution until it decouples from the thermal bath. After decoupling the DM comoving number density is constant and (for appropriate parameter values) can give the observed relic density. Models of freeze-in DM  provide a very different picture of the evolution of the DM abundance. In this setting it is supposed that the DM number density is initially negligible but over time an abundance suitable to match the relic density is produced due to interactions in the thermal bath involving a suppressed portal operator. This is illustrated in figure 1. For the DM abundance to be initially negligible, and subsequently set by the freeze-in mechanism (rather than freeze-out), the hidden sector must be thermally decoupled from the visible sector bath at all times, which implies that the portal operators must be extremely small. For instance, for TeV scale DM produced via 2 → 2 scattering of bath states involving a renormalisable portal interaction the coupling dressing this operator should be typically 10 −7 in order to avoid equilibration of the DM with the visible sector [2]. Such DM states are sometimes referred to as feebly interacting massive particles, or FIMPs. Freeze-in, as a general mechanism for DM production, was proposed only recently [1] 1 and thus many important aspects remain to be studied. In particular, a huge class of models has been largely neglected and the purpose of this paper is to rectify this. Freeze-in using renormalisable interactions has been considered in some detail [1][2][3][4][5][6][7][8][9][10]; here instead we examine the alternative possibility, that freeze-in production proceeds via non-renormalisable operators. A suitable DM abundance can potentially be generated by freeze-in via such effective operators, which we refer to as UltraViolet (UV) freeze-in, and in this case the DM abundance depends sensitively on the reheat temperature. Conversely, we use In-fraRed (IR) freeze-in to refer to the class of models in which the sectors are connected via renormalisable operators, in which case the DM abundance is set by IR physics and is independent of the reheat temperature. The different thermal histories associated to these DM frameworks are illustrated in figure 1.

JHEP03(2015)048
The two basic premises of the general freeze-in picture are that • The hidden and visible sectors are thermally disconnected, • The inflaton decays preferentially to the visible sector, not reheating the hidden sector.
Consequently, it is a model independent statement that, due to the out-of-equilibrium dynamics, DM production will proceed through freeze-in via any non-renormalisable operator which is not forbidden by symmetries. Further, the expectation from UV completions of the SM is that distinct sectors of the low energy theory are generically connected by UV physics.
On the other hand, IR freeze-in relies on a rather special construction in which the (renormalisable) portal operators have diminutive couplings, however the naïve expectation is that dimensionless parameters should be near unity. Whilst such feeble couplings are not JHEP03(2015)048 inconceivable (the electron Yukawa Y e ∼ 10 −6 is one example), such decoupling is readily achieved if the visible sector and hidden sector are only connected via high dimension operators.
The UV freeze-in scenario bears some resemblance to models of non-thermal DM [27,28]. The two frameworks both require the DM to be thermally decoupled from the visible sector, they also rely on particular realisations of inflation, and both lead to a DM abundance which is dependent on the reheat temperature. However, there are also important distinctions between these frameworks. In non-thermal DM, the DM has sufficiently small couplings with the visible sector such that energy exchange between the sectors is negligible and the DM relic density is set primarily by inflaton decay. In contrast, in UV freeze-in the DM is dominantly populated by energy transfer from the visible sector to the hidden sector.
In this paper we examine a range of motivated operators for UV freeze-in and discuss potential connections with other aspects of high scale physics. The paper is structured as follows: in section 2 we develop the physics behind UV freeze-in using a number of toy models which exemplify several interesting features. In particular, we discuss high dimension operators with many body final states. Further, we examine examples in which a field involved in the portal operator develops a vacuum expectation value (VEV). We consider the constraints which arise from avoiding sector equilibration in section 3 and find that this leads to bounds on the parameter space, but that large classes of viable models can be constructed. Subsequently, in section 4 we propose a number of simple models, motivated by beyond the Standard Model (BSM) physics, which realise the UV freeze-in picture. Specifically, we consider possible connections with axion models and Z portals. We also comment on the prospect of identifying the scale of UV physics given the DM mass, portal operator and the magnitude of the reheat temperature. In section 5 we provide a brief summary, alongside our closing remarks.

General possibilities for UV freeze-in
The possibility of freeze-in via non-renomalisable operators has been briefly discussed in [1,[11][12][13][14]26]. One of the distinguishing features of UV freeze-in is that DM production is dominated by high temperatures, and so the abundance is sensitive to the reheat temperature T RH . Whilst this possibility has been previously remarked upon as less aesthetic due to the dependence on the unknown value of T RH , it is nevertheless very well motivated as it is a generic expectation that sectors which are decoupled at low energy may communicate via high dimension operators. In this section we examine some general classes of toy models in which the hidden and visible sector are connected only by effective contact operators.

Dimension-five operators with two and three-body final states
We shall start by discussing the simple toy model of UV freeze-in outlined by Hall, March-Russell & West [1] (see also [25] for a similar model in the context of supersymmetry); this will provide a basis from which to examine more realistic scenarios in subsequent JHEP03(2015)048 sections. In this toy model a scalar DM state ϕ freezes-in due to a dimension five operator of the form 2 where φ is a boson in the thermal bath, ψ i are bath fermions, and Λ is the mass scale at which the effective operator is generated. Throughout this section we shall use ϕ and χ to denote scalar and fermion hidden sector states, respectively, and use φ and ψ to indicate scalars and fermions in the thermal bath. Let us suppose, for the time being, that φ does not develop a VEV. An abundance of ϕ can freeze-in via 2 → 2 scattering processes involving the bath states: The change in number density n can be described by the Boltzmann equation (see e.g. [29])ṅ 2E i and f i is the distribution function for a given state. We shall assume that the various states are in thermal equilibrium and thus Maxwell-Boltzmann distributed, f i ∼ e −E i /T i , with the visible sector states φ, ψ i distributed with respect to the temperature of the thermal bath T , whereas the DM ϕ is part of a cold hidden sector at initial temperature T ϕ 0. Correspondingly, this implies that f ϕ 0, and the initial number density of ϕ is negligible where g ϕ is the number of internal degrees of freedom of ϕ. Therefore, the latter term in the Boltzmann equation proportional to f ϕ (the back-reaction) can be neglected. This is the standard picture of the freeze-in scenario, which we shall adopt throughout. Further, we assume here that the portal operator is always sufficiently feeble that it does not bring the hidden sector into thermal equilibrium with the visible sector. We shall examine the specific requirement of this condition in section 3. It follows that the Boltzmann equation can be rewritten as an integral with respect to centre of mass energy as follows [1,30,31] where K 1 denotes a Bessel function of the second kind and

JHEP03(2015)048
In the limit that the particle masses involved in the scattering are negligible compared to the temperature, scattering via the dimension five operator 1 Λ ϕψ 1 ψ 2 φ is described by a matrix element the form where √ s is the centre of mass energy of the scattering at temperature T . Unless otherwise stated, throughout this paper we assume that the masses of the various states are substantially smaller than both Λ and the reheat temperature T RH .
It follows from eq. (2.4) and (2.6) that the Boltzmann equation reduces to the form [1,30,31] (2.7) Using the relationṪ = −HT [29], this can be re-expressing in terms of the yield Y ≡ n S (where S is the entropy density) to obtain in terms of the effective number of degrees of freedom in the bath g S,ρ * [29]. Using the definitions S = 2π 2 g S * T 3 45 , for M Pl the (non-reduced) Planck mass, then integrating with respect to temperature (between T = 0 and T = T RH ) gives [1] Y UV 180 The important thing to note is that the yield depends on the reheat temperature of the visible sector. This is in contrast to the case of freeze-in via renormalisable interactions, in which the yield only depends on the coupling λ and particle masses [1]. As we reproduce in appendix A, the DM yield due to IR freeze-in is parametrically (2.10) With the above example in mind, we extend this analysis to consider a range of effective operators of varying mass dimension and involving different combinations of fields. It is important to recognise that operators of large mass dimension typically lead to many-body final states. Indeed, as we examine below, even the simplest extension of the previous example to dimension five operators 1 Λ φ 1 φ 2 φ 3 φ 4 ϕ involving four bath scalars φ i and scalar DM ϕ, which allows freeze-in production via scattering φ 1 φ 2 → φ 3 φ 4 ϕ, leads to a 3-body phase space.
The Boltzmann equation describing DM production via 2 → 3 scattering is given bẏ

JHEP03(2015)048
where DLIPS 3 denotes the differential Lorentz invariant phase space for 3-body final states and the numerical prefactor accounts for permutations of initial and final states. An evaluation of the 3-body phase space (see appendix B) allows the Boltzmann equation to be rewritten in a form reminiscent of eq. (2.7) We have assumed here that the final state masses can be neglected. By dimensional analysis the associated matrix element is parametrically Substituting this into the Boltzmann equation and expressing our result in terms of the yield we obtain (2.14) Performing the integrals over s and, subsequently, temperature we find the form of the DM yield Y ϕ 135 Up to a numerical suppression of ∼ 10 −2 , this is similar in form to eq. (2.9) and, notably, also depends linearly on the reheat temperature. The DM yield may be related to the relic density as follows where ρ c denotes the critical density and S 0 is the present day entropy density, evaluated at T 0 = 2.75 K ∼ 10 −4 eV. In the latter expression we have approximated g ρ * g S * ∼ 10 3 . We can choose judicious parameter values such that the observed relic density (Ω ϕ h 2 ≈ 0.1) is obtained for a given value of the DM mass. For example, choosing a canonical DM mass of 1 TeV, eq. (2.15) can be rewritten It should be noted that in non-minimal models the DM produced via freeze-in might be able to subsequently annihilate. This would introduce further terms in the Boltzmann equation. If there are additional hidden sector interactions, and light hidden sector states into which the DM can annihilate, then this can give rise to a period of annihilation and freeze-out in the hidden sector, which may dilute the DM relic density. To maintain a degree of predictability, throughout we shall assume that there are no such additional hidden sector interactions which can lead to DM pair annihilation. On the other hand, if the DM only interacts via the suppressed portal operator then, as the DM never enters thermal equilibrium, the rate of annihilation back to the visible sector is always negligible compared to the rate of production. In some sense the DM is immediately frozen-out on production.

High dimension operators with many-body final states
We would like to understand how this generalises to operators of increasing mass dimension. For UV freeze-in involving an operator of mass dimension n + 4, the cut-off enters in the denominator of the yield as Λ 2n . Thus, still assuming that none of the fields involved acquire non-zero VEVs, the generic expectation is that the DM yield Y (n) due to this operator should scale as follows the factor of M Pl coming from the Hubble parameter. One important issue which arises, however, is that the phase space becomes increasingly large and complicated. Here we shall make certain assumptions and approximations to obtain an order-of-magnitude estimate of the yield. Consider the dimension-(n+4) operator 1 Λ n φ 1 φ 2 · · · φ n+3 ϕ which corresponds to (n+2)body final state phase space for scattering events φ 1 φ 2 → φ 3 · · · φ n+3 ϕ. Variant operators might be considered but this toy example should be illustrative of a more general issue regarding the interplay between mass dimension and phase space. For simplicity here we assume that there is only a single relevant high dimension operator; the converse scenario would imply a sum over the various operators in the Boltzmann equation.
The Boltzmann equation describing DM production via 2 → n + 2 scattering is given by 4ṅ The differential phase space grows like DLIPS (n+2) ∝ s n and we make the approximation (2.20) The square bracket provides a parametric estimate of the additional phase space suppression due to the (n + 2)-body final state. This is a somewhat crude approximation, but for low n should give an order-of-magnitude estimate. For increasing n the suppression to the cross section coming from the phase space should become more severe and strongly suppress operators of high (n 1) mass dimension. From dimensional analysis and thus the Boltzmann equation can be expresseḋ We will check the resulting estimate against the 3-body result calculated in section 2.1.

JHEP03(2015)048
The integral over s has a closed form expression which for n ∈ N is given by Using this result the Boltzmann equation can be rewritteṅ From the above we obtain an expression for the DM yield (2.25) The form of eq. (2.25) conforms with our expectations for the parametric scaling discussed in eq. (2.18). Moreover, this shows that a range of operators of varying mass dimension should be able to reproduce the observed relic density. We consider some examples below. Firstly, we can check this result against our 3-body calculation by setting n = 1 Comparing with eq. (2.15), we find that these two expressions agree up to an O(1) factor. For suitable parameter values the observed relic density can be reproduced for a large range of DM masses, comparing with eq. (2.16). As a concrete example, consider a model with 1 TeV DM, a yield of appropriate magnitude is found for (2.27) Observe, unlike typical models of freeze-out and IR freeze-in, the yield is independent of m ϕ , the DM mass, provided m ϕ T Λ. Thus in UV freeze-in one can find the observed DM relic density for different values of the DM mass by simply rescaling the yield.
Taking a few further examples, consider dimension-six (n = 2) and dimension-seven (n = 3) operators (with 4 and 5 body final states, respectively), the estimates for the freeze-in yield of ϕ for these cases are T RH 10 10 GeV T RH 8 × 10 10 GeV Thus for a given operator and fixed UV scale, the observed DM abundance can typically be obtained by adjusting the reheat temperature. We shall comment in section 4.3 on motivated values of T RH . Moreover, we learn that in the presence of multiple high dimension portal operators of varying mass dimension (but common Λ), generally the physics will be determined by the operator(s) with smallest mass dimension. This conforms with the standard intuition regarding effective field theory.

VEV expansions of high dimension operators
If any of the fields involved in the high dimension operator acquire VEVs then, by expanding the operator around the vacuum, one can construct a sequence of terms dressed by couplings of different mass dimension [1]. Let us again examine the simple example given in [1], involving the dimension-five operator in eq. (2.1). Suppose that the bath scalar field has a non-zero VEV φ = 0, expanding around this VEV gives To ensure validity of the effective field theory we will assume that φ Λ and therefore λ 1. In this example the yield will receive both an IR contribution from the first term (which appears as an operator with a dimensionless coupling after symmetry breaking) and a UV contribution from the latter term. The IR contribution is assumed to be generate by decays ψ 1 → ψ 2 ϕ, this is reproduced in eq. (A.6) of the appendix. By calculating the two contributions, it can be shown [1] that the yield is dominated by the IR contribution if It is notable that this condition does not depend on Λ. Let T * be the critical temperature associated to the spontaneous breaking of some symmetry, due to a scalar field involved in the UV freeze-in operator developing a VEV. It should be expected that T * ∼ φ . The VEV expansion of eq. (2.29) is only valid for T < T * , but as the yield due to IR operators is temperature independent, it will not depend on the temperature T * at which the IR operator is generated. Thus it is not required that T RH < T * in this case.
The situation is more complicated if the VEV expansion leads to additional UV freezein operators. For T * < T RH , one can find UV contributions which depend on T * , rather than T RH . This is because for temperatures T > T * there is no VEV expansion until thermal evolution (due to expansion) causes T to drop below T * . We shall illustrate this with an example below. Importantly, if the phase transition happens after the point at which DM freeze-in terminates, i.e. below the mass of the bath states involved in the freeze-in process (T * < m bath ), then no further (UV or IR) contributions will be generated.
Let us consider an example in which the VEV insertion does not lead to operators with dimensionless coefficients, but results in several UV contributions. One manner of realising this scenario is by dressing the SM Yukawa operators with a pair of DM states Note that the DM can be stabilised by the assumption of χ-parity, such that the Lagrangian is invariant under χ → −χ. After electroweak symmetry breaking (EWSB) one makes the appropriate expansion around the vacuum to obtain the four fermion operator. Henceforth we shall call O 1 the operator dressed by v/Λ 3 and the latter term O 2 .

JHEP03(2015)048
As one might anticipate the VEV expansion leads to two operators which provide UV contributions to the yield. The operator O 2 leads to freeze-in via 2 → 3 scattering processes such as hq → qχχ. From dimensional analysis the matrix element is of the form As previously, we can describe the production of DM through the Boltzmann equation given in eq. (2.12). It follows that the DM yield due to O 2 is of the form Similarly, the operator O 1 results in DM production via 2 → 2 scatteringqq →χχ The relevant Boltzmann equation is analogous to eq. (2.7) Hence the yield is given by The maximum temperature T max , which is the upper limit of the integral, depends on whether the reheating temperature is above or below the critical temperature T * at which the Higgs develops a VEV and can be expressed as T max = min(T * , T RH ). Thus there are two possible cases, which we examine below, depending on whether the reheat temperature is above or below the phase transition. As the operator in eq. (2.31) involves the Higgs VEV, the critical temperature is around T * ∼ 100 GeV, associated to the electroweak phase transition (EWPT). Provided T * > T RH the VEV expansion is valid, and thus the portal operator O 1 is active, for all physically relevant temperatures. Therefore T max = T RH , and identifying H = v, the ratio of the two contributions is (2.37) For VEV expansions which do not lead to dimensionless couplings, because More generally, if T * > T RH the expectation is that the contribution coming from the operator in the VEV expansion dressed by the coefficient with smallest (negative) mass dimension will dominate the yield. For T * < T RH the phase transition takes place during the cooling of the thermal bath, in this case production via the VEV expanded operator only occurs for T < T * and the dominant contribution will be generated at T ∼ T * . Further, thermal fluctuations may be JHEP03(2015)048 important in determining the field expectation value and it is expected 5 that H ∼ T . Evaluating eq. (2.36), the ratio of the contributions coming from eq. (2.31) in this case is instead (2.38) Since by assumption v < T RH , the operator O 2 is typically dominant. More generally, we expect that for alternative operators typically the term with no explicit VEVs in the expansion around the vacuum will provide the most significant contribution to the yield. In the case that the VEV expansion generates a dimensionless coupling, the associated yield is independent of T RH and T * , and the criteria under which this operator dominates will be described by an equation analogous to eq. (2.30).
For scattering processes hq → qχχ to occur in the thermal bath it is required that the SM states can be thermally produced. This implies that T RH 100 GeV if H is the SM Higgs. In the case that the Higgs is thermally produced, this pushes the model into the regime in which the operator O 2 always gives the dominate contribution. In the converse scenario that the Higgs is not thermally produced then O 2 will be exponentially suppressed, but O 1 can still potentially lead to DM production.
Further complications can arise if the VEV expansion gives several terms, and thus multiple contributions to the DM yield, or if there are multiple scalar fields, especially if the scalars develop VEVs due to spontaneous symmetry breaking taking place at different scales. We shall leave these more complicated possibilities until they arise in motivated examples.

Sector equilibration constraints
For the DM relic density to be established by freeze-in production (rather than freeze-out) it is imperative that the DM is not brought into thermal equilibrium with the visible sector due to its interactions through the portal operator. For the case of IR freeze-in via renormalisable interactions, for instance due to λψ 1 ψ 2 ϕ, it has been argued [2] that the hidden sector should not thermalise with the visible sector, provided λ 10 −6 m bath /100 GeV. 6 In this section we derive an analogous condition for the case that the portal is due to a non-renormalisable operator, leading to UV freeze-in.
It will be useful to introduce the freeze-out temperature T FO , the temperature at which a state in thermal equilibrium decouples from the thermal bath, defined such that at T = T FO n σv = H . (3.1) The requirement that the DM is always out of thermal equilibrium is equivalent to Y < Y eq . This implies two conditions: 5 A more careful study of these thermal effects would be of interest, but it is beyond the scope of this work. 6 In the case that a VEV expansion gives a renormalisable operator, then identifying λ = φ1 · · · φn /Λ n one should apply this IR constraint, in addition to requirements on high dimension operators in this section.

JHEP03(2015)048
• For T RH > T > m DM , the bath number density is n eq T 3 /π 2 and thus the DM number density is non-equilibrium provided n DM T 3 /π 2 .
• For m DM > T , the equilibrium number density is Boltzmann suppressed and thus to avoid n DM coinciding with n eq it is required that DM freeze-out occurs for T FO m DM .
The UV operator freezes-in a DM abundance at T RH , which is effectively frozen-out. Clearly, an upper bound is given by the scenario in which a near thermal abundance is generated in the early universe, subsequently decouples, and then evolves to the present modified only by entropy conservation. Including the entropy factor, this gives (for a real scalar DM particle) Further, we use that the DM relic density is given by Comparing with the observed value Ω DM ≈ 0.2, this gives a bound on the DM mass where we have used that g(T 0 )/g(T RH ) 3.36/100. Thus this places a model independent bound on the DM mass. Next we examine the second requirement: T FO m DM . The freeze-out temperature can be found by solving eq. (3.1), which we expand below Thus the requirement that freeze-out occurs before the mass threshold m DM is given by For a specific portal operator this can be re-expressed in terms of the UV scale Λ at which the operator is generated. Let us consider an example. Recall the Lagrangian term studied previously for UV freeze-in by 2 → 2 scattering: L ⊃ 1 Λ ϕψ 1 ψ 2 φ. In this case the matrix element is given by eq. (2.6) and the associated cross section is  Moreover, when combined with the constraint of eq. (3.2), this implies an extreme lower bound on the scale of new physics of 10 7 GeV, as indicated above.
In figure 2 we display the constraints on the parameter space of various models. It particular we look at the UV freeze-in portals considered in section 2.2 and 2.3. Using the derived forms of the yields, we present contours of Λ and T RH which give the observed relic density Ω DM h 2 ≈ 0.1 for a range of DM masses. We overlay these contour plots with the relevant constraints. Towards the lower right of the parameter space in each plot the yields are low and require increasingly larger DM masses in order to reproduce the observed relic density. In certain regions of parameter space the DM mass which would be required is larger than the T RH and thus such models can not give the correct DM abundance; this is indicated by the grey shaded regions in figure 2. The purple shaded regions indicate parameter values which lead to sector equilibration, as given in eq. (3.2), the DM enters thermal equilibrium Y = Y eq , thus the abundance will not be set via the freeze-in mechanism. In the red highlighted region T RH > Λ, the effective field theory breaks down, and the dark matter abundance is not set by UV freeze-in. Note also that because the DM is never in thermal equilibrium the unitarity bound of freeze-out DM [32] does not apply to freeze-in DM [1] and thus there is no upper bound on the DM mass.

UV freeze-in and BSM physics
Having discussed a range of possibilities in the context of toy models, we how turn to constructing explicit models based on extensions of the SM, motivated by outstanding problems. These BSM scenarios generally require additional states to be introduced at new physical scales above the weak scale. One interesting possibility is that the operator(s) responsible for UV freeze-in are generate at this scale of new physics. We shall also discuss JHEP03(2015)048 motivated BSM models which lead to high dimension operators with VEV expansions at temperatures above the weak scale. It was remarked in [1] that the high dimension operator responsible for UV freeze-in might arise from GUT scale physics, here we examine some alternative scenarios.

Z portal
First we consider the scenario in which the SM gauge group is extended by an additional U(1) gauge symmetry, which is broken at some high scale Λ (for a general review see e.g. [33]) If some visible sector states and the DM are both charged under this new gauge symmetry, then the associated massive gauge boson can provide a portal that links these two sectors. Additional U(1) gauge groups are a generic expectation of string theory compactifications, see e.g. [34][35][36][37], as supported by scans of vacua of Heterotic string theory [38,39]. Further, GUTs based on E 6 or SO (10) can introduce extra U(1) factors from the breaking of these larger groups [40] In type IIB theories extra U(1)'s can arise from isolated branes; moreover, brane stacks are associated to symmetry groups U(N ) ∼ SU(N ) × U(1), where the U(1) factor is (pseudo)-anomalous [34,35]. This U(1) anomaly is cancelled via the Green-Schwarz mechanism, and as a result the Z acquires a mass near the string scale. It is interesting to note that in type IIB theories the string scale can be lowered substantially compared to the Planck mass if the moduli are stabilised at LARGE volume [41].
Alternatively, from an IR perspective, it is conceivable that a global quantum number of the SM model might be gauged. In the SM baryon number B and lepton number L appear as accidental symmetries and are typically broken in extensions of the SM, for instance GUTs. However, it is possible that some global quantum number of the SM may arise from an exact gauged symmetry. An appealing possibility is that the combination B − L is gauged as this is anomaly free provided the spectrum includes right-handed neutrinos. If one assumes that DM is charged under B − L, then the U(1) B−L gauge boson can provide a portal operator which connects the DM and the SM fermions, see e.g. [42][43][44]. If B − L is gauged, typically it must be broken at a high scale to give masses to the neutrinos via the seesaw mechanism. Once the Z is integrated out this generates effective operators which connect the SM fermions and the DM, suppressed by the (intermediate) scale at which U(1) B−L is broken. This can potentially lead to UV freeze-in for appropriate parameter choices.
More generally, suppose that the DM χ,χ and the SM fermions are charged under some new group U(1) . It follows that the SM states can pair-annihilate and produce DM states viaqq → Z →χχ. In the UV, as usual, interactions mediated by the (heavy) Z appear in the Lagrangian through the covariant derivative in the gauge invariant kinetic terms

JHEP03(2015)048
where / D = / ∂ + iy / Z + · · · , the ellipsis denote the gauge fields of the SM and y is the U(1) charge. Once the Z is integrated out this leads to four-fermion interactions, suppressed by Λ = y q y χ /m Z , in the effective Lagrangian of the form This is similar to O 1 studied in eq. (2.31), however the prefactor is different, as is the Lorentz structure, as here we have integrated out a Lorentz vector. DM production will proceed fairly analogously and the DM yield is given by (cf. eq. (2.25)) It should be noted that the presence of an additional U(1) gauge group will generically lead to kinetic mixing via operators of the form F µν F µν . Such interactions can provide a renormalisable portal operator between the visible and hidden sectors. As we are primarily concerned here with UV freeze-in, we shall assume that such operators are negligible. This kinetic-mixing portal has been previously studied in the context of IR freeze-in [11]. Discussion on UV freeze-in via Z also appear in [13,14].

JHEP03(2015)048
In addition to the Yukawa couplings we can build the following Lorentz, gauge, and PQ invariant combinations of these fields: These field combinations allow us to construct the following Lagrangian terms involving the DM bilinear The scalar field S develops a non-vanishing VEV v a at scale f a = 2v a m Z , which preserves electroweak symmetry, but spontaneously breaks the PQ symmetry.
Let us assume the following mass hierarch: For T RH > f a , UV freeze-in can proceed via scattering SS † → χχ. The matrix element for this process is |M| 2 ∼ s/Λ 2 , which is similar in form to eq. (2.6). Thus the UV contribution to the yield is (up to O(1) factors) as eq. (2.9) There is a similar T RH -dependent contributions to the yield from the operators H † i H i χχ c for (i = u, d) which has the same form to that given above. The correct relic density is found for Λ 10 9 m χ T RH . A further model-dependent limit comes from the assumed mass hierarchy eq. (4.9) m χ Λ 2 10 18 T RH < Λ 2 10 18 m σ . (4.11) As we expect from the Lagrangian that m χ ∼ f 2 a /Λ and m σ ∼ f a , this implies the following consistency constraint on the hierarchy of scales: f a Λ/10 6 . The above requirements can be simultaneously satisfied with reasonable parameter values.
At energies above the EWPT, but after PQ breaking one can re-examine the operators appearing in eq. (4.8) following a VEV expansion around the vacuum of S. The singlet field can be decomposed into radial and axial components The axial field a is identified with the axion. Expanding around the VEV of S, the radial component σ provides the following op- For T m σ the state σ is part of the thermal bath, as it is kept in thermal contact via interactions involving products of the bilinear operators in eq. (4.7) involving H u , H d and S if these have O(1) coefficients.

JHEP03(2015)048
Once the temperature drops below the PQ breaking scale (but whilst still above m σ ), freeze-in can proceed via direct decay of heavy σ states to DM pairs, leading to an IR contribution to the yield. Comparing with the form of the IR freeze-in yield given in [1], and reproduced in eq. (A.6) of appendix A, one finds is the partial width of σ → χχ. The condition under which the UV contribution will be dominant is Note also that after EWSB there is a further VEV expansion involving the Higgs fields which leads to the portal operators of the form This leads to IR and UV freeze-in contributions to the yield, similar to eq. (4.10) and (4.14), in the case that m χ < m Z . It would be of interest to embed the above model into a supersymmetric extension of the SM, similar to e.g. [46]. The advantages of a supersymmetric implementation is that the type II structure of the Higgs sector (required to employ the DFSZ model and to avoid constrained flavour changing processes) is an automatic consequence of holomorphy of the superpotential. In addition this might also alleviate the naturalness problem [49] associated with destabilising the weak scale, and the DM can be stabilised by Z 2 R-parity if it is the lightest supersymmetric particle.

The reheat temperature
In section 4.1 and 4.2 we have discussed motivated scales of new physics which might generate the high dimension operators. An interesting alternative to this approach is to consider special values for the reheat temperature T RH and use this, in conjunction with the DM relic density, to identify the unknown UV scale. One drawback of this scenario is that very little is known about the reheat temperature. Precision measurements of primordial elements due to Big Bang nucleosynthesis are thought to imply that T RH few MeV. Models of inflation typically suggest an upper bound around T RH 10 16 GeV, see e.g. [51]. Moreover, if T RH is high then in principle this can lead to problems with long-lived exotic relics which can over-close the Universe, the classic example being the cosmological gravitino problem of supergravity [52]. This implies the upper bound on the reheat temperature in models of supergravity is typically T RH 10 10 GeV, with substantially stronger bounds if the gravitino is light. 8 8 It has been further argued [53] that, with some assumptions, once the contribution from axinos is also included this can lead to an even more stringent upper bound: TRH 10 5 GeV.

JHEP03(2015)048
The following extreme cases may be of particular interest: • The maximum expected from simple models of inflation: T RH ∼ 10 16 GeV.
• The upper bound from the cosmological gravitino problem: T RH ∼ 10 10 GeV.
• The lower bound from precision measurement of primordial elements: T RH ∼ 10 MeV.
The first two scenarios might be motivated through considerations of environmental selection if there is some anthropic pressure which favours high T RH , with the cosmological gravitino problem imposing a catastrophic boundary at T RH ∼ 10 10 GeV in supersymmetric models.
Let us consider a specific example involving dimension-(n + 4) operator 1 Λ n φ 1 φ 2 · · · φ n+3 ϕ of the scalar toy model studied in section 2.2. For DM with mass around 1 TeV the yield required to obtain the observed relic density is Y ∼ 10 −9 , as discussed in eq. (2.16). To obtain the correct relic abundance via freeze-in through the dimension-five operator with a reheat temperature of T RH ∼ 10 16 GeV, requires a UV scale as indicated below (4.17) The magnitude of the UV scale Λ in this example may be suggestive of a connection with the heterotic string scale [54]. It is not implausible that future observations might indicate the reheat temperature (given some assumptions regarding the model of inflation). Ultimately to test UV freeze-in DM, and disambiguate it from other frameworks, it will be necessary to determine the DM mass, the UV scale Λ and T RH . Recently the BICEP2 collaboration claimed they had observed primordial tensor modes [55], and thus could infer T RH , however this result is currently disputed [56]. If the BICEP2 signal survives further scrutiny we shall comment on this in a dedicated paper.

Conclusion
This work has provided an exploratory study of the model building opportunities which arise for the UV freeze-in mechanism. We considered general aspects of this scenario in the context of toy models and demonstrated that interesting and phenomenologically viable models can be constructed in motivated settings of BSM physics. In section 2, we examined various toy models which encapsulate the fundamental features of UV freeze-in. Typically high dimension operators lead to many body final states, and the case of DM production via 2 → 3 UV freeze-in was carefully studied. Subsequently, we attempted to quantify DM production associated to more complicated phase spaces. Further, we discussed the potential impact of spontaneous symmetry breaking on UV freeze-in, in particular, we identified a new case of interest in which the VEV expansion leads only to additional UV contributions, and does not generate an IR freeze-in portal. In this scenario the DM yield depends on both the reheat temperature and the critical temperature of symmetry breaking. Section 3 examined the constraints on UV freeze-in from the requirement that JHEP03(2015)048 the hidden sector and visible sector do not equilibrate and we argued that this can lead to bounds on the DM mass.
In section 4 we presented realistic models of UV freeze-in and related these to interesting BSM scenarios. We suggested that UV freeze-in might be connected to motivated solutions of prominent puzzles of the SM, specifically we consider an example involving the Peccei-Quinn mechanism. A further example was presented in which the UV freeze-in portal is generated by integrating out a heavy Z . This is appealing as Z arise in many extensions of the SM and additional U(1) gauge groups are common in realistic string compactifications. It should be evident from our discussions that UV freeze-in offers a large range of possibilities for DM model building and that there are many interesting aspects yet to be explored.
UV freeze-in presents a new manner of obtaining non-thermal DM, with a relic abundance directly related to the reheat temperature, and provides an interesting alternative to the conventional ideas regarding the DM thermal history. For the DM relic density to be set through UV freeze-in it is required that reheating of the hidden sector is negligible and that the DM is connected to the visible sector via non-renormalisable operators. Once one assumes that the abundance of DM is initially depleted, one might argue that freeze-in via high dimension contact operators presents a more generic mechanism than IR freeze-in portals, which require very small renormalisable couplings or complicated effective operators involving several scalar fields with non-vanishing VEVs. Moreover, from a UV perspective, it is a fairly general expectation that distinct sectors in the low energy theory may become coupled through the high scale physics, and we have presented some examples of this principle in the above.

JHEP03(2015)048
As previously, we re-express this as an integral with respect to centre of mass energẏ where P ij is defined in eq. (2.5). If the bath state masses can be neglected, this reduces tȯ Converting this to a yield one obtains the result [1] For the case of direct freeze-in due to decays of heavy bath states ψ H to a lighter state ψ L in the thermal bath and DM ϕ via the interaction λψ H ψ L ϕ this is insteaḋ