\section{Introduction} \label{sec:intro}
		The  $\Xi$ and $ \Omega$ hypernuclei offer a mixture of simplicity
		and fundamental interest in connection with the interactions between
		nucleons and strange particles. Hyperon interactions are not known sufficiently well due to the limited scattering data. However, more precise data on light hypernuclei (from high resolution gamma-ray experiments~\citep{Tamura2000}), advanced few-body theoretical methods~\cite{hiyama2000,hiyama2008,hiyama2014,hiyama2019,Moosavi2019,Sekihara2018,Garcilazo2018,Garcilazo2019}, quark delocalization, color screening~\cite{Ping199}, constituent quark models~\citep{huang2015}, lattice QCD calculations~\cite{etminanmpla2014,nemura2014,Sasaki2018} and femtoscopic analyses of pp, pA and AA collisions in the ALICE and STAR experiments~\citep{Acharya2019} have recently been providing us with valuable information. The KEK-E373 experiment reported the first evidence of a bound $\Xi$-hypernucleus $^{14}N+\Xi$, i.e., the so called KISO event~\citep{Nakazawa2015}. Recently, the first experimental observation of an attractive strong interaction between a proton and a hyperon $ \Xi $ has been reported by ALICE collaboration~\citep{Acharya2019}.

		The theoretical efforts of the lattice HAL QCD Collaboration have led to the derivation of baryon-baryon interactions near the physical pion mass~\citep{aoki2012}. The most recent results of these efforts hint to the existence of shallow bound states in $\Omega N$ systems~\citep{Iritani2019prb}. We study the $\Omega+\alpha$ system making use of this $\Omega N$ potential. The depth of the potential in $\Omega \alpha $ and $\Xi \alpha$ systems are uncertain, because there are not enough experimental data for their nuclear bound states. However, for the latter, some diverse phenomenological potentials are considered~\cite{dover1983,Filikhin2008,Garcilazo2016}.

		Motivated by the aforementioned description, and vast applications of $\Omega \alpha $ and $\Xi \alpha$ interactions in 3-,4- and 5-body cluster structures of $ \Omega$ hypernuclei (describing the response of 3$\alpha$ system to the addition of $\Omega$ particle)~\citep{hiyama2000}, $\alpha$ cluster model approach~\citep{hiyama2014} and variational four-body calculation~\citep{Garcilazo2020}, we present $\Omega \alpha $ and $\Xi \alpha$ interactions in a simple Woods-Saxon type form.

		In order to test and validate our method, we apply it to $\Xi \alpha$ systems and compare the results with the phenomenological Dover-Gal (DG) potential type~\citep{dover1983}.
		Additionally, we use ESC08c Nijmegen model $\Xi N$ potential in $^{3}S_{1}$ channel~\citep{Nagels}.

		Here, we consider a $Y +\alpha$ system because of the low compressibility of the $\alpha$-cluster and high reaction thresholds that enable us to employ a one-channel approximation over a wide energy range. The $Y+\alpha$ system is studied using a single-folding potential (SFP) method. In this model, the $Y+\alpha$ system consists of an alpha and a $ Y $ particle moving in the effective $Y \alpha$ potential. The effective nuclear potential is approximated by the single-folding of nucleon density $\rho\left(\vec{r}^{\prime}\right)$ in
		the $\alpha$-particle and hyperon-nucleon potential $V_{YN}\left(\left|\vec{r}-\vec{r}^{\prime}\right|\right)$ between the $Y$ particle at $\vec{r}^{\prime}$
		and the nucleon at $\vec{r}$~\cite{Satchler1979,Miyamoto2018}. Then, the resulting
		$Y\alpha$ potential is fitted to a separable form. Finally, we
		solve the Schr\"{o}dinger equation by using the fitted $Y \alpha$
		potential in the infinite volume and extract its scattering observables
		from the asymptotic behavior of the wave function. The model is expected
		to be accurate only for the low energy properties of the $Y+\alpha$ system
		since it is based on a $Y N$ potential, which is fitted to low
		energy $Y N$ scattering parameters.

		We should emphasize that the coupling of $ \Omega N $ to higher-mass ($ \Lambda\Xi^{*} $) and lower-mass ($ \Lambda\Xi $ and $ \Sigma\Xi $) channels is not taken into account since we assume that these contributions are of second order of smallness to the binding energy of few-body systems~\citep{Iritani2019prb}. To draw a definite conclusion about the binding energy of $\Omega +\alpha$ system, it is necessary to perform a coupled-channel analysis of the HAL QCD method~\citep{aoki2013}. Moreover, in the calculations, the Coulomb force is not taken into account.

		This paper is organized as follows. In Sec.~\textnormal{\ref{sec:folding-Model}}, after a
		brief discussion on the single-folding potential method, we introduce and parametrize all dominating set of input parameters, i.e., those of the two-body %
			potentials. In Sec.~\textnormal{\ref{sec:result}}, we present and discuss the results. Finally, in Sec.~\textnormal{\ref{sec:summery-conclusion}}, a summary and conclusions are presented.

		\section{Single-folding potential model} \label{sec:folding-Model}

		We obtain the effective potential of $Y+\alpha$ systems by
		using the single-folding potential model~\cite{Satchler1979,Miyamoto2018}. This method is briefly outlined in the following. The $Y \alpha$ potential is defined as

		\begin{equation}
		V_{Y\alpha}\left(\vec{r}\right)=\int\rho\left(\vec{r}^{\prime}\right)V_{Y N}\left(\left|\vec{r}-\vec{r}^{\prime}\right|\right)d\tau^{\prime},\label{eq:V_alfaOmega}
		\end{equation}

		where $\rho\left(\vec{r}^{\prime}\right)$ is the nucleon density in
		the $\alpha$-particle at a distance $\vec{r}^{\prime}$ from its
		center-of-mass, and is given by~\citep{Akaishi1986},

		\begin{equation}
		\rho\left(\vec{r}^{\prime}\right)=4\left(\frac{4\beta}{3\pi}\right)^{3/2}\exp\left(-\frac{4}{3}\beta r^{\prime2}\right).\label{eq:nucleon-density}
		\end{equation}

		The integration in Eq.~(\ref{eq:V_alfaOmega}) is over all space as permitted
		by $\rho\left(\vec{r}^{\prime}\right)$. The required normalization
		condition is satisfied by

		\begin{equation}
		\int\rho\left(\vec{r}^{\prime}\right)d\tau^{\prime}=4\left(\frac{4\beta}{3\pi}\right)^{3/2}\intop_{0}^{\infty}\mathop{\exp\left(-\frac{4}{3}\beta r^{\prime2}\right)4\pi r^{\prime2}dr^{\prime}}=4.
		\end{equation}

		The constant $\beta$ is determined from the rms radius of $\textrm{\ensuremath{^{4}}He}$~\citep{Akaishi1986},

		\begin{equation}
		\textrm{\ensuremath{r_{r.m.s}}}=\frac{3}{\sqrt{8\beta}}=1.47\:\textrm{fm}.
		\end{equation}

		In Eq.~(\ref{eq:V_alfaOmega}), $V_{YN}\left(\left|\vec{r}-\vec{r}^{\prime}\right|\right)$ is the potential in configuration space between the $Y$ particle at $\vec{r}$
		and the nucleon at $\vec{r}^{\prime}$.

		We take the $\Xi N$ potential
		in $^{3}S_{1}$ channel and simulate the ESC08c Nijmegen model, which consists
		of local central Yukawa-type potentials with attractive and repulsive terms~\cite{Nagels,Garcilazo2016},

		\begin{equation}
		V_{\Xi N}\left(\vec{r}\right)=-568\frac{\exp\left(-4.56r\right)}{r}+425\frac{\exp\left(-6.73r\right)}{r}. \label{eq:pot-XiN}
		\end{equation}

		The low-energy data of this potential is listed in Table~\ref{tab:ERE-ON-XiN}. In the case of $\Omega N$, we use S-wave and spin $2$ $\Omega N$
		potential, which is given by HAL QCD Collaboration with nearly physical
		quark masses~\citep{Iritani2019prb}. The lattice discrete potential is fitted by an analytic function composed of an attractive Gaussian core plus a long range
		\textrm{(Yukawa\ensuremath{)^{2}}} attraction with a form factor of~\citep{etminan2014},

		\begin{equation}
		V_{\Omega N}\left(r\right)=b_{1}e^{\left(-b_{2}r^{2}\right)}+b_{3}\left(1-e^{-b_{4}r^{2}}\right)\left(\frac{e^{-m_{\pi}r}}{r}\right)^{2},\label{eq:NOmega_pot}
		\end{equation}
		the pion mass in Eq.~(\ref{eq:NOmega_pot}), which is taken from the simulation, is $\textrm{\ensuremath{m_{\pi}}=146}$ MeV. The lattice results are fitted reasonably well, $\chi^{2}/d.o.f\simeq1$, with four different sets of parameters given in Table \ref{tab:Fit_para}.
		The low-energy data of this potential is also given in Table~\ref{tab:ERE-ON-XiN}.

		\begin{center}
			\tabcaption{ The low-energy parameters, scattering length, $a_{0}$, effective range, $r_{0}$,  and binding energy, $B_{Y N}$, of ESC08c Nijmegen $\Xi N$~\cite{Nagels,Garcilazo2016} given by Eq.~(\ref{eq:pot-XiN}), and the HAL QCD $\Omega N$ potential~\citep{Iritani2019prb} given by Eq.~(\ref{eq:NOmega_pot}). \label{tab:ERE-ON-XiN}}
			\footnotesize
			\begin{tabular*}{80mm}{c@{\extracolsep{\fill}}cccc}
				\toprule	System &Channel& $a_{0}$(fm)&$r_{0}$(fm)& $B_{Y N}$(MeV) \\
				\hline
				$\Omega N$ & $ ^5S_{2}$ & 5.30 & 1.26  & 1.54 \\
				$\Xi    N$ & $^3S_{1}$  & 4.91 & 0.527 & 1.67   \\
				\bottomrule
			\end{tabular*}
		\end{center}

		\begin{center}
			\tabcaption{Fitting parameters in Eq.~(\ref{eq:NOmega_pot}) for different models, $P_{i}$, for $^{5}S_{2}$ $\Omega N$ interaction~\citep{Iritani2019prb}.\label{tab:Fit_para} }
			\footnotesize
			\begin{tabular*}{80mm}{c@{\extracolsep{\fill}}cccc}
				\toprule	& $P_{1}$ & $P_{2}$ & $P_{3}$ & $P_{4}$\\
				\hline
				$b_{1}\left(\mathrm{MeV}
				\right)$ & -306.5  & -313.0 & -316.7 & -296\\
				$b_{2}\left(\mathrm{fm}^{-2}\right)$ & 73.9  & 81.7 & 81.9  & 64\\
				$b_{3}\left(\mathrm{MeV}.fm^{-2}\right)$ & -266  &  -252 & -237 & -272\\
				$b_{4}\left(\mathrm{fm}^{-2}\right)$ & 0.78  & 0.85  & 0.91  & 0.76\\
				\bottomrule
			\end{tabular*}
		\end{center}

		\section{Results} \label{sec:result}

		In order to test $\Xi \alpha$ potential obtained from SFP model,
		we employed the phenomenological potential of Woods-Saxon type for
		the $\Xi \alpha$ interaction by using Dover-Gal (DG) model given
		in~\cite{Filikhin2008,dover1983}

		\begin{equation}
		V_{\Xi\alpha}^{DG}\left(r\right)=-V_{0}\left[1+\exp\left(\frac{r-R}{c}\right)\right]^{-1} , \label{eq:DG}
		\end{equation}

		where $ V_{0}$ is the depth parameter, $ R=1.1A^{\left(1/3\right)}$ with $ A $ being the mass number of the nuclear core (here, $A=4$ for the alpha particle) and $c$ is the surface diffuseness. The values of these three parameters in DG model are given in Table~\ref{tab:ERE-DO-SFP}. The $\Xi+\alpha$ system is bound by this potential with an energy
		of $E_{B}=-2.1$ MeV. It is important to note that DG potential has no repulsive core.

	\end{multicols}

	\begin{center}
		\tabcaption{DG potential model parameters of Eq.~(\ref{eq:DG}) from~\cite{Filikhin2008,dover1983}. The values of fitting parameters of $\Xi \alpha$ are obtained by fitting SFP model to a function in the same form as Eq.~(\ref{eq:DG}). The corresponding low-energy parameters, scattering length, effective range and binding energy of both models are given. The obtained results, using experimental masses of $\alpha$ and $\Xi$, are $3727.38\:\textrm{MeV}/c^{2}$ and $1318.07\:\textrm{MeV}/c^{2}$, respectively.\label{tab:ERE-DO-SFP}}.
		\footnotesize
		\begin{tabular*}{95mm}{ccccccc}
			\toprule			Model & $ V_{0} $ (MeV)& $ R $(fm)& $ c $(fm)& $a_{0}$(fm)&$r_{0}$(fm)& $B_{\Xi \alpha}$(MeV) \\
			\hline
			DG  & 24   & 1.74 & 0.65 & -4.9 & 1.9 & -2.1\\
			SFP & 24.4 & 1.72 & 0.31 & -6.6 & 1.9 & -1.54\\
			\bottomrule
		\end{tabular*}
	\end{center}

	\begin{multicols}{2}

		To obtain observables, such as scattering phase shifts and binding
		energy, we fit $V_{Y\alpha}\left(\vec{r}\right)$ to the Wood-Saxon form using the function in Eq.~(\ref{eq:DG}) with three parameters $ V_{0}, R$ and $c$. The results of fitting for these parameters are presented in Table~\ref{tab:ERE-DO-SFP}.

		The obtained single-folding potential, $V_{\Xi\alpha}\left(\vec{r}\right)$, its corresponding fit function, $ V_{fit}(r) $, and  $V_{\Xi\alpha}^{DG}\left(r\right)$ (for comparison) are shown in Fig.~\ref{fig:vc-Xi-fit-DG}. We solve the Schr\"{o}dinger equation with the fitted potential in the infinite volume and extract its scattering observables from the asymptotic behavior of the wave function. For comparison, in Fig.~\ref{fig:phase-xialfa-DG-SFP}, the phase shifts from DG and SFP model potentials are shown.

		\begin{center}
			\includegraphics[width=8cm]{vc_xialfa_fit_dg.eps}
			\figcaption{The single-folding potential, $V_{\Xi\alpha}\left(\vec{r}\right)$, for the $\Xi N$  interaction in $^{3}S_{1}$ channel given in~\cite{Nagels,Garcilazo2016}. $V_{fit}(r)$ (Red line) shows the results of the fitting by using the same form as Eq.~(\ref{eq:DG}). For comparison, we also present Dover-Gal $\Xi\alpha$ potential (green dashed line)~\cite{Filikhin2008,dover1983}, i.e.,  $V_{\Xi\alpha}^{DG}\left(r\right)$ in Eq.~(\ref{eq:DG}).
				\label{fig:vc-Xi-fit-DG}}
		\end{center}

		\begin{center}
			\includegraphics[width=8cm]{phasexialfa_DO_SF.eps}
			\figcaption{ The extracted $\Xi\alpha$ phase shifts for two potential models, DG and SFP, are given in Fig.~\ref{fig:vc-Xi-fit-DG} for comparison. According to this figure, one can see a fairly good agreement between these two models. \label{fig:phase-xialfa-DG-SFP}}
		\end{center}

		The effective range expansion (ERE) of the phase shifts up to the
		next-leading-order (NLO) reads

		\begin{equation}
		k\cot\delta_{0}=-\frac{1}{a_{0}}+\frac{1}{2}r_{0}k^{2}+\mathcal{O}\left(k^{4}\right),\label{eq:ERE}
		\end{equation}

		with $\left(a_{0}\right)$ and $\left( r_{0}\right)$ being the scattering
		length and effective range, respectively. The results of the calculations for the binding energy and the ERE parameters $\left(a_{0}, r_{0}\right)$ are given in Table~\ref{tab:ERE-DO-SFP}. According to the results
		in Table~\ref{tab:ERE-DO-SFP}, a good agreement can be observed between the single-folding potential method and the phenomenological potential of Dover-Gal model.

		The single-folding potential $V_{\Omega\alpha}\left(\vec{r}\right)$
		for different models of $\Omega N$ interaction ($ P_{i} $, Table~\ref{tab:Fit_para}) are shown in Fig.~\ref{fig:vcF} (a). We summarize the results of fitting and the corresponding parameters in Fig.~\ref{fig:vcF} (b) and Table~\ref{tab:Fit_para-of-Oalfa-pot}, respectively. We adjust the depth $ V_{0} $ in such a way to produce the best fit for $R$ and $c$ parameters, i.e., $\chi^{2}/d.o.f\simeq1 $.

			Garcilazo and Valcarce~\cite{Garcilazo2018} showed that $\Omega N$  $^{5}S_{2}$ and $NN$ $^{3}S_{1}$ channels give rise to a $\Omega d$ bound state in the state with maximal spin $\left(I,J^{P}\right)=\left(0,5/2^{+}\right)$ with a binding energy of $\sim17$ MeV measured with respect to the $NN\Omega$ threshold by solving the three-body bound-state Faddeev equations. Here, we obtain an $\Omega\alpha$ binding energy of $\sim23$ MeV, as shown in Table~\ref{tab:ERE-oalfa}. Since the potential must be more attractive than the approximated single-folding potential, the resulting energy is only an upper bound for the $\Omega \alpha$ system, and the binding energy of $\Omega\alpha$ is greater than of $\Omega d$, which seems completely reasonable.

		\begin{center}
			\tabcaption{
				The values of parameters are obtained by fitting the $\Omega \alpha$ SFP potential to a function in the same form as Eq.~(\ref{eq:DG}), for different models, $P_{i}$ , of the $\Omega N$ interaction~\citep{Iritani2019prb}.\label{tab:Fit_para-of-Oalfa-pot}}
			\footnotesize
			\begin{tabular*}{55mm}{ccccc}
				\toprule
				& $P_{1}$ & $P_{2}$ & $P_{3}$ & $P_{4}$\\
				\hline
				$V_{0}\left(\textrm{MeV}\right)$ &-61  & -61  & -61  & -61\\
				$R\left(\textrm{fm}\right)$  & 0.47 & 0.47  & 0.47 & 0.47\\
				$c\left(\textrm{\textrm{fm}}\right)$ & 1.7 & 1.7 & 1.7 & 1.7 \\
				\bottomrule
			\end{tabular*}
		\end{center}

		\begin{center}
			\includegraphics[width=8cm]{vc_Oalfa_1234.eps}
			\includegraphics[width=8cm]{vc_fx_p1.eps}
			\figcaption{(a) The single-folding potential, $V_{\Omega\alpha}\left(\vec{r}\right)$,
				for different models of $\Omega N$ interaction ($ P_{i} $) given in Table~
				\ref{tab:Fit_para}. (b) The single-folding potential, $V_{\Omega\alpha}\left(\vec{r}\right)$,
				and $V_{fit}(r)$ fit function (Red line), which has the same form as Eq.~(\ref{eq:DG}) for the set $P_{1}$. \label{fig:vcF}}
		\end{center}

		Shown in Fig.~\ref{fig:Phase-shift.} is the S-wave scattering phase shift $\delta_{0}$ as a function of the kinetic energy. In Table~\ref{tab:ERE-oalfa}, we present the binding energies and ERE parameters $\left(a_{0},r_{0}\right)$ obtained from $\Omega \alpha$ phase shifts, for different models of $\Omega N$ interaction reported in
		Ref.~\citep{Iritani2019prb} and summarized in Table~\ref{tab:Fit_para}.

		\begin{center}
			\includegraphics[width=8cm]{phaseOalfa1234}
			\figcaption{The S-wave scattering phase shift $\delta_{0}$ as a function of the kinetic energy, $k^{2}/(2\mu)$.\label{fig:Phase-shift.}}
		\end{center}
	\end{multicols}

	\begin{center}
		\tabcaption{Scattering length, $a_{0}$, effective range, $r_{0}$, and binding energy, $B_{\Omega\alpha}$, of $\Omega\alpha$ for the different models of $\Omega N$ interaction given in Table~\ref{tab:Fit_para}~\citep{Iritani2019prb}. The results that are obtained using the experimental masses of $\alpha$ and $\Omega$ are $3727.38\:\textrm{MeV}/c^{2}$ and $1672.45\:\textrm{MeV}/c^{2}$, respectively. The values in parentheses correspond to the masses of $\alpha$ and $\Omega$ derived by the HAL QCD Collaboration, i.e., $3818.8\:\textrm{MeV}/c^{2}$and $1711.5\:\textrm{MeV}/c^{2}$, respectively~\citep{Iritani2019prb}.\label{tab:ERE-oalfa}}
		\footnotesize
		\begin{tabular*}{100mm}{@{\extracolsep{\fill}}ccccc}
			\toprule
			& $P_{1}$ & $P_{2}$ & $P_{3}$ & $P_{4}$\\
			\hline
			\footnotesize
			$a_{0}$ (fm) &-0.99(-0.93)& -1.01(-0.96) & -1.05(-1.00) & -0.98(-0.92)\\
			$r_{0}$ (fm) & 0.67(0.67) & 0.67(0.67) & 0.67(0.67) & 0.67(0.68) \\
			$B_{\Omega\alpha}$ (MeV) & -22.9 (-23.3) & -22.8(-23.2) & -22.4(-22.8) & -23.0(-23.4)\\
			\bottomrule
		\end{tabular*}
	\end{center}

	\begin{multicols}{2}

		\section{Summary and conclusions} \label{sec:summery-conclusion}

		We derived a simple Woods-Saxon type form for the potentials of $\Xi+\alpha$ and $\Omega+\alpha$ systems by making use of, respectively, the ESC08c Nijmegen $\Xi N$ potential in $^{3}S_{1}$ channel and the HAL QCD Collaboration $\Omega N$  potential in $^{5}S_{2}$ channel plus density function of the alpha particle in single-folding potential method.

			We showed that the effective central folding potential of $\Omega \alpha$ may have simple Wood-Saxon form and estimated the upper bound for binding energy of $\Omega$ particle on a $\alpha$.

		Our method was tested against the phenomenological potential of Woods-Saxon type for the $\Xi\alpha$ interaction by the Dover-Gal model and a fairly good agreement was found between the two methods.

		The scattering length and the effective range were obtained by solving the
		Schr\"{o}dinger equation using the resultant potential. The binding energies
		of $\Xi+\alpha$ and $\Omega+\alpha$ systems were about $-1.5$ and $-23$ MeV, respectively. These results indicate that $\Omega\alpha$ hypernuclei are deeply bound states or resonances, which may be experimentally observed in the real world.

		We should emphasize that the cases of coupling of $ Y N $ to higher- and lower-mass channels were not taken into account. The calculations do not also take into account the Coulomb force. To draw a definite conclusion about the binding energy of $Y  \alpha$, it is necessary to do a coupled-channel analysis.

		We hope that our results can be used as tests of various theoretical
		models for the exotic nuclei structures, especially in few $\alpha$ cluster structures of $ \Omega$ hypernuclei (describing the response of the few $\alpha$ systems to the addition of the $\Omega$ particle)~\citep{hiyama2000}, $\alpha$ cluster model approach~\citep{hiyama2014} and possible future experiments, where these lattice-QCD-based predictions may be tested.
		\vspace{3mm}


