\section{Introduction}

Since the introduction of Unruh-DeWitt particle detectors in quantum field theory \cite{Unruh1976db,dewitt1979quantum}, they have been employed to probe a plethora of features of quantum fields and of spacetimes, stretching from picking out properties of quantum states \cite{Ng2021enc} to identifying topological aspects of underlying backgrounds \cite{Ng2017iqh}. On a globally-hyperbolic spacetime with a timelike boundary, they are particularly sensible to the choice of an underlying boundary condition for the matter fields: phenomena such as the divergence of physical observables \cite{deSouzaCampos2021awm,Namasivayam2022bky} and the anti-Hawking effect \cite{Henderson2019uqo, DeSouzaCampos2020ddx} depend crucially on the boundary condition set. This is consonant with the fact that different boundary conditions are associated to different quantum states, to different dynamics, to different {\em physics} \cite{Ishibashi2003jd}.

In this arena we have recently shown that within free, scalar quantum field theories on static, curved spacetimes with a timelike boundary one may take into account generalized $(\gamma, v)$-Robin boundary conditions in the construction of physically meaningful states \cite{Campos2022byi}. These arise whenever the underlying dynamics can be reduced thanks to the background isometries to a second order differential equation on the half line and they are characterized by two data: a fixed parameter $\gamma\in\mathbb{R}$ and a so-called secondary solution $v$. For each admissible pair $(\gamma, v)$, one can construct a two-point function both for the ground state and also for arbitrary thermal states. All these correlation functions are physically sensible since, per construction, they enjoy the local {\em Hadamard} property, see \cite{Khavkine:2014mta}. In \cite{Campos2022byi}, we have discussed the example of a wave propagating on the two-dimensional half-Minkowski spacetime showing that an Unruh-DeWitt detector can discern not only the value of $\gamma$, but also the choice of secondary solution $v$. Here, we elaborate on this idea and we provide, in full detail, a full-fledged example that demonstrates how the hidden freedom in the mode expansion of the scalar field, which is associated to the choice of $v$, influences the rate of transition of an Unruh-DeWitt particle detector.

In this work, we consider a real, free, scalar field $\Psi$ with mass $m_0$ on $M\equiv\text{AdS}_2 \times \mathbb{S}^2$ spacetime. This background approximates the near horizon geometry of an extremal black hole with unit charge and it is known as the Bertotti-Robinson solution of Einstein-Maxwell equations, see e.g. \cite{Conroy2021aow} but also \cite{Bertotti, Robinson}. For $t\in\mathbb{R}$ and $r\in(0,\infty)$, its line-element can be written as
\begin{equation}
\label{eq: line-element}
  ds^2 = \frac{1}{r^{2}}\left(-dt^2 + dr^2 + r^2 d\mathbb{S}^{2}(\theta,\varphi)\right),
\end{equation}
where $d\mathbb{S}^{2}(\theta,\varphi)$ is the standard line-element on the unit $2$-sphere. Observe that the manifold possesses a timelike conformal boundary at $r=0$. On top of $M$ we consider a real scalar field $\Psi:\text{AdS}_2\times \mathbb{S}^2\rightarrow\mathbb{R}$ whose dynamics is ruled by the Klein-Gordon equation
 \begin{equation}
   \label{eq: KG}
 P\Psi : = (\Box - m_0^2)\Psi=0,
 \end{equation}
where $\Box$ is the D'Alembert wave operator built out of Equation \eqref{eq: line-element}, while $m_0$ is the mass parameter.

In Section \ref{sec: The Klein-Gordon field}, we recollect the main ingredients necessary to construct the two-point functions for ground and thermal states for $\Psi$. The procedure follows the prescription outlined in \cite{Campos2022byi} and it generalizes the results obtained in \cite{Dappiaggi2016fwc} by considering two different families of generalized $(\gamma, v)$-boundary conditions at the conformal boundary, dubbed $(\gamma, v_\kappa)$ with $\kappa=1,2$. The analysis for $\kappa=2$ coincides with that of \cite{Dappiaggi2016fwc}, while for $\kappa=1$ it yields novel two-point functions. Nonetheless, since also in this case the overall analysis is structurally identical to that for $\kappa=2$, we do not dwell into many details and we limit ourselves to sketching the main steps of the construction.

Subsequently, in Section \ref{sec: The transition rate} we obtain an explicit expression for the transition rate of a static Unruh-DeWitt detector interacting either with the ground or with a thermal state of a Klein-Gordon field. We scrutinize the response of the detector to grasp if and how it is affected by the choice of boundary condition. In particular, we focus on the consequences of the freedom in selecting a secondary solution in the construction of Section \ref{sec: The Klein-Gordon field}. In Section \ref{sec: Summary}, we summarize the main results of this work.

\section{Dynamics and Boundary Conditions}
\label{sec: The Klein-Gordon field}

Working with the coordinates of Equation \eqref{eq: line-element}, a real Klein-Gordon field $\Psi$ that solves Equation \eqref{eq: KG} can be decomposed as
\begin{equation*}
	\Psi(t,r,\theta,\varphi) = e^{-i\omega t} \psi(r) Y_\ell^m (\theta,\varphi),
\end{equation*}
where $Y_\ell (\theta,\varphi)$ are the spherical harmonics on the $2$-dimensional unit sphere. Setting $\lambda=-\ell(\ell+1)$, the function $\psi(r)$ solves the radial equation
\begin{align}
\label{eq: the radial eq}
\mathbf{L} \psi(r) :=
&\frac{d^2 \psi(r)}{dr^2} + \left[ \omega^2 + \frac{\lambda-m_0^2}{ r^2} \right]\psi(r)= 0.
\end{align}

\subsection{The radial equation}

\noindent In the following we assume that $m_0\geq 0$ and we introduce the auxiliary quantities
\begin{align}
\label{eq: nu and p}
  &\nu := \frac{1}{2}\sqrt{1-4\lambda  + 4m_0^2 } \quad \text{and}\quad
  p:= \sqrt{\omega^2},
\end{align}
subject to the constraint $\nu > 0, \,\forall \ell\geq0$. Equation \eqref{eq: the radial eq} is a Sturm-Liouville problem with eigenvalue $\omega^2$, which admits a basis of solutions written in terms of Bessel functions of first and second kind:
\begin{subequations}
\label{eq: basis y1 and y2}
\begin{align}
  y_1(r) &=\sqrt{r}  J_\nu (p r), \\
  y_2(r) &=\sqrt{r}  Y_\nu (p r).
\end{align}
\end{subequations}

Following \cite{Zettl:2005}, it is convenient to work on the Hilbert space $L^2((0,\infty), dr)$ and allow $p$ to take a priori any complex value. Hence, according to Weyl's endpoint classification $r\to\infty$ is a limit point, see also \cite{Dappiaggi2016fwc} for a succinct summary of the nomenclature. As a matter of fact, assuming that $\Imag p\neq 0$, the most general solution of Equation \eqref{eq: the radial eq} lying in $L^2((c,\infty),dr)$, $\forall\, c\in (0,\infty)$ can be written in terms of Hankel functions of first and second kind:
\begin{align*}
  \psi_{\scriptscriptstyle{\infty}}(r)  =\sqrt{r}  \left[ H^{\scriptscriptstyle{(1)}}_\nu(p r) \Theta(\Imag p)  +  H^{\scriptscriptstyle{(2)}}_\nu(p r) \Theta(-\Imag p) \right].
\end{align*}
On the other hand, it turns out that $r=0$ is a limit circle point if $\nu\in(0,1)$, since both $y_1,y_2\in L^2((0,c^\prime),dr)$, $\forall\, c^\prime\in(0,\infty)$, as one can infer from the following asymptotic expansion close to the origin:
\begin{align*}
  |y_1(r)|^2 &\overset{r\rightarrow 0}{\sim}r^{1+2\nu} , \\
  |y_2(r)|^2 &\overset{r\rightarrow 0}{\sim}r^{1-2\nu} .
\end{align*}
Observe that
\begin{equation}\label{Eq: mass range}
	\nu  \in[0,1) \iff \ell = 0 \text{ and } m_0^2\in\left[-\frac{1}{4}, \frac{3}{4}\right).
	\end{equation}
Still following the theory of Sturm-Liouville, this entails that only the $\ell=0$ mode calls for a boundary condition at $r=0$, provided that the mass lies in the range individuated in Equation \eqref{Eq: mass range}. When $\nu\geq1$, no boundary condition needs to be imposed at both ends and since this case is not of interest in this work, we shall not consider it further.

\subsection{Generalized $(\gamma, v)$-boundary conditions}

At the limit circle endpoint, $r=0$, we impose generalized $(\gamma,v)$-boundary conditions as introduced in \cite{Campos2022byi}. These identify self-adjoint extensions of the radial operator $\mathbf{L}$ in $L^2((0,\infty),dr)$ and they are fully characterized in terms of a parameter $\gamma\in\mathbb{R}$, of the principal solution, $u$, and of a secondary solution, $v$, which we introduce in the following. The principal solution at $r=0$ reads
\begin{align*}
&u := \sqrt{r} J_\nu (p r)  = y_1(r).
\end{align*}
This is unambiguously identified by the condition that, for any solution $v$ of Equation \eqref{eq: the radial eq} such that $v\neq \lambda u$, $\lambda\in\mathbb{C}$, then $\lim\limits_{r\to 0}\frac{u}{v}=0$.  On the contrary choosing a secondary solution at $r=0$, linearly independent from $u$, is fully arbitrary. In the following we consider two possible choices $v_\kappa$, $\kappa=1,2$:
\begin{subequations}
	\label{eq: sec sol v1 v2}
\begin{align}
v_1 :=& \,p \sqrt{r} Y_\nu (p r) \\  v_2 :=& - p^{2\nu}  \sqrt{r}  J_{-\nu} (p r).
\end{align}
\end{subequations}
In terms of the basis given in Equation \eqref{eq: basis y1 and y2}, they can be written as
\begin{align*}
v_1 :=& \, p \, y_2(r), \\  v_2 :=& -p^{2\nu}  \left[ \cos(\pi\nu ) y_1(r) - \sin(\pi\nu )y_2(r) \right].
\end{align*}
For $ \kappa\in\{1,2\}$ and $\gamma \in \mathbb{R}$, it follows that the solution
\begin{align*}
  \psi_{\kappa} := \cos(\gamma)u -\sin(\gamma) v_{\kappa} ,
