---
abstract: 'The Rényi and von Neumann entropies of the thermal state in the generalized uncertainty principle (GUP)-corrected single harmonic oscillator system are explicitly computed within the first order of the GUP parameter $\alpha$. While the von Neumann entropy with $\alpha = 0$ exhibits a monotonically increasing behavior in external temperature, the nonzero GUP parameter makes the decreasing behavior of the von Neumann entropy at the large temperature region. As a result, the von Neumann entropy is maximized at the finite temperature if $\alpha \neq 0$. The Rényi entropy $S_{\gamma}$ with nonzero $\alpha$ also exhibits similar behavior at the large temperature region. In this region the Rényi entropy exhibit decreasing behavior with increasing the temperature. The decreasing rate becomes larger when the order of the Rényi entropy $\gamma$ is smaller.'
author:
- 'MuSeong Kim$^1$, Mi-Ra Hwang$^1$, Eylee Jung$^1$, and DaeKil Park$^{1,2}$[^1]'
title: ' Rényi and von Neumann entropies of thermal state in Generalized Uncertainty Principle-corrected harmonic oscillator'
---

Introduction
============

The existence of a minimal length (ML) at Planck’s scale seems to be a universal characteristic of quantum gravity[@townsend76; @amati89; @garay94]. It appears in loop quantum gravity[@rovelli88; @rovelli90; @rovelli98; @carlip01], string theory[@konishi90; @kato90; @strominger91], path-integral quantum gravity[@padmanabhan85; @padmanabhan85-2; @padmanabhan86; @padmanabhan87; @greensite91], and black hole physics[@maggiore93]. ML also appeared in some microscope thought-experiment[@mead64]. From an aspect of quantum mechanics the existence of the ML modifies the uncertainty principle from Heisenberg uncertainty principle (HUP)[@uncertainty; @robertson1929] $\Delta P \Delta Q \geq \frac{\hbar}{2}$ to generalized uncertainty principle (GUP)[@kempf93; @kempf94]. This is because of the fact that the uncertainty of the position $\Delta Q$ should be larger than the ML. In order to examine the effect of the modified quantum mechanics the ML and GUP are recently used to explore several branches of physics such as quantum mechanics[@chung19], quantum electrodynamics[@bosso18], quantum cosmology[@bosso19], quantum gravity[@modak19], and black hole physics[@xiang18].

More recently, the Feynman propagators[@feynman; @kleinert] for the GUP-corrcted quantum mechanics were derived in $1d$ free particle system [@das2012; @gangop2019], $1d$ harmonic oscillator[@comment-1], and $d$-dimensional harmonic oscillator[@park-20]. Moreover, the singular potential problem in the GUP-corrected quantum mechanics was discussed in Ref.[@point-1].

The purpose of this paper is to discuss on the Rényi and von Neumann entropies of $1d$ harmonic oscillator thermal state in the presence of the ML. The specific form of GUP we will choose is a following form: $$\label{GUP-d-1}
\Delta P_i \Delta Q_i \geq \frac{\hbar}{2} \left[ 1 + \alpha \left(\Delta {\bf P}^2 + \langle {\bf P} \rangle^2 \right) + 2 \alpha \left( \Delta P_i^2 + \langle P_i \rangle^2 \right) \right]
\hspace{1.0cm}  (i = 1, 2, \cdots, D)$$ where $\alpha$ is a GUP parameter, which has a dimension $(\mbox{momentum})^{-2}$. Using $\Delta A \Delta B \geq \frac{1}{2} | \langle [A, B] \rangle |$, Eq. (\[GUP-d-1\]) induces the modification of the commutation relation as $$\begin{aligned}
\label{GUP-d-2}
&& \left[ Q_i, P_j \right] = i \hbar \left( \delta_{ij} + \alpha \delta_{ij} {\bf P}^2 + 2 \alpha P_i P_j  \right)    \\    \nonumber
&& \hspace{1.0cm} \left[ Q_i, Q_j \right] = \left[P_i, P_j \right] = 0.\end{aligned}$$

The existence of the minimal length is easily shown at $d = 1$. In this case Eq. (\[GUP-d-1\]) is expressed as $$\label{GUP-d-3}
\Delta P \Delta Q \geq \frac{\hbar}{2} \left( 1 + 3 \alpha \Delta P^2 \right)$$ if $\langle P \rangle = 0$. Then, the equality of Eq. (\[GUP-d-3\]) yields $$\label{minimal-length}
\Delta Q^2 \geq \Delta Q_{min}^2 = 3 \alpha \hbar^2.$$ In Fig. 1 the allowed region and minimal length of Eq. (\[GUP-d-3\]) is plotted when $\hbar = \alpha = 1$.

![(Color online) The minimal length and allowed region of one-dimensional GUP (\[GUP-d-3\]) when $\hbar = \alpha = 1$. This figure is copied from Ref. [@park-20]. ](fig1.pdf){height="6.0cm"}

If $\alpha$ is small, Eq. (\[GUP-d-2\]) can be solved as $$\label{GUP-d-4}
P_i = p_i \left(1 + \alpha {\bf p}^2 \right) + {\cal O} (\alpha^2)   \hspace{1.0cm} Q_i = q_i$$ where $p_i$ and $q_i$ obey the usual Heisenberg algebra $[q_i, p_j] = i \hbar \delta_{ij}$. Using Eq. (\[GUP-d-4\]) we will compute the Rényi and von Neumann entropies in this paper within the first order of $\alpha$ in the following.

The paper is organized as follows. In next section we will derive the thermal state $\rho_T$ of the GUP-corrected harmonic oscillator system. Since the thermal state is not Gaussian, it is highly difficult to solve the eigenvalue equation of $\rho_T$. However, it is shown in the following sections that the eigenvalue can be derived up to the first order of $\alpha$ without solving the eigenvalue equation explicitly. In section III we compute the ${\mbox Tr} \rho_T^n$ explicitly within ${\cal O} (\alpha)$. The purity function is also explicitly derived. It is shown in this section that while the thermal state with $\alpha = 0$ transforms from pure to completely mixed state with increasing the external temperature from $0$ to $\infty$, the state with nonzero $\alpha$ becomes the partially mixed state even at the limit of $T \rightarrow \infty$. By making use of ${\mbox Tr} \rho_T^n$ the Rényi and von Neumann entropies of $\rho_T$ are computed within ${\cal O} (\alpha)$ in section IV. The temperature-dependence of the entropies is analyzed in this section. In section V a brief conclusion is given. In the appendix $A$, $B$, and $C$ we derive some integral formulas, which are used in section II.

Thermal state
=============