\section{Introduction}
\label{sec:intro}

Cosmological observations provide detailed information about our Universe on the largest observable scales. Cosmic microwave background (CMB) measurements \cite{Gorski:1996cf,Bennett:2012zja, Aghanim:2019ame}, for instance, demonstrate that fluctuations in the matter and energy persist on cosmological scales. There is no causal explanation for the origin of these fluctuations in Standard Big Bang cosmology. Scenarios of early universe cosmology such as the Inflationary Universe \cite{Guth} provide a causal mechanism to generate these fluctuations. A key aspect of both inflationary cosmology and of other scenarios that provide an explanation for the origin of structure in the universe (see e.g. \cite{RHBrev} for a comparative review) is the existence of a phase in the early universe during which the Hubble horizon $H^{-1}(t)$, where $H(t)$ is the Hubble expansion rate as a function of time $t$, shrinks in comoving coordinates. The Hubble horizon provides the limiting length above which causal physics that is local in time cannot create fluctuations. In both inflationary cosmology and in the proposed alternatives, comoving length scales which are probed in current cosmological experiments were inside the Hubble horizon at early times.  It is postulated that fluctuations in both matter \cite{ChibMukh} and gravitational waves \cite{Starob} originate as quantum vacuum perturbations that exit the Hubble radius during the early phase, are squeezed and classicalize, and then re-enter the Hubble radius at late times to produce the CMB anisotropies and matter density perturbations that we observe today.

