## Minimal signatures of the Standard Model in non-Gaussianities

## Anson Hook, ∗

Maryland Center for Fundamental Physics, University of Maryland, College Park, MD 20742

Junwu Huang † and Davide Racco ‡ Waterloo, Ontario N2L 2Y5, Canada

Perimeter Institute for Theoretical Physics,

We show that the leading coupling between a shift symmetric inflaton and the Standard Model fermions leads to an induced electroweak symmetry breaking due to particle production during inflation, and as a result, a unique oscillating feature in non-Gaussianities. In this one parameter model, the enhanced production of Standard Model fermions dynamically generates a new electroweak symmetry breaking minimum, where the Higgs field classically rolls into. The production of fermions stops when the Higgs expectation value and hence the fermion masses become too large, suppressing fermion production. The balance between the above-mentioned effects gives the Standard Model fermions masses that are uniquely determined by their couplings to the inflaton. In particular, the heaviest Standard Model fermion, the top quark, can produce a distinct cosmological collider physics signature characterised by a one-to-one relation between amplitude and frequency of the oscillating signal, which is observable at future 21-cm surveys.

## I. INTRODUCTION AND SUMMARY

ton [11, 12]

Cosmological collider physics [1-8] provides an opportunity to search for new heavy particles that are not accessible at particle colliders. These heavy particles, produced through their interactions with the inflaton field, can accumulate to large enough densities and affect the bispectrum of density perturbations, leaving observable signatures in the cosmic microwave background and large scale structure of the universe [9, 10].

Given the exciting potential of this approach, it is worth asking 'What are the minimal signatures of the Standard Model (SM) in the context of cosmological collider physics?'. Of course the absolutely minimal signature is nothing if the SM does not couple directly to the inflaton. In this article, we explore the consequences of adding a single coupling to the SM fermions, the dimension five shift symmetric coupling between the inflaton and the SM fermions. Somewhat surprisingly, this single coupling by itself leads to an interplay between the dynamics of the SM fermions and the Higgs: The SM fermion production induces electroweak symmetry breaking during inflation, whereas a large expectation value of the Higgs field increases the masses of the SM fermions and suppresses their production. This interplay results in a very predictive scenario where the strength of the oscillating signature in non-Gaussianities is directly tied to the period of oscillation.

We assume that the fermions in the Standard Model couple to the inflaton through the lowest dimensional operator respecting the shift symmetry of the infla-

∗ hook@umd.edu

† jhuang@perimeterinstitute.ca

‡ dracco@perimeterinstitute.ca

$$\mathcal { L } \supset \frac { \partial _ { \mu } \phi } { \Lambda _ { f } } \left ( F ^ { \dagger } \overline { \sigma } ^ { \mu } F + f ^ { c \, \dagger } \overline { \sigma } ^ { \mu } f ^ { c } \right ) + y _ { f } \mathcal { H } F f ^ { c }, \quad ( 1 )$$

where φ is the inflaton, H = ( 0 , v + h √ 2 ) is the Higgs doublet, and F = Q, L and f c = u , d , e c c c are leftand right-handed fermions in the SM in two-component notation 1 . When 〈 ˙ φ 〉 glyph[negationslash] = 0, this coupling leads to the production of fermions during inflation whose effective number density is n f ∼ m λ f 2 f exp[ -πm /λ H 2 f f ], where λ f = ˙ φ Λ f , m f = y f v/ √ 2 and v/ √ 2 = 〈H〉 is the vacuum expectation value (vev) of the Higgs field.

It is well known that a high density of particles can change the properties of a scalar potential. Thermal effects are well known to favor symmetry restoration of the Higgs for temperatures above the electroweak phase transition [13]. On the other hand, chemical potentials favor symmetry breaking and in the context of the Higgs potential could prevent symmetry restoration even at high temperatures [14]. The coupling in Eq. (1) is very similar to a chemical potential. In fact, if the plus sign were a minus sign, it would be the familiar chemical potential for fermion number. As such, it is unsurprising to find that the effect of this coupling is to generate a correction to the Higgs potential that favors symmetry breaking of the form

$$\delta V _ { h } = - \frac { y _ { f } ^ { 2 } } { 2 \pi ^ { 2 } } \lambda _ { f } ^ { 2 } h ^ { 2 } \exp \left [ - \frac { \pi y _ { f } ^ { 2 } h ^ { 2 } } { 2 \lambda _ { f } H } \right ]. \quad \ \ ( 2 )$$