Now, let us consider a simple harmonic oscillator system, whose Hamiltonian is $$\label{hamil-single}
\hat{H} = \frac{1}{2 m}P^2 + \frac{1}{2} m \omega^2 Q^2 = \frac{p^2}{2 m} + \frac{\alpha}{m} p^4 + \frac{1}{2} m \omega^2 q^2 + {\cal O} (\alpha^2).$$ The Feynman propagator for this system was explicitly was explicitly derived up to the first order of $\alpha$ in Ref. [@comment-1], whose explicit expression is $$\label{GUP-kernel}
K[q_f, q_0: T] = \sqrt{\frac{m \omega}{2 \pi i \hbar \sin \omega T}} \left[ 1 + \alpha f(q_0, q_f: T) + {\cal O} (\alpha^2) \right] e^{\frac{i}{\hbar} (S_0 + \alpha S_1)},$$ where $$\begin{aligned}
\label{k-schrodinger-14}
&&S_0 = \frac{m \omega}{2 \sin \omega T} \left[ (q_0^2 + q_f^2) \cos \omega T - 2 q_0 q_f \right]                                                         \\     \nonumber
&&S_1 = - \frac{m^3 \omega^3}{32 \sin^4 \omega T}  \Bigg[ \left\{ 12 \omega T + 8 \sin 2 \omega T + \sin 4 \omega T \right\} (q_0^4 + q_f^4)               \\     \nonumber
&& \hspace{1.5cm} -4 \left\{ 12 \omega T \cos \omega T + 11 \sin \omega T + 3 \sin 3 \omega T \right\} q_0 q_f (q_0^2 + q_f^2)                                           \\     \nonumber
&& \hspace{4.5cm}   + 12 \left\{ 4 \omega T + 2 \omega T \cos 2\omega T + 5 \sin 2 \omega T \right\} q_0^2 q_f^2    \Bigg]                                               \\     \nonumber
&& f(q_0, q_f: T) = \frac{3 i \hbar m \omega}{8 \sin^2 \omega T}  \left( 2 \omega T + 5 \sin \omega T \cos \omega T + \omega T \cos 2 \omega T \right)                   \\     \nonumber
&& \hspace{2.8cm} - \frac{3 m^2 \omega^2}{8 \sin^3 \omega T} \Bigg[ 2 \omega T \left\{ 3 \cos \omega T (q_0^2 + q_f^2) -  2 (2 + \cos 2 \omega T )q_0 q_f  \right\}                     \\     \nonumber
&& \hspace{5.0cm} + 10 \sin \omega T (q_0^2 + q_f^2 - 2 q_0 q_f \cos \omega T) - 6 \sin^3 \omega T (q_0^2 + q_f^2)       \Bigg].\end{aligned}$$

The Brownian or Euclidean propagator is defined as $G[q_f, q_0: \tau] = K[q_f, q_0: T= - i \tau]$. Then, the thermal state of the system is given by $$\label{thermal-1}
\rho_T [q_f, q_0: \beta] = \frac{1}{{\cal Z}} G[q_f, q_0: \hbar \beta]$$ where $\beta = 1 / (k_B T)$. The parameters $k_B$ and $T$ are Boltzmann constant and external temperature. In Eq. (\[thermal-1\]) ${\cal Z}$ is a partition function defined as $$\label{partition-1}
{\cal Z} (\beta) \equiv \mbox{Tr} G = \int d q G[q, q: \beta].$$ Then, the explicit expression of the thermal state is $$\begin{aligned}
\label{thermal-2}
&&\rho_T[q_f, q_0: \beta]                                 \\    \nonumber
&&= \frac{1}{{\cal Z} (\beta)} \sqrt{\frac{m \omega}{2 \pi \hbar \sinh x}} \left[ 1 - 
\alpha \left(f_E (q_0, q_f: \beta)  + \frac{1}{\hbar} S_{1,E} (q_0, q_f: \beta) \right) + {\cal O} (\alpha^2) \right]   \\   \nonumber
&&\hspace{4.0cm} \times\exp \left[-\frac{m \omega}{2 \hbar \sinh x} \left\{ (q_0^2 + q_f^2) \cosh x - 2 q_0 q_f \right\} \right]\end{aligned}$$ where $x$ is a dimensionless parameter $x = \hbar \omega \beta$ and $$\begin{aligned}
\label{thermal-2-boso}
&& {\cal Z} (\beta) = \frac{1}{2 \sinh \frac{x}{2}} \left[ 1 - \frac{3 \alpha}{4} \left(\hbar m \omega x \right) \coth^2 \frac{x}{2} + {\cal O} (\alpha^2) \right]    \\    \nonumber
&&f_E (q_0, q_f: \beta) = A_1 - A_2 (q_0^2 + q_f^2) + 2 A_3 q_0 q_f                           \\   \nonumber
&&\frac{1}{\hbar} S_{1,E} (q_0, q_f: \beta) = B_1 (q_0^4 + q_f^4) + B_2 (q_0^3 q_f + q_0 q_f^3) + B_3 q_0^2 q_f^2\end{aligned}$$ with $$\begin{aligned}
\label{thermal-2-boso-2}
&&A_1 = \frac{3 \hbar m \omega}{8 \sinh^2 x} \left[ 2 x + 5 \sinh x \cosh x + x \cosh 2 x \right]             \\    \nonumber
&&A_2 = \frac{3 m^2 \omega^2}{8 \sinh^3 x} \left[ 6 x \cosh x + 10 \sinh x + 6 \sinh^3 x \right]                         \\    \nonumber
&&A_3 = \frac{3 m^2 \omega^2}{8 \sinh^3 x} \left[2 x (2 + \cosh 2 x) + 10 \sinh x \cosh x \right]                        \\    \nonumber
&&B_1 = \frac{m^3 \omega^3}{32 \hbar \sinh^4 x} \left[ 12 x + 8 \sinh 2 x + \sinh 4 x \right]                                       \\    \nonumber
&&B_2 = -\frac{m^3 \omega^3}{8 \hbar \sinh^4 x} \left[ 12 x \cosh x + 11 \sinh x + 3 \sinh 3 x  \right]                  \\    \nonumber
&&B_3 = \frac{3 m^3 \omega^3}{8 \hbar \sinh^4 x} \left[4 x + 2 x \cosh 2 x + 5 \sinh 2 x \right].\end{aligned}$$

In order to derive the Rényi and von Neumann entropies of $\rho_T[q_f, q_0: \beta]$ we need to solve the eigenvalue equation $$\label{eigen-1}
\int d q_0 \rho_T[q_f, q_0: \beta] f_n (q_0: \beta) = \lambda_n f_n (q_f: \beta).$$ Since $\mbox{Tr} \rho_T = 1$, the eigenvalue should obey $\sum_n \lambda_n = 1$. However, it is highly difficult to solve Eq. (\[eigen-1\]), because the thermal state is not a Gaussian. As we will show in the following, however, it is possible to derive the eigenvalue $\lambda_n$ exactly within the order of $\alpha$ without solving the eigenvalue equation (\[eigen-1\]) explicitly.