\end{align*}
satisfies
\begin{align}
\label{eq: robin bc}
  &\cos(\gamma)W_r[\psi_{\kappa},u] - \sin(\gamma)W_r[\psi_{\kappa},v_{\kappa}]  = 0 ,
\end{align}
where $W_r$ denotes the Wronskian $W_r[\psi_{\kappa},u]:=\psi_{\kappa}\partial_r u- u \partial_r\psi_\kappa$.
If Equation \eqref{eq: robin bc} holds true, we say that $\psi_\kappa$ satisfies a $(\gamma,v_\kappa)$-boundary condition at $r=0$.

\subsection{The radial Green function}

In the next step in our analysis we construct the Green function of the radial equation. To this end, it is convenient to rewrite the solution at infinity as a linear combination of $u$ and $v_\kappa$:
\begin{align*}
\psi_{\scriptscriptstyle{\infty}} &= a_\kappa u + b_\kappa v_{\kappa}.
\end{align*}
A direct computation yields
\begin{align*}
& a_1 = 1,\\
& a_2 = 1 +  i\sign(\Imag p) \cot(\pi\nu ),\\
& b_\kappa= \frac{i\sign(\Imag p)}{\Lambda_\kappa},\quad\kappa=1,2,
\end{align*}
where $\Lambda_1 := p$, and $\Lambda_2 := p^{2\nu}\sin(\pi\nu )$. Since $W_r[J_\nu(pr),Y_\nu(pr)]=\frac{2}{\pi r}$, it descends
\begin{align*}
W_r[\psi_{\kappa}, \psi_{\scriptscriptstyle{\infty}}] &= ( b_{\kappa}\cos(\gamma) + a_\kappa\sin(\gamma) )   \frac{2\Lambda_\kappa}{\pi }.
\end{align*}
Therefore, following \cite{greenBook}, if $\mathbf{L}$ is the operator as per Equation \eqref{eq: the radial eq}, the radial Green function, implicitly defined by
\begin{align*}
(\mathbf{L}\otimes \mathbbm{1}) G_\kappa(r,r') = ( \mathbbm{1}\otimes\mathbf{L}) G_\kappa(r,r') =\delta(r-r'),
\end{align*}
reads
\begin{align}
\label{eq: radial Green bessel kappa}
G_\kappa(r,r')
         &= \frac{\pi}{2 \Lambda_\kappa} \frac{ \psi_{\kappa}(r_<) \psi_{\scriptscriptstyle{\infty}}(r_>)}{ ( b_{\kappa}\cos(\gamma) + a_\kappa\sin(\gamma) )},
\end{align}
where $\psi_{\kappa}(r_<) \psi_{\scriptscriptstyle{\infty}}(r_>)$ is a shortcut notation for the integral kernel $\Theta(r-r^\prime)\psi_{\kappa}(r)\psi_{\scriptscriptstyle{\infty}}(r^\prime)+\Theta(r^\prime-r)\psi_{\scriptscriptstyle{\infty}}(r)\psi_{\kappa}(r^\prime)$.

\subsection{Symmetries of the radial Green function}

Seeing $G_\kappa$, defined in Equation \eqref{eq: radial Green bessel kappa}, as a function either of $\lambda = p^2$ or $\sqrt{\lambda} = p$, we find that
\begin{align}
&G_1(\overline{p}) = \overline{ G_1(p)}, \text{ but } G_1(\overline{\lambda} ) \neq \overline{ G_1(\lambda)}; \label{eq: radial Green bessel kappa complex conjugation G1} \\
&G_2(\overline{\lambda} ) = \overline{ G_2(\lambda)}, \text{ but } G_2(\overline{p}) \neq \overline{ G_2(p)}. \label{eq: radial Green bessel kappa complex conjugation G2}
\end{align}
\noindent We remark that in the subsequent analysis, particularly in Equations \eqref{eq: radial Green bessel kappa complex conjugation G1} and \eqref{eq: radial Green bessel kappa complex conjugation G2}, the map $\lambda \mapsto \overline{\lambda}$ corresponds to, respectively, $p\mapsto \overline{p}$ and $p\mapsto -\overline{p}$. In this last case, one chooses once and for all $\sign(\Imag p )\in\{-1,+1\}$.

\subsection{Bound States}

When choosing a boundary condition in Equation \eqref{eq: robin bc}, we consider as {\em admissible} the values of $\gamma$ for which there does not exist $\omega_b\in\mathbb{C}$ such that $\lim\limits_{\omega\rightarrow\omega_b}G_\kappa$ diverges. In other words, at a physical level, we rule out those boundary conditions for which $G_\kappa$ has {\em bound state} frequencies $\omega_b$. The case for $\kappa=1$ is analogous to the computation on Minkowski spacetime in spherical coordinates as in \cite[Pg.66]{DeSouzaCampos2022wsp}, while the case for $\kappa=2$ has been studied in detail in \cite[Pg.10]{Dappiaggi2016fwc}. Here, we report the final result: neither $G_1$ nor $G_2$ possesses bound state frequencies for $\gamma\in\left[0,\frac{\pi}{2}\right]$. Henceforth we restrict our attention to this interval. Observe that $\gamma=0$ corresponds to the usual Dirichlet boundary condition. It is tempting to claim that, if $\gamma=\frac{\pi}{2}$, we are instead considering a Neumann boundary condition. This is indeed an admissible nomenclature, but the reader should keep in mind that the solution at $r=0$ of the radial equation depends only on the choice of the secondary solution, contrary to the case $\gamma=0$ in which only $u$ plays a r\^{o}le. Hence there is no universal way to select a Neumann boundary condition, inasmuch as there does not exist a unique criterion to single out a secondary solution for Equation \eqref{eq: the radial eq}.

\subsection{The resolution of the identity}

In order to construct the two-point correlation function of the Klein-Gordon field $\Psi$ with arbitrary $(\gamma,v)$-boundary condition, a last necessary ingredient is the so-called resolution of the identity for the radial equation. Following \cite[Ch.7]{greenBook}, this reads
\begin{align}
\delta(r,r') & = -\frac{1}{2\pi i}\oint\limits_{C_\kappa^\infty} G_\kappa(x,x')dp^2  \nonumber \\ & = \int\limits_{-\infty}^{\infty} p\, dp \, \widetilde{\psi}_\kappa(x,x^\prime)\label{eq: resolution identity bessel example 2}
\end{align}
where the contours differ if $\kappa=1,2$ and they are illustrated in \cite[Fig.1]{Campos2022byi}, while the integrand $\widetilde{\psi}_\kappa(x,x^\prime)$ has been computed in \cite{DeSouzaCampos2022wsp}. Since its explicit form is not of relevance in this work, we shall not report it, although, for later convenience, we highlight that, whenever we restrict the attention to real, positive frequencies $p=\omega$,

\begin{align*}
\widetilde{\psi}_\kappa(r,r') &=\frac{ \psi_{\kappa}(r_<)\psi_{\kappa}(r_>) }{\mathcal{N}_\kappa},
\end{align*}
where
\begin{align*}
\mathcal{N}_1 &: =2(\cos^2(\gamma)+ p^2 \sin^2(\gamma)),\\
\mathcal{N}_2 &: =2(\sin^2(\gamma) p^{4 \nu } + \sin(2\gamma) \cos (\pi  \nu ) p^{2 \nu }+\cos^2(\gamma)).
\end{align*}

\subsection{Two-point correlation functions}

In the preceding analysis we have individuated all the ingredients necessary to construct the two-point correlation function both of a ground and of a thermal state for a real Klein-Gordon field whose mass is such that one can endow it with physically admissible boundary conditions, {\it i.e.}, $m_0^2\in\left[-\frac{1}{4}, \frac{3}{4}\right)$. As in the previous sections, we shall not dwell into a detailed construction which has been already outlined in \cite{DeSouzaCampos2022wsp}, but we limit ourselves to recalling the final result

\begin{widetext}
\begin{align}
\label{eq: two point function ground and thermal psi}
\omega_{2,\kappa}(t,r,\theta,\varphi,t',r',\theta',\varphi') =\lim_{\varepsilon\to 0^+}\sum\limits_{\ell=0}^\infty \sum\limits_{m=-\ell}^\ell \int_0^\infty d\omega\, T_{\beta,\varepsilon}(t-t')\widetilde{\psi}_\kappa(r,r') Y_\ell^m  (\theta,\varphi)\overline{Y_\ell^m  (\theta',\varphi')},
\end{align}
\end{widetext}
\noindent where
\begin{equation}
\label{eq: two point function ground and thermal time function}
T_{\beta,\varepsilon}(t-t') := \left[ \frac{e^{-i \omega (t-t' - i \varepsilon)}}{1-e^{-\beta \omega}} + \frac{e^{i \omega (t-t' + i \varepsilon)}}{e^{\beta\omega}-1}\right].
\end{equation}
For finite $\beta$, Equation \eqref{eq: two point function ground and thermal psi} together with Equation \eqref{eq: two point function ground and thermal time function} characterize a thermal state at temperature $1/\beta$. In the limit $\beta\rightarrow\infty$, we obtain instead the two-point correlation function of the ground state:
\begin{equation}
T_{\infty,\varepsilon}(t-t') := e^{-i \omega (t-t' - i \varepsilon)}.
\end{equation}
In the next section we shall use this class of two-point correlation functions to show that the choice of different secondary solutions associated to Equation \eqref{eq: the radial eq} is not a mere mathematical exercise, but it bears strong physical consequences.

\section{The transition rate}
\label{sec: The transition rate}

In this section we consider an Unruh-deWitt detector, namely a spatially localized two-level system that interacts with the underlying scalar field via a monopole type interaction. Here we are interested in computing and analyzing the transition rate when the detector follows static trajectories on $\text{AdS}_2 \times \mathbb{S}^2$. The detector can be initially either in the ground-state $|0\rangle$ or in the excited state $|\Omega\rangle$, see \cite{Unruh1976db, dewitt1979quantum}. This information is codified in the sign of its energy gap: $\Omega>0$ corresponds to excitations, $\Omega<0$ corresponds to de-excitations.

We let the detector interact for an infinite proper time with the Klein-Gordon field, whose initial state is either the ground state or a thermal state at inverse-temperature $\beta$ with $(\gamma,v_\kappa)$-boundary conditions, see Equation \eqref{eq: two point function ground and thermal psi}. In addition, we pa\-ram\-e\-triz\-e the trajectory at fixed spatial coordinates $(r,\theta,\varphi)$ by its proper time $\tau=t/r$. In this setting, letting $s := \tau - \tau'$, the transition rate reads \cite{Fewster2016ewy,Louko2007mu}
\begin{align}
	\label{eq:transition rate}
	\dot{\mathcal{F}}_{\kappa} &= \int_\mathbb{R}ds\, e^{-i\Omega s} \omega_{2,\kappa}\left(r\, \tau,r,\theta,\varphi,r\,\tau',r,\theta,\varphi\right).
\end{align}

Using the two-point functions as per Equation \eqref{eq: two point function ground and thermal time function}, as well as the addition formula between spherical harmonics, Equation \eqref{eq:transition rate} yields
\begin{align*}
	\dot{\mathcal{F}}_{\kappa} &= C_\beta(\Omega) \sum\limits_{\ell=0}^{\infty}(2\ell+1) \frac{ \psi_{\kappa}(r)^2 }{2r\mathcal{N}_\kappa}\bigg|_{p = \frac{|\Omega|}{r}},
\end{align*}
where, being $\Theta$ the Heaviside step function,
\begin{align}
	\label{eq: def Cbeta}
C_\beta(\Omega)  :=
\begin{cases}
\Theta(-\Omega), & \text{for } \beta=\infty,\\
\frac{\sign\Omega}{e^{\beta  \frac{\Omega }{r} } - 1} & \text{for } \beta<\infty.
\end{cases}
\end{align}
\noindent Observe that the variable $ p$ must be evaluated at $\frac{|\Omega|}{r}$, but, for simplicity, we shall not indicate it explicitly in the following.

Recalling that only the $\ell=0$ term calls for a boundary condition, let us decompose the transition rate as
\begin{align*}
\dot{\mathcal{F}}_{\kappa} = \dot{\mathcal{F}}_\kappa^{(\gamma)} + \dot{\mathcal{F}}^{\text{(c)}}.
\end{align*}
The component
\begin{align}
	\label{eq: def trans kappa gamma}
\dot{\mathcal{F}}_\kappa^{(\gamma)} := C_\beta(\Omega)  \frac{\psi_{\kappa}(r)^2}{2r\mathcal{N}_\kappa}
\end{align}
corresponds to the $\ell=0$ contribution, while the remainder
\begin{align}
	\label{eq: trans common term}
\dot{\mathcal{F}}^{\text{(c)}}:=  C_\beta(\Omega)  \sum_{\ell=1}^\infty \frac{2\ell+1}{4r} u ^2
\end{align}
is independent of both $\kappa$ and $\gamma$.

There are two scenarios, i) and ii) as described in the following, where we can take the sum over $\ell$ giving an analytic expression for $\dot{\mathcal{F}}_{\kappa}$.
\begin{itemize}
\item[i)]\label{page item i}
		If $m_0=0$, then $\nu=\ell+1/2$ and thus $\nu=1/2$ for $\ell=0$. Consequently $\dot{\mathcal{F}}_1 \equiv \dot{\mathcal{F}}_2$ and, using the identity
\begin{align*}
  \sum_{\ell=0}^\infty (2\ell+1)j_\ell(z)^2 = 1,
\end{align*}
we obtain
\begin{align}
	\label{eq: m0 gammaNOTnecessarily0 transition rate}
	\dot{\mathcal{F}}^{\text{(c)}} &= C_\beta(\Omega)   r\left[\frac{p}{2\pi}  - \frac{1}{2\pi}  \frac{\sin^2(pr)}{p r^2}\right].
\end{align}

\item[ii)]
		If $\gamma=0$, the secondary solution plays no r\^{o}le and it holds true that $\dot{\mathcal{F}}_1 \equiv \dot{\mathcal{F}}_2$.