In \cite{Jerome} (see also \cite{Niemeyer}) it was realized that if the inflationary phase lasts somewhat longer than the minimal period, then the length scales we observe today originate from modes that are smaller than the Planck length during inflation. This was called the {\it trans-Planckian problem} for cosmological fluctuations (see also \cite{Greene}).  This problem was viewed not so much as an issue with a particular model, but more as a question of how to treat trans-Planckian modes in such a situation.  It has been conjectured \cite{Bedroya:2019snp}, however, that this trans-Planckian problem can never arise in a consistent theory of quantum gravity and that all the models which would lead to such issues are inconsistent and belong to the Swampland.  This is called the {\it Trans-Planckian Censorship Conjecture} (TCC).

According to the TCC no length scales which exit the Hubble horizon could ever have had a wavelength smaller than the Planck length.  In Standard Big Bang cosmology no modes ever exit the Hubble horizon, and the TCC has no implications (indeed the TCC is automatically satisfied for all models with a $w\geq {-1/3}$). However, in all early universe scenarios which can provide an explanation for the origin of structure, modes exit the Hubble horizon in an early phase. If $a_i$ is the value of the cosmological scale factor at the beginning of the new early universe phase, and $a_f$ is the value at the time of the transition from the early phase to the phase of Standard Big Bang expansion, the TCC reads
\be \label{TCC}
\frac{a_f}{a_i} \, < \,  \frac{M_{pl}}{H_f} \, ,
\ee
where $H_f$ is the radius of the Hubble horizon at the final time $t_f$ and $M_{pl}$ is the reduced Planck mass.