calculation of $\mbox{Tr} \rho_T^n$
===================================

In this section we will compute $\mbox{Tr} \rho_T^n$ up to order of $\alpha$. One can show $\mbox{Tr} \rho_T^n$ reduces to $$\begin{aligned}
\label{trace-1}
&&\mbox{Tr} \rho_T^n                                                                                                                                                            
\equiv \int dx_1 \cdots dx_n \rho_T[x_1, x_2: \beta] \rho_T[x_2, x_3: \beta] \cdots \rho_T[x_{n-1}, x_n: \beta] \rho_T[x_n, x_1: \beta]                                         \\    \nonumber
&&=\frac{1}{{\cal Z}^n} \left(\frac{m \omega}{2 \pi \hbar \sinh x} \right)^{n/2}                                                                                     \\    \nonumber
&& \times \int dx_1 \cdots dx_n 
\Bigg[ 1 - \alpha \bigg\{f_E (x_1, x_2: \beta) + \cdots +f_E(x_{n-1}, x_n: \beta) + f_E(x_n, x_1: \beta) \bigg\}                                                                  \\    \nonumber
&&\hspace{2.0cm} -\frac{\alpha}{\hbar} \bigg\{S_{1, E}(x_1, x_2: \beta) + \cdots + S_{1, E}(x_{n-1}, x_n: \beta) + S_{1, E} (x_n, x_1: \beta) \bigg\} + {\cal O} (\alpha^2)     \Bigg]         \\    \nonumber
&&\hspace{8.0cm}  \times \exp \left[- {\bm X} G_n {\bm X}^{\dagger} \right],\end{aligned}$$ where ${\bm X} = (x_1, \cdots, x_n)$ and $G_n$ is a $n \times n$ matrix in a form $$\begin{aligned}
\label{matrix-Gn}
G_n = \left(                  \begin{array}{cccccc}
                  2 a & \hspace{0.2cm}    -b    & \hspace{0.2cm}         & \hspace{0.2cm}         & \hspace{0.2cm}         & \hspace{0.2cm}  -b            \\
                  -b  & \hspace{0.2cm} 2 a      & \hspace{0.2cm} \bullet & \hspace{0.2cm}         & \hspace{0.2cm}         & \hspace{0.2cm}                \\
                      & \hspace{0.2cm} \bullet  & \hspace{0.2cm} \bullet & \hspace{0.2cm} \bullet & \hspace{0.2cm}         & \hspace{0.2cm}                \\
                      & \hspace{0.2cm}          & \hspace{0.2cm} \bullet & \hspace{0.2cm} \bullet & \hspace{0.2cm} \bullet & \hspace{0.2cm}                \\
                      & \hspace{0.2cm}          & \hspace{0.2cm}         & \hspace{0.2cm} \bullet & \hspace{0.2cm} \bullet & \hspace{0.2cm}   -b           \\
                   -b   & \hspace{0.2cm}          & \hspace{0.2cm}         & \hspace{0.2cm}         & \hspace{0.2cm}     -b  & \hspace{0.2cm}  2a 
                              \end{array}                        \right).\end{aligned}$$ with $$\label{Gn-boso-1}
a =\frac{m \omega \cosh x}{2 \hbar \sinh x}    \hspace{1.0cm}  b = \frac{m \omega}{2 \hbar \sinh x}.$$ In Eq. (\[matrix-Gn\]) matrix elements in empty space are zero. Inserting Eq. (\[thermal-2-boso\]) into Eq. (\[trace-1\]) $\mbox{Tr} \rho_T^n$ becomes $$\begin{aligned}
\label{trace-2}
&&\mbox{Tr} \rho_T^n                                                                                                                                                            \\    \nonumber                                                                                                                                                          
&&=\frac{1}{{\cal Z}^n} \left(\frac{m \omega}{2 \pi \hbar \sinh x} \right)^{n/2}                                                                                     \\    \nonumber
&& \times \int dx_1 \cdots dx_n 
\Bigg[ 1 - \alpha \left\{n A_1 - 2 A_2 (x_1^2 + \cdots + x_n^2 ) + 2 A_3 (x_1 x_2 + \cdots + x_{n-1} x_n + x_n x_1) \right\}                                                                  \\    \nonumber
&&\hspace{3.0cm} -\alpha \bigg[2 B_1 (x_1^4 + \cdots + x_n^4)                                                                                                                               \\   \nonumber                                                                                                                            
&&\hspace{4.0cm}+ B_2 \left\{ x_1 x_2 (x_1^2 + x_2^2) + \cdots + x_{n-1} x_n (x_{n-1}^2 + x_n^2) + x_n x_1 (x_n^2 + x_1^2) \right\}                                                                       \\    \nonumber                                                                            
&& \hspace{4.0cm}                            + B_3 (x_1^2 x_2^2 + \cdots + x_{n-1}^2 x_n^2 + x_n^2 x_1^2 )  \bigg] + {\cal O} (\alpha^2)     \Bigg]                                                        \\    \nonumber       
&&\hspace{8.0cm}  \times \exp \left[- {\bm X} G_n {\bm X}^{\dagger} \right],\end{aligned}$$ It is easy to show $\mbox{det} G_2 = 4 (a^2 - b^2)$, $\mbox{det} G_3 = 2 (a - b) (2 a + b)^2$, and $\mbox{det} G_4 = 16 a^2 (a^2 - b^2)$. Generalizing to $n$ one can show $$\label{det-Gn}
\mbox{det} G_n = \frac{1}{2^n} \left[ \left( \sqrt{a + b} + \sqrt{a - b} \right)^n -  \left( \sqrt{a + b} - \sqrt{a - b} \right)^n \right]^2.$$ Thus, it is possible to derive $$\begin{aligned}
\label{int-formula-1}
&&\int dx_1 \cdots dx_n  \exp \left[- {\bm X} G_n {\bm X}^{\dagger} \right] = \frac{\pi^{n/2}}{\sqrt{\mbox{det} G_n}} \equiv g_n                                 \\     \nonumber
&&\int dx_1 \cdots dx_n (x_1^2 + \cdots + x_n^2) \exp \left[- {\bm X} G_n {\bm X}^{\dagger} \right]                                                              \\     \nonumber
&&= - \frac{1}{2} \frac{\partial g_n}{\partial a} = g_n \frac{n}{4 \sqrt{a^2 - b^2}} \frac{(\sqrt{a + b} + \sqrt{a - b})^n + (\sqrt{a + b} - \sqrt{a - b})^n}
                                                                                          {(\sqrt{a + b} + \sqrt{a - b})^n - (\sqrt{a + b} - \sqrt{a - b})^n}    \\     \nonumber