\end{itemize}
Note that if either i) or ii) holds true, then $ \dot{\mathcal{F}}_1 \equiv \dot{\mathcal{F}}_2 =: \dot{\mathcal{F}}$. Moreover, if both hold true, then using the expression obtained in item i) above, we get
$
	\dot{\mathcal{F}} = C_\beta(\Omega)  \frac{|\Omega|}{2\pi},
$ which is consistent with the fact that $\text{AdS}_2 \times \mathbb{S}^2$ is conformal to Minkowski spacetime

\subsection{The $\ell=0$ contribution}

In this section we focus on the $\ell=0$ term: $\dot{\mathcal{F}}_\kappa^{(\gamma)}$.
First, note that the secondary solutions $v_1$ and $v_2$, as per Equation \eqref{eq: sec sol v1 v2}, are markedly different when $\nu$ is close to $1$, that is for $m_0\sim \sqrt{\frac{3}{4}}$. Let us thus focus in this section on this mass range. We shall show that, for $\gamma\in [0,\frac{\pi}{2}]$ and $\nu_0\sim 1$, the detector behaves in a qualitatively different way for the various choices of secondary solution. We can verify it analytically and visualize it numerically, as shown next. As a matter of fact expanding  $\dot{\mathcal{F}}_\kappa^{(\gamma)}$ as $\nu_0\sim 1$ yields
\begin{widetext}
\begin{subequations}
\label{eq: an exp F1gamma F2gamma nu0 1}
\begin{align}
	\dot{\mathcal{F}}_1^{(\gamma)} = & C_\beta(\Omega) \bigg\{ \nonumber \\
	& (\nu_0 -1)^0\cdot \frac{1}{4}\frac{\left(r \cos (\gamma ) J_1(|\Omega| ) - |\Omega|  \sin (\gamma ) Y_1(|\Omega|)\right)^2}{r^2\cos ^2(\gamma )+|\Omega|^2 \sin ^2(\gamma )} + \nonumber \\
        + &(\nu_0 -1)^1\cdot \frac{1}{2} \frac{ (r \cos (\gamma ) J_1(| \Omega | )-| \Omega |  \sin (\gamma ) Y_1(| \Omega | )) \left(r \cos (\gamma ) \left[\partial_{\nu_0} J_{\nu_0} (| \Omega | )|_{\nu_0=1}\right] - |\Omega |  \sin (\gamma )  \left[\partial_{\nu_0} Y_{\nu_0} (| \Omega | )|_{\nu_0=1}\right]\right)}{ r^2 \cos ^2(\gamma ) + |\Omega|^2  \sin ^2(\gamma ) } + \nonumber \\
				+& o(\nu_0-1)^2\bigg\},\label{eq: an exp F1gamma F2gamma nu0 1a} \\
	\dot{\mathcal{F}}_2^{(\gamma)} = & C_\beta(\Omega) \bigg\{ \nonumber\\& (\nu_0 -1)^0\cdot 	\frac{1}{4}  J_1(|\Omega| )^2 + \nonumber \\
	+& (\nu_0 -1)^1 \cdot  \frac{1}{2} \frac{  J_1(| \Omega | ) \left(r^2 \cos (\gamma \left[\partial_{\nu_0} J_{\nu_0} (| \Omega | )|_{\nu_0=1}\right]-| \Omega | ^2 \sin (\gamma )  \left[\partial_{\nu_0} J_{\nu_0} (| \Omega | )|_{\nu_0=-1}\right] \right)}{ r^2 \cos ^2(\gamma ) - |\Omega|^2  \sin ^2(\gamma ) } + \nonumber\\
	+& o(\nu_0-1)^2\bigg\}.\label{eq: an exp F1gamma F2gamma nu0 1b}
\end{align}
\end{subequations}
\end{widetext}

\noindent The consequences of Equation \eqref{eq: an exp F1gamma F2gamma nu0 1} are two-fold:
\begin{itemize}
	\item[1)] at zeroth order, the zeros of $\dot{\mathcal{F}}_2^{(\gamma)}$ are those of $J_1(|\Omega|)$ and they do not depend on $\gamma$, while those $\dot{\mathcal{F}}_1^{(\gamma)}$ do;
	\item[2)] by direct inspection of Equation \eqref{eq: an exp F1gamma F2gamma nu0 1b}, $|\Omega| = r \cot^{\frac{1}{2\nu_0}}(\gamma) $ is a pole of $\dot{\mathcal{F}}_2^{(\gamma)}$, while Equation \eqref{eq: an exp F1gamma F2gamma nu0 1a} entails that no pole occurs if $\kappa=1$.
\end{itemize}
These features are manifest in Figure \ref{fig: l0 term} which is displayed at the end of the paper for layout purposes. It contains line and contour plots of  $\dot{\mathcal{F}}_1^{(\gamma)}$ and of $\dot{\mathcal{F}}_2^{(\gamma)}$, respectively in the left and right columns. This numerical evaluation was performed using the software Mathematica and it is presented in a notebook available online \cite{github_lsc}. For this analysis, we considered a field at inverse-temperature $\beta=2\pi$ and of mass $m_0=\sqrt{\frac{3}{4}}-10^{-2}$, which corresponds to $\nu_0 \approx 0.99$. Consistently with the asymptotic expansion in Equation \eqref{eq: an exp F1gamma F2gamma nu0 1}, Figure \ref{fig: l0 term} displays the following notable features:
\begin{itemize}
	\item The dashed, vertical white lines of subfigures (e) and (f) in Figure \ref{fig: l0 term} correspond to the zeros of $J_1$. These lines align with the zeros of $\dot{\mathcal{F}}_2^{(\gamma)}$ for all $\gamma\in\left[0,\frac{\pi}{2}\right]$ as shown in subfigure (f), but this is not the case for $\kappa=1$, as shown in subfigure (e). Accordingly, we refer to this difference in behavior between $\dot{\mathcal{F}}_1^{(\gamma)}$ and $\dot{\mathcal{F}}_2^{(\gamma)}$ as a {\em zeroth-order high-mass effect}.
	\item The dashed, peaked, red line of subfigure (f) of Figure \ref{fig: l0 term} codifies the identity $|\Omega| = r \cot^{\frac{1}{2\nu_0}}(\gamma) $. The restriction of this curve to subfigure (b) is given by the three highlighted points: the circle, the triangle and the square. We can see that these points indicate a {\em first-order high-mass effect} manifest only for $\kappa=2$, and that brings about an extra zero for $\dot{\mathcal{F}}_2^{(\gamma)}$ that depends on the choice of $\gamma$.
\end{itemize}