In \cite{Bedroya:2019snp}, the relationship between the TCC and other {\it swampland conjectures} \cite{swamp} which have recently attracted a lot of attention (see e.g. \cite{swamprevs} for reviews) was discussed. Here, we focus on the consequences of the TCC for inflationary cosmology.

It is clear from the form of (\ref{TCC}) that the TCC will have strong implications for inflationary cosmology. In the case of de Sitter expansion, $a_f$ is exponentially larger than $a_i$, and hence (\ref{TCC}) strictly limits the time duration of any inflationary phase. The implications for some alternative early universe scenarios are weaker. For example, in {\it String Gas Cosmology} \cite{BV}, the early phase is postulated to be quasi-static. Hence, the condition (\ref{TCC}) is satisfied: no modes which were larger than the Hubble scale at the beginning of the Standard Cosmology phase ever had a wavelength smaller than the Planck length. The same is true for the various bouncing scenarios (e.g. the  {\it matter bounce} \cite{FB}, and the {\it Pre-Big-Bang} \cite{PBB} and {Ekpyrotic} \cite{Ekp} scenarios)\footnote{Of course it is possible that these models have issues related to other Swampland conditions.}, where the initial phase is one of contraction. This is true as long as the energy scale at the bounce point is smaller than the Planck scale. In the following we will study the consequences of the TCC for inflationary cosmology.

The outline of this paper is as follows: In the following section we discuss general constraints imposed by the TCC on the energy scale of inflation and the resulting consequences for the amplitude of gravitational waves. These conclusions do not depend on what drives inflation, only that it occurred. In Section 3, we then specialize to slow-roll inflation models, and show that consistency with the TCC leads requires fine-tuning of the initial conditions.

We will work in the context of homogeneous and isotropic cosmologies with $4$ space-time dimensions. For simplicity, we assume spatially flatness so that the metric can be written in the form
\be \label{FRW}
ds^2 \, = \, dt^2 - a(t)^2 d{\bf x}^2 \, ,
\ee
where ${\bf x}$ are the spatial comoving coordinates and $a(t)$ is the scale factor (which can be normalized such that $a(t_0) = 1$, where $t_0$ is the present time). The Hubble expansion rate is
\be
H(t) \, \equiv \, \frac{{\dot{a}}}{a} \, ,
\ee
and its inverse is the Hubble radius. As is well known (see \cite{MFB} for an in-depth review of the theory of cosmological fluctuations and \cite{RHBfluctrev} for an overview), quantum-mechanical fluctuations oscillate on sub-Hubble scales, whereas they freeze out and become squeezed when the wavelength is larger than the Hubble radius. During a phase of accelerated expansion, the proper wavelengths of fluctuations initially smaller than the Hubble scale can be stretched to super-Hubble scales. This transition from sub-Hubble to super-Hubble is referred to as horizon-crossing. The TCC prohibits horizon-crossing for modes with initial wavelengths smaller than the Planck length.