&&\int dx_1 \cdots dx_n (x_1 x_2 + \cdots + x_{n-1} x_n + x_n x_1)  \exp \left[- {\bm X} G_n {\bm X}^{\dagger} \right]                                           \\     \nonumber
&&= \frac{1}{2} \frac{\partial g_n}{\partial b} = g_n \frac{n b}{2 \sqrt{a^2 - b^2}} \frac{(\sqrt{a + b} + \sqrt{a - b})^{n-2} + (\sqrt{a + b} - \sqrt{a - b})^{n-2}}
                                                                                           {(\sqrt{a + b} + \sqrt{a - b})^n - (\sqrt{a + b} - \sqrt{a - b})^n}.\end{aligned}$$ Furthermore, in appendix the following integral formula are derived:

\[int-formula-2\] $$\label{int-formula-2-1}
\int dx_1 \cdots dx_n (x_1^4 + \cdots + x_n^4)  \exp \left[- {\bm X} G_n {\bm X}^{\dagger} \right] = g_n \frac{3 n}{4} \frac{(\mbox{det} H_{n-1})^2}{(\det G_n)^2}$$ $$\begin{aligned}
\label{int-formula-2-2}
&&\int dx_1 \cdots dx_n (x_1^2 x_2^2 + \cdots + x_{n-1}^2 x_n^2 + x_n^2 x_1^2) \exp \left[- {\bm X} G_n {\bm X}^{\dagger} \right]                                       \\    \nonumber
&&= g_n \frac{n}{4 (\det G_n)^2} \Bigg[ 12 a^2 (\det H_{n-2})^2 - 12 a b^2 (\det H_{n-2}) (\det H_{n-3})                                                                 \\    \nonumber
&& \hspace{5.0cm}  + 3 b^4 (\det H_{n-3})^2 - 2 (\det H_{n-2}) (\det G_n) \bigg]                                                                                         \end{aligned}$$ $$\begin{aligned}
\label{int-formula-2-3}
&&\int dx_1 \cdots dx_n \left[ x_1 x_2 (x_1^2 + x_2^2) + \cdots + x_{n-1} x_n (x_{n-1}^2 + x_n^2) + x_n x_1 (x_n^2 + x_1^2) \right]                                       \\    \nonumber
&& \hspace{8.0cm}  \times \exp \left[- {\bm X} G_n {\bm X}^{\dagger} \right]    \\    \nonumber
&&= g_n \frac{3 n}{2} \frac{ \left[ b^{n-1} + b (\det H_{n-2}) \right] \left[ 2 a (\det H_{n - 2}) - b^2 (\det H_{n - 3}) \right]}{(\det G_n)^2}\end{aligned}$$

where $$\begin{aligned}
\label{det-Hn}
&&\det H_n = \frac{1}{2^n \sqrt{a^2 - b^2}} \Bigg[ a \left[ \left( \sqrt{a + b} + \sqrt{a - b} \right)^{2 n} - \left( \sqrt{a + b} - \sqrt{a - b} \right)^{2 n} \right]                 \\    \nonumber
&&      \hspace{4.0cm}                                      - b^2 \left[ \left( \sqrt{a + b} + \sqrt{a - b} \right)^{2 n - 2} - \left( \sqrt{a + b} - \sqrt{a - b} \right)^{2 n - 2} \right]  \Bigg].\end{aligned}$$ Inserting Eqs. (\[thermal-2-boso-2\]), (\[Gn-boso-1\]), (\[int-formula-1\]), and (\[int-formula-2\]) into Eq. (\[trace-2\]), one can derive $$\label{trace-rho-n}
\mbox{Tr} \rho_T^n = \frac{2^{n-1} \sinh^n \frac{x}{2}} {\sinh \frac{n x}{2}}
\left[ 1 + \frac{3 \alpha n}{4} (\hbar m \omega x) \left( \coth^2 \frac{x}{2} - \coth^2 \frac{n x}{2} \right) + {\cal O} (\alpha^2) \right].$$

From Eq. (\[trace-rho-n\]) the purity function is given by $$\label{purity}
P (T) \equiv \mbox {Tr} \rho_T^2 = \tanh \frac{x}{2} \left[ 1 + \frac{3 \alpha}{2} (\hbar m \omega x) \left( \coth^2 \frac{x}{2} - \coth^2 x \right) + {\cal O} (\alpha^2) \right].$$

![(Color online) The external temperature dependence of the purity function when $\hbar = m = \omega = k_B = 1$. The black solid, red dashed, and blue dotted lines correspond to $\alpha = 0$, $\alpha = 0.04$, and $\alpha = 0.08$, respectively. As expected, the thermal state $\rho_T$ becomes pure and completely mixed state at $T = 0$ and $T = \infty$ when $\alpha = 0$. When, however, $\alpha \neq 0$, the state reduces to the partially mixed state even in the limit of $T \rightarrow \infty$. ](fig2.pdf){height="6.0cm"}

The external temperature dependence of the purity function is plotted in Fig. 2 with varying the GUP parameter $\alpha$. For convenience we fix other constants as $\hbar = m = \omega = k_B = 1$. In Fig. 1 the black solid, red dashed, and blue dotted lines correspond to $\alpha = 0$, $\alpha = 0.04$, and $\alpha = 0.08$, respectively. As expected, the thermal state $\rho_T$ becomes pure and completely mixed state at $T = 0$ and $T = \infty$ when $\alpha = 0$. When, however, $\alpha \neq 0$, the state reduces to the partially mixed state even in the limit of $T \rightarrow \infty$. This can be easily seen from $\lim_{T \rightarrow \infty} P(T) = \frac{9 \alpha}{4}$.

Rényi and von Neumann entropies
===============================

![(Color online) The external temperature dependence of the von Neumann entropy when $\hbar = m = \omega = k_B = 1$. The black solid, red dashed, and blue dotted lines correspond to $\alpha = 0$, $\alpha = 0.01$, and $\alpha = 0.02$, respectively. When $\alpha = 0$, the von Neumann entropy monotonically increases with increasing temperature. The remarkable fact is that the von Neumann entropy with nonzero $\alpha$ exhibits a decreasing behavior at the high-temperature region. ](fig3.pdf){height="6.0cm"}