\begin{widetext}

	\begin{figure}\begin{minipage}{\textwidth}
			\caption{\label{fig: l0 term} The $\ell=0$ contribution to the transition rate at $r=3$ for a massive field with $m_0=\sqrt{\frac{3}{4}}-10^{-2}$ at inverse-temperature $\beta=2\pi$ admitting the generalized $(\gamma,v_1)$- and $(\gamma,v_2)$-boundary conditions, respectively, on the left and on the right columns. On top, seen as a function of $\Omega$, each curve corresponds to a fixed value of $\gamma$. In the middle, density plots for varying both $\gamma$ and $\Omega$. On the bottom, the same density plots of the middle, but with dashed vertical (white) lines corresponding to the first four zeros $j_i$, $i\in\{1,2,3,4\}$, of $J_1(\Omega)$. At the bottom right, the wiggling, dashed, red line is given by $|\Omega| = r \cot^{\frac{1}{2\nu_0}}(\gamma) $. The three points marked by a circle, a triangle and a square on image (b) map to those of image (f). All images share the $\Omega$-axis at the bottom.}
			%
			\scriptsize

			\begin{tabular}{c c c}
				%
				\hspace{-10pt}\includegraphics[align=c,width=.41\textwidth]{diffbcs1}    			 &
				\hspace{-10pt}\includegraphics[align=c,width=.41\textwidth]{diffbcs2}    			 &
				\includegraphics[align=c,width=.08\textwidth]{plegends2}  			 \\
				%

				%
				\hspace{.65cm}(a) & \hspace{.65cm}(b) &  \\
				%

				\includegraphics[align=c,width=.4\textwidth]{p1}        			   &
				\includegraphics[align=c,width=.4\textwidth]{p2}        			   &
				\includegraphics[align=c,width=.08\textwidth]{plegends1}  			 \\
				%

				%
				\hspace{.65cm}(c) & \hspace{.65cm}(d) &  \\
				%

				\includegraphics[align=c,width=.4\textwidth]{p1withlinesnolabel} &
				\includegraphics[align=c,width=.4\textwidth]{p2withlinesnolabel} &
				\begin{minipage}{.08\textwidth}
					\setlength{\extrarowheight}{7pt}
					\scriptsize	\vspace{-.4cm}
					\begin{tabular}{c}
						Zeros of $J_1$\\
						$j_1\sim 3.83$ \\
						$j_2\sim 7.01$ \\
						$j_3\sim 10.2$ \\
						$j_4\sim 13.3$ \\ \\ \\

					\end{tabular}
					\includegraphics[width=1.6\textwidth]{plegendsredline}
				\end{minipage}\\
				%

				%
				\hspace{.65cm}(e) & \hspace{.65cm}(f) &  \\
				%
				& & \\
			\end{tabular}
			%
	\end{minipage}\end{figure}
\end{widetext}
\normalsize

\subsection{The sum over $\ell$}

The sum over $\ell$ conceals the high-mass effect discussed in the previous section, see Figure \ref{fig:summed up to 100} that shows the transition rate summed up to $\ell_{\scriptscriptstyle{\text{max}}}=100$ with the same parameters of Figure \ref{fig: l0 term}. Nonetheless we can still disambiguate between the choices of secondary solution for the radial equation if, instead of checking out its zeros, we focus on the behaviour of the transition rate at  $\Omega=0$. Figure \ref{fig:summed up to 100} corroborates this statement, but it does not clarify what is happening at such locus. Yet, we can analytically compute the limit of $\dot{\mathcal{F}}_\kappa$ for vanishing energy gap, as we show in the following.

\begin{widetext}
\begin{figure}[H]
  \centering
\includegraphics[align=t,width=.12\textwidth]{p100_0}
	 \includegraphics[align=t,width=.35\textwidth]{p100_1}    \includegraphics[align=t,width=.35\textwidth]{p100_2} \includegraphics[align=t,width=.12\textwidth]{p100_6}
  \caption{The transition rate summed up to $\ell_{\scriptscriptstyle{\text{max}}}=100$, with the same parameters and legends as those of Figure \ref{fig: l0 term}. On the left, for $\kappa=1$. On the right, for $\kappa=2$. The three small plots on the right side provide a zoom-in to a neighborhood of the points highlighted by a circle, a triangle and a square. The three small plots on the left correspond to a zoom-in to the same region of those on the right. The contrast observed between the two sets of three small plots is due to the {\em high-mass effects}.}
  \label{fig:summed up to 100}
    \end{figure}
\end{widetext}

For $\ell\geq1$, we have that $\nu> \frac{1}{2}$ regardless of the value of the field mass. Thus, using the following two identities

\begin{align*}
	&\lim\limits_{\Omega\rightarrow 0} J_\nu(\Omega) = 0, \text{ for }\nu>0,\\

  & \lim\limits_{\Omega\rightarrow 0} \frac{J_\nu^2(\Omega)}{e^{\Omega}-1} = 0, \text{ for }\nu>\frac{1}{2},
\end{align*}
one can show that
	\begin{equation*}
	 \lim\limits_{\Omega\rightarrow 0} \dot{\mathcal{F}}^{\text{(c)}} = 0, \text{ for }\beta\leq \infty.
	\end{equation*}
Consequently, the behavior of $ \dot{\mathcal{F}}_\kappa$ at $\Omega=0$ depends only on the contribution from the $\ell=0$ term: $\dot{\mathcal{F}}_\kappa^\gamma$, defined in Equation \eqref{eq: def trans kappa gamma}.

For vanishing mass, $m^2_0=0$, $\dot{\mathcal{F}}:=\dot{\mathcal{F}}_1\equiv \dot{\mathcal{F}}_2$. Since
\begin{align*}
  & \lim\limits_{\Omega\rightarrow 0} \frac{J_{\nu}^2(\Omega)}{e^{c\Omega}-1} = \frac{2}{\pi c}, \text{ for }c>0, \, \nu=\frac{1}{2},
\end{align*}
we obtain that, for $\gamma\in\left[0,\frac{\pi}{2}\right)$,
				\begin{equation}
					\label{eq: lim transition m0 gamma0}
			 		\lim\limits_{\Omega\rightarrow 0} \dot{\mathcal{F}} =
					\begin{cases}
						0, &\text{ for }\beta=\infty;\\
						\frac{(r+\tan(\gamma))^2}{2\pi r\beta}, &\text{ for }\beta<\infty.
					\end{cases}
				\end{equation}

Observe that if $\gamma=\pi/2$ then $\lim\limits_{\Omega\to 0}\dot{\mathcal{F}}$ diverges for $\beta\leq\infty$.
For positive masses, $\dot{\mathcal{F}}_1 \not\equiv \dot{\mathcal{F}}_2$ and the results are summarized in Table \ref{tab: limits as Omega 0 for m greater than 0}.

\renewcommand{\arraystretch}{1.5}
\begin{table}[H]
   \centering
	\caption{The limits for $m_0\in\left(0,\sqrt{\frac{3}{4}}\right)$.}
   \begin{tabular}{c c c c}
   \toprule
	 State &
   $\gamma$&
   $\lim\limits_{\Omega\rightarrow 0}\dot{\mathcal{F}}_1 $&
   $\lim\limits_{\Omega\rightarrow 0}\dot{\mathcal{F}}_2$ \\

	 \hline\hline
      \multirow{ 2}{*}{$ \beta=\infty$ } &  $\left[0,\frac{\pi}{2}\right)$ & $0$      &  $0$       \\
			                          &  $\frac{\pi}{2}$                & $\infty$ &  $\infty$  \\\hline
	  \multirow{ 3}{*}{$ \beta<\infty$ } &  $0$                            & $0$      &  $0$       \\
		\rowcolor{cyan} $ \beta<\infty$     &  $\left(0,\frac{\pi}{2}\right)$ & $\infty$ &  $0$       \\
		                           &  $\frac{\pi}{2}$                & $\infty$ &  $\infty$  \\
 \bottomrule
   \end{tabular}
\label{tab: limits as Omega 0 for m greater than 0}
\end{table}

All in all, the boundary condition chosen for the mode $\ell=0 $ affects the behavior of the transition rate in its full form, both quantitatively and qualitatively. In fact, its limit at $\Omega\rightarrow 0$ receives no contribution from the terms with $\ell>0$, as shown above. Markedly, the transition rate of an Unruh-DeWitt detector interacting with a Klein-Gordon field at a finite temperature state, with $m_0\in\left(0,\sqrt{\frac{3}{4}}\right)$, and admitting a generalized $(\gamma,v)$-boundary condition with $\gamma\in\left(0,\frac{\pi}{2}\right)$ explicitly discriminates between $\kappa=1$ and $\kappa=2$ in Equation \eqref{eq: sec sol v1 v2}, as highlighted in blue in Table \ref{tab: limits as Omega 0 for m greater than 0}.

\section{Outlook}
\label{sec: Summary}

We have shown that choosing one among the generalized $(\gamma,v)$-boundary conditions is not a mere mathematical exercise, rather it has notable physical consequences. For definiteness we have considered a Klein-Gordon field on $\text{AdS}_2 \times \mathbb{S}^2$ spacetime and we have constructed the two-point correlation function of the ground and thermal states for two families of boundary conditions, characterized by different choices of secondary solution, $v_1$ and $v_2$, see Equation \eqref{eq: sec sol v1 v2}. We probed the underlying system with an Unruh-DeWitt detector and we showed that its transition rate differentiates between the two families of boundary conditions.

There are two notable facts that highlight the hidden freedom in the mode expansion of the Klein-Gordon field and connect our results to \cite{Campos2022byi}. First, the boundary conditions with $(\gamma,v_1)$ and $(\gamma,v_2)$ are equivalent at the level of the radial differential equation, but cease to be such at the level of the whole Klein-Gordon equation. This is hinted by the map that relates the parameter $\gamma$ in the two settings. It depends on the Fourier frequency associated to the time coordinate and, therefore, we can map a solution satisfying Equation \eqref{eq: robin bc} for $\kappa =1$, into one satisfying it for $\kappa=2$ via
	\begin{equation*}
		\tan(\gamma) \mapsto  \frac{p^{2\nu}\sin(\pi\nu)}{p}\frac{\tan(\gamma)}{1+\tan(\gamma)p^{2\nu}\cos(\pi\nu)}.
	\end{equation*}