We will be considering models of inflation in which a canonically normalized scalar field $\phi$ with potential energy $V(\phi)$ constitutes the matter field driving the accelerated expansion of space. We will be using
units in which the speed of sound, Boltzmann's constant and $\hbar$ are set to 1.

\section{Implications of the TCC for the Energy Scale of Inflation}

In this section we work in the approximation that the Hubble expansion rate during the period
of inflation is constant. In order for inflation to provide a solution to the structure formation problem of Standard Big Bang cosmology, the current comoving Hubble radius must originate inside
the Hubble radius at the beginning of the period of inflation (see Fig. 1). This condition reads
\be \label{success}
\frac{1}{H}\cdot e^{N_+} \cdot \frac{a_{R}}{a_{end}}  \cdot \frac{T_R\  g_*(T_R)^{1/3}}{T_0 \ g_*(T_0)^{1/3}}\simeq \, \frac{1}{H_0} ,
\ee
where $H$ is the Hubble scale during inflation, $a_{\rm{end}}$ and $a_R$ are the values of the scale factor at the
end of inflation and when reheating is completed, and $g_{*}$ indicates
the number of spin degrees of freedom in the thermal bath.
Here $N_{+}$ is the number of e-foldings accrued during the inflation after the CMB-scale modes exit the horizon, $T_0$ is the temperature of the CMB at the present
time, and $T_R$ is the corresponding temperature after reheating.   Equation \eqref{success} can be summarized as follows. We start with a Hubble horizon scale which at the time of the inflation is $1/H$, by the end of inflation it is magnified by $e^{N_+}$, by reheating it has grown again by $a_R/a_{\rm end}$, and between reheating and the present day it grows by the ratio of the $T_R/T_0$  (corrected by the number of degrees of freedom). To solve the horizon problem, this scale should then be larger than the Hubble scale of the current universe $1/H_0$.
To obtain the order of magnitude
of the constraints, we consider
rapid reheating and take the reheating temperature to be
given by the potential energy at the end of inflation, and
hence set $a_R \sim a_{\rm{end}}$.
For simplicity we also assume that the ratio of $g_*$'s is 1.

Under the assumption that the
period of reheating lasts less than one Hubble time $T_R$ is given by the potential energy $V$
during inflation via $T_R \approx V^{1/4}$. Using the Friedmann equation, the Hubble scale $1/H_0$ is given by the
current energy density $\rho_0$ via
\be
\frac{1}{H_0} \, = \, \sqrt{3}\rho_0^{-1/2} M_{pl} \, .
\ee
In turn, $\rho_0$ can be re-expressed in terms of the temperature $T_0$:
\be
\rho_0 \, \approx \, T_0^4 \frac{T_{eq}}{T_0}{1\over \Omega_m}  \,
\ee
where $\Omega_m$ is the fraction of energy density in matter and
 $T_{eq}$ is the temperature at the time of equal matter and radiation, and we
have used the fact that the matter energy density today is larger than the radiation energy
density $T_0^4$ by the factor $T_{eq} / T_0$. Using $H = V^{1/2}/(\sqrt{3} M_{pl})$,
the condition (\ref{success}) then becomes
\be \label{success2}
e^{N_+} \, \simeq \, \frac{V^{1/4}}{(T_0 T_{eq})^{1/2}}\sqrt{\Omega_m} \sim \frac{V^{1/4}}{(T_0 T_{eq})^{1/2}}  \, .
\ee

\begin{figure}[t]
\begin{center}

\begin{tikzpicture}[x=0.5pt,y=0.5pt,yscale=-1,xscale=1]

\draw [color={rgb, 255:red, 8; green, 59; blue, 245 }  ,draw opacity=1 ]   (230.5,423) .. controls (230,367) and (247.5,311) .. (297.5,280) ;

\draw [line width=1.5]    (494.5,71) -- (297.5,280) ;

\draw [color={rgb, 255:red, 0; green, 0; blue, 0 }  ,draw opacity=1 ][line width=1.5]    (297.5,280) -- (296.5,424) ;

\draw [color={rgb, 255:red, 155; green, 155; blue, 155 }  ,draw opacity=1 ] [dash pattern={on 4.5pt off 4.5pt}]  (201.5,71) -- (494.5,71) ;

\draw [color={rgb, 255:red, 15; green, 17; blue, 243 }  ,draw opacity=1 ]   (296.5,424) .. controls (297.5,331) and (376.5,301) .. (428.5,281) ;

\draw [color={rgb, 255:red, 13; green, 14; blue, 226 }  ,draw opacity=1 ][line width=0.75]    (494.5,71) -- (428.5,281) ;

\draw [color={rgb, 255:red, 155; green, 155; blue, 155 }  ,draw opacity=1 ] [dash pattern={on 4.5pt off 4.5pt}]  (201.5,424) -- (296.5,424) ;

\draw  (187.5,281.22) -- (560.5,281.22)(201.56,44) -- (201.56,436) (553.5,276.22) -- (560.5,281.22) -- (553.5,286.22) (196.56,51) -- (201.56,44) -- (206.56,51)  ;

\draw (202,23) node   {$t$};
\draw (499,168) node   {$\lambda _{0}( t)$};
\draw (578,282) node   {$X$};
\draw (227,331) node   {$\lambda _{1}( t)$};
\draw (340,201) node   {$H^{-1}( t)$};
\draw (154,276) node   {$t_{R}$};
\draw (149,415) node   {$t_{i}$};
\draw (159,74) node   {$t_{0} \sim \frac{1}{H_{0}}$};

\end{tikzpicture}

\caption{Space-time sketch of inflationary cosmology. The vertical axis is time, the horizontal axis
represents physical distance. The inflationary period lasts from $t_i$ to $t_R$. Shown are the
Hubble radius $H^{-1}(t)$ and two length scales $\lambda_0(t)$ and $\lambda_1(t)$ (fixed wavelength
in comoving coordinates). For inflation to provide a possible explanation for the observed fluctuations
on large scales, the scale $\lambda_0(t)$ corresponding to the current Hubble horizon must originate inside
of the Hubble radius at the beginning of inflation. This leads to the condition (\ref{success}). The TCC, on
the other hand, demands that the length scale $\lambda_1(t)$which equals the Hubble radius at the end of inflation was never trans-Planckian. In the sketch, both conditions are marginally satisfied.}
\end{center}
\end{figure}

In the approximation of constant value of $H$ during inflation, the TCC condition (\ref{TCC}) can be
written in the form
\be \label{TCC2}
e^{N_+} \, < \frac{ M_{pl}}{H} \, .
\ee
The equation (\ref{success2}) for $N_+$ and the upper bound (\ref{TCC2}) on $N_+$ coming from the TCC are compatible only provided that the condition
\be
V^{3/4} \, < \sqrt{3}\, M_{pl}^2 (T_0 T_{eq})^{1/2}
\ee
is satisfied. Inserting the values of $T_0$, $T_{eq}$ and $M_{pl}$ we obtain
\be \label{Vbound}
V^{1/4} \, < 6\times\, 10^8 {\rm{GeV}}\sim 3\times 10^{-10}M_{pl} \, .
\ee
Note that this conclusion is independent of the assumption that quantum fluctuations during inflation are
the seeds for primordial structure formation. While we have used a potential $V$ to describe the energy density during inflation, our analysis holds for more general scenarios and Eq.~(\ref{Vbound}) can be interpreted as a bound on the energy density during the inflationary epoch.