Using Eq. (\[trace-rho-n\]) it is possible to derive the eigenvalue $\lambda_n$ in the eigenvalue equation (\[eigen-1\]) as following. Now, we take a trial solution of $\lambda_n$ in a form $$\label{value-1}
\lambda_n (\beta) = (1 - e^{-x}) e^{-n x} \left[ 1 + \alpha (\hbar m \omega x) h_n (\beta) + {\cal O} (\alpha^2) \right].$$ Since $$\label{value-2}
\sum_{n=0}^{\infty} \left( \lambda_n \right)^k = \mbox{Tr} \rho_T^k$$ for all positive integer $k$, one can derive a condition of $h_n (\beta)$ as $$\label{value-3}
\sum_{n=0}^{\infty} e^{-n k x} h_n (\beta) = \frac{3}{4 \sinh^2 \frac{x}{2}} \frac{1}{1 - e^{-k x}} + \frac{3}{(1 - e^{- k x})^2} - \frac{3}{(1 - e^{- k x})^3}.$$ Eq. (\[value-3\]) can be solved as $$\label{value-4}
h_n (\beta) = 3 \left[ \frac{e^{- x}}{(1 - e^{- x})^2} - \frac{n (n + 1)}{2} + \zeta_n (\beta) \right],$$ where $\zeta_n$ should satisfy $\sum_{n=0}^{\infty} e^{-n k x} \zeta_n (\beta) = 0$ for all positive integer $k$. Inserting Eq. (\[value-4\]) into Eq. (\[value-1\]), one can show $$\label{value-5}
\lambda_n (\beta) = (1 - e^{-x}) e^{-n x} \left[ 1 + 3 \alpha (\hbar m \omega x) \left\{\frac{e^{- x}}{(1 - e^{- x})^2} - \frac{n (n + 1)}{2} + \zeta_n (\beta) \right\} + {\cal O} (\alpha^2) \right].$$ It is straightforward to show $\sum_{n=0}^{\infty} \lambda_n = 1$, which is consistent with $\mbox{Tr} \rho_T = 1$.

Now, we will show $\zeta_n (\beta) = 0$ as following. From Eq. (\[trace-rho-n\]) the Rényi entropy with integer order is $$\label{renyi-1}
S_{\gamma = n} = \frac{1}{1 - n} \ln {\mbox Tr} \rho_T^n.$$ If we take $n \rightarrow 1$ limit in $S_{\gamma = n}$, one can derive the von Neumann entropy explicitly in a fome: $$\label{von-1}
S_{von} = - \ln (1 - e^{-x}) +  \frac{ x e^{-x}}{1 - e^{-x}} 
- 3 \alpha (\hbar m \omega x^2) \frac{e^{-x} (1 + e^{-x})}{(1 - e^{-x})^3} + {\cal O} (\alpha^2).$$ The von Neumann entropy should coincide with $-\sum_{n=0}^{\infty} \lambda_n \ln \lambda_n$, which gives $\zeta_n (\beta) = 0$.

The temperature dependence of the von Neumann entropy is plotted in Fig. $3$ when $\hbar = m = \omega = k_B = 1$. The black solid, red dashed, and blue dotted lines correspond to $\alpha = 0$, $\alpha = 0.01$, and $\alpha = 0.02$, respectively. Globally, the nonzero GUP parameter $\alpha$ decreases the von Neumann entropy. This is due to the fact that the last term of Eq. (\[von-1\]) has a minus sign. The remarkable fact is that while the von Neumann entropy with $\alpha = 0$ monotonically increases with increasing $T$, the von Neumann entropy with nonzero $\alpha$ has a maximum at finite temperature $T_*$. This temperature can be computed by solving $$(1 - e^{-\beta_* \omega})^2 - 3 \alpha (\hbar m \omega) \left[ (-2 + \beta_* \omega) + 4 \beta_* \omega e^{-\beta_* \omega} + (2 + \beta_* \omega) e^{-2 \beta_* \omega} \right] + {\cal O} (\alpha^2) = 0,$$ where $\beta_* = 1 / T_*$. Mathematically, this extraordinary behavior of the von Neumann entropy with nonzero $\alpha$ can be understood as a competition of the last term in Eq. (\[von-1\]) with the remaining terms. However, it is unclear, at least for us, to understand this behavior from a point of physics.

![(Color online)The external temperature dependence of the Rényi entropy when $\hbar = m = \omega = k_B = 1$ when (a) $\alpha = 0$ (b) $\alpha = 0.01$. In both figures the black solid, red dashed, and blue dotted lines correspond to $\gamma = 0.8$, $\gamma = 1.8$, and $\gamma = \infty$, respectively. When $\alpha = 0$, the Rényi entropy decreases with increasing the order $\gamma$. This is well-known in the HUP quantum mechanics. The remarkable fact is that the Rényi entropy with nonzero $\alpha$ exhibits a decreasing behavior at the high-temperature region. ](fig4a.pdf "fig:"){height="5.0cm"} ![(Color online)The external temperature dependence of the Rényi entropy when $\hbar = m = \omega = k_B = 1$ when (a) $\alpha = 0$ (b) $\alpha = 0.01$. In both figures the black solid, red dashed, and blue dotted lines correspond to $\gamma = 0.8$, $\gamma = 1.8$, and $\gamma = \infty$, respectively. When $\alpha = 0$, the Rényi entropy decreases with increasing the order $\gamma$. This is well-known in the HUP quantum mechanics. The remarkable fact is that the Rényi entropy with nonzero $\alpha$ exhibits a decreasing behavior at the high-temperature region. ](fig4b.pdf "fig:"){height="5.0cm"}

The Rényi entropy with arbitrary real order $\gamma$ is given by $$\begin{aligned}
\label{renyi-final}
S_{\gamma} &\equiv& \frac{1}{1 - \gamma} \ln \sum_{n=0}^{\infty} \left( \lambda_n \right)^{\gamma}                         \\    \nonumber
&=& \frac{1}{1 - \gamma} \Bigg[\gamma \ln (1 - e^{-x}) - \ln (1 - e^{-\gamma x})                     \\    \nonumber
&&       \hspace{2.0cm}    + 3 \alpha \gamma (\hbar m \omega x) \left\{ \frac{e^{-x}}{(1 - e^{-x})^2} - \frac{e^{-\gamma x}}{(1 - e^{-\gamma x})^2} \right\} + {\cal O} (\alpha^2) \Bigg].\end{aligned}$$ In Fig. $4$ the temperature dependence of the Rényi entropy is plotted when $\alpha = 0$ (Fig. $4$(a)) and $\alpha = 0.01$ (Fig. $4$(b)) with choosing when $\hbar = m = \omega = k_B = 1$. When $\alpha = 0$, the the Rényi entropies with various $\gamma$ exhibit increasing behavior with increasing temperature. With increasing the order $\gamma$ the entropy decreases gradually and eventually approaches to blue dotted line, which corresponds to $\gamma = \infty$. When $\alpha = 0.01$, however, the behavior of the Rényi entropy is drastically changed at large temperature region. In this region the entropies exhibit decreasing behavior in temperature. Furthermore, with increasing the order $\gamma$ it increases and approaches to the blue dotted line. We do not exactly know how to interpret this extraordinary behavior physically.

Conclusion
==========