There is an explicit dependence of the Fourier parameter $p$ which in turn related to the underlying frequency. This entails that, while, at the level of radial equation, this is a mere algebraic relation, from a fully covariant viewpoint, this is no longer true as one can realize by an inverse Fourier transform. Finally, our work sets the ground for multiple future investigations that are worth mentioning:
\begin{itemize}
	\item comparisons between the physical phenomena associated to different choices of secondary solution, possibly considering other observables such as the stress-energy tensor;
	\item proving the universality of the features highlighted in this paper by studying their occurrence in other globally hyperbolic spacetimes with a timelike boundary, such as for example global monopoles and cosmic strings;
	\item studying the influence of the $(\gamma,v)$-boundary conditions on specific phenomena such as the anti-Hawking effect, see \cite{DeSouzaCampos2020ddx, Brenna:2015fga, Robbins:2021ion, deSouzaCampos:2020bnj};
	\item a further investigation of the behaviour of the transition rate at $\Omega\rightarrow 0$ particularly in connection with the notion of a structureless scalar source \cite{Cozzella2020gci}.
\end{itemize}

\section{Acknowledgments}
The work of L.C. is supported by a postdoctoral fellowship of the Department of Physics of the University of Pavia, while that of L.S. by a PhD fellowship of the University of Pavia.


\section{Introduction}\label{sec.I}

\begin{figure}[!htb]
\includegraphics[width=.8\hsize]{pic-1.png}
\caption{Comparison of the two decay modes: 1)  a regular three-body decay mode of the neutron into protons, electrons and neutrinos, with the continuum energy spectra for all three; 2) an exotic decay mode of the neutron into a hydrogen atom and an anti-neutrino, which has a extremely small cross section and the final state particles are both monoenergetic.}
\label{fig1}
\end{figure}

Neutron decay is closely related to fundamental physics. The fact that, no new particles other than the Higgs boson have been found at the LHC, has prompted a turn to the search for signatures of new physics in the low and medium-energy regions~\cite{1-BoB-th4}. Experiments related to $\beta$ decay are an important avenue to explore the standard model of electroweak interactions, e.g.,  a neutron decaying into a proton, an electron, and an electron antineutrino, which is a concise model for studying weak interactions. In 1947 Daudel, Jean and Lecoin~\cite{1-BoB-th5,1-BoB-th6} have predicted the existence of a two-body $\beta$-decay mode where the daughter nucleus and the electron remain bound. In the case of the free neutron, the neutron two-body decay mode is $\rm{n \rightarrow H+\bar {\nu}_e}$, which is also referred to as "bound $\beta$-decay" (BoB). As this decay results in a two-body final state, it is theoretically quite straightforward to determine the states of individual particles. Many theoretical studies have been carried out on this two-body decay~\cite{1-BoB-th6,1-BoB-th1}, although it has never been observed in experiments so far. The successful measurement of this process is of great physical importance~\cite{1-BoB-th1}. As the BoB leads to a two-body final state, the spin state of the anti-neutrino is thus mirrored by the outgoing hydrogen atom. As such, measuring the hyperfine spin state of the hydrogen atom would therefore contain full information relating to the momentum direction of the antineutrino~\cite{McAndrew:2014iia}. Moreover, the measurement of the BoB branching ratio can also be a solution to the recently proposed neutron lifetime puzzle~\cite{1-BoB-th8,Czarnecki:2018okw}.

 The major challenges in observing the BoB decay and investigating its properties lie in the small predicted branching ratio $\sim4\times10^{-6}$  of the dominated three-body decay mode~\cite{1-BoB-th9,1-BoB-th8,1-BoB-th10}, together with the detection of low-energy electrically neutral hydrogen atoms in the final state. According to the theory, the hydrogen atom in the final state has a kinetic energy of 325.7 eV, corresponding to a velocity $\sim 10^5 {\rm m/s}$, and is expected to be populated with zero angular momentum, specifically, with $83.2\%$ of atoms in the 1s-state and $10.4\%$ in the 2s-state and the remainder in an $n$s-state where $n\geq3$~\cite{1-BoB-th2}.

In this study, we propose a scheme to detect the BoB that will rely on the  measurement of H atoms through its energy spectrum structure, in order to observe for the first time this hitherto unobserved decay mode. As known, the standard three-body decay mode is $n \rightarrow H^{+}+e^{-}+\bar {\nu}_{e}$, where the proton and electron are charged particles that are easily captured by solid materials, while the antineutrino is extremely difficult to capture. The total energy in the three-body decay is conserved, leading to continuum energy spectra for each of particles. Since the mass of the proton is much larger than that of the electron and neutrino, the decay energy is obtained mainly by the electron and neutrino, both of which have kinetic energies in the range of 0--782 keV, while the energy spectrum of the proton lies in the range of 0--750 eV~\cite{1-BoB-th3}.

For the BoB mode, the conservation of momentum and energy during the decay leads to the energy carried by the antineutrino is constant, while the kinetic energy of the hydrogen atom is also constant (325.7 eV),
making it easy to be captured by solid materials. Since the hydrogen atom in the $n$s ($\geq 2$) state after the two-body decay can be captured by the detector, it will de-excite to the 1s state very quickly, transferring its thermal energy greater than 10.2 eV to the detector. Thus, the kinetic energy spectra of the hydrogen atoms on the detector would be a multiple line structure dominated by 325.7 eV and 335.9 eV. The lower right panel of Fig.~\ref{fig2} shows only these two spectral lines, and the other spectral lines from $n$s states are ignored for the moment. Moreover, three-body-decay protons will introduce a strong background continuum spectrum contaminating the measurement of the hydrogen atoms, which needs to be removed to obtain a high signal-to-noise ratio. Since the charged protons and electrons are susceptible to an external electric field, one can set up positive and negative electrodes at about 1 kV in a low-temperature vacuum, resulting in a large difference in signal amplitudes between electrons with  $E>1$ keV and hydrogen atoms at 325.7 eV, which therefore can be easily identified. In Fig.~\ref{fig3}, the expected shapes of energy spectrum before and after electric field screening is shown. As seen, combine with the 1kV screening electric field, the central energy of the electrons is in the energy region of several hundred keV,  which is very different from the line signature of the two-body-decay hydrogen atoms.

A microcalorimeter is a sensitive detector for heat signals, with extremely high energy resolution at eV level and no dead layer, as well as a wide range of available absorption materials~\cite{2-uCal-theory,3-basic-theory}. The energy of the particle deposition is converted into heat and the resulting temperature rise is measured. Since hydrogen atoms, electrons and protons can all generate heat signals at the detector, and the energy resolution at the eV level fully satisfies the requirement for resolution of a single energy peak at 325.7 eV, the microcalorimeter detector is thus an ideal detector for measuring the BoB according to energy spectrometry. Although electrons and protons below 1 keV can be completely eliminated by the electric field, there are still some background processes that would produce a continuum at the detector. As seen from the lower right panel of Fig.~\ref{fig2}, if the initial kinetic energy of neutrons is much less than 1 eV and the resolution of the detector is around 1 eV, not only can the hydrogen atom be clearly seen in the 1s or 2s or even a higher state, but also a high signal-to-background can be obtained. For the transition edge sensor (TES)-based microcalorimeter, in the 1 keV energy range, the best reported energy resolution is about 0.75 eV~\cite{4-TES0.75eV}, and the best one we obtained so far is about 1.4 eV~\cite{5-TES1.4eV}, so that the detector performance is basically appropriate for the BoB measurement. More importantly, as the energy range of interest is around 325.7 eV, a better energy resolution can be obtained by further reducing the heat capacity of the microcalorimeter~\cite{3-basic-theory}.

This paper is organized as follows. In Sect.~\ref{sec.II}, we will present the theoretical calculations for the BoB and summarize various neutron sources, then discuss the the requirements for the BoB measurements, in addition we briefly review the working principle of microcalorimeter. In Sect.~\ref{sec.III}, we outline our experimental concept and investigate suitable neutron sources, as well as discuss measurement related issues, such as the $^3$He problem in refrigerators, the elimination of proton effects, effects of high-energy particles, the method of readout electronics, etc. Finally, we draw our summary and conclusions in Sect.~\ref{sec.IV}.

\begin{figure*}[!htb]
\includegraphics[width=0.8\hsize]{pic-2.png}
\caption{Comparison of the probability densities  of electrons (upper left), neutrinos (upper right) and protons (lower left) produced by neutron three-body decay and hydrogen atoms (lower right) by two-body decay~\cite{1-BoB-th3}. The broadening of the hydrogen atom energy spectrum is determined by the initial kinetic energy of the parent neutron. With varying the initial kinetic energy, the hydrogen atom energy spectra broadened to 1 eV, 2 eV, 5 eV, and 10 eV are shown. Note that, the predicted probablity density of hydrogen atoms is highly suppressed by the two-body-decay branching ratio that is assumed to be $4\times 10^{-6}$.}
\label{fig2}
\end{figure*}

\begin{figure*}[!htb]
\includegraphics[width=0.8\hsize]{pic-3.png}
\caption{Same as Fig.~\ref{fig2}, but for comparison of the BOB-induced spectral structures before (left) and after (right) electric field screening. Adding a voltage of 1 kV or more can significantly eliminate the effects of electrons and protons. The protons can be completely eliminated and are thus invisible on the plot, since their kinetic energies are less than 1 keV.}
\label{fig3}
\end{figure*}

\section{basic principles}\label{sec.II}

\subsection{Theoretical Calculations of Neutron Two-Body Decay}\label{A}

