\section{\@startsection{section}{1
  \z@{.7\linespacing\@plus\linespacing}{.5\linespacing
  {\normalfont\scshape\centering}}
\def\subsection{\@startsection{subsection}{2
  \z@{.5\linespacing\@plus.7\linespacing}{-.5em
  {\normalfont\bfseries}}
\makeatother
\titleformat*{\section}{\large\bfseries}
\titleformat*{\subsection}{\large\bfseries}
\newtheorem{theorem}{Theorem}[section]
\newtheorem{definition}{Definition}[section]
\newtheorem{lemma}{Lemma}[section]
\newtheorem{prop}{Proposition}[section]
\newtheorem{cor}{Corollary}[section]
\theoremstyle{remark}
\newtheorem{rem}{Remark}
\renewcommand{\labelenumi}{{\upshape(\roman{enumi})}}
\title
{On a generalisation of the Riemann $\xi$-function}
\author{Hirotaka Kobayashi}
\date{}
\address{Graduate School of Mathematics, Nagoya University, Furocho, Chikusaku, Nagoya 464-8602, Japan}
\email{m17011z@math.nagoya-u.ac.jp}
\begin{document}

\maketitle

\begin{abstract}
It is known that we can construct the meromorphic function $Z_k(s)$ associated with the higher derivative of Hardy's $Z$-function.
In this paper, we introduce the entire function derived from $Z_k(s)$, a generalisation of the Riemann $\xi$-function and prove some properties.
\end{abstract}

\section{Introduction}
Hardy's $Z$-function is defined by
\begin{equation*}
Z(t)=\chi\left(\frac{1}{2}+it \right)^{-\frac{1}{2}}\zeta\left(\frac{1}{2}+it\right) \quad (t \in \mathbb{R}),
\end{equation*}
where $\chi(s)=2^s\pi^{s-1}\sin(\pi s/2)\Gamma(1-s)$ which comes from the functional equation for $\zeta(s)$.
When we take the higher derivative of $Z(t)$, it is possible to construct a meromorphic function $Z_k(s)$ which satisfies
\begin{equation*}
\left|Z^{(k)}(t) \right|=\left|Z_k\left(\frac{1}{2}+it \right) \right|.
\end{equation*}

This function has some similar properties to those of $\zeta(s)$.
Matsumoto and Tanigawa \cite{M-T} studied the zeros of $Z_k(s)$.
Later, Matsuoka \cite{M} considered the relation between the zero of $Z^{(k)}(t)$ and that of $Z^{(k+1)}(t)$.
They obtained some properties of $Z_k(s)$.
In later sections we will present their results. Especially, we reproduce the proofs of Matsuoka's results because his preprint \cite{M} is unpublished.

In this article, we consider $Z_k(s)$ rather than $Z^{(k)}(t)$. Actually, we generalise the Riemann $\xi$-function to define
\begin{equation*}
\xi_k(s):=\pi^{-\frac{s}{2}}s(s-1)\frac{Z_k(s)}{\Gamma(\frac{s}{2})^{k-1}\Gamma(\frac{1-s}{2})^k},
\end{equation*}
and give some properties of $\xi_k(s)$.
When $k=0$, this coincides with the Riemann $\xi$-function. We may expect that $\xi_k(s)$ is a natural generalisation of the Riemann $\xi$-function because the facts on $\xi_k(s)$ proved in this note are valid for $k=0$.
We will prove that this function is entire in Section 3.
For entire functions, the concept of the order is defined. The order of $\xi_k(s)$ is $1$. Thus we can factorize this function by Hadamard's factorization theorem and we can determine some constants appear in the factorized form.
We show these facts in Section 5.
In Section 6, we give the special values of $\xi_k(s)$ at positive integers. As we see in Section 2, $\xi_k(s)$ has the functional equation. It implies that we can obtain the special values at all integers.
Finally, we show Matsuoka's theorem, or more precisely, Mozer's formula by a different way from Matsuoka's proof. The factorization of $\xi_k(s)$ leads to this different way. This method is a generalisation of the method by which Edwards (\cite{Ed}) and Mozer (\cite{Mo}) (maybe independently) proved the special case of Matsuoka's result.

\begin{rem}
Actually, Matsumoto and Tanigawa considered the function denoted by $\eta_k(s)$ in their notation. However, this function is essentially the same as $Z_k(s)$.
In fact, we can show that $-\frac{1}{2}\omega(s)\eta_k(s)=Z_k(s)$.
\end{rem}

\section{The definition and basic properties of $Z_k(s)$}
Let $Z_0(s)=\zeta(s)$, and according to Matsuoka \cite{M} we define $Z_k(s)$ for $k\geq 1$ by
\begin{equation}\label{def}
Z_{k+1}(s)=Z_k'(s)-\frac{1}{2}\omega(s)Z_k(s) \quad (k\geq 0),
\end{equation}
where
\begin{equation*}
\omega(s):=\frac{\chi'}{\chi}(s)=\log 2\pi-\frac{\Gamma'}{\Gamma}(s)+\frac{\pi}{2}\tan\frac{\pi s}{2}.
\end{equation*}

\begin{rem}
The function $Z_k(s)$ was already introduced in 1990 by Y\i ld\i r\i m \cite{Y} in the following alternative way:
\begin{equation}\label{Y's-def}
Z_k(s):=(\chi(s))^{\frac{1}{2}}\frac{d^k}{ds^k}\left((\chi(s))^{-\frac{1}{2}}\zeta(s) \right).
\end{equation}
We can see that this definition coincides with Matsuoka's.
\end{rem}

We can see that
\begin{equation*}
\omega(s)=\omega(1-s).
\end{equation*}

When we put
\begin{equation*}
\theta(t):=\arg \Gamma \left(\frac{1}{4}+\frac{it}{2} \right)-\frac{t}{2}\log \pi
\end{equation*}

then we see that $\chi(1/2+it)=e^{-2i\theta(t)}$, hence
\begin{equation}\label{omegatheta}
\omega \left(\frac{1}{2}+it \right)=-2\theta'(t).
\end{equation}

\begin{prop}[Proposition 2.1 in \cite{M}]
For any non-negative $k$, we have
\begin{equation}\label{Z^k}
Z^{(k)}(t)=i^{k}Z_k\left(\frac{1}{2}+it \right)e^{i\theta(t)}.
\end{equation}
\end{prop}

\begin{proof}
The case $k=0$ is the definition of $Z(t)$.
If we assume that the equation is true for $k$, then
\begin{equation}
Z^{(k+1)}(t)=i^{k+1}e^{i\theta(t)}\left(Z_k'\left(\frac{1}{2}+it\right)+\theta'(t)Z_k\left(\frac{1}{2}+it \right)\right).
\end{equation}
By $(\ref{def})$ and $(\ref{omegatheta})$, we find that the equation is true for $k+1$.
\end{proof}

\begin{prop}[The Functional Equation, Proposition 2.2 in \cite{M}]
For any non-negative $k$, we have
\begin{equation}\label{f-e}
\chi(s)Z_k(1-s)=(-1)^k Z_k(s).
\end{equation}
\end{prop}
\begin{proof}
The case $k=0$ is the functional equation of the Riemann $\zeta$-function.
If we assume that the equation is true for $k$, then by the definition,
\begin{align*}
\chi(s)Z_{k+1}(1-s) &= \chi(s)\left(Z_{k}'(1-s)-\frac{1}{2}\omega(1-s)Z_k(1-s) \right) \\
&= \chi'(s)Z_k(1-s)-(-1)^{k}Z_k'(s)-\frac{(-1)^k}{2}\omega(s)Z_k(s) \\
&= (-1)^{k+1}Z_k'(s)+(-1)^k\omega(s)Z_k(s)-\frac{(-1)^k}{2}\omega(s)Z_k(s) \\
&= (-1)^{k+1}\left(Z_k'(s)-\frac{1}{2}\omega(s)Z_k(s) \right) \\
&= (-1)^{k+1}Z_{k+1}(s).
\end{align*}
The proof is completed.
\end{proof}

Therefore we can obtain the following theorem.
\begin{theorem}
For all $s$, we have
\begin{equation}\label{fcte-xi}
\xi_k(s)=(-1)^k\xi_k(1-s).
\end{equation}
\end{theorem}
\begin{proof}
We can transform $\chi(s)$ to
\begin{equation*}
\chi(s)=\frac{\Gamma (\frac{1-s}{2})}{\Gamma(\frac{s}{2})}\pi^{s-\frac{1}{2}}.
\end{equation*}
Thus
\begin{align*}
\xi_k(s)
&=
\pi^{-\frac{s}{2}}s(s-1)\frac{Z_k(s)}{\Gamma(\frac{s}{2})^{k-1}\Gamma(\frac{1-s}{2})^k} \\
&=
\pi^{-\frac{s}{2}}s(s-1)\frac{(-1)^k\chi(s)Z_k(1-s)}{\Gamma(\frac{s}{2})^{k-1}\Gamma(\frac{1-s}{2})^k} \\
&=
(-1)^k\pi^{-\frac{1-s}{2}}s(s-1)\frac{Z_k(1-s)}{\Gamma(\frac{1-s}{2})^{k-1}\Gamma(\frac{s}{2})^k} \\
&=
(-1)^k\xi_k(1-s).
\end{align*}
\end{proof}

Let $f_0(s)=1$, and define $f_k(s)$ for $k\geq 1$ by
\begin{equation}
f_{k+1}(s)=f_k'(s)-\frac{1}{2}\omega(s)f_k(s) \quad (k\geq 1).
\end{equation}

The following proposition is inspired by  Proposition 3 in \cite{M-T}.
\begin{prop}\label{exZ_k}
For any non-negative $k$, we have
\begin{equation}
Z_k(s)=\sum_{j=0}^{k}\binom{k}{j}f_{k-j}(s)\zeta^{(j)}(s).
\end{equation}
\end{prop}

\begin{proof}
The case $k=0$ is clear. We assume that this is valid for $k$. By the definition,
\begin{align*}
Z_{k+1}(s)&=Z_k'(s)-\frac{1}{2}\omega(s)Z_k(s)\\
&=\sum_{j=0}^{k}\binom{k}{j}f_{k-j}'(s)\zeta^{(j)}(s)+\sum_{j=0}^{k}\binom{k}{j}f_{k-j}(s)\zeta^{(j+1)}(s) \\
&\quad-\frac{1}{2}\omega(s)\sum_{j=0}^{k}\binom{k}{j}f_{k-j}(s)\zeta^{(j)}(s)\\
&=\sum_{j=0}^{k}\binom{k}{j}f_{k+1-j}(s)\zeta^{(j)}(s)+\sum_{j=0}^{k}\binom{k}{j}f_{k-j}(s)\zeta^{(j+1)}(s)\\
&=\sum_{j=0}^{k+1}\left\{\binom{k}{j}+\binom{k}{j-1}\right\}f_{k+1-j}(s)\zeta^{(j)}(s)\\
&=\sum_{j=0}^{k+1}\binom{k+1}{j}f_{k-j}(s)\zeta^{(j)}(s).
\end{align*}
Here to obtain the last equality, we use the relation
\begin{equation*}
\binom{k}{j}+\binom{k}{j-1}=\binom{k+1}{j}.
\end{equation*}
\end{proof}

The function $f_k(s)$ can be expressed explicitly as follows.

\begin{prop}\label{exf_k}
For $k\geq 1$, we have
\begin{equation*}
f_k(s)=k!\sum_{a_1+2a_2+\dots+ka_k=k}\left(-\frac{1}{2} \right)^{a_1+\dots+a_k}\prod_{l=1}^{k}\frac{1}{a_{l}!}\left(\frac{\omega^{(l-1)}(s)}{l!} \right)^{a_l}.
\end{equation*}
\end{prop}
We can guess this expression by (\ref{Y's-def}) and Fa\`{a} di Bruno's formula (see \cite{J}, \cite{R}).

\begin{proof}
We will prove this statement by induction on $k$. When $k=1$, this formula is trivial.
We assume that this is true for $k=j$. By the definition, we have
\begin{align*}
f_{j+1}(s)
&=
f_j'(s)-\frac{1}{2}\omega(s)f_j(s) \\
&=
j!\sum_{a_1+2a_2+\dots+ja_j=j}\left(-\frac{1}{2} \right)^{a_1+\dots+a_j}\sum_{\substack{m=1 \\ a_m\geq 1}}^{j}\prod_{\substack{l=1 \\ l\neq m}}^{k}\frac{1}{a_{l}!}\left(\frac{\omega^{(l-1)}(s)}{l!} \right)^{a_l} \\
&\quad
\times \frac{1}{(a_m-1)!}\left(\frac{1}{m!} \right)^{a_m}\omega^{(m-1)}(s)^{a_m-1}\omega^{(m)}(s) \\
&\quad
-\frac{1}{2}\omega(s)j!\sum_{a_1+2a_2+\dots+ja_j=j}\left(-\frac{1}{2} \right)^{a_1+\dots+a_j}\prod_{l=1}^{j}\frac{1}{a_{l}!}\left(\frac{\omega^{(l-1)}(s)}{l!} \right)^{a_l}.
\end{align*}
On the first term, let $b_m=a_m-1, b_{m+1}=a_{m+1}+1$ and $b_l=a_l \ (l\neq m, m+1)$ for each $m \ (1\leq m\leq j-1)$. When $m=j$, we put $b_j=a_j-1$ and $b_{j+1}=1$. Then we can see that the first term is
\begin{align*}
&\quad
j!\sum_{b_1+2b_2+\dots+jb_j=j+1}\left(-\frac{1}{2} \right)^{b_1+\dots+b_j}\sum_{m=1}^{j-1}(m+1)b_{m+1}\prod_{l=1}^{j}\frac{1}{b_{l}!}\left(\frac{\omega^{(l-1)}(s)}{l!} \right)^{b_l} \\
&\quad
+(j+1)!\sum_{\substack{b_1+2b_2+\dots+(j+1)b_{j+1}=j+1 \\ b_{j+1}=1}}\left(-\frac{1}{2} \right)^{b_1+\dots+b_{j+1}}\prod_{l=1}^{j+1}\frac{1}{b_{l}!}\left(\frac{\omega^{(l-1)}(s)}{l!} \right)^{b_l} \\
&=
(j+1)!\sum_{b_1+2b_2+\dots+(j+1)b_{j+1}=j+1}\left(-\frac{1}{2} \right)^{b_1+\dots+b_{j+1}}\prod_{l=1}^{j+1}\frac{1}{b_{l}!}\left(\frac{\omega^{(l-1)}(s)}{l!} \right)^{b_l} \\
&\quad
-j!\sum_{b_1+2b_2+\dots+(j+1)b_{j+1}=j+1}\left(-\frac{1}{2} \right)^{b_1+\dots+b_{j+1}}b_1\prod_{l=1}^{j+1}\frac{1}{b_{l}!}\left(\frac{\omega^{(l-1)}(s)}{l!} \right)^{b_l}.
\end{align*}
As for the second term, when we put $c_1=a_1+1$ and $c_l=a_l \quad (2\leq l\leq j)$ then
\begin{align*}
j!\sum_{c_1+2c_2+\dots+jc_j=j+1}\left(-\frac{1}{2} \right)^{c_1+\dots+c_j}c_1\prod_{l=1}^{j}\frac{1}{c_{l}!}\left(\frac{\omega^{(l-1)}(s)}{l!} \right)^{c_l}.
\end{align*}
Therefore
\begin{align*}
&\quad
f_{j+1}(s) \\
&=
(j+1)!\sum_{a_1+2a_2+\dots+(j+1)a_{j+1}=j+1}\left(-\frac{1}{2} \right)^{a_1+\dots+a_{j+1}}\prod_{l=1}^{j+1}\frac{1}{a_{l}!}\left(\frac{\omega^{(l-1)}(s)}{l!} \right)^{a_l}.
\end{align*}
This completes the proof.
\end{proof}

\section{The poles of $f_k(s)$ and $Z_k(s)$}
We investigate the poles of $f_k(s)$ and $Z_k(s)$.

\begin{lemma}[Lemma 1 in \cite{M-T}]\label{pole-omega}
The poles of $\omega(s)$ are all simple, and located at positive odd integers with residue $-1$ and at non-positive even integers with residue $1$.
\end{lemma}

\begin{proof}
Since $\chi(s)=2^s\pi^{s-1}\sin (\pi s/2)\Gamma(1-s)$, the zeros of $\chi(s)$ are non-positive even integers and the poles are positive odd integers. Therefore we obtain the lemma.
\end{proof}

\begin{lemma}[Lemma 4.2 in \cite{M}]\label{f_poles}
For $k\geq 0$, the function $f_k(s)$ has poles of order $k$ which are located only at $\dots, -4,-2,0,1,3,5, \dots$.
\end{lemma}

\begin{proof}
The case $k=1$ is the previous lemma.
We assume that the lemma is valid for $k\geq 1$.
Let $a$ be a pole of $f_k(s)$. Then by the Laurent expansion at centre $a$,
we have
\begin{equation*}
f_k(s)=\frac{c_k}{(s-a)^k}+\cdots,
\end{equation*}
where $c_k$ does not vanish.
By the definition and the previous lemma, we have
\begin{equation*}
f_{k+1}(s)=\frac{-kc_k\pm\frac{c_k}{2}}{(s-a)^{k+1}}+\cdots.
\end{equation*}
Since $-kc_k\pm c_k/2\neq 0$, the lemma is true for $k+1$. This completes the proof of the lemma.
\end{proof}

	This lemma and Proposition \ref{exZ_k} immediately lead to the following lemma.
\begin{lemma}[Lemma 4.4 in \cite{M}]\label{poles}
For $k\geq 0$, the function $Z_k(s)$ has poles of order $k$ located at $0,3,5,7, \cdots$,
that of order $k+1$ located at $1$, and those of order $k-1$ located at $-2,-4,-6, \cdots$.
\end{lemma}
We understand that ``poles of order $-1$'' means zeros of order $1$.

This lemma implies $\xi_k(s)$ is entire, because $\Gamma(s/2)$ and $\Gamma((1-s)/2)$ have poles at non-positive even integers and positive odd integers respectively.

\section{Some estimates on $Z_k(s)$}
In this section, we prove some estimates on $Z_k(s)$.
Proposition \ref{exZ_k} implies that we need the estimates of $f_k(s)$ and $\zeta^{(k)}(s)$ for any non-negative integer $k$.

First we consider $f_k(s)$.
Stirling's formula implies that
\begin{lemma}\label{estgam}
For $\sigma \geq \frac{1}{4}$, we have

\begin{equation*}
\Gamma(s)=\sqrt{2\pi}s^{s-\frac{1}{2}}e^{-s}(1+O(|s|^{-1})),
\end{equation*}

\begin{equation*}
\frac{\Gamma'}{\Gamma}(s)=\log s+O(|s|^{-1}) \label{degamma}
\end{equation*}
and
\begin{equation*}
\frac{d^k}{ds^k}\frac{\Gamma'}{\Gamma}(s)=O_k(|s|^{-k}) \quad (k\geq 1).
\end{equation*}
\end{lemma}

Define the set $\mathscr{D}$ by removing all small circles whose centres are odd positive integers and even non-positive integers with radii depending on $k$ from the complex plane.

\begin{lemma}[Lemma 3.1 in \cite{M}]\label{esttan}
In the region $\{s \in \mathscr{D}|\sigma>1/4 \}$, we have
\begin{equation}\label{tan}
\tan \frac{\pi s}{2}=\begin{cases}
           i+O(e^{-2t}) & \text{$(t\geq 0)$}, \\
           -i+O(e^{2t}) & \text{$(t\leq 0)$},
           \end{cases}
\end{equation}
and
\begin{equation}\label{hightan}
\frac{d^k}{ds^k}\tan \frac{\pi s}{2}=O_k(e^{-2|t|}) \quad(k\geq 1).
\end{equation}
\end{lemma}

By this lemma and Lemma $\ref{estgam}$, we have
\begin{lemma}[Lemma 3.2 in \cite{M}]\label{estome}
For $s \in \mathscr{D}$, we have
\begin{equation}\label{omega}
\omega(s)=-\log |s|+O(1),
\end{equation}
and
\begin{equation}\label{highomega}
\omega^{(k)}(s)=O_k(1) \quad(k\geq 1).
\end{equation}
\end{lemma}

Therefore we can see that $f_k(s)\ll_k (\log (|s|+2))^k$ for $s \in \mathscr{D}$ by Proposition \ref{exf_k}.

By Lemma \ref{f_poles},
\begin{equation*}
\left(\cos \frac{\pi s}{2}\right)^kf_k(s)
\end{equation*}
has no pole in the half plane $\Re s\geq \frac{1}{2}$.

From the above, for $\Re s\geq \frac{1}{2}$ there is a constant $C_1=C_1(k)$ such that
\begin{equation}\label{f-est}
\left(\cos \frac{\pi s}{2}\right)^kf_k(s)\ll_k e^{C_1|s|}.
\end{equation}

Next we consider $\zeta^{(k)}(s)$.
Let $k=0$.
It is known that
\begin{equation*}
\zeta(s)=\frac{1}{s-1}-s\int_{1}^{\infty}\frac{x^2-[x]-\frac{1}{2}}{x^{s+1}}dx+\frac{1}{2}.
\end{equation*}
Therefore we can see that when $s$ is not close to $1$,
\begin{equation*}
\zeta(s)\ll_{\eta} |s| \quad (\Re s >\eta)
\end{equation*}
for any $\eta>0$.
Therefore, by Cauchy's integral formula applied to a circle which has centre $s$ and radius $1/\log (|s|+2)$,

\begin{equation*}
\zeta^{(k)}(s)\ll_{\eta'} |s|\log^k (|s|+2) \quad (\Re s>\eta')
\end{equation*}
for $k\geq 1$ and any $\eta'=\eta'(k)>0$ unless $s$ is close to $1$.
Thus for all $\Re s>\eta'$, there is a constant $C_2=C_2(k)$ such that

\begin{equation}\label{zeta-est}
\left(\cos \frac{\pi s}{2} \right)^k(s-1)\zeta^{(k)}(s)\ll_k e^{C_2|s|}.
\end{equation}

\section{The factorization of $\xi_k(s)$}
By the estimates showed in the previous section, we can prove the following theorem.
\begin{theorem}\label{xi-facto}
For $k\geq 0$, there are constants $a_k$ and $B_k$ such that

\begin{equation}\label{factral-xi}
\xi_k(s)=e^{A_k+B_ks}\prod_{\rho_k}\left(1-\frac{s}{\rho_k} \right)e^{\frac{s}{\rho_k}}
\end{equation}

for all $s$. Here the product is extended over all zeros $\rho_k$ of $\xi_k(s)$.

\end{theorem}

We note that it is possible that some of the $\rho_k$'s are real. Actually, Anderson (\cite{An}) proved the existence of real zeros for $k=1$. This fact is reproduced by Hall (\cite{H}).

\begin{proof}
It is sufficient to prove that the order of $\xi_k(s)$ is $1$, for then Hadamard's factorization theorem can be applied.

Let $\Re s \geq \frac{1}{2}$.
By the formulas
\begin{equation*}
\Gamma(s)\Gamma(1-s)=\frac{\pi}{\sin \pi s} \quad \text{and} \quad \Gamma(s)\Gamma \left(s+\frac{1}{2} \right)=\sqrt{\pi}2^{1-2s}\Gamma(2s),
\end{equation*}
we see that
\begin{equation*}
\Gamma \left(\frac{s}{2}\right)\Gamma \left(\frac{1-s}{2} \right)=2^{s-1}\sqrt{\pi}\frac{\Gamma(\frac{s}{2})^2}{\Gamma(s)\cos \frac{\pi s}{2}}.
\end{equation*}

Then, by the definition of $\xi_k(s)$ and Proposition \ref{exZ_k},
\begin{align*}
\xi_k(s)
&=
2^{-k(s-1)}\pi^{-\frac{s+k}{2}}s(s-1)\Gamma \left(\frac{s}{2} \right)\left(\frac{\Gamma(s)}{\Gamma \left(\frac{s}{2} \right)^2} \right)^k\left(\cos \frac{\pi s}{2}\right)^k \\
&\quad
\times\sum_{j=0}^{k}\binom{k}{j}f_{k-j}(s)\zeta^{(j)}(s).
\end{align*}

Lemma \ref{estgam} implies that
\begin{equation}\label{xi_app}
\begin{split}
\xi_k(s)
&=
2^{-\frac{s+k}{2}+1}\pi^{-\frac{s-1}{2}-k}s^{\frac{s+k+1}{2}}e^{-\frac{s}{2}} \\
&\quad
\times\sum_{j=0}^{k}\binom{k}{j}\left(\cos \frac{\pi s}{2} \right)^kf_{k-j}(s)(s-1)\zeta^{(j)}(s)(1+O_k(|s|^{-1})).
\end{split}
\end{equation}
Finally, by the estimates in the previous section (\ref{zeta-est}) and (\ref{f-est}), we have
\begin{equation*}
\xi_k(s) \ll_k \exp(C|s|\log (|s|+2)),
\end{equation*}
where the constant $C$ depends only on $k$. The functional equation for $\xi_k(s)$ (\ref{fcte-xi}) implies that this estimate is valid for $\Re s\leq \frac{1}{2}$ too.
Thus the order of $\xi_k(s)$ is at most $1$.
By (\ref{xi_app}), we have $\log \xi_k(\sigma)\sim \frac{1}{2}\sigma\log \sigma$ for $\sigma\rightarrow \infty$.
This completes the proof.
\end{proof}

We can determine constants $e^{A_k}$ and $B_k$ in (\ref{factral-xi}).
\begin{theorem}
In the previous theorem, for $k\geq 0$ we have
\begin{equation*}
e^{A_k}=\xi_k(0)=\frac{(-1)^k(2k-1)!!}{(4\sqrt{\pi})^k}
\end{equation*}
and
\begin{equation*}
B_k=-\frac{2k(k-1)}{2k-1}\log 2-\frac{1}{2(2k-1)}\log 4\pi+\frac{\gamma}{2(2k-1)}-1,
\end{equation*}
where $\gamma$ is the Euler constant.
\end{theorem}

Here, we note that
\begin{equation*}
(2k-1)!!=\prod_{l=1}^k(2l-1)=(2k-1)(2k-3)\cdots3\cdot2\cdot1,
\end{equation*}
and we define $(-1)!!=1$ and $(-3)!!=-1$.

\begin{proof}
The logarithmic derivative leads to
\begin{equation*}
\frac{\xi_k'}{\xi_k}(0)=B_k,
\end{equation*}
by the previous theorem. Thus, by the definition of $\xi_k(s)$,
\begin{equation*}
B_k=
-\frac{1}{2}\log \pi-1+\frac{k}{2}\frac{\Gamma'}{\Gamma}\left(\frac{1}{2} \right)+\lim_{s\rightarrow 0}\left(\frac{1}{s}-\frac{k-1}{2}\frac{\Gamma'}{\Gamma}\left(\frac{s}{2} \right)+\frac{Z_k'}{Z_k}(s)\right).
\end{equation*}

Since
\begin{equation*}
\frac{1}{s}+\frac{\Gamma'}{\Gamma}(s)=\frac{\Gamma'}{\Gamma}(s+1),
\end{equation*}
we have
\begin{align*}
\frac{1}{s}-\frac{k-1}{2}\frac{\Gamma'}{\Gamma}\left(\frac{s}{2} \right)+\frac{Z_k'}{Z_k}(s)
&=
\frac{k}{s}+\frac{Z_k'}{Z_k}(s)-\frac{k-1}{2}\left(\frac{2}{s}+\frac{\Gamma'}{\Gamma}\left(\frac{s}{2} \right)\right) \\
&=
\frac{k}{s}+\frac{Z_k'}{Z_k}(s)-\frac{k-1}{2}\frac{\Gamma'}{\Gamma}\left(\frac{s}{2}+1 \right).
\end{align*}

We consider the logarithmic derivative of $Z_k(s)$.
We recall $Z_k(s)$ has a pole with order $k$ at $s=0$.
Then we can express
\begin{equation*}
Z_k(s)=\sum_{l=-k}^{\infty}a_{k,l}s^l.
\end{equation*}

Then we obtain
\begin{align*}
\frac{k}{s}+\frac{Z_k'}{Z_k}(s)
&=
\frac{k}{s}+\frac{-ka_{k,-k}s^{-k-1}-(k-1)a_{k,-k+1}s^{-k}+O(|s|^{-k+1})}{a_{k,-k}s^{-k}+a_{k,-k+1}s^{-k+1}+O(|s|^{-k+2})} \\
&=
\frac{ka_{k,-k}s^{-k-1}+ka_{k,-k+1}s^{-k}+O(|s|^{-k+1})}{a_{k,-k}s^{-k}+a_{k,-k+1}s^{-k+1}+O(|s|^{-k+2})} \\
&\quad
+\frac{-ka_{k,-k}s^{-k-1}-(k-1)a_{k,-k+1}s^{-k}+O(|s|^{-k+1})}{a_{k,-k}s^{-k}+a_{k,-k+1}s^{-k+1}+O(|s|^{-k+2})} \\
&=
\frac{a_{k,-k+1}+O(|s|)}{a_{k,-k}+a_{k,-k+1}s+O(|s|^{2})}.
\end{align*}

Therefore
\begin{equation*}
\lim_{s\rightarrow 0}\left(\frac{k}{s}+\frac{Z_k'}{Z_k}(s)\right)=\frac{a_{k,-k+1}}{a_{k,-k}}.
\end{equation*}
By the definition (\ref{def}) and Lemma \ref{pole-omega}, we can see that for $k\geq 1$
\begin{equation*}
a_{k,-k}=-(k-1)a_{k-1,-k+1}-\frac{1}{2}a_{k-1,-k+1}=\left(-k+\frac{1}{2} \right)a_{k-1,-k+1},
\end{equation*}
(see the argument in the proof of Lemma \ref{f_poles}) and
\begin{align*}
\frac{a_{k,-k}}{a_{0,0}}
&=
\frac{a_{k,-k}}{a_{k-1,-k+1}}\cdot\frac{a_{k-1,-k+1}}{a_{k-2,-k+2}}\cdots\frac{a_{1,-1}}{a_{0,0}} \\
&=
\left(-k+\frac{1}{2} \right)\left(-k+\frac{3}{2} \right)\cdots \left(-\frac{1}{2}\right)=\left(-\frac{1}{2}\right)^{k}(2k-1)!!.
\end{align*}
Since $a_{0,0}=\zeta(0)=-1/2$, we have
\begin{equation*}
a_{k,-k}=\left(-\frac{1}{2}\right)^{k+1}(2k-1)!!.
\end{equation*}

On the other hand, by the definition (\ref{def}) and Lemma \ref{pole-omega}, we can obtain
\begin{equation*}
a_{k,-k+1}=-\frac{1}{2}(2k-3)a_{k-1,-k+2}-\frac{C_0}{2}a_{k-1,-k+1},
\end{equation*}
where $C_0$ is the constant term of the Laurent expansion of $\omega(s)$ around $0$ and it is $\log 2\pi+\gamma$.
When we put
\begin{equation*}
a'_k=\frac{a_{k,-k+1}}{-\frac{C_0}{2}a_{k-1,-k+1}}
\end{equation*}
we have
\begin{equation*}
a'_k=a'_{k-1}+1=a'_1+k-1.
\end{equation*}
Thus we can see that
\begin{equation*}
a_{k,-k+1}=\left(\frac{\gamma}{\log 2\pi+\gamma}+k-1 \right)(\log 2\pi+\gamma)\left(-\frac{1}{2}\right)^{k+1}\frac{(2k-1)!!}{2k-1}
\end{equation*}
by the fact that $a_{1,0}=(\zeta'(0)- \zeta(0)\log 2\pi- \zeta(0)\gamma)/2=\gamma/4$.

Hence we have
\begin{align*}
B_k
&=
-\frac{1}{2}\log \pi-1+\frac{k}{2}\frac{\Gamma'}{\Gamma}\left(\frac{1}{2} \right)-\frac{k-1}{2}\frac{\Gamma'}{\Gamma}(1) \\
&\quad
+\left(\frac{\gamma}{\log 2\pi+\gamma}+k-1 \right)\frac{\log 2\pi+\gamma}{2k-1} \\
&=
-\frac{2k(k-1)}{2k-1}\log 2-\frac{1}{2(2k-1)}\log 4\pi+\frac{\gamma}{2(2k-1)}-1.
\end{align*}
Here, to show the last equality, we use
\begin{equation*}
\frac{\Gamma'}{\Gamma}(1)=-\gamma \quad \text{and} \quad \frac{\Gamma'}{\Gamma}\left(\frac{1}{2} \right)=-\gamma-2\log 2.
\end{equation*}

Lastly, we calculate $\xi_k(0)=e^{A_k}$.
It is known that
\begin{equation*}
\Gamma(s+1)=s\Gamma(s).
\end{equation*}
Hence we obtain
\begin{align*}
\xi_k(0)
&=
-\Gamma \left(\frac{1}{2} \right)^{-k}\lim_{s\rightarrow 0}\frac{sZ_k(s)}{\Gamma \left(\frac{s}{2} \right)^{k-1}} \\
&=
-\Gamma \left(\frac{1}{2} \right)^{-k}\lim_{s\rightarrow 0}\left(\frac{s}{2}\right)^{k-1}\frac{sZ_k(s)}{\Gamma \left(\frac{s}{2}+1 \right)^{k-1}} \\
&=
-\frac{2^{-k+1}}{\Gamma(\frac{1}{2})^k}\lim_{s\rightarrow 0}s^kZ_k(s)=-2^{-k+1}\pi^{-\frac{k}{2}}a_{k,-k}=\frac{(-1)^k(2k-1)!!}{(4\sqrt{\pi})^k}.
\end{align*}

\end{proof}

\section{The special values of $\xi_k(s)$ at integer points}
In this section, we give an explicit expression of special values of $\xi_k(s)$ at positive integer points.
\begin{theorem}\label{sp-value}
Let $k\geq 0$. We have
\begin{equation*}
\xi_k(1)=\frac{(2k-1)!!}{(4\sqrt{\pi})^k},
\end{equation*}
and for $n\geq 1$,
\begin{align*}
&\quad
\xi_k(2n+1) \\
&=
(-1)^{kn+1}\pi^{-\frac{2n+k}{2}}\frac{(2n+1)!2n}{4^n\cdot n!}\left(\frac{4^{n-1}(n!)^2}{(2n)!}\right)^k(2k-3)!!\zeta(2n+1)
\end{align*}
and
\begin{align*}
\xi_k(2n)
&=
(-1)^{kn}\pi^{-\frac{2n+k}{2}}2n(2n-1)(n-1)!\left(\frac{(2n)!}{4^nn!(n-1)!}\right)^kZ_k(2n).
\end{align*}
\end{theorem}

Before starting the proof, we note that Proposition \ref{exZ_k} and \ref{exf_k} imply
\begin{align*}
Z_k(2n)
&=
\sum_{j=0}^{k}\frac{k!}{j!}\zeta^{(j)}(2n)\sum_{a_1+2a_2+\dots+(k-j)a_{k-j}=k-j}\left(-\frac{1}{2} \right)^{a_1+\dots+a_{k-j}} \\
&\quad
\times \prod_{l=1}^{k-j}\frac{1}{a_{l}!}\left(\frac{\omega^{(l-1)}(2n)}{l!} \right)^{a_l}.
\end{align*}
On the $\omega^{(l-1)}(2n)$, since
\begin{equation*}
\left.\frac{d^{l-1}}{ds^{l-1}}\tan \frac{\pi s}{2}\right|_{s=2n}
=
\begin{cases}
\displaystyle
\frac{(-1)^{m-1}B_{2m}(4^m-1)\pi^{2m-1}}{2m} & l=2m, \\
0 & l=2m+1,
\end{cases}
\end{equation*}
and for $l\geq 2$
\begin{equation*}
\left.\frac{d^{l-1}}{ds^{l-1}}\frac{\Gamma'}{\Gamma}(s)\right|_{s=2n}=(-1)^l(l-1)!\left(\zeta(l)-\sum_{j=1}^{2n-1}\frac{1}{j^l} \right),
\end{equation*}
we have
\begin{equation*}
\omega^{(l-1)}(2n)
=
\begin{cases}
\displaystyle
(-1)^m\frac{B_{2m}(2-4^m)\pi^{2m}}{4m}+(2m-1)!\sum_{j=1}^{2n-1}\frac{1}{j^{2m}} & l=2m, \\
\displaystyle
-(2m)!\left(\sum_{j=1}^{2n-1}\frac{1}{j^{2m+1}}-\zeta(2m+1) \right) & l=2m+1,
\end{cases}
\end{equation*}
where $B_{2m}$ is the $2m$-th Bernoulli number and therefore $\omega^{(l-1)}(2n)>0$.
When $l=1$, we can obtain
\begin{equation*}
\omega(2n)=\log 2\pi+\gamma-\sum_{j=1}^{2n-1}\frac{1}{j}
\end{equation*}
because
\begin{equation*}
\frac{\Gamma'}{\Gamma}(n)=-\gamma+\sum_{j=1}^{n-1}\frac{1}{j}.
\end{equation*}

\begin{proof}[Proof of Theorem \ref{sp-value}]
The case of $s=1$ is trivial.

Let $n\geq 1$. First we consider the case $s=2n$. When $s=2n$ we have
\begin{align*}
\xi_k(2n)
&=
\pi^{-n}2n(2n-1)\frac{Z_k(2n)}{\Gamma(n)^{k-1}\Gamma(\frac{1}{2}-n)^k} \\
&=
(-1)^{kn}\pi^{-\frac{2n+k}{2}}2n(2n-1)(n-1)!\left(\frac{(2n)!}{4^nn!(n-1)!}\right)^kZ_k(2n).
\end{align*}

Next we treat the case $s=2n+1$. In this case, we see that
\begin{align*}
\xi_k(2n+1)
&=
\pi^{-\frac{2n+1}{2}}(2n+1)2n\frac{1}{\Gamma(\frac{2n+1}{2})^{k-1}}\lim_{s\rightarrow 2n+1}\frac{Z_k(s)}{\Gamma(\frac{1-s}{2})^k} \\
&=
\pi^{-\frac{2n+1}{2}}(2n+1)2n\frac{(4^nn!)^{k-1}}{((2n)!)^{k-1}\pi^{\frac{k-1}{2}}} \\
&\quad
\times \lim_{s\rightarrow 2n+1}\frac{\prod_{m=0}^{n-1}(2m+1-s)^k}{2^{kn+k}\Gamma(\frac{2n+3-s}{2})^k}(2n+1-s)^kZ_k(s) \\
&=
(-1)^{kn+k}\pi^{-\frac{2n+k}{2}}(2n+1)2n\frac{(n!)^k(4^nn!)^{k-1}}{2^k((2n)!)^{k-1}} \\
&\quad
\times \lim_{s\rightarrow 2n+1}(s-2n-1)^kZ_k(s).
\end{align*}

We consider the Laurent expansion at $s=2n+1$ i.e.
\begin{equation*}
Z_k(s)=\sum_{l=-k}^{\infty}b_{k,l}(s-2n-1)^l.
\end{equation*}
By the definition (\ref{def}) and Lemma \ref{pole-omega}, we can see that
\begin{equation*}
b_{k,-k}=-(k-1)b_{k-1,-k+1}+\frac{1}{2}b_{k-1,-k+1}=-\left(k-\frac{3}{2} \right)b_{k-1,-k+1}.
\end{equation*}
Thus we have
\begin{align*}
\frac{b_{k,-k}}{b_{0,0}}
&=
\frac{b_{k,-k}}{b_{k-1,-k+1}}\cdot \cdots \cdot \frac{b_{1,-1}}{b_{0,0}} \\
&=
\left(-k+\frac{3}{2} \right)\cdots \left(-\frac{1}{2} \right)\frac{1}{2}=-\left(-\frac{1}{2} \right)^{k}(2k-3)!!.
\end{align*}
Since $b_{0,0}=\zeta(2n+1)$, we obtain
\begin{equation*}
b_{k,-k}=-\left(-\frac{1}{2}\right)^k(2k-3)!!\zeta(2n+1).
\end{equation*}
This leads to
\begin{align*}
&\quad
\xi_k(2n+1) \\
&=
(-1)^{kn+1}\pi^{-\frac{2n+k}{2}}\frac{(2n+1)!2n}{4^n\cdot n!}\left(\frac{4^{n-1}(n!)^2}{(2n)!}\right)^k(2k-3)!!\zeta(2n+1).
\end{align*}

\end{proof}

\section{Necessary Additional statements}
To show Matsuoka's theorem (Theorem \ref{M-th}), we need some additional statements.

By Proposition \ref{exf_k} and Lemma \ref{estome},
\begin{equation}\label{f_kest}
f_k(s)=\left(-\frac{\omega(s)}{2}\right)^k(1+O_k((\log |s|)^{-2}))
\end{equation}
for $s\in \mathscr{D}(2m)=\{ s\in \mathscr{D} \mid \sigma \geq 2m \}$ with sufficiently large integer $m=m(k)$.
This implies $\Re f_k(s)>0$ in this region.
Thus, by the argument principle, we have

\begin{lemma}[Lemma 4.3 in \cite{M}]\label{f_zeros-outside}
For $s\in \{ s \mid \sigma<1-2m \} \cup \{ s \mid \sigma>2m \}$, the zeros of $f_k(s)$ are in small circles centred negative even integers and positive odd integers,
and the number of those is $k$ in each circles.
\end{lemma}

By Proposition \ref{exZ_k}, (\ref{f_kest}) and the facts that
\begin{equation*}
\zeta(s)=1+O(2^{-\sigma}), \quad \zeta^{(k)}=O_k(2^{-\sigma}),
\end{equation*}
we see that
\begin{align*}
Z_k(s)
&=
f_k(s)\left \{\zeta(s)+\sum_{j=1}^{k}\binom{k}{j}\frac{f_{k-j}(s)}{f_k(s)}\zeta^{(j)}(s) \right \} \\
&=
f_k(s)\{1+O(2^{-\sigma})\}
\end{align*}
for $s\in \mathscr{D}(2m)=\{ s\in \mathscr{D} \mid \sigma \geq 2m \}$ and $k\geq 1$.

Hence, by (\ref{f_kest})
\begin{equation}\label{Zkest}
Z_k(s)=\left(-\frac{\omega(s)}{2}\right)^k\{1+O_k((\log |s|)^{-2}) \} \quad(k\geq 1)
\end{equation}
for $s\in \mathscr{D}(2m)$. This leads to $\Re Z_k(s)>0$ in this region.
Hence, by the same argument as in Lemma \ref{f_zeros-outside}, we obtain

\begin{lemma}[Lemma 4.5 in \cite{M}]\label{zeros-outside}
For $s\in \{ s \mid \sigma<1-2m \} \cup \{ s \mid \sigma>2m \}$, the zeros of $Z_k(s)$ are all located in small circles centred negative even integers and positive odd integers,
and the number of those is $k$ in each circles.
\end{lemma}

Finally, we present a result of Matsumoto and Tanigawa on the zeros of $Z_k(s)$ in the ``critical strip''.
\begin{prop}[see the proof of Theorem 2 in \cite{M-T}]\label{zeros-Z_k}
If the Riemann Hypothesis is true, then the number of zeros of $Z_k(s)$ in $\{ s \mid 1-2m\leq \sigma \leq 2m, \sigma \neq1/2 \}$ is $O_k(1)$.
\end{prop}
They proved that the difference of the number of zeros of $Z^{(k)}(t)$ and that of $Z_k(s)$ is $O_k(1)$. From this fact, they obtained the following approximate formula;
\begin{prop}[Theorem 2 in \cite{M-T}]\label{number-zero}
Let $N_k(T)$ be the number of the zero of $Z^{(k)}(t)$ in the interval $(0,T)$. If the Riemann Hypothesis is true, then for any $k\geq 1$,
\begin{equation*}
N_k(T)=\frac{T}{2\pi}\log \frac{T}{2\pi}-\frac{T}{2\pi}+O_k(\log T).
\end{equation*}
\end{prop}

\section{An Alternative Proof of Matsuoka's result}
In this section, we give an alternative proof of Matsuoka's result;
\begin{theorem}[Theorem 1.1 in \cite{M}]\label{M-th}
If the Riemann Hypothesis is true, then for any non-negative integer $k$ there exists a $T=T(k)>0$ such that, for $t\geq T$, $Z^{(k+1)}(t)$ has exactly one zero between consecutive zeros of $Z^{(k)}(t)$.
\end{theorem}
More precisely, we prove Mozer's formula
\begin{equation}\label{Mozer}
\frac{d}{dt}\frac{Z^{(k+1)}}{Z^{(k)}}(t)=-\sum_{\gamma_k}\frac{1}{(t-\gamma_k)^2}+O_k(t^{-1})
\end{equation}
in the different way from Matsuoka's proof. Here the sum is taken over zeros of $Z^{(k)}(t)$.
If we can prove this formula, then we see that
\begin{align*}
\frac{d}{dt}\frac{Z^{(k+1)}}{Z^{(k)}}(t)
&<
-\sum_{0<\gamma_k<t}\frac{1}{(t-\gamma_k)^2}+At^{-1} \\
&<
t^{-1}(A-N_k(t)t^{-1}).
\end{align*}
By Proposition \ref{number-zero}, this is negative for large $t$.
Matsuoka's proof is inspired by Anderson's method (\cite{An}), which implies the above theorem for $k=1$. Matsuoka considered the integral
\begin{equation*}
\int_{C}\frac{G_k'}{G_k}(w)\frac{s}{w(s-w)}dw,
\end{equation*}
where $C$ is an appropriate rectangle and
\begin{equation*}
G_k(w)=h(w)\frac{Z_k(w)}{f_k(w)}
\end{equation*}
with $h(w)=\pi^{-s/2}\Gamma(s/2)$. However, in view of Theorem \ref{xi-facto}, we can prove the formula more easily.

\begin{proof}[Proof of (\ref{Mozer})]
By the definition of $\xi_k(s)$ and (\ref{Z^k}), we have
\begin{equation*}
\xi_k\left(\frac{1}{2}+it \right)=-i^{-k}\pi^{-\frac{1}{4}}\left(\frac{1}{4}+t^2 \right)\left|\Gamma \left(\frac{1}{4}+\frac{it}{2} \right) \right|^{-2k+1}Z^{(k)}(t).
\end{equation*}
Hence when we put
\begin{equation*}
g_k(t)=i^{-k}\pi^{-\frac{1}{4}}\left(\frac{1}{4}+t^2 \right)\left|\Gamma \left(\frac{1}{4}+\frac{it}{2} \right) \right|^{-2k+1}
\end{equation*}
then, by the logarithmic derivative with respect to $t$, we can obtain
\begin{equation}
i\frac{\xi_k'}{\xi_k}\left(\frac{1}{2}+it \right)=\frac{g_k'}{g_k}(t)+\frac{Z^{(k+1)}}{Z^{(k)}}(t).
\end{equation}
As for the function $g_k'/g_k(t)$, we can see that
\begin{equation*}
\frac{g_k'}{g_k}(t)
=
-(2k-1)\frac{d}{dt} \Re \log \Gamma \left(\frac{1}{4}+\frac{it}{2} \right)+\frac{2t}{\frac{1}{2}+t^2}
\end{equation*}
and hence
\begin{equation*}
\frac{d}{dt}\frac{g_k'}{g_k}(t)\ll_k t^{-1}
\end{equation*}
by Lemma \ref{estgam}.

On the other hand, Theorem \ref{xi-facto} implies that
\begin{equation*}
\frac{\xi_k'}{\xi_k}\left(\frac{1}{2}+it \right)=B_k+\sum_{\rho_k}\left(\frac{1}{\frac{1}{2}+it-\rho_k}+\frac{1}{\rho_k} \right).
\end{equation*}
Therefore we have
\begin{align*}
\frac{d}{dt}\frac{\xi_k'}{\xi_k}\left(\frac{1}{2}+it \right)
&=
\sum_{\rho_k}\frac{-i}{(\frac{1}{2}+it-\rho_k)^2} \\
&=
i\sum_{\gamma_k}\frac{1}{(t-\gamma_k)^2}+\sum_{\substack{\rho_k \\ \beta_k\neq \frac{1}{2}}}\frac{-i}{(\frac{1}{2}+it-\rho_k)^2} \\
&=
i\sum_{\gamma_k}\frac{1}{(t-\gamma_k)^2} \\
&\quad
+\sum_{\substack{\rho_k \\ \beta_k<1-2m, \ 2m<\beta_k}}\frac{-i}{(\frac{1}{2}+it-\rho_k)^2}+O(t^{-2}) \\
&=
i\sum_{\gamma_k}\frac{1}{(t-\gamma_k)^2}+O(t^{-1})
\end{align*}
by Proposition \ref{zeros-Z_k} and Lemma \ref{zeros-outside}. To show the last equality, we use the same argument as Matsuoka's (see p.15 in \cite{M});
\begin{equation*}
\sum_{\substack{\rho_k \\ \beta_k<1-2m, \ 2m<\beta_k}}\frac{1}{(\frac{1}{2}+it-\rho_k)^2}\ll_k \sum_{n=m}^{\infty}\frac{k}{(t+2n+1)^2}\ll_k\int_{2m+1}^{\infty}\frac{dx}{(t+x)^2}.
\end{equation*}
Thus we have
\begin{equation*}
\frac{d}{dt}\frac{Z^{(k+1)}}{Z^{(k)}}(t)=-\sum_{\gamma_k}\frac{1}{(t-\gamma_k)^2}+O_k(t^{-1}).
\end{equation*}
\end{proof}

\section*{Acknowledgements}
The author would like to thank my supervisor Professor Kohji Matsumoto for useful advice.
The author is grateful to the seminar members for some helpful remarks and discussions.


\section{Introduction}

General probability theory (GPT)  was conceived in part to have a common framework for both classical and quantum probability theory and in part to explore ``all possibilities that nature could have chosen'' instead of quantum
probability theory. Perhaps the root of this kind of intention was the study of possible correlations that can arise between two parts of a quantum system. Without any formal GPT structure, one may just consider the set of possible correlations subject to the {\it non-signalling} principle to observe
\cite{tsirelson,popescu_rohrlich} that the ones realizable by a quantum bipartite system form a strict subset. Consequently, it is natural to look for some further physical principle that somehow singles out quantum correlations \cite{infcaus} and to investigate how much better certain tasks could be performed using more general correlations (see e.g.\! \cite{zeroerror,sajat}).

Of course, when we consider a system {\it without} a division into parts, the structure becomes less rich. The concepts we can still talk about are the {\it state space} of the system and the set of possible {\it measurements}, which -- for our purposes -- are simply affine functions from the state space to a simplex (i.e.\! to the state space of a classical system). Since the first of these concepts completely determines\footnote{We assume the {\it no restriction principle} \cite{norestr}: all measurements that are theoretically possible (i.e.\! not ruled out by the convex structure of states) are feasible.} the latter one, ultimately every property of the undivided system is determined by the geometrical shape of the state space; e.g.\! its ``signalling dimension'' \cite{sigdim} or its ``information storability'', which in particular is known to depend on the amount of asymmetry \cite{geoinfcap} of the state space.

In both classical and quantum probability theory, if any two states of a certain collection can be perfectly distinguished, then in fact the whole collection is jointly perfectly distinguishable. However, in general -- when the state space $S$ of the system is not assumed to be neither a simplex (the classical case) or
the set of density operators on $\mathbb C^n$ (the quantum case), it can happen that there are three states $s_1,s_2,s_3\in S$ such that any two are perfectly distinguishable by {\it some} binary measurement, yet there is no measurement with three possible outcomes that would distinguish all of them with zero error.
In fact, in order for the states $s_1,s_2, \ldots s_{n+1} \in S$ to be jointly perfectly distinguishable (or henceforth in short: to be a ``j.p.d.\! collection''), at least they need to form a convex independent set, requiring $S$ to be at least $d=n$ dimensional. This is to be confronted with the fact that for every $d=1,2,\ldots$ one can give a $d$ dimensional state space $S$ and $n=2^d$ elements of $S$ such that any {\it two} can be perfectly distinguished; see \cite{rej}. So clearly, in general the condition of pairwise perfect distinguishability is much weaker than joint perfect distinguishability.

At this point one might wonder what can we say about the structure of perfect distinguishability in general. Let $S$ be a state space of a GPT (i.e.\! any convex set) and fix some elements $s_1,s_2,\ldots s_n\in S$ and consider the set of all subsets of indices
$$
\mathcal A \equiv \{H\subset [n]\, | \, (s_j)_{j\in H} \; \textrm{is a j.p.d. collection} \}
$$
where we have used the notation $[n]\equiv \{1,2,\ldots n\}$. It is clear that all one element subsets are in $\mathcal A$ and as smaller collections are easier to distinguish, if $H\in \mathcal A$ and $L \subset H$ then $L\in \mathcal A$; so $\mathcal A$ is a so-called {\it independence system} on $[n]$. But is there anything else that can be established at this level of generality? The answer is negative, as we shall seen this by the explicit construction of Section
\ref{sec:constr} which, given an independence system $\mathcal A$ on $[n]$ produces a convex set $S$ and elements $s_1,s_2, \ldots s_{n} \in S$ such that for any $H\subset [n]$ we have that $(s_j)_{j\in H}$  is a j.p.d.\! collection if and only if $H \in \mathcal A$.

This is somewhat similar to the situation of joint measurability, except that for an arbitrary structure of joint measurability it suffices to consider the quantum case (so in the case of joint measurability there is no need to invoke a generic GPT model). A collection of measurements is said to be jointly measurable, if they have a common refinement. Given any independence system $\mathcal A$ on $[n]$, one can produce some measurements $M_1,M_2,\ldots M_n$ on some quantum system such that for any $H\subset [n]$ we have that $(M_j)_{j\in H}$  is a jointly measurable collection if and only if $H \in \mathcal A$; see \cite{joint_measurability} for both the precise definition of joint measurability and the actual construction.

we shall conclude this introduction by pointing out that though perfect distinguishability is a nice mathematical concept, physically it is more relevant to discuss the {\it amount} by which some states can be distinguished. Given some states $s_1,s_2, \ldots s_{n} \in S$, regardless whether they can be jointly perfectly distinguished or not, we might try to find a measurement $M$ with possible outcomes numbered from $1$ to $n$ such that for each $k$, in case the system is in state $s_k$, the measurement $M$ ends up with a ``good enough'' probability with its $k^{\rm th}$ outcome, rightly indicating that the system was in state $s_k$. The usually considered ``figure of merit'', showing how well we can distinguish these states by performing measurement $M$ is the so-called {\it symmetric error probability}, which is the sum (or in some other convention: the average) of the error probabilities; that is, the quantity
$$
\sum_{k=1}^n P({\rm outcome}\neq k|M\, {\textrm{is measured when the system is in state}\; s_k} ).
$$
Following the notations of \cite{milanek}, let us denote by $P_e^*(s_1,s_2\ldots s_n)$ the minimum (or infimum, if there is no minimum) of
the above quantity when all measurements are considered. By its definition, $P_e^*$ is symmetrical in its arguments, takes values in the interval $[0,1]$ and equals to zero if its arguments form a j.p.d.\! collection. So for a fixed collection of states $s_1,s_2,\ldots s_n$
let us consider the function
$$
[n]\supset H\mapsto F(H):=P_e^*((s_j)_{j\in H}).
$$
Since less states are easier to distinguish, if $L\subset H$, then $F(L)\leq F(H)$. Moreover, given a collection and one more state, we can still perform the optimal measurement for the initial collection (although it will result in a sure error, if the state of system is the added new one), showing that $F(H\cup \{j\})\leq F(H)+1$. However, apart from these trivial properties, what else can be established about $F$, in general? To my knowledge, there is very little we know about this, though for the case when the state space is assumed to be that of a quantum system, \cite{milanek} has a bound on $P_e^*(s_1,s_2\ldots s_n)$ in terms of the pairwise error probabilities $P_e^*(s_j,s_k)$ ($j\neq k$). So both the quantum case and the general one should be topic of further researches.

\section{Preliminaries}
\subsection{The simplex and the affine hyperplane embedding it}

Let $\mathbb R^n_1\subset \mathbb R^n$ denote the affine hyperplane
\begin{equation}
\mathbb R^n_1\equiv \{(t_1,\ldots t_n)\in \mathbb R^n \, | \, t_1+\ldots +t_n=1\}.
\end{equation}
Let further $\Delta_n\subset \mathbb R^n_1$ be the $n$-simplex; i.e.\! the convex set of classical probability distributions with $n$ terms:
\begin{equation}
\Delta_n \equiv \{(p_1,\ldots p_n)\in \mathbb R^n_1 \, |\, \forall j: p_j\geq 0 \}.
\end{equation}
For $m<n$, one might view $\Delta_m$ as a subset of $\Delta_n$ (a probability distribution with $m$ terms can be also viewed as a probability distribution with $n$ terms by appending $n-m$ zero terms to it) and similarly, $\mathbb R^m_1$ as a subset of
$\mathbb R^n_1$. Then it is meaningful to consider the orthogonal projection $\pi_{n,m}:\mathbb R^n_1 \to \mathbb R^m_1$; we have
that
\begin{equation}
\label{defofpi}
\pi_{n,m}((t_1,t_2,\ldots, t_k,t_{k+1},\ldots t_n)) = (t_1+r,t_2+r,\ldots t_k+r)
\end{equation}
where $r=(t_{m+1}+\ldots +t_{n})/m$. Note that $\pi_{n,m}:\mathbb R^n_1 \to \mathbb R^m_1$ is an {\it affine} (i.e.\! convex combination preserving) map and that $\pi_{n,m}(\Delta_n)=\Delta_m$.

\subsection{GPT in a nutshell}

To describe a physical system, we want to give a mathematical meaning to the physical notions of {\it state} and {\it measurement} (with finite many possible outcomes). Correspondingly, we shall have a set $S$ which we will refer to as the {\it state space} of the system and whose elements will represent  the possible states of the system, and for each $k=1,2,\ldots$, a set $\mathcal M_k$ whose elements will represent the {\it measurements} with $k$ (numbered) outcomes that we can perform. Given a measurement $M\in \mathcal M_k$ and a state $s\in S$, the model should specify a classical probability distribution $p_{M,s}\in \Delta_k$ which is interpreted as the probability distribution of the outcomes for $M$ if the measurement is performed on the system when it is in state $s$.

From a physical point of view, two states should be regarded the same if for all possible measurements they lead to the same outcome probability distribution. Similarly, two measurement should be regarded the same if in every state they always produce the same outcome probability distribution. Thus, through the evaluation $p$, elements of $S$ and $\mathcal M=\cup_{k=1}^\infty \mathcal M_k$ should {\it separate} each other. In particular, we may identify a state $s$ with the function $p_{\cdot,s}:\mathcal M \to \cup_{k=1}^\infty  \Delta_k$ which in fact, for each $k=1,2,\ldots$, maps $\mathcal M_k$ into $\Delta_k$. As such, it is intrinsically meaningful to consider convex combinations of states; after all, convex combinations of the {\it functions} $p_{\cdot,s_1}, p_{\cdot,s_2}$ are well-defined. Importantly, such a combination is again regarded as a state of the system. Note that in this way, the outcome probability distribution for a fixed measurement, as a function of the state, is -- by definition -- affine. Note also that the convex combination of states has a clear operational meaning: to obtain outcome statistics of the measurement $M$ in the ``mixed state'' $\lambda s_1 +(1-\lambda)s_2$ where $s_1,s_2\in S$ and the coefficients $\lambda,1-\lambda \in [0,1]$, one prepares the system in either state $s_1$ or $s_2$ with corresponding probabilities $\lambda$ and $1-\lambda$, then performs $M$.

Thus by what was explained, we assume that the state space is endowed with a structure making it a {\it convex set} and every measurement is given as an affine map from $S$ to a certain simplex. In what follows, we shall assume the {\it no restriction principle} \cite{norestr}: every ``theoretically possible measurement'', that is, every affine map from $S$ to a simplex, is actually a realizable measurement. Hence the whole probability model is completely determined by a single input: by the choice of the convex set $S$ playing the role of the state space of the physical system.

This principle has some physical motivations and it certainly holds in both the classical and quantum case. In the (finite) classical case $S=\Delta_n$ for some $n$, while in the (finite) quantum case $S$ is the set of density operators on $\mathbb C^n$ for some $n$. As is well-known, $\Phi$ is an affine function from the set of density operators on $\mathbb C^n$ to $\Delta_k$ if and only if there exist $k$ positive operators $E_1,E_2,\ldots E_k\geq 0$ on $\mathbb C^n$ whose sum is the identity, and $\Phi$ is the map
$$
\rho \mapsto ({\rm Tr}(\rho E_1), {\rm Tr}(\rho E_2),\ldots {\rm Tr}(\rho E_k)),
$$
giving back the usual picture where measurements on a quantum system are represented as {\it positive operator valued measures}.

\subsection{Joint perfect distinguishability}
Let $S$ be the state space of a GPT model. A collection of states $s_1,\ldots s_n\in S$ is said to be {\it jointly perfectly distinguishable} or in short: j.p.d., if there exists a measurement with outcomes $1,\ldots n$ such that for all $k=1,\ldots n$, in case the system is in state $s_k$, the measurement ends up with its $k^{\rm th}$ outcome with probability one (i.e.\! rightly indicates that out of the listed states, the system must had been in state $s_k$). Thus, by what was explained in words, $s_1,\ldots s_n\in S$ are j.p.d.\! if and only if there exists an affine map $\Phi:S\to \Delta_n$ such that for all $k=1,\ldots n:$ $(\Phi(s_k))_k=1$. In the quantum case, as is well-known, this requirement is equivalent to asking that the $n$ density operators have (pairwise) diagonal images.

Let $k\leq n$ and $s_1,\ldots s_n\in S$ a j.p.d.\! collection. This means that there is an affine map $\Phi:S\to \Delta_n$ such that $\Phi(j)$ is precisely the $j^{\rm th}$ vertex of $\Delta_n$ for each $j=1,\ldots n$. Then $\pi_{n,k}\circ\Phi$ is an affine map from $S$ to $\Delta_k$ still having the property that $\Phi(j)$ is precisely the $j^{\rm th}$ vertex of $\Delta_k$ for each $j=1,\ldots k$. Hence the smaller collection of states $s_1,\ldots s_k$ is j.p.d., as we naturally expected.

Finally note that the j.p.d.\! requirement for two states $s_1,s_2\in S$ is geometrically equivalent to asking $s_1,s_2$ to be a so-called {\it antipodal} pair of points of the convex set $S$, see \cite{rej}. Thus, the j.p.d.\! property may be viewed as a natural generalization of the mentioned concept giving rise to the notion of ``jointly antipodal collection'' of points.

\section{The construction}
\label{sec:constr}

Let $n>1$ be an integer and $\mathcal A$ be an independence system on $[n]$. We shall say that a subset
$H\subset [n]$ is {\it minimally dependent} (with respect to $\mathcal A$), if $H\notin \mathcal A$, but dropping any element of $H$ makes it belong to $\mathcal A$. Note that any minimally dependent set must have at least $2$ elements. Let $\mathcal O_\mathcal A$ be the set of all minimally dependent (with respect to $\mathcal A$) subsets of $[n]$.

Our construction will start with the vertices of the $n$-simplex; so we shall set
\begin{equation}
s_1=(1,0,0,\ldots), \; s_2=(0,1,0,\ldots), \, \ldots
\end{equation}
and so on.
The problem is that if the state space $S$ was simply the convex hull of these points (i.e.\! the simplex $\Delta_n$),
then $s_1,\ldots s_n$ would be a j.p.d.\! collection. In order to take account of the given independence system $\mathcal A$, for each
$H\in \mathcal O_{\mathcal A}$ we shall add a number of new points to ``ruin'' the j.p.d.\! property for $(s_j)_{j\in H}$ and hence also for all larger collections, while leaving it ``intact'' for any $(s_j)_{j\in K}$ such that $K$ does not contain $H$.

For a subset $H\subset [n]$ we shall denote by $1_H$ the {\it indicator} of $H$; i.e.\! $1_H$ is the point of $\mathbb R^n$ whose $j^{\rm th}$ coordinate is $1$ if $j\in H$ and zero otherwise. Note that in particular $1_{\{j\}}=s_j$.

Now suppose $H\subset [n]$ has at least $m:=|H|\geq 2$ elements.
If $H$ is a proper subset (i.e.\! $m<n$) then for every $j\in H$ we shall set
\begin{equation}
\label{defofq}
q_H^{(j)}:= -2\epsilon\, 1_{\{j\}} +\left(\frac{1-\epsilon(m-2)}{m-1} \right)1_{H\setminus\{j\}}+
\left(\frac{\epsilon m}{n-m}\right)1_{[n]\setminus H}
\end{equation}
where the constant $\epsilon=1/(3n^2)$. For the full set $[n]$ instead, we shall set
\begin{equation}
q_{[n]}^{(j)}:= -\epsilon 1_{\{j\}} + \left(\frac{1+\epsilon}{n-1}\right) 1_{[n]\setminus \{j\}}.
\end{equation}
In both cases, one can easily check that in fact, the thus defined point $q_H^{(j)}\in \mathbb R^n$  belongs to $\mathbb R^n_1$; i.e.\!
the sum of its coordinates is $1$. However, because its $j^{\rm th}$ coordinate is negative, $q_H^{(j)}\notin \Delta_n$ (but note that all of is other coordinates are positive). Note also that in both cases $-2\epsilon$ is a lower bound on all coordinates and as by a straightforward check
\begin{equation}
\min\{\left(\frac{1-\epsilon(m-2)}{m-1} \right), \left(\frac{1+\epsilon}{n-1}\right)\} \geq \frac{2}{3n}=2n\epsilon,
\end{equation}
that $2n\epsilon$ is a lower bound on the value of the coordinates whose index belongs to $H\setminus j$. Finally, note that for $1<m<n$, we have that $\frac{1-\epsilon(m-2)}{m-1}> \frac{1}{3n}=\epsilon n > \frac{\epsilon m}{n-m}$, so in case $1<|H|<n$, we have that the coordinate value of $q_H^{(j)}$ associated to an index in $H\setminus \{j\}$ is always greater than that one associated to an index in $[n]\setminus H$.

We shall now introduce our state space of choice. Let $S$ be the convex hull of the set points formed by the vertices $s_j$ $(j\in [n])$ of the simplex $\Delta_n$ together with the newly introduced points $q^{(j)}_H$ for all minimally dependent set $H\in \mathcal O_{\mathcal A}$ and $j\in H$:
\begin{equation}
S:= \textrm{Conv}
    \big\{ \,
    \{s_j|j\in [n]\}
    \cup
    \{q_H^{(j)} | \,
	H\in \mathcal O_{\mathcal A},\, j\in H \}
	\,
	\big\}.
\end{equation}

\begin{Pro} For any $H\in \mathcal O_{\mathcal A}$, we have that $(s_j)_{j\in H}$ is {\bf not} a j.p.d. collection of the state space $S$.
\end{Pro}
\begin{proof}
Without loss of generality, we may assume that $H=[m]$ is precisely the first $m$ elements of $[n]$, where $m=|H|\geq 2$. We have to show that there exists no affine function $\Phi:S\to \Delta_m$ that would map each of the points $s_j$ $(j=1,\ldots m)$ into a separate vertex of the simplex $\Delta_m$ (in fact, without loss of generality: maps $s_j$ into the $j^{\rm th}$ vertex of $\Delta_m$, which -- with an abuse of notations -- we shall still denote by $s_j$).

Since $\Delta_m\subset S$, an affine function $\Phi:S\to\Delta_m$ extends in a unique manner to be a $\hat{\Phi}:\mathbb R^n\to\mathbb R^m$ linear map which preserves the sum of the coordinates (as it must map $\mathbb R^n_1$ into $\mathbb R^m_1$).

So we shall prove our statement by showing that if $\hat{\Phi}:\mathbb R^n\to\mathbb R^m$ is a linear map such that it preserves the sum of coordinates and $\hat{\phi}(s_j)=s_j$ for all $j\in [m]$, then there exists a $j\in [m]$ such that $\hat{\Phi}(q_{[m]}^{(j)})\notin \Delta_m$ and hence that $\hat{\Phi(S)}\not\subset \Delta_m$. We shall treat separately the cases when $m=n$ and when $m<n$.

\smallskip

{\it Case 1: $m=n$}.  The points $s_j$ ($j\in [m]=[n]$) form a basis of $\mathbb R^n$, so the map $\hat{\Phi}$ must be the identity and hence the proof is finished since, as was noted, $q_{[n]}^{(j)}\notin \Delta_n$.

\smallskip

{\it Case 2: $2\leq m<n$}. The sum of the coordinates of $\xi:= \hat{\Phi}(1_{[n]\setminus H})$ is equal to that of $1_{[n]\setminus H}$: it is $n-m$. Hence at least one coordinate of $\xi\in \mathbb R^m$ -- say the $j^{\rm th}$ one must be smaller than or equal to $\frac{n-m}{m}$:
\begin{equation}
c:=(\xi)_j\leq \frac{n-m}{m}.
\end{equation}
Since $1_{\{j\}}=s_j$ and
\begin{equation}
1_{H\setminus \{j\}}=1_{[m]\setminus \{j\}}=\sum_{k\leq m, k\neq j }s_k,
\end{equation}
taking account of the fact that $\hat{\phi}(s_j)=s_j$ for all $j\in [m]$, a straightforward
substitution of the point $q_H^{(j)}$ given by formula (\ref{defofq}) into $\hat\Phi$ gives
\begin{equation}
\hat{\Phi}(q_{[m]}^{(j)})= -2\epsilon\, s_j +\left(\frac{1-\epsilon(m-2)}{m-1} \right)\sum_{k\leq m, k\neq j} s_k +
\left(\frac{\epsilon m}{n-m}\right)\xi,
\end{equation}
whose $j^{\rm th}$ coordinate is
\begin{equation}
(\hat{\Phi}(q_{[m]}^{(j)}))_j=-2\epsilon + \left(\frac{\epsilon m}{n-m}\right) c \leq -2\epsilon + \epsilon < 0,
\end{equation}
showing that $\hat{\Phi}(q_{[m]}^{(j)})\notin \Delta_m$.
\end{proof}

\begin{Pro} For any $H\in \mathcal A$, we have that $(s_j)_{j\in H}$ is a j.p.d.\! collection of the state space $S$.
\end{Pro}

\begin{proof}
As in the previous proof, we may assume that $H=[m]$ is precisely the first $m$ elements of $[n]$
The statement is trivial if $m<2$ or $m=n$ (in this latter case, $\mathcal O_{\mathcal A}$ is empty and $S=\Delta_n$).
So suppose that $2\leq m<n$. We know that $\pi_{n,m}$ is affine and $\pi_{n,m}(s_j)=s_j$ for all $j\in [m]$; so if
$\pi_{n,m}(S)\subset \Delta_m$ also holds, then this shows that $(s_j)_{j\in [m]}$ is a j.p.d.\! collection of the
state space $S$.

In what follows we shall show that this is indeed the case. We will need to check that the coordinates of
$\pi_{n,m}(q_K^{(j)})$ are all non-negative whenever $K\in \mathcal O_{\mathcal A}$ and $j\in K$. Since
$\pi_{n,m}(q_K^{(j)}) \in \mathbb R^m_1$, this non-negativity will then show that $\pi_{n,m}(q_K^{(j)}) \in \Delta_m$
whenever $K\in \mathcal O_{\mathcal A}$ and $j\in K$, which in turn will imply that $\pi_{n,m}(S)\subset \Delta_m$.

What remains is a computation of coordinate values of $\pi_{n,m}(q_K^{(j)})$. We will need a case-by-case
treatment depending on the possible position of $K$ and $j$ relative to the subset $H=[m]$. In what follows,
we shall use the notation $\ell:=|K\setminus [m]|$. Note that $\ell$ must be at least one:
since $[m]=H\in \mathcal A$ and $K$ is a minimal dependent set, $K$ cannot be a subset of $[m]$.

\smallskip

{\it Case 1: $j\in [m]$}. With the exception of its $j^{\rm th}$ coordinate, all coordinates of $q_K^{(j)}$ are positive. Taking account of the formula (\ref{defofpi}), the definition $q_K^{(j)}$ and the estimates we had on its coordinates and the fact that at least one element of $K\setminus\{j\}$ is outside of $[m]$, it follows that the smallest coordinate of $\pi_{n,m}(q_K^{(j)})$ is the $j^{\rm th}$ one, whose actual value is
\begin{equation}
(q_K^{(j)})_j + \frac{1}{m}\sum_{b=m+1}^n(q_K^{(j)})_b \geq -2\epsilon + \frac{1}{n} 2n\epsilon = 0.
\end{equation}

{\it Case 2.a: $j \notin [m]$ and $\ell>1$}.
Apart from $j$, there must be at least one more element of $K$ outside of $[m]$. Hence the sum of coordinates of $q_K^{(j)}$ corresponding to indices in $K\setminus [m]$ is at least $-2\epsilon + 2n \epsilon\geq 0$. Taking account of the fact that coordinates of $q_K^{(j)}$ corresponding to indices in $[m]$ are positive, this implies that all coordinates of $\pi_{n,m}(q_K^{(j)})$ are positive, too.
\smallskip

{\it Case 2.b1: $j \notin [m], \ell=1$ and $[m]\subset K$}.
It is straightforward to check that in this case all coordinates of $\pi_{n,m}(q_K^{(j)})$ are equal and hence they are all of magnitude
$\frac{1}{m}$, which is positive.
\smallskip

{\it Case 2.b2: $j \notin [m], \ell=1$ and $[m]\not\subset K$}. Let $k:=|K|$; we have that $|K\cup [m]|=k-\ell=k-1<m$,
since $[m]\not\subset K$. The value of every coordinate of $q_K^{(j)}$ associated to an index in $[n]\setminus ([m]\cup \{j\})$
is $\frac{\epsilon k}{n-k}$, so the sum of the coordinates of $q_K^{(j)}$ associated to indices in $[n]\setminus [m]$ is
\begin{equation}
-2 \epsilon + (n-m-1) \frac{\epsilon k}{n-k}.
\end{equation}
On the other hand, the minimal value of the coordinates of $q_K^{(j)}$ associated to indices in $[m]$ is also $\frac{\epsilon k}{n-k}$.
Thus the minimal coordinate value of $\pi_{n,m}(q_K^{(j)})$ is
\begin{equation}
\frac{\epsilon k}{n-k} + \frac{1}{m}\left(-2 \epsilon + (n-m-1) \frac{\epsilon k}{n-k}\right).
\end{equation}
After some rearrangement, multiplying the above expression by $\frac{m(n-k)}{\epsilon}$ makes $\epsilon$ and $m$
disappear from it and yields
\begin{equation}
k(n-1)-2(n-k)
\end{equation}
which is evidently positive, since $n>1$ and $k\geq 2$.
\smallskip

We exhausted all possibilities, so the proof is complete.
\end{proof}

\begin{Cor}
For any $H\subset [n]$, we have that $(s_j)_{j\in H}$ is a j.p.d.\! collection of the state space $S$
if and only if $H\in \mathcal A$.
\end{Cor}