For λ f &gt; H glyph[greatermuch] v ew ( v ew being the electroweak vacuum), this contribution to the potential induces spontaneous

1 We neglect the typically smaller dimension five anomalous couplings to the gauge bosons.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00019v1-figure-a8fbd0d1-49d4-4f8a-8f6f-352f753d7d05.jpg)

glyph[negationslash]

FIG. 1. Higgs potential at zero temperature and λ f ∼ 0 (blue line) and in presence of a top quark condensate induced by λ t H &gt; v ew (purple line). The Higgs field sits in the dynamically generated minimum v during inflation.

breaking of the electroweak symmetry, as shown in Fig. 1. However, as the Higgs vev increases, the fermion masses increase and their particle production is exponentially suppressed as m &gt; λ H 2 f ∼ f . Therefore, quite insensitively to any other term of the Higgs potential, e.g. the value of the quartic λ h , the Higgs gets a vev during inflation

$$v \sim \frac { 1 } { y _ { f } } \sqrt { \lambda _ { f } H } \,. \quad \quad \quad ( 3 ) \quad \text{ by } \\ \text{th}.$$

In the SM, due to the large hierarchy between the top quark and the lighter leptons and quarks, this effect is determined entirely by the top quark. Incidentally, such a scenario can only occur for fermions with an O (1) Yukawa coupling, since ˜ λ f = λ /H f cannot be arbitrarily large. Thus, in what follows, we focus on the coupling of the top quark with the inflaton.

The observational signature associated to this coupling is the generation of a large non-Gaussian oscillating pattern in the squeezed limit. The Feynman diagram to be calculated is shown in Fig. 2 [11, 12].

1 FIG. 2. Contribution to the inflaton bispectrum from a loop of SM fermions. Two fermions f are produced at a time τ 3 by the interaction with a soft inflaton leg δφ , and annihilate later at τ 1 ∼ τ 2 into two hard inflaton legs with k , k 1 2 glyph[greatermuch] k 3 . The non-analytic contribution to the bispectrum is due to the time propagation of the fermions from τ 3 to τ 1 , τ 2 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00019v1-figure-a29884a6-e388-4818-b9a1-62cb048f9486.jpg)

The interesting feature in this case is that the dynamics of the Higgs potential ensures that the mass of the fermion f that produces the largest observable non-Gaussianity is solely determined by the inflaton coupling λ f to that fermion. As a result, the non-Gaussian signal in the squeezed limit simplifies into

$$f _ { N L } ^ { ( \text{clock} ) } & \approx \frac { 4 \sqrt { 2 } N _ { e } \mathcal { P } _ { \zeta } } { 3 e } \widetilde { \lambda } _ { f } ^ { 1 3 / 2 }, \quad \quad ( 4 )$$

depending only on the size of the inflaton fermion coupling in Hubble units ˜ λ f = λ f H = ˙ φ Λ f H . This signal oscillates in ln( k /k 3 1 ) with a frequency ∼ ˜ λ f in the squeezed limit. This relation between the amplitude ( f (clock) NL ) and the frequency of the oscillating signal offers a simple cross check of this mechanism. Such a signature offers a direct probe of the induced electroweak symmetry breaking during inflation, and could help us to shed light on the inflationary sector.

The rest of paper is organized as follows. In Sec. II, we discuss the fermion particle production during inflation and calculate the subsequent effect on the Higgs potential. In Sec. III, we study the dynamical equilibrium reached during inflation and how the fermion mass is set by its coupling to the inflaton. In Sec. IV, we discuss the signature in the bispectrum of density perturbations ultimately sourced by the dynamically generated Higgs vev, and the implications for the detection of this signal.

## II. PARTICLE PRODUCTION AND THE HIGGS POTENTIAL

In this section, we will discuss particle production and how the fermion density affects the Higgs potential. To calculate the correction to the Higgs mass term, we calculate the standard diagram shown in Fig. 3 that corrects the energy density of the state.

p

p

←-

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00019v1-figure-89f185f5-62a3-4c44-9f3c-fd6b89517972.jpg)

←-

1 FIG. 3. Feynman diagrams for the contribution of the top fermion condensate to the Higgs potential. For notations, see the main text and references [12, 15, 16].