Let us consider the BoB process, ${\rm n} \rightarrow {\rm H}+\bar{\nu}_{\rm e}$, where the expression for the kinetic energy of the hydrogen atom can be derived from the conservation of energy and momentum. In the following, the subscripts of n, H, $\bar{\nu}_e$ represent neutron, hydrogen atom and electron antineutrino, respectively. $E$, $P$ and $M$ denote the energy, the momentum and the rest mass, respectively. We will use natural units in which $c = 1$ throughout. First, the conservation of energy leads to
\begin{equation}\label{eq:E}
E_{\rm n} = E_{\rm{H}}+E_{\bar{\nu}_e}\,,
\end{equation}
According to the energy–momentum relation and the momentum conservation, $P_{\rm{H}}+P_{\bar{\nu}_e}=0$, Eq.~\ref{eq:E} is then expressed by
\begin{equation}\label{eq:pm}
M_{n} = \sqrt{M_{\rm{H}}^2+P_{\rm{H}}^2}+\sqrt{M_{\bar{\nu}_e}^2+P_{\bar{\nu}_e}^2}\,,
\end{equation}
which yields
\begin{equation}\label{eq:pm1}
\left[M_n -\sqrt{M_{\rm{H}}^2+P_{\rm{H}}^2}\right]^2 = M_{\bar{\nu}_e}^2+P_{\rm H}^2\,,
\end{equation}
and by expanding it, one obtain
\begin{equation}\label{eq:pm2}
\frac{M_{\rm n}^2+M_{\rm H}^2-M_{\bar{\nu}_e}^2}{2M_n}= \sqrt{M_{\rm H}^2+P_{\rm H}^2}
\end{equation}
In the non-relativistic limit for H,
\begin{equation}\label{eq:lim}
\sqrt{M_{\rm H}^2+P_{\rm H}^2}\approx M_{\rm H} + T_{\rm H}\,,\quad {\rm with}~ T_{\rm H} = \frac{P_{\rm H}^2}{2M_{\rm H}}
\end{equation}
Inserting it into Eq.~\ref{eq:pm2}, the expression for the kinetic energy of H, $T_{\rm H}$, then reads
\begin{equation}\label{eq:TH}
T_{\rm H} = \frac{(M_{\rm n}-M_{\rm H})^2-M_{\bar{\nu}_e}^2}{2M_{\rm n}}\,.
\end{equation}
As known, $M_{\rm n}-M_{\rm H}\approx 782.15$ keV, $M_{\rm n}=939.56542$ MeV, and the neutrino energy $M_{\bar{\nu}_e}\ll 1$ eV. Based on these values, the recoil energy of H is about $T_{\rm H}=325.56$ keV. Inversely, from Eq.~\ref{eq:TH}, the mass of a neutrino can be obtained by precisely measuring the kinetic energy of the recoiling nucleus,
\begin{equation}\label{eq:Mnu}
M_{\bar{\nu}_e}\simeq \left[(M_{\rm n}-M_{\rm H})^2 - 2M_{\rm n}T_{\rm H}\right]^{1/2}\,.
\end{equation}
Taking derivative of both sides of Eq.~\ref{eq:TH}, we can find
\begin{equation}\label{eq:dT}
dT_{\rm H} = -\frac{M_{\bar{\nu}_e}}{M_{\rm n}}dM_{\bar{\nu}_e}\,,
\end{equation}
It follows from the prefactor of $M_{\bar{\nu}_e}/M_{\rm n}$ that accurate measurement of $M_{\bar{\nu}_e}$ requires a very high measurement accuracy for $T_{\rm H}$, and, e.g., constraining the neutrino mass down to 0.1 eV requires an accuracy of $T_{\rm H}$ to the level of  $10^{-11}$ eV.

\subsection{Energy Loss of High-Energy Particles in Microcalorimeter Absorbers}\label{B}

Let us first investigate the various possible background signals, which will involve interactions between high-energy particles and experimental materials. The types of particles include mainly high-energy muons, electrons , $\alpha$ particles and $\gamma$ rays, as well as low-energy hydrogen atoms and protons. Neutrinos are also the main secondary particles during the neutron decays, whereas the interactions with the material can be completely neglected due to their extremely high penetration capacity. High-energy muons come from cosmic rays with a typical energy of 4 GeV, and high-energy electrons and lower-energy protons primarily originate from the three-body decay of neutrons, with electron energies $\lesssim782$ keV and proton energy $\lesssim750$ eV. Due to the use of a 1 kV electrode device, proton and electron signals with energies $\lesssim1$ keV are fully screened out and can therefore be safely ignored. Moreover, $\gamma$ rays from neutron sources have a typical energy around MeV, and $\alpha$ particles are produced by the $\alpha$ decay of radioactive nuclei, also with a typical energy of MeV. As mentioned before, the two-body-decay induced hydrogen atoms typically have the kinetic energy of 325.7 eV. For the MeV $\gamma$ rays, the penetration power is known to be extremely high, resulting in an extremely low cross section for direct interactions with the microcalator, so we can also ignore their effects completely. For $\alpha$ particles, due to their relatively weak penetration power, if it comes from the natural radioactive nuclei in the microcalorimeter itself, all their energies will be deposited on the microcalorimeter and thus a detailed analysis of their energy loss in the microcalorimeter absorber is not necessary. Similar to the $\alpha$ particles, the low-energy hydrogen atoms with kinetic energy of 325.7 eV can only penetrate metallic materials of a few nanometers, and all their energies will be deposited in the microcalorimeter, according to the simulation of the~\texttt{SRIM} software~\footnote{http://www.srim.org/}.

Therefore, in the following we will focus on the analysis of the interactions of high-energy muons and electrons with microcalorimeter absorbers. Both muons and electrons are charged particles, and when they interact with materials, the mean rates of energy loss can be given by the Bethe-Bloch formula~\cite{8-Detector-theory}, which reads
\begin{equation}\label{eq:BB}
-\left<\frac{dE}{dx}\right> = \frac{4\pi z^2e^4NZ}{m_0v^2}\left[\ln\left(\frac{2m_0v^2}{I(1-\beta^2)}\right)-\beta^2-\frac{C}{Z}\right]\,,
\end{equation}
which is for a particle with speed $v$, charge $z$ (in multiples of the electron charge $e$) and energy $E$, traveling a distance $x$ into a target with the mean excitation potential $I$ (usually measured through experiments). Here $\beta=v/c$ is the velocity ratio, $N$ is Avogadro constant, $Z$ is the atomic number of the atoms of the absorber material, and $C/Z$ is the shell correction.  $m_0$ is the mass of the incident particle.

According to Eq.~\ref{eq:BB}, the energy loss rate is about 2 MeV per $\rm g/{cm}^2$ for muons with $\beta\rightarrow 1$~\cite{8-Detector-theory}. For 200 nm thick copper and 80 nm thick molybdenum, the energy loss for vertical penetration is about 535 eV. Due to the broad distribution of incident directions of the muons, the energy left in the TES per penetration will be $\ge 535$ eV. Since the lowest signal produced by the muons on the microcalorimeter is about 200 eV higher than the target signal of 325.7 eV, it is easily excluded from signal amplitude and the high-energy muons do not  contribute a continuous background on the target signal.

\begin{figure*}[!htb]
\includegraphics[width=0.8\hsize]{pic-4.png}
\caption{{\it Left}: energy loss rate(Stopping Power) of electrons in copper due to collision and radiation at different kinetic energies. {\it Right}: total energy loss rates of electrons in different metallic materials at different kinetic energies.}
\label{fig4}
\end{figure*}

For high-energy relativistic electrons, the energy loss in the interaction with the material varies for its kinetic energy. Since the mass of electrons is much smaller than that of heavy ions, electromagnetic waves are radiated during the acceleration of electrons. Therefore, when the high-energy electrons collide with a material, there is a certain amount of radiated energy loss along with the collision energy loss.

The left panel of Fig.~\ref{fig4} shows the collision and radiated energy loss rates at different electron kinetic energies. As seen, when $E<10$ MeV, the energy loss is dominated by the collisions. The electron kinetic energy from neutron three-body decay is $E <0.78$ MeV, so that we mainly account for the collision energy loss. It can be seen that in the energy range of interest, at the same material thickness, the energy loss of high-energy electrons is rather low, because if the electron energy is high enough to penetrate the material, the high-energy electrons leave only a small fraction of their energies in the detector. In the right panel, the energy loss rates for different metallic materials are shown. As will be mentioned later, molybdenum and copper are the two main elements in our microcalorimeter, and copper and aluminum are the most used elements in the main structure of the refrigerator.

\subsection{Classification and characteristics of neutron sources}\label{C}

Nuclear reactors, radioisotopic neutron emitters and accelerators are the main sources of neutrons, with typical intensities of $10^2$--$10^8$ and $10^8$--$10^{12}$ n$\rm s^{-1}$ for the first two, respectively, and the flux of $10^{10}$--$10^{15}$ n$\rm s^{-1} cm^{-2}$ for the last one.

The neutron emissions are mainly from $\alpha$ neutron, $\gamma$ neutron, and spontaneous fission neutron sources. The nuclear potential barrier increases with increasing the nuclear charge, so the $alpha$ neutron sources generally use light nuclei as targets. Po-Be, Ra-Be and Am-Be sources are the most frequently used neutron $\alpha$ sources. The $\gamma$ neutron source produces neutrons by irradiating $\gamma$ with an energy higher than the binding energy of a target, e.g., Be. The neutron kinetic energy is equal to the photon energy minus the neutron binding energy. Some heavy nuclei can make intense neutron sources, such as californium-252 that emits $2.3\times 10^{12}$ neutrons per gram per second through spontaneous fission. When a few very heavy nuclei absorb a neutron, they become so destabilized that they usually split and eject a few neutrons, which is called nuclear fission. Nuclear reactors thus produce neutron sources through fission nuclei and producing a self-sustaining nuclear chain reaction. Moreover, accelerator-based neutron sources produce neutrons by bombarding specific nuclei with high-energy charged particles produced by accelerators, which has the advantage of high intensity and can obtain monochromatic pulsed neutron beams in a wide energy interval. Neutron sources from accelerators are very intensive, since bombarding the target nucleus with GeV particles will emit multiple neutrons simultaneously. Since accelerator-based source is also pulsed and the target position is well determined, the decay position of the neutrons can be roughly estimated from the neutron velocity. In turn, measurements can be made at specific times and specific locations, allowing for background subtraction for neutron decay measurements.

\begin{figure}[!htb]
\includegraphics[width=0.8\hsize]{pic-5.png}
\caption{Schematic diagram of the structure of microcalorimeter, which directly measures heat signals, insensitive to the type of incident particles.}
\label{fig5}
\end{figure}

\subsection{Working principle and characteristics of microcalorimeter}\label{D}

Microcalorimeter is a new type of detector proposed by D. McCammon et al. in the 1980s~\cite{2-uCal-theory,4-Si-uCal,4-uCal-Rev}, which measures heat changes, with very low energy thresholds and excellent energy resolution.
In principle, any type of particles that is deposited with the same energy and without quenching effects will produce the same signal. In principle, the signal produced by any type of particles remains the same as long as they are deposited with the same energy and without quenching effects. Then, the energy resolution of microcalorimeter, $\delta E_{\rm FWHM}$, is independent of the incident particle energy ~\cite{3-basic-theory} in the linear regime,
\begin{equation}\label{eq:deltaE}
\delta E_{\mathrm{FWHM}} \propto \sqrt{4 k_{\mathrm{B}} T_0^2 C / \alpha_I}\,,
\end{equation}
where $T_0$ is the temperature of absorber, $C$ is the heat capacity of absorber, $\alpha_I$ is the temperature sensitivity of  thermometer.

Moreover, the choice of absorber material for the detector is very flexible, allowing the material to be changed according to the type of particles being detected. Also, since the absorber surface of this detector has the same structure as the inner one, there is no problem of dead layer that occurs in semiconductor detectors, so that it can be used especially for measuring the energies of particles with very weak penetration power.

A microcalorimeter is conceptually a simple device made of four components: an absorber, a thermometer, a weak thermal link and a heat sink, which is shown in Fig.~\ref{fig5}. The absorber is thick enough to stop the particles of interest, but must have the minimum possible specific heat capacity, at cryogenic temperature, to allow for a large temperature rise proportional to the energy deposition of the particles. A weak thermal link to a heat sink is employed to cool a pixel back to its base temperature. The electrical resistance of the thermometer changes inversely as temperature rise, resulting in a voltage signal. As such, the energy of incident particles can be derived from the received electrical signals, along with an appropriate calibration. To balance parameters such as cooling power of the thermostat and energy resolution, a microcalorimeter typically operates at 50 mK or lower. Three classes of microcalorimeters by thermistor technology are frequently used, including semiconductor, transition edge sensor (TES)~\cite{5-TES-uCal}, metallic magnetic calorimeter (MMC)~\cite{6-MMC-uCal}. A TES thermistor has a very high responsivity, with an excellent energy resolution~\cite{4-TES0.75eV}, and a MMC is a dissipation-less system
 that greatly improve some thermal management issues on large detector arrays, enabling excellent repeatability of measurements.

\section{Design Scheme}\label{sec.III}
\subsection{Neutron source requirements}\label{A}
\textbf{Initial kinetic energy of neutrons:}
According to two factors, the width of the energy spectrum and the decay distance of the neutrons from the initial kinetic energy, lower energy neutrons are required for the purpose of this study. Since the neutrino mass in the BoB is completely negligible, almost all of the initial kinetic energy of the neutron is transferred to the hydrogen atom. Since emitted neutrons have an energy broadening, and their energy will be further broadened uniformly upon their interactions with moderator. Therefore, the initial neutrons are not monochromatic, leading to hydrogen atoms with broadening energies. The only way to control the energy broadening of hydrogen atoms is to reduce the initial kinetic energy of neutrons. Since the energy resolution of the detector is intended to be better than 1 eV, it is thus required that the initial kinetic energy of the neutron should be lower than 1 eV. In addition, since the probability of neutron decay near the microcalorimeter is proportional to its side length and inversely proportional to the neutron velocity, the coverage of the microcalorimeter's pixels for neutron trajectories at lower kinetic energies will be better.

For a typical neutron of 4Å wavelength from the Chinese Spallation Neutron Source with a kinetic energy of about 5 meV ($v\simeq 976 {\rm m/s}$), the travel time through a microcalorimeter with the side length of 500 $\mu$m is about 0.512 $\mu$s, and the probability of a single neutron decaying in the path is about $5.8\times10^{-10}$. If an ultracold neutron with an initial kinetic energy of $10^{-7}$ eV were used instead, the probability would increase by about two orders of magnitude.

\textbf{Neutron flux requirements:} with a sufficiently high neutron flux, the single-pixel count rate can be kept around 10 cps (count per second), according to the optimal resolution requirement of the microcalorimeter. According to the predicted two-body-decay branching ratio, the count rate is about $4\times10^{-5}$ cps. In order to accumulate hundreds of counts, one pixel needs to accumulate $2.5\times10^6$ s (about 29 days). This time requirement is relatively high, so when the neutron flux is high enough, it is necessary to consider increasing the single-pixel count rate to obtain more counts.

At present, the number of pixels of microcalorimeters can reach a scale of 100 or more, which means that it takes about 3 days to accumulate a few thousand two-body decay events, so the measurement can be completed in an affordable time, while the flux requirement of ultracold neutrons is about the order of $10^9$--$10^{10}/\rm cm^2$. In addition, neutron sources are often accompanied by $\gamma$ rays, which will be briefly discussed later.

\begin{figure*}[!htb]
\includegraphics[width=0.8\hsize]{pic-6.png}
\caption{Schematic diagram of the $^3$He/$^4$He dilution refrigerator. The refrigerator provides continuous cooling to sufficient low temperatures for the BoB measurement, while a shielding layer of neutrons to protect the $^3$He is necessarily required.}
\label{fig6}
\end{figure*}

\subsection{$^3$He related issues in refrigerator}\label{B}

The microcalorimeter needs to reach a temperature of 50 mK or even lower to obtain a sufficiently high energy resolution. At present, adiabatic demagnetization refrigeration (ADR) or dilution refrigeration (DR) is commonly used to obtain such low temperatures.

The ADR uses the magnetocaloric effect caused by changing the external magnetic fields to change the state of the magnetic material, with the minimum operating temperature of 40 mK and without the influence of gravity. The ADR has high refrigeration efficiency, allowing for miniaturization and weight reduction. However, the ADR cannot continue cooling and prevent temperature rise during the magnetization process as it have to operate with a repeated magnetization and demagnetization process. Since the ADR must be terminated when the temperature rises to about 100 mK in the magnetization process, it is not available for experiments that require long periods of time, such as the BoB measurements.

The DRs utilize the properties of $^3$He and $^4$He mixtures to achieve cooling, which have the advantages of lowest temperature down to several mK, continuous cooling, and high cooling power. The operation of a DR depends on gravity to keep $^3$He and $^4$He in their correct chambers and the size is typically quit big but manageable. In particular, it can operate continuously for a long time in the range of 10--300 mK, which is suitable for our needs in this study. However, $^3$He is a strong absorber of neutrons, reacting by absorbing neutrons and producing a proton and a triton. Such reaction will lead to the failure of DR. Therefore, in order to shield the neutrons to avoid this problem, the $^3$He loop needs to be covered with a protective layer of resin for a sustainable operation of DR.

\begin{figure*}[!htb]
\includegraphics[width=0.8\hsize]{pic-7.png}
\caption{Electrode layout scheme for voltages of $\pm 1$ and 0 kV.}
\label{fig7}
\end{figure*}

\subsection{Eliminating impacts of protons}\label{C}
The three-body-decay produced protons will contribute a strong continuum background of about $E<750$ eV against the signal of hydrogen atoms from the BoB process. Therefore, such background have to be subtracted completely.
It is suggested to add the high-voltage of $\pm1$ kV at the upper and lower ends parallel to the neutron beam, respectively, to remove all charged particles with $E<1$ keV. Also, the front of the microcalorimeter needs to be equipped with a shield to screen the electric field by the high voltages.

\subsection{Impacts of high-energy particles}\label{D}

In measurements, the radiation background of the environment has some impacts on the determination of the hydrogen signal, which will be briefly discussed as follows. In the energy spectrum, especially near the target energy of 325.7 eV, there are mainly two dominated components: the continuum energy spectrum and the elemental spectral lines.

The sources of the continuum spectrum around 325.7 eV originate mainly from the following scenarios, including 1) energetic charged particles penetrating the microcalorimeter, 2) scintillation lights shining on the microcalorimeter, where the lights are produced from the interactions of high-energy particles with the material surrounding the microcalorimeter, and 3) secondary particles with energy levels of a few hundred eV, from energetic particles interacting  with the material around the microcalorimeter.