We now add the assumption that quantum fluctuations of the inflaton are responsible for the origin of structure. In this case,
the power spectrum ${\cal{P}}$ of the curvature fluctuation ${\cal{R}}$ (see \cite{MFB}) is given by

\be \label{spec}
{\cal{P}}_{\cal{R}}(k) \,  = \, \frac{1}{ 8\pi^2\epsilon} \bigl( \frac{H(k)}{M_{pl}} \bigr)^2 \, ,
\ee
where $k$ is the comoving wavenumber of the fluctuation mode and $H(k)$ is the value of $H$
at the time when the mode $k$ exits the Hubble radius. The parameter $\epsilon$ determines
the deviation of the equation of state in the inflationary phase compared to pure de Sitter:
\be
\epsilon \, \equiv \, \frac{3}{2} \bigl( \frac{p}{\rho} + 1 \bigr) \, ,
\ee
where $p$ and $\rho$ are pressure and energy densities, respectively. For inflation to provide the source of structure in the Universe, we need \cite{Planck}
\be \label{value}
{\cal{P}}_{\cal{R}}(k) \, \sim \, 10^{-9} \, .
\ee
Combining (\ref{Vbound}) (\ref{spec}) and (\ref{value}) leads to an upper bound on $\epsilon$
\begin{align}\label{epsilon}
\epsilon \, \sim \,10^9 \frac{1}{8\pi^2 }\bigl( \frac{H(k)}{M_{pl}} \bigr)^2 \, \sim \, 10^9 \frac{V}{24\pi^2M_{pl}^4} \, < \, 10^{-31} \, .
\end{align}

Since the power spectrum of gravitational waves is given by
\be \label{valued}
{\cal{P}}_{h}(k) \, \sim \,   \bigl( \frac{H(k)}{M_{pl}} \bigr)^2 \, ,
\ee
the tensor to scalar ratio $r$ is given by
\be \label{result}
r \, = \, 16 \epsilon \, < \, 10^{-30} \, ,
\ee
where the factor $16$ comes from the different normalization conventions for the scalar and tensor
spectra. While the discussion above assumed that the inflaton dominated the scalar perturbations it is important to note that the TCC constrains the absolute amplitude of the primordial gravitational waves. The bound on $r$ therefore relies only on the TCC bound on the energy in Eq.~(\ref{Vbound}) and the observed amplitude of $\mathcal{P}_{\mathcal{R}}$. Allowing scalar perturbations from additional fields or a modified sound speed for the inflaton, for example, will not relax Eq.~(\ref{result}).

From (\ref{result}) we draw the conclusion that any detection of primordial gravitational waves on
cosmological scales would provide evidence for a different origin of the primordial gravitational wave
spectrum than any inflationary model consistent with the TCC. Note that a
number of cosmological scenarios alternative to inflation do predict  significant primordial tensor
modes on cosmological scales. One example is String Gas Cosmology which predicts both a
scale-invariant spectrum of cosmological perturbations with a slight red tilt \cite{NBV} and a
roughly scale-invariant spectrum of gravitational waves with a slight blue tilt \cite{BNPV2}.

Note that the above analysis applies not only to single field inflation, but also to multi-field
inflation. The conclusions only depend on the fact that the parameter $\epsilon$ is $\ll 1$ which is self-consistent with what we found. The constraint also
applies to {\it warm inflation} \cite{warm} models, models which can be consistent with
the de Sitter swampland conjecture \cite{Kamali}.

In this section, we have been general and have not assumed a slow-roll inflation.
In the following section we will study the consequences of the TCC for slow-roll inflation.

\section{Application to Slow-Roll Inflation}

The equation of motion of a canonically normalized scalar field $\phi$ in a homogeneous and isotropic metric of the form (\ref{FRW}) is
\be
{\ddot{\phi}} + 3 H {\dot{\phi}} + V' \, = \, 0 \, ,
\ee
where the prime indicates the derivative with respect to $\phi$. Here, we have assumed that there is no important coupling of $\phi$ to other matter fields during inflation. Thus, the analysis in this section applies to cold inflation but not to warm inflation. The Friedmann equation takes the form
\be
3 H^2 M_{pl}^2 \, = \frac{1}{2} {\dot{\phi}}^2 + V \, ,
\ee