We compute in the paper the Rényi and von Neumann entropies of the thermal state $\rho_T$ within the first order of $\alpha$. Both entropies with $\alpha = 0$ exhibit monotonically increasing behavior with increasing the external temperature. It is also found that the nonzero $\alpha$ reduces the entropies compared with $\alpha = 0$ case. Furthermore, there is a striking difference between $\alpha = 0$ and $\alpha \neq 0$. For the case of von Neumann entropy it has a maximum at finite temperature $T_*$. At $T > T_*$ the von Neumann entropy exhibits a decreasing behavior with increasing the temperature.

For the case of the Rényi entropy $S_{\gamma} (T)$ it reduces with increasing the order $\gamma$ when $\alpha = 0$. If, however, $\alpha \neq 0$, the temperature dependence of the Rényi entropy is drastically changed at the large temperature region. In this region the Rényi entropy exhibits a decreasing behavior with increasing the external temperature. The decreasing rate becomes larger with decreasing the order $\gamma$. As a result, the Rényi entropy reduces with decreasing the order $\gamma$ at the large temperature region.

It is of interest to extend our case to one-dimensional two-coupled harmonic oscillator system, whose Hamiltonian is $$\begin{aligned}
\label{two-coupled-1}
H_2 &=& \frac{1}{2 m} (P_1^2 + P_2^2 ) + \frac{1}{2} \left[ k_0 (X_1^2 + X_2^2) + J (X_1 - X_2)^2  \right]    \\     \nonumber  
&=& h_1 + h_2 + \frac{1}{2} J (x_1 - x_2)^2 + {\cal O} (\alpha^2),\end{aligned}$$ where $$\label{two-coupled-2}
h_j = \frac{1}{2 m} \left( p_j^2 + 2 \alpha p_j^4 \right) + \frac{1}{2} k_0 x_j^2$$ with $j = 1, 2$. It is possible to derive the vacuum and thermal states of this coupled system. It is of interest to compute the entanglement of formation[@benn96] of these states and compare with the HUP cases[@park18; @park19]. We guess the GUP parameter protects the entanglement against the external temperature. However, it should be checked by explicit calculation. We hope to address this issue in the future.