Furthermore, the characteristic spectral lines of elements near 325.7 eV are in principle negligible. This is because the elements that make up the microcalorimeter itself and its surrounding mainly include copper, molybdenum, carbon, aluminum, gold, and silicon, which have no characteristic lines in the range of 10--325.7 eV. Therefore, only the effects of the continuum energy spectrum needs to be considered. Cosmic rays and naturally occurring radioactive radiations are the main sources of the radiation background, and in addition, if spallation neutron sources are used to provide neutrons, the neutron beam is often accompanied by MeV $\gamma$ rays.

\textbf{Impacts of cosmic rays:}
at ground level, muons are the most abundant particles produced from the interactions of primary cosmic rays at the top of the atmosphere, with energies of the order of GeV and the typical flux of about 1 $\rm s^{-1}cm^{-2}$. When mouns passing through the absorber, due to its $\mu$m scale and the extremely low ionization cross section of muons, the energy residuals in the absorber are expected to be the order of $<1$ ev. The residuals is thus so different from the target energy that their effects are fully negligible. However, the muons will hit the interior of the thermostat to produce scintillation lights, and the microcalorimeter will be able to respond to such lights. Therefore, the surrounding components of TES should be minimized in designing the structure. Also, these lights can be effectively reduced by light absorbing coating.

\textbf{Impacts of naturally-occurring radionuclides:}
the absorber contains about $10^{17}$ atoms, and if the main component is copper, assuming 1ppb of natural radioactive contaminants, there are about $10^{8}$ radionuclides. According to their typical half-life of about $10^8$  years, there is about one decay event per year, thus with no contamination on the BoB measurement.

\textbf{Impacts of  $\gamma$ rays:} with a very low possibility of interactions with the absorber ($\mu$m thickness), $\gamma$ rays have an extremely high penetration power and they do not interact directly with microcalorimeters. Similar to the cosmic muons, the effects of the $\gamma$-ray induced scintillation lights are also expected to be negligibly small.

\subsection{Heat release from absorption of hydrogen atoms}\label{E}
When the hydrogen atoms interact with the absorber, some chemical energies are released at the eV level, leading to a final microcalorimeter-measured heat signal slightly higher than 325.7 eV. The shift of the spectral lines at the eV level does not affect the confirmation of the BoB mode, although the spectral line shifts of the order of eV have a significant impact on inferring the physical properties of neutrinos. In principle, this shift would be well corrected by theoretical calculations and experimental calibrations, and we will leave this for future investigations

\subsection{Design of the Microcalorimeter and Test System}\label{F}
\subsubsection{Parameters requirements of the Microcalorimeter}
\textbf{Energy resolution:}
In Sect.~\ref{A}, the use of ultracold neutrons with kinetic energies below    $10^{-7}$ eV, would hardly lead to broadening of the energy spectrum of hydrogen atoms. When hydrogen atoms interact with materials, a certain amount of chemical energy will be released, and the hydrogen atoms will eventually bind to the absorber in the form of covalent or ionic bonds. The binding process will certainly yield energy fluctuations, and thus lead to an intrinsic broadening of the amplitude of the total heat signal, with width of $\ll 1$ eV . In this sense, the energy resolution of the microcalorimeter needs to be of the order of eV or even sub-eV. Meanwhile, the higher the energy resolution of the detector, the more the number of background signals can be subtracted from the continuum energy background radiation, thus also requiring the energy resolution to be achieved as high as possible. From Eq.\ref{eq:deltaE}, the high energy resolution corresponds to a small heat capacity $C$, thus demanding a smaller microcalorimeter pixel size.  However, this conflicts with the requirement of a large absorption surface for the BoB measurement. We will perform a more detailed investigation on the effective area in a moment.