In a companion paper [12] (see in particular Sec. 3 and Appendix A), we show in detail how to estimate and calculate these diagrams. For a single fermion flavor and color, the diagrams in Fig. 3 contribute as

$$N _ { a b } & = - y _ { f } ^ { 2 } \epsilon ^ { \alpha \beta } \epsilon ^ { \dot { \alpha } \dot { \beta } } \text{ab} \iint _ { - \infty } ^ { 0 } \frac { d \tau _ { 1 } } { H \tau _ { 1 } } \frac { d \tau _ { 2 } } { H \tau _ { 2 } } G _ { a } ( \vec { k } _ { 1 }, \tau _ { 1 } ) G _ { b } ( - \vec { k } _ { 1 }, \tau _ { 2 } ) \\ & \quad \times \int \frac { d ^ { 3 } q } { ( 2 \pi ) ^ { 3 } } \left ( D _ { a b \alpha \dot { \beta } } ( \vec { p } _ { 1 2 }, \tau _ { 1 }, \tau _ { 2 } ) D _ { a b \beta \alpha } ( - \vec { p } _ { 2 1 }, \tau _ { 1 }, \tau _ { 2 } ) + D _ { b \alpha \alpha \dot { \beta } } ( - \vec { p } _ { 1 2 }, \tau _ { 2 }, \tau _ { 1 } ) D _ { b \alpha \dot { \alpha } } ( \vec { p } _ { 2 1 }, \tau _ { 2 }, \tau _ { 1 } ) \right ),$$

where τ denotes conformal time, glyph[vector] p 12 = glyph[vector] k 1 + glyph[vector] p 21 = glyph[vector] q and | glyph[vector] q | glyph[greatermuch] | glyph[vector] k 1 | is the internal momentum. The indices α, β, α, β ˙ ˙ ∈ { 1 2 , } are spinor indices. The antisymmetric glyph[epsilon1] tensors are defined by glyph[epsilon1] 12 = -glyph[epsilon1] 12 = 1. The in-in indices a , b take values in { +1 , - } 1 , denoted respectively by {⊕ glyph[circleminus]} , to distinguish them from spinor indices. The functions G a ( glyph[vector] k, τ ) = √ H / k 2 2 3 (1 -i a kτ ) e i a kτ and D ab αβ ˙ ( glyph[vector] k, τ 1 , τ 2 ) are the propagator of the Higgs and the fermion fields respectively, and can be found in Sec. A.2 of [12]. By the same logic of the Electroweak Hierarchy problem, the momentum integral in Eq. (5) is quadratically divergent. The leading quadratic divergence is absorbed into the definition of the physical Higgs mass that we observe today when λ f = 0 (this operation automatically removes some of the subleading corrections in the ˜ λ f expansion). In the massless limit, the integral in Eq. (5) can be simplified (by exploiting properties of the Whittaker functions appearing in the fermion mode functions) to be

$$N _ { \text{aa} } = y _ { f } ^ { 2 } \int \frac { d ^ { 3 } q } { ( 2 \pi ) ^ { 3 } } \iint \frac { d \tau _ { 1 } } { H \tau _ { 1 } } \frac { d \tau _ { 2 } } { H \tau _ { 2 } } 2 \theta ( \tau _ { 2 } - \tau _ { 1 } ) e ^ { - i \alpha ( p _ { 1 2 } + p _ { 2 1 } ) ( \tau _ { 2 } - \tau _ { 1 } ) } \left ( \left ( \frac { \tau _ { 1 } } { \tau _ { 2 } } \right ) ^ { 2 i \tilde { \lambda } _ { f } } + \left ( \frac { \tau _ { 2 } } { \tau _ { 1 } } \right ) ^ { 2 i \tilde { \lambda } _ { f } } \right ) G _ { \text{a} } ( \vec { k } _ { 1 }, \tau _ { 1 } ) G _ { \text{a} } ( - \vec { k } _ { 1 }, \tau _ { 2 } ) + \cdots$$

glyph[negationslash]

There is not a clean analytic formula interpolating between the small and large mass regions and a full result would need to be obtained numerically. For simplicity, we interpolate between the two expressions using the following potential in the SU (2) L symmetric form:

$$\text{result}, \quad & \text{$|$ from} \\ \text{tween} \\ \text{ctions} & & \text{$V_{h}=-\mu_{h}^{2}|\mathcal{H}|^{2}+\lambda_{h}|\mathcal{H}|^{4}-\frac{N_{c}y_{f}^{2}}{\pi^{2}}\lambda_{f}^{2}|\mathcal{H}|^{2}$} \exp \left [ - \frac { \pi y _ { f } ^ { 2 } | \mathcal{H}|^{2} } { \lambda_{f} H } \right ], \\ & \text{$\min to corrections that are independent of $\lambda_{e}$ where $u$}.$$

up to corrections that are independent of λ f , where µ h and λ h are the coefficients of the quadratic and quartic terms of the Higgs potential, and m h = √ 2 µ h is the physical mass of the Higgs boson.

## III. NON-GAUSSIAN SIGNATURE

The induced Higgs potential in Eq. (9) has a new minimum at

$$\mathfrak { i } _ { \mathfrak { l } } \quad v = \frac { 1 } { y _ { f } } \sqrt { \frac { 2 } { \pi } \lambda _ { f } H } \left ( 1 - \frac { e \pi \lambda _ { h } / y _ { f } ^ { 4 } } { N _ { C } \widetilde { \lambda } _ { f } } + \mathcal { O } ( \lambda _ { h } ^ { 2 } ) \right ). \quad ( 1 0 )$$

for v glyph[greatermuch] v ew . In the SM, the top quark Yukawa is y t = √ 2 m /v t EW ≈ 1 at the weak scale and runs to about 0 6 at very high energies, while the Higgs quartic is . λ h = m / v 2 h 2 2 EW ≈ 0 13 at the weak scale and runs to . λ h &lt; ∼ 0 01 . at very high energies. This ensures that the small quartic expansion is a very good approximation for λ /H f glyph[greatermuch] 1.

In [12], we discuss in detail how to compute the non analytic signal induced in the bispectrum by fermions interacting with the inflaton. The result in the squeezed where we only consider the terms where a = b . The terms with a = b are not divergent, not unrelated to the fact that the counterterm δm h / 2 2 2 for the Higgs mass can only cancel the divergences in δm 2 ⊕⊕ ( N ⊕⊕ ) and δm 2 glyph[circleminus]glyph[circleminus] ( N glyph[circleminus]glyph[circleminus] ). Moreover, these two counterterms are related by a minus sign (due to complex conjugation). As a result, we can find the correction to the Higgs potential from the diagrams in Fig. 3 by taking the difference between the two diagrams, removing in this way the corrections that are linearly divergent [15]

$$\frac { \delta m ^ { 2 } } { 3 k ^ { 3 } } = \frac { 1 } { 3 k ^ { 3 } } \frac { \delta m _ { \oplus \oplus } ^ { 2 } - \delta m _ { \ominus \ominus } ^ { 2 } } { 2 } = - \frac { y _ { f } ^ { 2 } \lambda _ { f } ^ { 2 } } { 3 \pi ^ { 2 } k ^ { 3 } }. \quad ( 7 ) \quad \text{term} \\ \text{phys}$$

This gives the correction from a single massless fermion flavor/color to the Higgs potential

$$\delta V _ { h } = - \frac { y _ { f } ^ { 2 } } { 2 \pi ^ { 2 } } \lambda _ { f } ^ { 2 } h ^ { 2 }$$

at leading order in the large ˜ λ f expansion 2 and to leading order in the Higgs field.

In the large m f limit, we obtain an exponential suppression with exponent ∼ -y 2 f h /λ H 2 f . We cannot calculate the coefficient of the exponent. As an approximation, we match this exponential suppression for the correction to the Higgs potential to the exponential suppression of the fermion density 3 .

2 In the massless limit, the potential can be computed to next-to-leading-order in large ˜ λ f expansion to be δV h = -y 2 f 2 π 2 ( λ 2 f -H 2 2 ) h 2 .

3 We acknowledge the potential existence of subleading terms in the small m /H f limit due to our inability to resum and reproduce the exact exponential suppression in the fermion density.

limit and at leading order in the large ˜ λ f expansion is

$$S ( k _ { 1 }, k _ { 2 }, k _ { 3 } ) \stackrel { k _ { 3 } \ll k _ { 1 } \sim k _ { 2 } } { \simeq } f _ { N L } ^ { ( \text{clock} ) } \left ( \frac { k _ { 3 } } { k _ { 1 } } \right ) ^ { 2 - 2 i \widetilde { \lambda } _ { f } }, \quad ( 1 1 ) \quad \text{struct} \quad \text{$\stackrel{ \sim \, \cdots } { \text{mode} } }$$