[99]{} P. K. Townsend, [*Small-scale structure of spacetime as the origin of the gravitational constant*]{}, Phys. Rev. [**D 15**]{} (1977) 2795. D. Amati, M. Ciafaloni, and G. Veneziano, [*Can spacetime be probed below the string size?*]{}, Phys. Lett. [**B 216**]{} (1989) 41. L. J. Garay, [*Quantum gravity and minimum length*]{}, Int. J. Mod. Phys. [**A 10**]{} (1995) 145 \[gr-qc/9403008\]. C. Rovelli and L. Smolin, [*Knot Theory and Quantum Gravity*]{}, Phys. Rev. Lett. [**61**]{} 1988 1155. C. Rovelli and L. Smolin, [*Loop space representation of quantum general relativity*]{}, Nucl. Phys. [**B 331**]{} (1990) 80. C. Rovelli, [*Loop Quantum Gravity*]{}, Living Rev. Relativity, [**1**]{} (1998) 1 \[gr-qc/9710008\]. S. Carlip, [*Quantum Gravity: a Progress Report*]{}, Rep. Prog. Phys. [**64**]{} (2001) 885 \[gr-qc/0108040\]. K. Konishi, G. Paffuti, and P. Provero, [it Minimum physical length and the generalized uncertainty principle in string theory]{}, Phys. Lett. [**B 234**]{} (1990) 276. M. Kato, [*Particle theories with minimum observable length and open string theory*]{}, Phys. Lett. [**B 245**]{} (1990) 43. A. Strominger, [*Quantum Gravity and String Theory, What Have We Learned?*]{}, hep-th/9110011. T. Padmanabhan, [*Physical significance of planck length*]{}, Ann. Phys. [**165**]{} (1985) 38. T. Padmanabhan, [*Planck length as the lower bound to all physical length scales*]{}, Gen. Rel. Grav. [**17**]{} (1985) 215. T. Padmanabhan, [*The role of general relativity in the uncertainty principle*]{}, Class. Q. Grav. [**3**]{} (1986) 911. T. Padmanabhan, [*Limitations on the operational definition of spacetime events and quantum gravity*]{}, Class. Quant. Grav. [**4**]{} (1987) L107. J. Greensite, [Is there a minimum length in D=4 lattice quantum gravity?]{}, Phys. Lett. [**B 255**]{} (1991) 375. M. Maggiore, [*A Generalized Uncertainty Principle in Quantum Gravity*]{}, Phys. Lett. [**B 304**]{} (1993) 65 \[hep-th/9301067\]. C. A. Mead, [*Possible Connection Between Gravitation and Fundamental Length*]{}, Phys. Rev. [**135**]{} (1964) B849. W. Heisenberg, Über den anschaulichen Inhalt der quantentheoretischen Kinematik und Mechanik, Z. Phys. [**43**]{} (1927) 172. H. P. Robertson, [*The Uncertainty Principle*]{}, Phys. Rev. [**34**]{} (1929) 163. A. Kempf, [*Uncertainty Relation in Quantum Mechanics with Quantum Group Symmetry*]{}, J. Math. Phys. [**35**]{} (1994) 4483 \[hep-th/9311147\]. A. Kempf, G. Mangano, and R. B. Mann, [*Hilbert Space Representation of the Minimal Length Uncertainty Relation*]{}, Phys. Rev. [**D 52**]{} (1995) 1108 \[hep-th/9412167\]. W.S. Chung and H. Hassanabadi, [*A new higher order GUP: one dimensional quantum system*]{}, Eur. Phys. J. [**C 79**]{} (2019) 213. P. Bosso, [*Rigorous Hamiltonian and Lagrangian analysis of classical and quantum theories with minimal length*]{}, Phys. Rev. [**D 97**]{} (2018) 126010 \[arXiv:1804.08202 (hep-th)\]. P. Bosso, O. Obregón, [*Minimal Length Effects on Quantum Cosmology and Quantum Black Hole Models*]{}, Class.Quant.Grav. [**37**]{} (2020) 045003 \[arXiv:1904.06343 (hep-th)\]. C. Villalpando and S.K. Modak, [*Minimal length effect on the broadening of free wave packets and its physical implications*]{}, Phys. Rev. [**D 100**]{} (2019) 024054 \[arXiv:1812.06112 (qr-qc)\]. L. Xiang, Y. Ling, Y.-G. Shen, C.-Z. Liu, H.-S. He, and L.-F. Xu, [*Generalized uncertainty principles, effective Newton constant and the regular black hole*]{}, Ann. Phys. [**396**]{} (2018) 334 \[arXiv:1611.09016 (gr-qc)\]. R. P. Feynman and A. R. Hibbs, Quantum Mechanics and Path Integrals (McGraw-Hill, 1965, New York). H. Kleinert, Path integrals in Quantum Mechanics, Statistics, and Polymer Physics (World Scientific,1995, Singapore). S. Das and S. Pramanik, [*Path Integral for non-relativistic Generalized Uncertainty Principle corrected Hamiltonian*]{}, Phys. Rev. [**D 86**]{} (2012) 085004 \[arXiv:1205.3919 (hep-th)\]. S. Gangopadhyay and S. Bhattacharyya, [*Path-integral action of a particle with the generalized uncertainty principle and correspondence with noncommutativity*]{}, Phys. Rev. [**D 99**]{} (2019) 104010 \[arXiv:1901.03411 (quant-ph)\]. DaeKil Park and Eylee Jung, [*Comment on “Path-integral action of a particle with the generalized uncertainty principle and correspondence with noncommutativity”*]{}, Phys. Rev. [**D 101**]{} (2020) 068501 \[arXiv:2002.07954 (quant-ph)\]. DaeKil Park, [Generalized uncertainty principle and $d$-dimensional quantum mechanics]{}, Phys. Rev. [**D 101**]{} (2020) 106013 \[arXiv:/3109031(quant-ph\]. DaeKil Park and Eylee Jung, [*Generalized uncertainty principle and point interaction*]{}, Phys. Rev. [**D 101**]{} (2020) 066007 \[arXiv:2001.02850 (quant-ph)\]. C. H. Bennett, D. P. DiVincenzo, J. A. Smokin and W. K. Wootters, [*Mixed-state entanglement and quantum error correction*]{}, Phys. Rev. [**A 54**]{} (1996) 3824 \[quant-ph/9604024\]. DaeKil Park, [*Dynamics of entanglement and uncertainty relation in coupled harmonic oscillator system: exact results*]{}, Quant. Inf. Proc. [**17**]{} (2018) 147 \[arXiv:1801.07070 (quant-ph)\]. DaeKil Park, [*Thermal Entanglement Phase Transition in Coupled Harmonic Oscillators with Arbitrary Time-Dependent Frequencies*]{}, Quant. Inf. Proc. [**19**]{} (2020) 129 \[arXiv:1903.03297 (quant-ph)\].

**Appendix A: Derivation of Eq. (\[int-formula-2-1\])**

In this appendix we derive Eq. (\[int-formula-2-1\]). First, let us consider a $n \times n$ matrix $$\begin{aligned}
\label{matrix-tilde-Gn}
\tilde{G}_n = \left(                  \begin{array}{cccccc}
                  2 a' & \hspace{0.2cm}    -b    & \hspace{0.2cm}         & \hspace{0.2cm}         & \hspace{0.2cm}         & \hspace{0.2cm}  -b              \\
                  -b  & \hspace{0.2cm} 2 a      & \hspace{0.2cm} \bullet & \hspace{0.2cm}         & \hspace{0.2cm}         & \hspace{0.2cm}                \\
                      & \hspace{0.2cm} \bullet  & \hspace{0.2cm} \bullet & \hspace{0.2cm} \bullet & \hspace{0.2cm}         & \hspace{0.2cm}                \\
                      & \hspace{0.2cm}          & \hspace{0.2cm} \bullet & \hspace{0.2cm} \bullet & \hspace{0.2cm} \bullet & \hspace{0.2cm}                \\
                      & \hspace{0.2cm}          & \hspace{0.2cm}         & \hspace{0.2cm} \bullet & \hspace{0.2cm} \bullet & \hspace{0.2cm}   -b           \\
                   -b   & \hspace{0.2cm}          & \hspace{0.2cm}         & \hspace{0.2cm}         & \hspace{0.2cm}     -b  & \hspace{0.2cm}  2a 
                              \end{array}                        \right).\end{aligned}$$ Then one can show $$\label{app-a-1}
\det \tilde{G}_n = \det G_n + 2 (a' - a) \det H_{n-1}$$ where $\det G_n$ is given by Eq. (\[det-Gn\]) and $H_n$ is a following $n \times n$ tridiagonal matrix: $$\begin{aligned}
\label{matrix-Hn}
H_n = \left(                  \begin{array}{cccccc}
                  2 a & \hspace{0.2cm}    -b    & \hspace{0.2cm}         & \hspace{0.2cm}         & \hspace{0.2cm}         & \hspace{0.2cm}                \\
                  -b  & \hspace{0.2cm} 2 a      & \hspace{0.2cm} \bullet & \hspace{0.2cm}         & \hspace{0.2cm}         & \hspace{0.2cm}                \\
                      & \hspace{0.2cm} \bullet  & \hspace{0.2cm} \bullet & \hspace{0.2cm} \bullet & \hspace{0.2cm}         & \hspace{0.2cm}                \\
                      & \hspace{0.2cm}          & \hspace{0.2cm} \bullet & \hspace{0.2cm} \bullet & \hspace{0.2cm} \bullet & \hspace{0.2cm}                \\
                      & \hspace{0.2cm}          & \hspace{0.2cm}         & \hspace{0.2cm} \bullet & \hspace{0.2cm} \bullet & \hspace{0.2cm}   -b           \\
                      & \hspace{0.2cm}          & \hspace{0.2cm}         & \hspace{0.2cm}         & \hspace{0.2cm}     -b  & \hspace{0.2cm}  2a 
                              \end{array}                        \right).\end{aligned}$$ It is easy to show that $\det H_n$ satisfies the recursion relation $$\label{app-a-2}
\det H_n = 2 a \det H_{n-1} - b^2 \det H_{n-2}$$ with $\det H_0 = 1$ and $\det H_1 = 2 a$. Using Eq. (\[app-a-2\]) it is possible to show $$\label{app-a-3}
\det H_n = \frac{1}{\sqrt{a^2 - b^2}} \left[ a \sqrt{\det G_{2 n}} - \frac{b^2}{2} \sqrt{\det G_{2 n - 2}} \right]$$ where $\det G_n$ is given in Eq. (\[det-Gn\]). If Eq. (\[det-Gn\]) is inserted into Eq. (\[app-a-3\]), $\mbox{det} H_n$ reduces to Eq. (\[det-Hn\]).

Thus, it is easy to show $$\begin{aligned}
\label{app-a-4}
&&\int dx_1 \cdots dx_n x_1^4 \exp \left[- {\bm X} \tilde{G}_n {\bm X}^{\dagger} \right]                                   \\    \nonumber
&&= \frac{1}{4} \frac{\partial^2}{\partial a'^2}  \frac{\pi^{n/2}}{\sqrt{\det G_n + 2 (a' - a) \det H_{n-1}}}\Bigg|_{a' = a}    \\    \nonumber
&&= g_n \frac{3 \left(\det H_{n-1} \right)^2}{4 \left( \det G_n \right)^2}.\end{aligned}$$

Next, by making use of the redefinition of the integral variables it is possible to show that $$\int dx_1 \cdots dx_n x_j^4 \exp \left[- {\bm X} \tilde{G}_n {\bm X}^{\dagger} \right]$$ is independent of $j$. Thus, $$\label{app-a-5}
\int dx_1 \cdots dx_n \left(x_1^4 + \cdots + x_n^4 \right) \exp \left[- {\bm X} \tilde{G}_n {\bm X}^{\dagger} \right] 
= n \int dx_1 \cdots dx_n x_1^4 \exp \left[- {\bm X} \tilde{G}_n {\bm X}^{\dagger} \right],$$ which exactly coincides with Eq. (\[int-formula-2-1\]).

**Appendix B: Derivation of Eq. (\[int-formula-2-2\])**

In this appendix we derive Eq. (\[int-formula-2-2\]). Let us consider a $n \times n$ matrix $$\begin{aligned}
\label{matrix-dasi-Gn}
G_n' = \left(                  \begin{array}{cccccc}
                  2 a_1 & \hspace{0.2cm}    -b    & \hspace{0.2cm}         & \hspace{0.2cm}         & \hspace{0.2cm}         & \hspace{0.2cm}  -b              \\
                  -b  & \hspace{0.2cm} 2 a_2     & \hspace{0.2cm} \bullet & \hspace{0.2cm}         & \hspace{0.2cm}         & \hspace{0.2cm}                \\
                      & \hspace{0.2cm} \bullet  & \hspace{0.2cm} 2 a & \hspace{0.2cm} \bullet & \hspace{0.2cm}         & \hspace{0.2cm}                \\
                      & \hspace{0.2cm}          & \hspace{0.2cm} \bullet & \hspace{0.2cm} \bullet & \hspace{0.2cm} \bullet & \hspace{0.2cm}                \\
                      & \hspace{0.2cm}          & \hspace{0.2cm}         & \hspace{0.2cm} \bullet & \hspace{0.2cm} \bullet & \hspace{0.2cm}   -b           \\
                   -b   & \hspace{0.2cm}          & \hspace{0.2cm}         & \hspace{0.2cm}         & \hspace{0.2cm}     -b  & \hspace{0.2cm}  2a 
                              \end{array}                        \right).\end{aligned}$$ Then, one can show straightforwardly $$\label{app-b-1}