\textbf{Count rate:}
provided that the neutron flux of the neutron source is sufficiently high, the higher the single-pixel count rate of the microcalometer will of course lead to the higher the statistics achievement, thus improving the measurement accuracy and reducing measurement time. For conventional microcalometer covering energy below 10 keV, typical pulse recovery times are a few ms. As neutron decay is a quasi-random signal in time domain, there is no significant signal accumulation for single-pixel count rates of 50-100 cps. Since the target energy is only 325.7 eV, it is more appropriate to choose a smaller heat capacity $C$ of the microcalorimeter. According to the time estimate $\tau=C/G$, where $G$ is the thermal conductivity and generally a relatively constant, the typical pulse recovery time of the microcalorimeter can be on the order of sub-milliseconds, i.e., in principle the design allowing for single-pixel count rates in the range of 500-1000 cps. However, limited by the ultracold neutron flux from the current spallation neutron source, a count rate of 50 cps is sufficient for the measurement.

\textbf{Effective area:}
As mentioned above, the low flux of the current ultracold neutron source requires increasing the effective area to improve the collection efficiency of hydrogen atoms. The heat capacity $C$ of the absorber is proportional to its area for a fixed absorber thickness. Eq.~\ref{eq:deltaE} indicates that $C$ is inversely correlated with the energy resolution. As such, a large effective area help for achieving high energy resolution. Currently, the effective area of microcalorimeters with energy resolution below 1 eV is about $200\times200 \mu{\rm m}^2$~\cite{4-TES0.75eV}. To compromise the effective area and energy resolution, a larger effective area can be obtained by increasing the number of pixels, but this demands to multiplex the readout system .

\subsubsection{Parameter requirements for the microcalorimeter}

As mentioned in Sect.~\ref{D}, the TES microcalorimeter can achieve the highest energy resolution~\cite{4-TES0.75eV}, which is the most important parameter for the BoB measurement. Compared with the MMC which has a better linearity, the TES can achieve a higher multiplexing factor for SQUID, which gives a great advantage in expanding the total effective area. In short, TES microcalorimeter would be the best choice for measuring neutron two-body decay. At present, the achievement of the highest energy resolution is based on a molybdenum-copper bilayer, and therefore we will focus our investigation on it in the following.

\begin{figure}[!htb]
\includegraphics[width=1\hsize]{pic-8.png}
\caption{Dependence of the absorbed energy of two-layer films consisting of 200 nm thick copper and 80 nm thick molybdenum on the incident electron energy.}
\label{fig8}
\end{figure}

\subsubsection{Structure design of the Microcalorimeter}
\textbf{Absorber Design:}
At present, TES for particle detection mostly uses molybdenum-copper, molybdenum-gold and titanium bilayer films, or an aluminum-manganese alloy film, etc. Due to the poor penetration power of the hydrogen atoms, the thermometer itself can capture it quite well, the TES-microcalorimeter for BoB does not require the use of an additional absorber. According to the calculation, the initial kinetic energy of hydrogen atoms is so small that a material with a thickness of a few nanometers can dissipate all its energy.
Taking the molybdenum-copper bilayer TES as an example, the transition temperature below 100 mK would be obtained with a thickness of about 80 nm for molybdenum and about 200 nm for copper, which itself, in other words, can meet the measurement requirements as an absorber.

In the reality, the protons from three-body neutron decay  will have been eliminated by the 1 kV electrode, but the electrons with energy higher than 1  keV will hit the microcalorimeter and produce a signal. For electrons with energies below 6 keV, their kinetic energy is not sufficient to penetrate the molybdenum-copper bilayer, so electrons with energies below this are all deposited into the microcalorimeter. For electrons with energies higher than 6 keV, their kinetic energies are so high that they will penetrate the microcalorimeter, leaving some energy being absorbed. According to the energy loss rates of high-energy electrons (see the right panel of Fig.\ref{fig4} based on Eq.~\ref{eq:BB}), the resulting energy absorption below 782 keV in the molybdenum-copper bilayer film are derived and shown in Fig.~\ref{fig8}.
The residual energy of electrons of 782 keV is about 346 eV, about 20 eV shift from the target spectral line of 325.7 eV, thus without biasing our BoB measurement.

\begin{figure*}[!htb]
\includegraphics[width=1\hsize]{pic-9.png}
\caption{Schematic diagram of {\bf XXXX}. {\it Upper right}: the microcalorimeter chip is placed parallel to the neutron source beam and each pixel is distributed in one dimension on the chip. According to the design, each pixel contains molybdenum wires, molybdenum-copper TES, and copper for heat conduction.}
\label{fig9}
\end{figure*}

\textbf{Structure design of the microcalorimeter:}
To efficiently utilize the neutron source beam, the microcalorimeter chip needs to be placed parallel to the beam direction, and each pixel needs to be arranged in one dimension on the chip. In doing so, it helps to improve the process uniformity between pixels. For individual pixels, molybdenum-copper bilayer films are used as TES with a thickness of about 80 nm for molybdenum and about 200 nm for copper. The transition temperature of single-layer molybdenum is about 900 mK, well above the operating temperature, so molybdenum can be used as a wire, reducing micro-machining steps and thus improving process consistency. A single layer of copper is used as a thermal conductor to ensure that the heat induced by the high-energy particles is conducted out timely. In order to weld the heat sink of the refrigerator with the heat conduction wire, a solder pad is arranged at each end of the chip.

\subsection{Multiplexing method of readout electronics for TES-type microcalorimeters}\label{G}
Increasing the number of pixels in a microcalorimeter can directly increase the total effective area and significantly reduce the measurement time. The readout electronics currently in TES microcalorimeters commonly use the  superconducting quantum interference device (SQUID). However, the SQUID is costly, so a SQUID needs to read out multiple TES to increase the number of pixels. Four SQUID multiplexing techniques currently read out large TES arrays in common use: time-division multiplexing (TDM), frequency-division multiplexing (FDM), a code-division multiplexing (CDM) and microwave SQUID multiplexing (uMUX). The first three all use a DC-SQUID amplifier, converting a current change caused by the incident signal into a magnetic flux change and coupling it into a SQUID current amplifier to amplify the output. TDM is the most frequent and established multiplexed readout technology for TES array detectors at present. The TDM is fully compatible with a single-pixel TES, and it is easy to switch between it and TES array through a back-end control circuit. uMUX couples the SQUID current to SQUID loaded microwave resonators, with each resonator tuned to a different frequency, allowing the low power TES current to modulate a higher power RF signal. The uMUX technology will potentially enable ultra large-scale TES array detector readout circuits with 100,000 pixel level, playing a key role in the next-generation multiplexing readout.

\section{Summary and Outlook}\label{sec.IV}
In this paper, a microcalorimeter-based scheme for measuring the two-body neutron decay with extremely high energy resolution is proposed. The signal of the two-body neutron decay is expected to be two single-energy spectral lines at 325.7 eV and 335.9 eV, while it will be accompanied by continuum energy spectra of electrons and protons arising from the neutron three-body decay. From the estimation, this background contamination can be significantly suppressed by adding high-voltage electrodes of $\pm1$ kV. Also, muons and electrons and other particles have potential impacts on the target spectral line measurements. However, we find that the signals produced by these particles are significantly different from the target spectral lines and thus can be easily excluded. Then, the requirements of the BoB measurement for the neutron source are discussed. This experiment would ideally use an ultracold neutron source. It also requires a neutron flux of the order of $10^9$--$10^{10}$ so that the count rate collected by the absorber can reach about 10 cps, equivalent to the accumulation of thousands of neutron two-body decay events within three days. Finally, the energy resolution, count rate and other information required for the BoB are provided. Based on the microcalorimeter with the TES of molybdenum-copper bilayer film, the choice of absorber material and structural design are investigated. A readout scheme for electronic multiplexing of TES is also given, which can be used to improve the measurement sensitivity and reduce the cost.

1). {\it The primary objective} of this measurement scheme is to efficiently determine the existence of the bound state of neutron decay using the energy spectrum-based microcalorimeter, thus promising to confirm for the first time the theoretical hypothesis never observed so far.

2). {\it The next objective} is to measure the branching ratio of the two-body and the three-body decays as well as the population ratios of the ground state (1s) and the excited one ($n$s), and check whether these ratios are consistent with the theoretical expectations. The difficulty of completing such measurements does not increase significantly.

3). {\it An ambitious and long-term objective} is the direct measurement of the neutrino mass, $M_{\bar{\nu}_e}$. If the kinetic energy of $H$ can eventually be measured with accuracy to the order of $10^{-11}$ eV, the mass of the electron neutrino can be obtained. To obtain a high confidence level, this requires a tremendous amount of measurements. Therefore, it demands a massive expansion of the number of pixels in the microcalorimeter, but which is very costly.

In addition to ground-based particle physics experiments, cosmological observations currently provide important complementary information on neutrino masses, placing the strongest upper limit on neutrino mass. The large number of Big Bang neutrinos affects the evolution of the universe and provides an excellent playground for studying their physical properties over the cosmic scale of tens of billions of light years. Theoretically, the observable effects of cosmology depend mainly on the total mass of neutrinos. In the standard $\Lambda$CDM model, massive neutrinos change the expansion rate of the cosmic background, affecting the anisotropy of the cosmic microwave background (CMB) and its angular power spectrum. Meanwhile, as the universe expands, mass neutrinos, after becoming non-relativistic particles, still have large thermal velocities and strong diffusivity, which will smooth out the small-scale perturbations in the matter density, and cause a suppression on the small-scale matter power spectrum.

Thus, the total neutrino mass affects the large-scale structure of the universe and changes the gravitational lensing effect of the CMB. Currently, combining the Planck CMB data and the large-scale galaxy surveys, astronomical observations give an upper limit of 0.12 eV (95\% C.L.) for the total neutrino mass~\cite{Planck:2018vyg}. The next generation of sky observations (e.g., CMB-S4~\cite{CMB-S4:2016ple}, DESI~\cite{Font-Ribera:2013rwa}) will further improve the neutrino mass constraint, and the measurement error is expected to be narrowed to 0.03 eV (95\% C.L.). Thus, in the next decade, the combination of ground-based particle physics experiments and astronomical observations is expected to accurately determine the absolute masses of neutrinos and answer major physics questions such as mass ordering.