where m f = √ λ H/π f glyph[lessmuch] λ f , and

$$f _ { N L } ^ { ( \text{clock} ) } \approx & \ \frac { N _ { c } } { 6 \pi } \mathcal { P } _ { \zeta } ^ { - 1 / 2 } \left ( \frac { m _ { f } } { \Lambda _ { f } } \right ) ^ { 3 } \widetilde { \lambda } _ { f } ^ { 2 } & \quad \text{where} \\ & \quad \cdot \frac { e ^ { \pi \widetilde { \lambda } _ { f } } \widetilde { \mu } _ { f } \Gamma ( - i \widetilde { \mu } _ { f } ) ^ { 2 } \Gamma ( 2 i \widetilde { \mu } _ { f } ) ^ { 3 } } { 2 \pi \Gamma ( i ( \widetilde { \lambda } _ { f } + \widetilde { \mu } _ { f } ) ) ^ { 3 } \Gamma ( i ( \widetilde { \mu } _ { f } - \widetilde { \lambda } _ { f } ) + 1 ) } \,, \quad ( 1 2 ) \quad \text{so this} \\ \quad \cdot \frac { \Lambda _ { f } } { 2 \pi \, \mathcal { D } } \\ \dots \dots \,. \quad & \quad \sqrt { \sum _ { 2 } + \dots ^ { 2 } } \dots \alpha \, \approx \quad. \ \, \, \,$$

where µ f = √ λ 2 f + m 2 f and ˜ µ f = µ /H f . In the scenario discussed in this paper, the contribution to the Higgs potential from the top quark condensate dynamically adjusts the Higgs vev so that m f = √ λ H/π f , and the nonanalytic contribution to the bispectrum is just a function of ˜ λ f :

$$f _ { N L } ^ { ( \text{clock} ) } \approx \frac { 4 \sqrt { 2 } N _ { c } \mathcal { P } _ { \zeta } } { 3 e } \widetilde { \lambda } _ { f } ^ { 1 3 / 2 }. \quad \ \ \ ( 1 3 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

We show the contribution to f NL from the top quark in our scenario as a function of λ /H t in Fig. 4.

FIG. 4. Size of f NL generated by the top quark condensate in the scenario described in this paper (Eq. 13), as a function of λ /H t = c φ/ t ˙ Λ f H (or ˙ φ/ Λ , assuming 2 f c t = 1).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00019v1-figure-65ff4715-b96f-46ca-9489-af684e61b582.jpg)

## IV. REMARKS

In this article, we described the signatures of shift symmetric couplings between the inflaton and the SM, in particular the one to the top quark. This coupling dynamically sets the masses of the SM fermions and simultaneously determines the strength of the non-Gaussianity giving rise to the unique feature that the frequency of the oscillating pattern in ln( k /k 3 1 ) fixes the amplitude of the non-Gaussian signal. This provides a way to distinguish the scenario that we describe in this paper from other models generating this feature.

Surprisingly, high energy colliders probing the flavor structure of the SM [17] provide an alternative route for testing our scenario. For a typical slow-roll inflation model [18-21], the inflaton mass m inf is 4

$$m _ { \inf } & \ll H _ { \inf } \lesssim 1 0 ^ { - 2 } \Lambda _ { f } \,, \quad \quad ( 1 4 ) \\ \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot$$

where the last inequality comes from the requirement that ˜ λ f = P -1 / 2 ζ 2 π H Λ f = P -1 / 4 ζ √ 2 π ( ˙ φ Λ 2 f ) 1 / 2 satisfies ˙ φ &lt; ∼ Λ 2 f so that we can rely on the EFT expansion in powers of Λ . Experimental constraints from BaBar and exotic mef son decays imply a bound Λ f &gt; ∼ 10 TeV for scalars that couple with 'Yukawa-like couplings' to the three generations of SM fermions (see Secs. 4.1 and 4.2 in [17]). A detailed study of the flavor signatures of this scenario is beyond the scope of this paper.

A potentially observable signal in future measurements of the bispectrum (say f NL &gt; ∼ 10) requires a large ˜ λ f , which together with Λ f &gt; ∼ 10 TeV translates to a lower bound on the Higgs vev during inflation:

$$\int _ { 0 } ^ { 1 3 } v \approx \frac { 1 } { y _ { t } } \sqrt { \frac { 2 \lambda _ { f } H } { \pi } } \approx \frac { \sqrt { 8 \pi } } { y _ { t } } \widetilde { \lambda } _ { f } ^ { 3 / 2 } \mathcal { P } _ { \zeta } ^ { 1 / 2 } \Lambda _ { f } \gtrsim 3 0 0 \, \text{GeV}. \quad ( 1 5 )$$

Similarly to the correction to the Higgs vev during inflation from λ h , there is also a correction from v EW ≡ √ m / λ 2 h 2 h when λ f ∼ v EW . Such a correction amounts to an expansion in v 2 ew v 2 eπλ h y 4 t N λ c ˜ f , which is also a small correction as long as v &gt; v ew .

∼ A potentially significant correction to the scenario we have described can arise if there is a significant coupling between the inflaton and the Higgs boson in the form of

$$c _ { 2 } \left ( \frac { \partial \phi } { \Lambda _ { \mathcal { H } } } \right ) ^ { 2 } \mathcal { H } ^ { \dagger } \mathcal { H }.$$

This coupling can be generated from UV dynamics or by integrating out SM fermions that couple directly to the inflaton. In order for the scenario described in this paper to follow through, we needed to be able to suppress any large corrections to the Higgs mass. During inflation, Eq. (16) leads to a Higgs mass term of order c 2 ( ˙ φ/ Λ H ) 2 , potentially comparable to the mass term shown in Eq. (10). This mass term (like the term -µ 2 h |H| 2 in Eq. (9)) can lead to corrections to the Higgs vev v during inflation. Requiring that Eq. (16) is a negligible correction to v implies c 2 &lt; ∼ N y c 2 t π e 2 ≈ 0 1 if Λ . f = Λ H . If c 2 is large and positive, the minima in the Higgs potential would disappear (together with the signal), while if c 2 is large and negative, then we go back to the cases

4 We assume for simplicity a scenario where the inflaton mass today does not differ too much from the inflaton mass m inf during inflation.

studied in [6, 12], where the Higgs vev is not dynamically adjusted by the mechanism discussed here.

The signature that we discuss can also arise whenever there is a fermion whose mass is generated from an O (1) Yukawa coupling. For example, if the dark sector contains dark fermions whose mass is generated (also partially) from Yukawa couplings to a dark Higgs boson, then the leading order shift symmetric couplings between the dark fermions and the inflaton would result in a signal similar to the one we discuss (with a different N c ).

Soon, the limits on non-Gaussanities and the scalarto-tensor ratio will improve dramatically [9, 10, 22, 23]. It is thus extremely interesting to ask what sort of signals might be discovered in the upcoming years. In this article, we described the most minimal way in which the SM itself may reappear in the sky. In the near future, a combination of flavor measurements at particle colliders and bispectrum measurements at the cosmological collider may combine to uncover deep mysteries of the universe.

- [1] X. Chen and Y. Wang, JCAP 1004 , 027 (2010), arXiv:0911.3380 [hep-th].
- [2] D. Baumann and D. Green, Phys. Rev. D85 , 103520 (2012), arXiv:1109.0292 [hep-th].
- [3] N. Arkani-Hamed and J. Maldacena, (2015), arXiv:1503.08043 [hep-th].
- [4] H. Lee, D. Baumann, and G. L. Pimentel, JHEP 12 , 040 (2016), arXiv:1607.03735 [hep-th].
- [5] X. Chen, Y. Wang, and Z.-Z. Xianyu, JHEP 04 , 058 (2017), arXiv:1612.08122 [hep-th].
- [6] S. Kumar and R. Sundrum, JHEP 05 , 011 (2018), arXiv:1711.03988 [hep-ph].
- [7] S. Alexander, S. J. Gates, L. Jenks, K. Koutrolikos, and E. McDonough, (2019), arXiv:1907.05829 [hep-th].
- [8] S. Lu, Y. Wang, and Z.-Z. Xianyu, (2019), arXiv:1907.07390 [hep-th].
- [9] J. B. Mu˜ noz, Y. Ali-Ha¨ ımoud, and M. Kamionkowski, Phys. Rev. D92 , 083508 (2015), arXiv:1506.04152 [astroph.CO].
- [10] P. D. Meerburg, M. M¨ unchmeyer, J. B. Mu˜ noz, and X. Chen, JCAP 1703 , 050 (2017), arXiv:1610.06559 [astro-ph.CO].
- [11] X. Chen, Y. Wang, and Z.-Z. Xianyu, JHEP 09 , 022 (2018), arXiv:1805.02656 [hep-ph].
- [12] A. Hook, J. Huang, and D. Racco, arXiv:1907.10624.