\det G_n' = 4 (a_1 - a) (a_2 - a) \det H_{n-2} + 4 a (a_1 + a_2 - 2 a) \det H_{n-2} - 2 b^2 (a_1 + a_2 - 2 a) \det H_{n-3} + \det G_n.$$ Thus, Eq. (\[int-formula-2-2\]) is verified as following: $$\begin{aligned}
\label{app-b-2}
&&\int dx_1 \cdots dx_n \left(x_1^2 x_2^2 + \cdots + x_{n-1}^2 x_n^2 + x_n^2 x_1^2 \right) \exp \left[- {\bm X} G_n {\bm X}^{\dagger} \right]           \\   \nonumber
&&= \frac{n}{4} \frac{\partial^2}{\partial a_1 \partial a_2} \frac{\pi^{n/2}}{\sqrt{\det G_n'}} \Bigg|_{a_1 = a_2 = a}                                          \\   \nonumber
&&= g_n \frac{n}{4 \left( \det G_n \right)^2} \Bigg[ 12 a^2 (\det H_{n-2})^2 - 12 a b^2 (\det H_{n-2}) (\det H_{n-3})                                           \\   \nonumber
&& \hspace{5.0cm} + 3 b^4 (\det H_{n-3})^2 - 2 (\det H_{n-2}) (\det G_n ) \Bigg],\end{aligned}$$ which exactly coincides with Eq. (\[int-formula-2-2\]).

**Appendix C: Derivation of Eq. (\[int-formula-2-3\])**

In this appendix we derive Eq. (\[int-formula-2-3\]). Let us consider a $n \times n$ matrix

$$\begin{aligned}
\label{matrix-twodasi-Gn}
G_n'' = \left(                  \begin{array}{cccccc}
                  2 a_1 & \hspace{0.2cm}    -b_1    & \hspace{0.2cm}         & \hspace{0.2cm}         & \hspace{0.2cm}         & \hspace{0.2cm}  -b              \\
                  -b_1  & \hspace{0.2cm} 2 a_1     & \hspace{0.2cm} -b & \hspace{0.2cm}         & \hspace{0.2cm}         & \hspace{0.2cm}                \\
                      & \hspace{0.2cm} -b  & \hspace{0.2cm} 2 a & \hspace{0.2cm} \bullet & \hspace{0.2cm}         & \hspace{0.2cm}                \\
                      & \hspace{0.2cm}          & \hspace{0.2cm} \bullet & \hspace{0.2cm} \bullet & \hspace{0.2cm} \bullet & \hspace{0.2cm}                \\
                      & \hspace{0.2cm}          & \hspace{0.2cm}         & \hspace{0.2cm} \bullet & \hspace{0.2cm} \bullet & \hspace{0.2cm}   -b           \\
                   -b   & \hspace{0.2cm}          & \hspace{0.2cm}         & \hspace{0.2cm}         & \hspace{0.2cm}     -b  & \hspace{0.2cm}  2a 
                              \end{array}                        \right).\end{aligned}$$

Then, one can show straightforwardly $$\begin{aligned}
\label{app-c-1}
&&\det G_n'' = \left[4 (a_1 - a) (a_1 + a) - (b_1 - b) (b_1 + b) \right] \det H_{n-2}                       \\    \nonumber
&& \hspace{2.0cm}- 4 b^2 (a_1 - a) \det H_{n-3} - 2 (b_1 - b) b^{n-1} + \det G_n.\end{aligned}$$ Thus, Eq. (\[int-formula-2-3\]) is proved as following: $$\begin{aligned}
\label{app-c-2}
&&\int dx_1 \cdots dx_n \left\{x_1 x_2 (x_1^2 + x_2^2) + \cdots + x_{n-1} x_n (x_{n-1}^2 + x_n^2) + x_n x_1 (x_n^2 + x_1^2) \right\}                             \\   \nonumber   
&&  \hspace{9.0cm}       \times \exp \left[- {\bm X} G_n {\bm X}^{\dagger} \right]                                                         \\   \nonumber
&&= -\frac{n}{4} \frac{\partial^2}{\partial a_1 \partial b_1} \frac{\pi^{n/2}}{\sqrt{\det G_n''}} \Bigg|_{a_1 = a, b_1 = b}                                          \\   \nonumber
&&= g_n \frac{3 n}{2} \frac{ \left[ b^{n-1} + b (\det H_{n-2}) \right] \left[ 2 a (\det H_{n - 2}) - b^2 (\det H_{n - 3}) \right]}{(\det G_n)^2}\end{aligned}$$ which exactly coincides with Eq. (\[int-formula-2-3\]).

[^1]: dkpark@kyungnam.ac.kr