In the case of single field slow-roll inflation, the slow-roll parameter $\epsilon$ is
\be \label{SRepsilon}
\epsilon \, \simeq \, \frac{M_{pl}^2}{2} \bigl( \frac{V^{'}}{V} \bigr)^2 \, .
\ee
The slow-roll equation of motion is
\be
3 H {\dot{\phi}} \, = \, - V^{'} \, .
\ee

The field range $\Delta \phi$ which the inflaton field $\phi$ moves during the period of
inflation is given by
\be
|\Delta \phi| \, \simeq \, |{\dot{\phi}} \Delta t| \, ,
\ee
where $\Delta t$ is the time period of inflation. We will show that $\Delta \phi$ is very small
compared to the Planck mass. In this case, it is self-consistent to assume that $H$ and ${\dot{\phi}}$
are constant. In this case using the TCC
\be
\Delta t \, = \, H^{-1} N \, < \, H^{-1} {\rm{ln}} \bigl( \frac{ M_{pl}}{H} \bigr) \, ,
\ee
and, making use of (\ref{SRepsilon}), the field range becomes
\begin{align} \label{range}
 |\Delta \phi| \, < &\sqrt{2}\epsilon^{1/2}  {\rm{ln}} \bigl( \frac{M_{pl}}{H} \bigr) M_{pl}\nonumber \\
 <&{10^{9/2}V^{1/2} \over M_{pl}}\ln( {M^2_{pl}\over {\sqrt V}})\nonumber \\
 < &10^{-13}M_{pl}\, ,
 \end{align}

where in the last inequality we used the monotonicity of $[x\ln(1/x)]$ for $ x<e$ to substitute the upper bound on $V $ from \eqref{Vbound}.
As first studied in \cite{Kung}, in the case of {\it large field inflation}, i.e. $|\Delta \phi| \gg M_{pl}$, the
inflationary slow-roll trajectory is a local attractor in initial condition space, even taking into account metric
fluctuations \cite{Hume,East}. Small field inflation as is the case here, on the other hand, is not an attractor in initial
condition space, as reviewed in \cite{Goldwirth}. If the field range for slow-roll inflation is constrained
by the TCC conjecture to obey (\ref{range}), then the inflationary scenario is faced with an initial
condition problem. The expected initial field velocity is
\be
{\dot{\phi}}_i^2 \, \sim \, V \,
\ee
and hence
\be
\frac{{\dot{\phi}}_{SR}}{{\dot{\phi}}_{i}} \, \sim \, \epsilon^{1/2}< 10^{-15}\, ,
\ee
and it takes fine tuning of the initial velocity in order to be sufficiently close to the slow-roll trajectory.

In the following we propose a model which can consistently explain observations, including the observational value of the tilt. We consider an inverted parabola potential $V(\phi)=V_0-{|V''|\phi^2/2}$ over a small field range $[\phi_i,\phi_f]$ such that $\delta V/V\ll1$ over the field range. Given the smallness of $\epsilon$ from \eqref{epsilon}, we have
\begin{align}\label{M1}
M_{pl}^2\frac{V''}{V}\simeq \frac{n_s-1}{2},
\end{align}
where $n_s=1+2\eta-6\epsilon$ is the tilt parameter and $\eta=M_{pl}^2V''/V$ is the second slow-roll parameter. This fixes $V''$ from observation. From the equations \eqref{spec} and \eqref{SRepsilon} we find
\begin{align}\label{M2}
\frac{V_0}{12\pi^2M_{pl}^2\phi_{CMB}^2}\simeq\mathcal{P}(k)(\frac{n_s-1}{2})^2.
\end{align}
where $\phi_{CMB}$ is the value of the field when the modes on CMB scales exited the Hubble horizon. Assuming $H$ remains almost constant $H\simeq \sqrt{V_0/3M_{pl}^2}$ during the slow-roll inflation, one can show
\begin{align}\label{e-fold2}
 \ln(\frac{\phi_2}{\phi_1})&\simeq|\eta|N(\phi_1\rightarrow \phi_2),
\end{align}
where $N(\phi_1\rightarrow \phi_2)$ is the number of e-folds accrued as $\phi$ goes from $\phi_1$ to $\phi_2$. If we plug $\phi_1=\phi_{CMB}$ and $\phi_2=\phi_{f}$ into the above identity and use the equation \eqref{success2}, we find
\begin{align}\label{efold}
    \frac{\phi_{f}}{\phi_{CMB}}=e^{|\eta|N_+}\simeq\frac{V_0^\frac{|\eta|}{4}\Omega_m^\frac{|\eta|}{2}}{(T_0T_{eq})^\frac{|\eta|}{2}}.
\end{align}
Plugging $\phi_{CMB}$ from \eqref{M2} into \eqref{efold} leads to
\begin{align}\label{end}
    \phi_{f}\simeq \frac{V_0^{\frac{1}{2}+\frac{|\eta|}{4}}\Omega_m^\frac{|\eta|}{2}}{M_{pl}^2(T_0T_{eq})^\frac{|\eta|}{2}12^\frac{1}{2}\mathcal{P}^\frac{1}{2}\pi|\eta|}\simeq 3.9 \times10^5\cdot \left(\frac{V_0}{M_{pl}}\right)^{0.505},
\end{align}
where in the last step we substituted $|\eta|\simeq0.02$, $\mathcal{P}\simeq2\cdot10^{-9}$, $T_0\simeq3K$, $T_{eq}\sim10^4K$, and $\Omega_m\simeq0.3$. This fixes the end of the field range $\phi_f$ in terms of the energy scale $V_0$. The only free parameters left are $V_0$ and $\phi_i$.

Now we impost the TCC for the slow-roll trajectory to find a constraint in terms of $\phi_i$ and $V_0$. Plugging $\phi_1=\phi_i$ and $\phi_2=\phi_{CMB}$ in \eqref{e-fold2} gives
\begin{align}\label{fieldefold}
\phi_{CMB}\simeq \phi_i e^{\frac{1-n_s}{2}N_-},
\end{align}
where $N_{-}$ is the number of e-folds accrued before the modes on CMB scales exit the horizon. The total number of e-foldings is $N_{total}=N_{-}+N_{+}$. From \eqref{success2} we find
\begin{align}
e^N\simeq e^{N_-}\frac{V^{1/4}}{(T_0 T_{eq})^{1/2}}.
\end{align}
On the other hand, from the TCC, we know that the total number of e-folds is bounded by $e^N<M_{pl}/H$. Using \eqref{fieldefold}, this can be expressed as
\begin{align}\label{ineqfield}
(\frac{\phi_{CMB}}{\phi_i})^\frac{2}{1-n_s}<\frac{(3M_{pl}^4T_0T_{eq})^\frac{1}{2}}{V_0^\frac{3}{4}}.
\end{align}
Plugging $\phi_{CMB}$ from \eqref{M2} in \eqref{ineqfield} with $\mathcal{P}\simeq 2\cdot10^{-9}$, $n_s\simeq 0.96$, $T_0\simeq 3K$, and $T_{eq}\sim 10^4 K$ leads to
\begin{align}\label{TCCM}
(\frac{V_0}{M_{pl}^4})^{1.03}<6.6\times10^{-12}\cdot(\frac{\phi_i}{M_{pl}})^2.
\end{align}

The above inequality is necessary for the potential to be consistent with the TCC, but it is not sufficient. This is because a potential is consistent with the TCC if the inequality \eqref{TCC} is satisfied for every expansionary trajectory, not just one particular trajectory.

For energy scale $V_0^{1/4}=10^{-10}M_{pl}$ the potential $V_0(1-0.02\phi^2)$ defined over the field range $[\phi_i,\phi_f]=[9.7\times 10^{-16}M_{pl},2.4\times 10^{-15}M_{pl}]$ satisfies all the criteria \eqref{Vbound}, \eqref{end}, and \eqref{TCCM}. These criteria were imposed by observation and consistency with \eqref{TCC} for the slow-roll trajectory. By numerical analysis, we further verified the consistency of this potential with the inequality \eqref{TCC} for every expansionary trajectory. This is an example of a simple potential that can explain the observation and be consistent with the TCC at the same time, however, due to its short field range, it suffers from the fine-tuning problem.

\section{Conclusions and Discussion} \label{conclusion}

We have studied the implications of the recently proposed the TCC for inflationary cosmology. Demanding that the TCC holds and that the largest scales that we currently probe in cosmology are sub-Hubble at the beginning of the inflationary phase (a necessary condition for the causal generation mechanism of fluctuations of inflationary cosmology to work) leads to an upper bound on the energy scale of inflation which is of the order of $10^9 {\rm{GeV}}$. Demanding that the amplitude of the cosmological perturbations agree with observations then leads to an upper bound on the generalized slow-roll parameter $\epsilon$ of the order of $\epsilon < 10^{-31}$. As a consequence, the tensor to scalar ratio is predicted to be smaller than $10^{-30}$. A detection of primordial gravitational waves via B-mode polarization, pulsar timing measurements or direct detection, assuming the TCC, would then imply that the source of these gravitational waves is {\it not} due to quantum fluctuations during inflation.

The above conclusions are independent of any assumptions on the possible single-field nature of inflation. If we then specialize the discussion to the case of single field slow-roll inflation with a canonically normalized inflaton field, we find that the range $\Delta \phi$ which the inflaton field traverses during the inflationary phase is of the order of $\epsilon^{1/2} M_{pl}$. This raises an initial condition problem for most of the models since the expected field velocity is much larger than the field velocity along the slow-roll trajectory.

We proposed an inverted parabola potential as a simple example that is consistent with the TCC and can explain the observation at the same time.

\section*{Acknowledgement}

\noindent We thank the Simons Center for Geometry and Physics for hospitality during the 2019 Summer Workshop on String Theory during which the ideas for this project were developed. RB thanks the Institute for Theoretical Physics of the ETH for hospitality. The research at McGill is supported in part by funds from NSERC and from the Canada Research Chair program. ML is supported by DOE DE-SC0017848. The research of CV is supported in part by the NSF grant PHY-1719924 and by a grant from the Simons Foundation (602883, CV).


\subsection*{{\sc \large Supplementary Material}}
\twocolumngrid

\vspace{-3mm}
\section*{A.~Size of the PBH}
\vspace{-2mm}

The Schwarzschild radius of a black hole is given by
\begin{equation}\begin{aligned}
r_{\rm BH}=\frac{2GM_{\rm BH}}{c^2}\simeq 4.5{\rm cm}\left(\frac{M_{\rm BH}}{5M_\oplus}\right)~.
\end{aligned}\end{equation}
In Figure \ref{fig1} we provide an exact scale image of a $5M_\oplus$ PBH. The associated DM halo however extends to the stripping radius $r_{t,\odot}\sim8$AU, this would imply a DM halo which extends roughly the distance from Earth to Saturn (both in real life and relative to the image).

\onecolumngrid

\begin{figure}[h!]
\vspace{1mm}
\includegraphics[height=9.35cm]{Fig9cmBH}
	\caption{Exact scale (1:1) illustration of a $5M_\oplus$ PBH. Note that a  $10M_\oplus$ PBH is roughly the size of a ten pin bowling~ball.}\label{fig1}
\end{figure}

\newpage
\twocolumngrid

\section*{B.~Dark Matter Freeze-in via a $Z'$}

To demonstrate a viable freeze-in model with potentially detectable indirect detection signals we consider a specific implementation in which the Standard Model is supplemented by two new ``hidden sector'' states, the DM particle $\chi$ with mass  $m$ and a mediator $Z'$ with mass $M_{Z'}$. We assume a sizeable interaction $gZ'\overline{\chi}\chi$ for $g\sim\mathcal{O}(10^{-2})$, with mass ordering $M_{Z'}>2m$, and that the hidden sector only couples to the Standard Model fermions $f$ (with mass $m_f$) via $\lambda Z'\overline{f}f$, which involves a tiny coupling  $\lambda\sim10^{-12}$.

Here we relabel the mediator $Z'$ (from $\phi$ in the main text) and highlight that it can naturally arise from a  kinetic mixing portal \cite{Holdom:1985ag} involving a new U(1)$'$ gauge symmetry via the Lagrangian term $\Delta{\mathcal L}=\delta B_{\mu\nu}X^{\mu\nu}$
where $X$ is the $Z'$ field strength and $B$ is the hypercharge field strength, where $\delta$  parameterises the mixing induced by states  (with masses $M$ and $M'$)  charged under both U(1)$'$ and hypercharge $\delta\propto \frac{g_Yg_X}{32\pi^2}{\rm log}\left(M^2/M'{}^2\right)$. In motivated models this can lead to couplings appropriate for freeze-in \cite{Dienes:1996zr,Gherghetta:2019coi} and one can identify the freeze-in coupling $\lambda\simeq\delta g_{Y}$ where $g_Y$ are Standard Model fermion hypercharge couplings. For further details on $Z'$ freeze-in portals see e.g.~\cite{Chu:2013jja,Chu:2011be,Elahi:2014fsa,Heeba:2019jho,Foldenauer:2018zrz}.

The initial  $Z'$ and $\chi$ abundances are negligible but populations are generated via the portal interaction $\lambda Z'\overline{f}f$.
 The $Z'$ abundance evolves according to \cite{Hall:2009bx}
\begin{equation}\begin{aligned}
\dot n_{Z'}+3n_{Z'} H\simeq  \int & {\rm d}\Pi_{Z'}  {\rm d}\Pi_{\overline{f} }{\rm d}\Pi_{f}  (2\pi)^4\\
 &\delta^4(p_{Z'}-p_f-p_{\overline{f}})|\mathcal{M}|^2_{Z'\rightarrow \overline{f}f}f_{Z'}^{\rm eq}
\end{aligned}\end{equation}
where $f_{Z'}^{\rm eq}={\rm exp}(-E_\chi/T)$ is the equilibrium distribution.

The generated abundance of $Z'$ is  \cite{Hall:2009bx}
\begin{equation}\begin{aligned}
Y_{Z'}=\frac{135}{8\pi^3g_*^{S}\sqrt{g_*^{\rho}}}\left(\frac{M_{\rm Pl}\Gamma_{Z'}}{M_{Z'}^2}\right)
\end{aligned}\end{equation}
where $\Gamma_{Z'}$ is the partial width involving bath states.

Since we take $M_{Z'}>2m$ and $g\sim\mathcal{O}(10^{-2})$ any $Z'$ produced decays promptly to DM pairs and thus the resulting DM relic density is $Y_\chi\approx2Y_{Z'}$ for $\lambda\ll g$. This can be related to the observed quantity
\begin{equation}\begin{aligned}
\Omega_{\rm DM}&\simeq
 0.2 \left(\frac{m}{100~{\rm GeV}}\right)
\left(\frac{\lambda}{6\times10^{-12}}\right)^2
\left(\frac{10~{\rm TeV}}{M_{Z'}}\right)~,
\end{aligned}\end{equation}
where we have parameterised $\Gamma_{Z'}\simeq \lambda^2M_{Z'}/8\pi$.

The feeble coupling $\lambda\ll1$ is characteristic of freeze-in models and is actually required to satisfy the consistency condition $Y_{\rm DM}(T)\ll Y_{\rm DM}^{\rm eq}$(T). If $\lambda$ is too large the hidden sector enters equilibrium ($Y_{\rm DM}\sim Y_{\rm DM}^{\rm eq}$) and then subsequently freeze-out dynamics sets the DM relic density.
However, since $\lambda$ is very small, direct and indirect detection is typically challenging. For the $Z'$ model the s-wave annihilation cross-section for $m_f\ll m<M_{Z'}$ is \cite{Alves:2015pea}
\begin{equation}\begin{aligned}
\sigma_0\simeq\sum_f \frac{N_c \lambda^2}{2\pi }
 \left(\frac{2g_{\chi v}^2   m^2}{(M_{Z'}^2-4m^2)^2 } + g_{\chi a}^2 \frac{m_f^2}{M_{Z'}^4}\right),
\end{aligned}\end{equation}
where the sum runs over Standard Model fermions and $g_{\chi a}$ and $g_{\chi v}$ are the axial and vector couplings to $\chi$. Then eq.~(\ref{xc}) is obtained from the above in the limit $ m\ll M_{Z'}$. In the main text we highlight that if $f_{\rm PBH}$ is relatively large, indirect detection may be a viable route to search for this class of models in certain parameter ranges.

\onecolumngrid
\vspace{3mm}


\section*{Acnowledgments}
We thank the Collider Accelerator Department and the RHIC/AGS Operation Groups. We also would like to thank A.~Bazilevsky, B.Z.~Kopeliovich, and M.~Krelina for useful discussions. B.Z.~Kopeliovich read the manuscript and made valuable comments. This work was supported by Brookhaven Science Associates, LLC under Contract No.~DE-AC02-98CH\,10886 with the U.S. Department of Energy. Funding was also provided from the RIKEN BNL Research Center. N.H.B. is grateful for partial support from the School of Mathematics.
\end{acknowledgments}