## ACKNOWLEDGMENTS

We would like to acknowledge Asimina Arvanitaki for relentlessly inviting us to think about this scenario. We also thank Prateek Agrawal, David Curtin, Savas Dimopoulos, Jiji Fan, Daniel Ega˜ na-Ugrinovic, Matthew Johnson, Soubhik Kumar, Gustavo Marques-Tavares, Maxim Pospelov, Antonio Riotto, Raman Sundrum, Zhong-Zhi Xianyu for many useful discussions. JH and DR thank the Stanford Institute for Theoretical Physics, and AH and JH thank MIAPP and KITP Santa Barbara, for generous hospitality during the completion of this work. This research was supported in part by Perimeter Institute for Theoretical Physics. Research at Perimeter Institute is supported in part by the Government of Canada through the Department of Innovation, Science and Economic Development Canada and by the Province of Ontario through the Ministry of Economic Development, Job Creation and Trade. AH is supported in part by the NSF under Grant No. PHY-1620074 and by the Maryland Center for Fundamental Physics (MCFP). This research was supported in part by the National Science Foundation under Grant No. NSF PHY-1748958.

- [13] A. D. Linde, Contemp. Concepts Phys. 5 , 1 (1990), arXiv:hep-th/0503203 [hep-th].
- [14] K. M. Benson, J. Bernstein, and S. Dodelson, Phys. Rev. D44 , 2480 (1991).
- [15] X. Chen, Y. Wang, and Z.-Z. Xianyu, JHEP 08 , 051 (2016), arXiv:1604.07841 [hep-th].
- [16] X. Chen, Y. Wang, and Z.-Z. Xianyu, JCAP 1712 , 006 (2017), arXiv:1703.10166 [hep-th].
- [17] M. J. Dolan, F. Kahlhoefer, C. McCabe, and K. Schmidt-Hoberg, JHEP 03 , 171 (2015), [Erratum: JHEP07,103(2015)], arXiv:1412.5174 [hep-ph].
- [18] A. H. Guth, Phys. Rev. D23 , 347 (1981), [Adv. Ser. Astrophys. Cosmol.3,139(1987)].
- [19] A. D. Linde, Quantum Cosmology , Phys. Lett. 108B , 389 (1982), [Adv. Ser. Astrophys. Cosmol.3,149(1987)].
- [20] D. H. Lyth and A. Riotto, Phys. Rept. 314 , 1 (1999), arXiv:hep-ph/9807278 [hep-ph].
- [21] D. Baumann, in Physics of the large and the small, TASI 09, proceedings of the Theoretical Advanced Study Institute in Elementary Particle Physics, Boulder, Colorado, USA, 1-26 June 2009 (2011) pp. 523-686, arXiv:0907.5424 [hep-th].
- [22] K. N. Abazajian et al. (CMB-S4), (2016), arXiv:1610.02743 [astro-ph.CO].
- [23] S. Shandera, N. Dalal, and D. Huterer, JCAP 1103 , 017 (2011), arXiv:1010.3722 [astro-ph.CO].
Anson Hook Junwu Huang Davide Racco
Minimal signatures of the Standard Model in non-Gaussianities
We show that the leading coupling between a shift symmetric inflaton and the
Standard Model fermions leads to an induced electroweak symmetry breaking due
to particle production during inflation, and as a result, a unique oscillating
feature in non-Gaussianities. In this one parameter model, the enhanced
production of Standard Model fermions dynamically generates a new electroweak
symmetry breaking minimum, where the Higgs field classically rolls into. The
production of fermions stops when the Higgs expectation value and hence the
fermion masses become too large, suppressing fermion production. The balance
between the above-mentioned effects gives the Standard Model fermions masses
that are uniquely determined by their couplings to the inflaton. In particular,
the heaviest Standard Model fermion, the top quark, can produce a distinct
cosmological collider physics signature characterised by a one-to-one relation
between amplitude and frequency of the oscillating signal, which is observable
at future 21-cm surveys.