Holographic entanglement entropy in the fourth-order scale-invariant gravity
الموضوعات : Journal of Theoretical and Applied PhysicsTahereh Hamedi 1 , M. Reza Tanhayi 2
1 - Department of Physics, Central Tehran Branch, Islamic Azad University (IAUCTB), Tehran, Iran
2 - Department of Physics, Central Tehran Branch, Islamic Azad University (IAUCTB), Tehran, Iran
الکلمات المفتاحية: Entanglement Entropy, Holographic method, Log gravity, mutual information,
ملخص المقالة :
In this paper, we use the holographic method to compute the entanglement entropy for different regions in the fourth-order scale-invariant theory of gravity. In four dimensions, the action of scale-invariant gravity contains a parameter where the moduli space has some distinguished points, precisely, the space of solutions contains the Log-gravity and for this specific solution, we compute the entanglement entropy. We also use the numerical method to investigate mutual information and tripartite information. Moreover, we make a comment on the sign of these quantities.
\documentclass[a4paper,12pt]{article}
\newcommand{\head}[1]{\textnormal{\textbf{#1}}}
\usepackage[onehalfspacing]{setspace}
\usepackage{amsmath}
\usepackage{amssymb}
\usepackage{hyperref}
\usepackage{fontenc}
\usepackage{inputenc}
\usepackage{authblk}
\usepackage{color}
\usepackage{graphicx}
\usepackage{cite}
\usepackage{slashed}
\hypersetup{colorlinks=true,linkcolor=magenta,linktoc=page,citecolor=blue}
\usepackage{tikz}
\usepackage{multirow}
\usepackage{hhline}
\usepackage[top=1.5cm,bottom=2.5cm,right=2.54cm,left=2.54cm]{geometry}
\definecolor{c1}{rgb}{0.5,0,1}
\colorlet{aqua}{-red!75}
\numberwithin{equation}{section}
%\begin{document}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\newcommand{\be}{\begin{equation}}
\newcommand{\bea}{\begin{eqnarray}}
\newcommand{\eea}{\end{eqnarray}}
\newcommand{\ba}{\begin{align}}
\newcommand{\ea}{\end{align}}
\newcommand{\ee}{\end{equation}}
\newcommand{\nn}{\nonumber}
\newcommand{\half}{\frac{1}{2}}
\newcommand{\Lt}{\tilde{L}}
\newcommand{\A}{{\cal{A}}}
\newcommand{\K}{{\cal{K}}}
\newcommand{\CO}{{\cal{O}}}
%%%% page setup for A4 portrait normal version %%%%
%\setlength{\topmargin}{-0.5cm}
%\setlength{\oddsidemargin}{-0.3cm}
%\setlength{\evensidemargin}{-0.3cm}
%\setlength{\textheight}{23cm}
%\setlength{\textwidth}{16.5cm}
\begin{document}
\begin{titlepage}
\vspace{10mm}
%\begin{flushright}
%IPM/P-2015/nnn \\
%FPAUO-1/10\\
%\end{flushright}
\vspace*{20mm}
\begin{center}
{\Large {\bf Holographic Entanglement Entropy in the Fourth-order Scale-invariant Gravity }\\
}
\vspace*{10mm} {T. Hamedi and M. Reza Tanhayi$^\star$}
\vspace*{.5cm}
{Department of Physics, Central Tehran Branch, Islamic Azad University (IAUCTB), P.O. Box 14676-86831,
Tehran, Iran}
\vspace*{0.5cm}
{E-mail: {\tt $^\star$mtanhayi@ipm.ir }}%
\vspace*{1cm}
%%\maketitle
\end{center}
\begin{abstract}
In this paper, we use the holographic method to compute the entanglement entropy for different regions in the fourth-order scale-invariant theory of gravity. In four dimensions, the action of scale-invariant gravity contains a parameter where the moduli space has some distinguished points, precisely, the space of solutions contains the Log-gravity and for this specific solution, we compute the entanglement entropy. We also use the numerical method to investigate mutual information and tripartite information. Moreover, we make a comment on the sign of these quantities.
{\tt Keywords}: Scale invariance; Conformal gravity; Log-gravity
\end{abstract}
%\tableofcontents
\end{titlepage}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Introduction and Motivation}
In principle, scale-invariance is an important feature of a theory implies that if the variables, such as energy, length or other scales are multiplied by a certain factor, then the content of the theory does not change, and in this sense, it may present a universality of the theory. The underlying transformation of scale-invariance is identified by the dilatation, which dilatations can be considered as a subclass of a larger family of transformations named by conformal symmetry. The scale-invariant theories have many applications in most ares of physics: For example, the theory of phase transitions in statistical mechanics; the strength of particle interactions in elementary particle physics are described by the scale-invariant theory. Moreover, it is argued that a scale-invariant function shall be used to describe the power spectrum of the spatial distribution of the cosmic microwave background. Additionally, this theory has its own application in gravity theory in describing the issue of dark matter \cite{Nojiri:2006ri, Tanhayi:2011dh, Mohsenzadeh:2013bba}. In our previous paper \cite{Hamedi:2022pfu,me1}, we investigated the scale-invariant gravity theory in four dimensions, and showed that the theory has a critical point, in which, the theory admits several non-trivial solutions such as a logarithmic solution. In this paper, in order to investigate the other features of the theory, we use the holographic methods which based on Anti-de Sitter/Conformal Field Theory (AdS/CFT) correspondence, and consider some non-local probes of entanglement. More precisely, we compute the holographic entanglement entropy (HEE) and mutual information for the logarithmic solution in the scale-invariant theory in four dimensions. First, let us review the theory. \\
In four dimensions the action of the scale-invariant gravity is given by
\bea\label{action-RS1}
{\cal I}=
-\frac{\kappa}{32\pi}\int \sqrt{-g}\,\,d^{4}x\,\left[(\sigma_0-6)R_{\mu\nu\rho\sigma}^2-2(\sigma_0-12)
R_{\mu\nu}^2+(\frac{\sigma_0}{3}-5) R^2\right],
\eea
in this notation, one receives conformal gravity at
point $\sigma_0\rightarrow \infty$ on the other hand, at $\sigma_0=0$ the $R^2$ gravity is obtained. We also argued there is a critical point at $\sigma_0=6$ where the theory reduces to the log gravity . The corresponding equations of motion are given by
\be\label{EOMT}
\left(\nabla^\sigma\nabla^\rho-\frac{1}{2}R^{\sigma\rho}\right)C_{\mu\sigma\nu\rho}=\frac{1}{2\sigma_0}
\left(R R_{\mu\nu}-g_{\mu\nu}\frac{R^2}{4}-\nabla_\mu\nabla_\nu R+g_{\mu\nu}\Box R\right),
\ee
where $\Box=\nabla^\mu\nabla_\mu$. The model has some black hole solutions \cite{Kehagias:2015ata}
\be\label{Sol1}
ds^2=\frac{L^2}{r^2}\left(-F(r)\;dt^2+\frac{dr^2}{F(r)}+d\Sigma^2_{2,k}\right),\;\;\;\;\;\;\;
F(r)=\lambda+k r^2+c_3{r^3},
\ee
where $L$ stands for the radius of curvature and $k=1,-1$ and $0$ corresponds to $\Sigma_{2,k}= S^2,H_2$ and $ R^2$, respectively. It is worth mentioning that
$$R_{\mu\nu}
=\frac{3\lambda}{L^2}g_{\mu\nu},\,\,\,\,\,\,\lambda=\pm 1,0$$
One can find that an AdS wave solution as follows (see
also \cite{Gullu:2011sj,Lu:2012xu})
\be\label{wave}
ds^2=\frac{L^2}{r^2}\bigg({dr^2}+dy^2-2dx_-dx_++k(x_+,r)\;dx_+^2\bigg),\;\;\;\;\;\;\;
k(x_+,r)=c_0(x_+)+c_3(x_+)r^3.
\ee
In order to explore the space of solutions, the AdS wave Ansatz can be utilized as follows
\be
ds^2=\frac{L^2}{r^2}\bigg({dr^2}+dy^2-2dx_-dx_++k(x_+,r)\;dx_+^2\bigg),
\ee
where this Ansatz together with \eqref{EOMT} leads one to write
\be
k(x_+,r)=c_0(x_+)+c_3(x_+) r^3+b_1(x_+) r^{\frac{3}{2}-\frac{1}{2} \sqrt{\frac{\sigma_0 +48}{\sigma_0 }}}+b_2(x_+) r^{\frac{3}{2}+\frac{1}{2} \sqrt{\frac{\sigma_0 +48}{\sigma_0 }}}.
\ee
Noting that the case $\sigma_0=6$, a logarithmic solution is obtained as follows
\bea\label{LOG}
k(x_+,r)=c_0(x_+)+b_0(x_+)\log r+\bigg(c_3(x_+)+b_3(x_+)\log r\bigg)r^3.
%&&\sigma_0=-48:\;\;\;\;\;\;k(x_+,r)=c_0(x_+)+c_3(x_+)r^3+\bigg(b_1(x_+)+b_2(x_+)\log r\bigg)
%r^{\frac{3}{2}}
\eea
By making use of the Ansatz \be
ds^2={dr^2}+dy^2-2dx_-dx_++k(x_+,r)\;dx_+^2,\;\;\;\;\;\;\;
\ee
the equations of motion reduces to
\be
\sigma_0\frac{\partial^4 k}{\partial^4r}=0,
\ee
and for $\sigma_0\neq 0$ the above equation has a non-trivial solution given by
\be
k(x_+,r)=c_0(x_+)+c_1(x_+)r+c_2(x_+)r^2+c_3(x_+)r^3.
\ee
Moreover, the theory admits the Lifshitz solution which is given by
\be
ds^2=\frac{L^2}{r^2}\left(-\frac{dt^2}{r^{2z}}+dr^2+dx_1^2+dx_2^2\right),\;\;\;\;\;\;\;
z=\frac{\sigma_0-6 +\sqrt{(\sigma_0 -6) (4 \sigma_0 +3)}}{\sigma_0 +3}
\ee
Therefore, at different points of moduli
space, the model has some non-trivial solutions. It is worth noting that for $\sigma_0=6$ the action \eqref{action-RS1} reads
\bea\label{actc}
{\cal I}=
-\frac{3 \kappa}{32\pi}\int d^{4}x\,\sqrt{-g}\,\left[4
R_{\mu\nu}^2- R^2\right].
\eea
We note, however, that although the above action has degenerate equations
of motion giving rise to a logarithmic solution, modifying the IR limit by adding a linear scalar
curvature term that would remove this solution. It is worth to mention that at $\sigma_0=0$, one obtains the $R^2$ gravity. In order
to make the finite action, one might use the Gauss-Bonnet term, where
the corresponding equations of motion become
\be\label{actR2}
\left( R_{\mu\nu}-\frac{1}{4}g_{\mu\nu}R\right)R+\bigg(g_{\mu\nu}\Box -
\nabla_\mu\nabla_\nu\bigg) R=0.
\ee
One can show that all the Einstein solutions, previously mentioned, are the
solutions to the above equations, as well. However, $R=0$ while $R_{\mu\nu}\neq 0$ can be considered as a new class of solutions. Namely, one can say that the theory has a larger class of the solutions, such as the Ricci scalar flat which is given by
\be\label{Sol3}
ds^2=\frac{L^2}{r^2}\left(-F(r)\;dt^2+\frac{dr^2}{F(r)}+d\Omega^2_{2,k}\right),\;\;\;\;\;\;\;
F(r)=k r^2+c_3 r^3+c_4 r^4,
\ee
where $c_3$ and $c_4$ are two constant, noting that the Ricci scalar becomes zero whereas, due to the $c_4\neq 0$, the Ricci tensor
is non-zero. This can be stated as follows
\be
R_{\mu\nu}=-\frac{c_4 r^4}{2L^2}\; (\eta_{\mu}^{\rho} g_{\rho\nu}+\eta_{\nu}^{\rho} g_{\rho\mu})
,\;\;\;\;\;\;{\rm with}\;\;
\eta_{\mu}^{\nu}={\rm diag}(-1,-1,1,1).
\ee
It is important to mention that for $k=1,c_4=c^2$ and $c_3=-2c$ where
$F(r)=r^2(1-c r)^2$ the solution reduces to an extremal black hole solution and the near horizon geometry
is given by ${\rm AdS}_2\times S^2$. By making use of $r=\frac{1}{\xi}$ one obtains
\be
ds^2=L^2\left[-\left(1-\frac{c}{\xi}\right)^2dt^2+\frac{d\xi^2}{\left(1-\frac{c}{\xi}\right)^2}+\xi^2
d\Omega_2^2\right],
\ee
with the following geometry
\be
ds^2=L^2c^2\left(-\rho^2dt^2+\frac{ d\rho^2}{\rho^2}+d\Omega_2^2\right),
\ee
where we have used $\xi-c=c^2 \rho$.
By setting $k=0$ one can obtain another class of solutions. Let us set $c_4=1$ and suppose $c_3$ is a free parameter. Therefore by setting $r=\frac{1}{\rho}$ the corresponding metric is given by
\be\label{hs}
ds^2=L^2\rho^2\left(-(1+c_3\rho) \frac{dt^2}{\rho^4}+\frac{d\rho^2}{1+c_3\rho}+dx_1^2+dx_2^2\right).
\ee
The above solution can be considered as the hyperscaling violating metric with $\theta=4, z=3$, noting that a $d+2$ dimensional
hyperscaling violating geometry is
\be
ds^2=\frac{L^2}{r^2} r^{2\frac{\theta}{d}}\left(-f(r)\frac{dt^2}{r^{2(z-1)}}+\frac{dr^2}{r}+d\vec{x}_d^2\right),\;\;\;\;\;\;\;\;{\rm with}\;\;f(r)=1-m\;r^{d-\theta+z}.
\ee
The rest of the paper is organized as follows. In the next section, we use holographic methods and study the entanglement entropy for three different entangling regions: Strip, circle and cusp. We also will explore different
features of the model, especially at the critical point, we compute the mutual information and also the tripartite information. We use numerical analysis, to explore the sign of these quantities, which become important the monogamy of mutual information.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Holographic Entanglement Entropy}
In this section in order to further explore the scale-invariant theory, let us use the holographic method and study the HEE in this theory. As mentioned in the previous section the theory has critical points leading to some non-trivial solutions. It is worth to mention that due to higher order terms, to compute the HEE one should use the generalized Ryu-Takayanagi prescription, \cite{{Fursaev:2013fta},{Dong:2013qoa},
{Camps:2013zua}}. First let us review the proposal briefly. For the following action
\be\label{ACT2}
{\cal I}=-\frac{\kappa}{32\pi}\int \sqrt{-g} \,d^{4}x\bigg(a R^2+b R_{\mu\nu}
R^{\mu\nu}+c R_{\mu\nu\rho\sigma}R^{\mu\nu\rho\sigma}\bigg),
\ee
thus it is argued that in order to obtain the HEE one should minimize the following functional
\begin{eqnarray}\label{EE}
S\!\!=\!\!\frac{ \kappa}{8}\! \int \;\sqrt{\gamma}\; \!d^2\zeta \bigg[2a R+b \left({ R}_{\mu\nu}n^\mu_i n^\nu_i-\frac{1}{2}\mathcal{K}^i\mathcal{K}_i\right)\!+\!2c\bigg( R_{\mu\nu\rho\sigma}n^\mu_i n^\nu_j n^\rho_i n^\sigma_j-\mathcal{K}^i_{\mu\nu}\mathcal{K}_i^{\mu\nu}\bigg)\!\bigg],
\end{eqnarray}
where $i=1,2$ refers to the two transverse directions to a co-dimension two
hypersurfaces in the bulk,
$n_i^\mu$ are two mutually orthogonal unit vectors to the hypersurface. ${\cal K}^{(i)}$ are the traces of two extrinsic curvature tensors defined by
\be \label{ECUR}
{\cal K}_{\mu\nu}^{(i)}=\pi^\sigma_{\ \mu} \pi^\rho_{\ \nu}\nabla_\rho(n_i)_\sigma,\;\;\;\;\;\;\;
{\rm with}\;\;\;\; \pi^\sigma_{\ \mu}=\epsilon^\sigma_{\ \mu}+\xi\sum_{i=1,2}(n_i)^\sigma(n_i)_\mu\ .
\ee
In the above relation, $\gamma$ is the
induced metric on the hypersurface whose coordinates are denoted by $\zeta$, and, $\xi=1$ and $\xi=-1$ are used for time-like and space-like vectors, receptively. Now we use this method and study the HEE for three different entangling regions: circle, strip and
cone.
\subsection{HEE of a Strip}
In order to find HEE for a strip entangling region which is given by ($t={\rm constant}$,)
\be
-\frac{\ell}{2}\leq y\leq \frac{\ell}{2}, \;\;\;\;\;\;\;\;0\leq x\leq H,
\ee
it is useful
to parametrize the AdS metric as follows
\be
ds^2=\frac{L^2}{r^2}(-dt^2+dr^2+dx^2+dy^2).
\ee
The co-dimension two hypersurface in the bulk can be parametrized by setting $y=f(r)$.
After doing some algebra, one obtains
\begin{eqnarray}
S=\kappa H\left[(\sigma_0-6)\left(\frac{1}{\epsilon }-\frac{2 \pi \Gamma \left(\frac{3}{4}\right)^2}{\Gamma \left(\frac{1}{4}\right)^2}\;\frac{1}{\ell}\right)+(\sigma_0-6)\left(-\frac{1}{\epsilon}\right)\right].
\end{eqnarray}
The above equation can be divided in two terms: The first term comes form the dynamical part of the action, whereas, second term is due to the Gauss-Bonnet term. The second one plays the role of regulator in the theory. Now let us compare the HEE with that of Einstein gravity which is
\begin{eqnarray}
S_{\rm Ein}=\frac{L^2 H}{2G}\left(\frac{1}{\epsilon }-\frac{2 \pi \Gamma \left(\frac{3}{4}\right)^2}{\Gamma \left(\frac{1}{4}\right)^2}\;\frac{1}{\ell}\right),
\end{eqnarray}
which by setting
$\kappa=\frac{L^2}{2(\sigma_0-6)G}$, the dynamical part is reproduced.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\subsection{HEE of a Circle}
To compute entanglement entropy for a circle it is useful to parametrize the
metric as follows
\be\label{AdSc}
ds^2=\frac{L^2}{r^2}(-dt^2+dr^2+d\rho^2+\rho^2d\phi),
\ee
and the entangling region is given by
\be
t={\rm contant},\;\;\;\;\;\;\;\;\;\;\;\;\rho=H={\rm constant},
\ee
by which, using the symmetry of the system, the corresponding co-dimension two hypersurface
in the bulk may be parametrized by $\rho=f(r)$. Therefore two mutually orthogonal unit vectors
to the hypersurface are given by
\be
n_1=\frac{L}{r}\bigg(1,0,0,0\bigg),\;\;\;\;\;\;\;\;\;\;\;\;\;\;n_2=\frac{L }{r \sqrt{1+f'(r)^2+1}}
\bigg(0,-f'(r),1,0\bigg).
\ee
Using these vectors and going through the minimization procedure one gets following expression
for the holographic entanglement entropy (see \cite{Bhattacharyya:2013gra} for five dimensional
case)
\be
S=\frac{\pi \kappa H}{2}(3 (4 \lambda _1+\lambda _2)+2 \lambda _3) \int_\epsilon^H dr \frac{ 1}{r^2}=\frac{\pi \kappa }{2}(3 (4 \lambda _1+\lambda _2)+2 \lambda _3) \left(
\frac{H}{\epsilon}-1\right).
\ee
Now going back to the scale-invariant case where $\lambda_1=\frac{\sigma_0}{3}-5,\lambda_2=
-2(\sigma_0-12)$ and $\lambda_3=\sigma_0-6$ one finds that the corresponding entanglement
entropy is identically zero. Indeed this may be understood as follows. Let us write the action as follows
\bea
{\cal I}=
-\frac{\kappa}{32\pi}\int d^{4}x\,\sqrt{-g}\,\left(2\sigma_0 (R_{\mu\nu}R^{\mu\nu}-\frac{1}{3}R^2)+
R^2+(\sigma-6) {\rm GB}_4\right).
\eea
Then one can find the contribution of each terms separately. The result is as follows
\bea
S&=&-\frac{\pi \kappa}{4}\bigg[4\sigma_0\left(-\frac{H}{\epsilon}+1\right)-24\left(-\frac{H}{\epsilon}+1\right)
+4(\sigma_0-6) \left(\frac{H}{\epsilon}-1\right)\bigg]\cr &&\cr
&=&-\pi \kappa\bigg[(\sigma_0-6)\left(-\frac{H}{\epsilon}+1\right)+(\sigma_0-6) \left(\frac{H}{\epsilon}-1\right)\bigg]=0.
\eea
Therefore it is zero due to the contribution of the Gauss-Bonnet term which plays the role of regulator.
Indeed the situation is the same as that in four-dimensional conformal gravity
\cite{Alishahiha:2013dca}. Actually since the Gauss-Bonnet term does not contribute to the
equations of motion the whole dynamics must be encoded in the two first terms in the
above action. Therefore if we would like to find the contributions of the dynamical degrees of freedom
to the entanglement entropy, we should only consider the dynamical parts of the action in which
the corresponding entanglement entropy is found
\be\label{SC4-Sph}
S=\pi \kappa(\sigma_0-6)\left(\frac{H}{\epsilon}-1\right),
\ee
which clearly has UV divergent term. \\
It is also illustrative to compare the contribution of dynamical part to that obtained from just
Einstein gravity. Of course for Einstein gravity the holographic entanglement entropy is given by
minimizing the area of a co-dimension two hypersurface. The corresponding entanglement entropy
is \cite{{RT:2006PRL},{RT:2006}}
\be
S_{\rm Ein}=\frac{\pi L^2}{2G} \left(\frac{H}{\epsilon}-1\right).
\ee
which is the same as \eqref{SC4-Sph} if one identifies $\kappa=\frac{L^2}{2(\sigma_0-6)G}$.
It is obvious that this identification is valid as long as $\sigma_0\neq 6$.
\subsection{HEE of a Cusp}
Finally, it is also straightforward to compute the universal term of the holographic entanglement
entropy for the case where the entangling region has a cusp. To proceed we will parametrize the
AdS metric as \eqref{AdSc}.
In this case, the entangling region is given by $t={\rm constant},\;\;-\Omega\leq \phi\leq \Omega$. Indeed such computation has already been
done for general quadratic action in \cite{Bueno:2015xda} where it was shown that the Gauss-Bonnet
term does not contribute to the universal term and therefore one arrives at
\be
S_{\rm univ}={\kappa}(\sigma_0-6)a(\Omega)\log\frac{R}{\epsilon},
\ee
where
\be
a(\Omega)=\int_0^\infty dy\left[1-\sqrt{\frac{1+h_0^2(1+y^2)}{2+h_0^2(1+y^2)}}\right],
\ee
and $h_0$ can be found from the following equation
\be
\Omega=\int_0^{h_0} dh\frac{2h^2\sqrt{1+h_0^2}}{\sqrt{1+h^2}\sqrt{(h_0^2-h^2).(h_0^2+(1+h_0^2)h^2)}}
\ee
Again the result reduces to that of the Einstein gravity if one identifies $\kappa=\frac{L^2}{2(\sigma_0-6)G}$. \\
An interesting feature of this new universal term is that in the smooth limit where $\Omega
\rightarrow \frac{\pi}{2}$ the coefficient of the universal term vanishes, though with certain form.
More precisely one has
\be
\kappa (\sigma_0-6) a(\Omega)=
\frac{\kappa(\sigma_0-6)}{4\pi}\left(\Omega-\frac{\pi}{2}\right)^2\equiv C
\left(\Omega-\frac{\pi}{2}\right)^2
,\;\;\;\;\;\;\;{\rm at}\;\Omega
\rightarrow \frac{\pi}{2}.
\ee
One can then show that $C$ is related to $C_T$ which is the central charge appearing in the
two point function of stress tensor: $C=\frac{\pi^2}{24} C_T$\cite{Bueno:2015xda}.
\section{HEE for the logarithmic solution }
The gauge/gravity duality tells us that the log term in the
AdS solution might be identified by deformed the dual CFT by an irrelevant operator. This means that such term destroy the conformal symmetry at UV
and therefore, applying the AdS/CFT correspondence should be clear. However, following \cite{Skenderis:2009nt}, let us assume that the deformation is sufficiently small in a way that this term may be treated perturbatively. Holographic entanglement entropy for logarithmic solution in Log-gravity
has been studied in \cite{Alishahiha:2013zta}. In what follows we will use the procedure
of \cite{Alishahiha:2013zta} for the action \eqref{actc}. To proceed let us set $c_0=c_3=b_3=0$ and $b_0=\beta\neq 0$ in the solution \eqref{LOG} and therefore one arrives at
\be
ds^2=\frac{L^2}{r^2}\bigg(\beta \log\left(\frac{r}{L}\right) dx_+^2-2dx_-dx_++dy^2+dr^2\bigg).
\ee
Although the constant $\beta$ can be removed by a rescaling, in order to trace the effect of logarithmic term we will keep it in the metric. Obviously setting $\beta=0$ the solution reduces
to an AdS metric. In what follows we will compute holographic entanglement entropy
for the above metric at leading order in $\beta\rightarrow 0$ limit.
We will consider an entangling region given by $x_++x_-=0$ and $\frac{\ell}{2}\leq y\leq\frac{\ell}{2}$ which is a strip with width $\ell$ along the $y$ direction. Then the co-dimension two hypersurafce
in the bulk may be parametrized by $y=f(r)$. In this case the induced metric reads
\begin{eqnarray}
\label{inducedd-dim}
ds^2=\frac{L^2}{r^2}\left((1+f'^2)dr^2+(2+\beta\log \frac{r}{L})dx_+^2 \right)\ .
\end{eqnarray}
On the other hand two unit vectors normal to the co-dimension two hypersurfaces are
\begin{eqnarray}\label{sigma1}
&&x_++x_-=\text{const.}\hspace*{1.36cm}n_1=\frac{L}{r \sqrt{2+\beta \log \frac{r}{L}}}\;(0,1,1,0)\ ,\nonumber\\
&&\rho -f(r)=\text{const.}\hspace*{1.33cm}n_2=\frac{L}{r \sqrt{1+f'^2}}\;(-f',0,0,1)\ .
\end{eqnarray}
It is then straightforward to compute the entropy functional \label{EE} for the case of
$\lambda_1=-3, \lambda_2=12$ and $\lambda_3=0$ according to the action \eqref{actc}.
Setting $F=2+\beta \log\frac{r}{L}$ one finds
\be
S_A=-\frac{\kappa}{4}\;\frac{3 H }{\sqrt{2}}\int dr \;
\frac{\sqrt{\left(f'(r)^2+1\right)F} }{2 F^2 r^2}\left(\frac{\left(2 r f''(r)F+(\beta -4 F) \left(f'(r)^3+f'(r)\right)\right)^2}{\left(f'(r)^2+1\right)^3}+12 \beta F\right),
\ee
where $H$ comes from the integration over $x_+$ and $\sqrt{2}$ is due to our notation of
light cone coordinates, $x_\pm=\frac{t\pm x}{\sqrt{2}}$. Now the aim is to minimize the above
entropy function and find the corresponding profile of $y=f(r)$ which we will do it for in
$\beta=0$ limit. Indeed at leading order one finds
\be
f'(r)=\frac{r^2}{r_t^4-r^4}+{\cal O}(\beta),
\ee
where $r_t$ is the turning point of the hypersurface. Plugging this expression into the
above entropy functional one gets
\bea
S_A&=&\frac{9 H \kappa\beta }{4}\;r_t^2\int_\epsilon^{r_t} dr
\frac{ 1}{r^2\sqrt{r_t^4-r^4}}+{\cal O}(\beta^2)\cr &&\cr
&=&\frac{9 H \kappa\beta }{4}\left(\frac{1}{\epsilon}-\frac{v_0}{r_t}\right)+{\cal O}(\beta^2).
\eea
where $v_0=\frac{\sqrt{2} \pi ^{3/2} }{\Gamma \left(\frac{1}{4}\right)^2}$. On the other hand
using the fact that $\ell=2\int_0^{r_t} dr f'+{\cal O}(\beta)=2 r_t v_0+{\cal O}(\beta)$, one arrives at
\be
S_A=\frac{9 \kappa\beta}{4}H\left(\frac{1}{\epsilon }-\frac{2 \pi \Gamma \left(\frac{3}{4}\right)^2}{\Gamma \left(\frac{1}{4}\right)^2}\;\frac{1}{\ell}\right)+{\cal O}(\beta^2).
\ee
Interestingly enough if one choose the coupling constant $\kappa$ properly,
the above expression at leading order exactly reduces to that of Einstein gravity, even though
the logarithmic solution is not an Einstein solution.
\section{Concluding Remarks: Holographic Mutual information }
Entanglement entropy as a probe of entanglement has a UV divergence, and due to the UV cut-off, the entanglement entropy can not be considered as an universal quantity. However, for a system that has two disjoint parts, one can define the mutual information. The mutual information
quantifies the amount of information which is
shared between two systems. This quantity is finite given by
\begin{align}\label{HMI}
I(A_1,A_2)=S(A_1)+S(A_2)-S(A_1\cup A_2),
\end{align}
where for a region $A$, the $S(A)$ is the entanglement entropy whereas, the entanglement entropy for the union
of the two entangling regions is $S(A_1\cup A_2)$. For three regions, the corresponding quantity is given by the tripartite information defined as follows
\begin{align}\label{3par}
{I^{[3]}}\left( {{A_1},{A_2},{A_3}} \right) &= S\left( {{A_1}} \right) + S\left( {{A_2}} \right) + S\left( {{A_3}} \right) - S\left( {{A_1} \cup {A_2}} \right) - S\left( {{A_1} \cup {A_3}} \right) \notag\\
&- S\left( {{A_2} \cup {A_3}} \right) + S\left( {{A_1} \cup {A_2} \cup {A_3}} \right).
\end{align}
Based on holographic method, the union parts of the regions play a key role, and in order to investigate this point, it is important to mention that the tripartite information can be rewritten in terms of the mutual information as follows
\begin{align}\label{3par1}
{I^{[3]}}({A_1},{A_2},{A_3}) &= I({A_1} , {A_2}) + I({A_1}, {A_3}) - I({A_1} ,{A_2} \cup {A_3}),
\end{align}
and by this modification, the aim is to compute the mutual information.
At the leading order for the logarithmic solution, the holographic entanglement entropy for a strip entanglement region one has
\be
S(\ell)=\frac{9 \kappa\beta}{4}H\left(\frac{1}{\epsilon }-\frac{2 \pi \Gamma \left(\frac{3}{4}\right)^2}{\Gamma \left(\frac{1}{4}\right)^2}\;\frac{1}{\ell}\right).
\ee
In order to compute the mutual information, we should consider two entangling regions; we assume two strips with the same length of $\ell$ separated by distance $h$. After doing some calculation, one obtains
\begin{equation}\label{mutualcon}
I =\frac{9\kappa \beta}{4}H
\Big(\frac{1}{2\ell+h}+\frac{1}{h}-\frac{2}{\ell}\Big).
\end{equation}
\begin{figure}[h!]
\centering \includegraphics[width=15cm]{a2.png}
\caption{ For three subsystem, in reading the holographic entanglement entropy of the union part, one must use the minimal configuration. In this figure, we have plotted the configurations of the two and three regions. }
\label{imutual1}
\end{figure}
Similarly, for three strips with the same length $\ell$ separated by distance $h$, as mentioned, the key point is to identify the union part of entanglement entropies, noting that the holographic principle forces us to find a minimal configuration in the bulk. For three strips, we have plotted all the union parts in Fig.\ref{imutual1}. In each case, in computing the $S(A_i\cup A_j)$ and $S\left( {{A_1}\cup{A_2}\cup{A_3}} \right)$ one must find the minimum among the possible diagrams. This needs a numerical computation and the holographic tripartite information is given by
\begin{equation}
{I^{\left[ 3 \right]}}\left( {A,B,C} \right) = 3S\left( \ell \right)- \min \left\{ {{S_1},{S_3}} \right\} - 2\min \left\{ {{S_1},{S_2}} \right\} + \min \left\{ {{S_4},{S_5},{S_6},{S_7}} \right\},
\end{equation}
noting that $\min \left\{ {{S_1},{S_2}} \right\}$ means that the minimum configuration between $S_1$ and $S_2$ must be used. Note that we have used the following simplifications: \[\begin{array}{l}
S\left( {{A_i} \cup {A_j}} \right):\left\{ \begin{array}{l}
2S\left( \ell \right) \equiv {S_1} \\
S\left( {2\ell + h} \right) + S\left( h \right) \equiv {S_2} \\
S\left( {3\ell + 2h} \right) + S\left( {\ell + 2h} \right) \equiv {S_3} \\
\end{array} \right.
\end{array}\]
and for three entangling regions, one has
\[\begin{array}{l}
S\left( {{A_1} \cup{A_2} \cup{A_3}} \right):\left\{ \begin{array}{l}
3S\left( \ell \right) \equiv {S_4} \\
S\left( {3\ell + 2h} \right) + S\left( {\ell + 2h} \right) + S\left( \ell \right) \equiv {S_5} \\
S\left( {2\ell + h} \right) + S\left( \ell \right) + S\left( h \right) \equiv {S_6} \\
S\left( {3\ell + 2h} \right) + 2S\left( h \right) \equiv {S_7} \\
\end{array} \right. \\
\end{array}\]
We have used Mathematica software to compute the numeric analysis of finding the minimal configuration and figure \ref{WDWP} shows the numerical results for the mutual and tripartite information for some certain values of $\ell$ and $h$.
Our numerical analysis shows that the mutual information is always positive, whereas the tripartite information becomes negative, at least this happens for the range of parameters that we have used. Previously. we have observed this feature for some specific gravity theories as well \cite{Mirabi:2016elb, Tanhayi:2016uui, Tanhayi:2017wcd, Ng }.
\begin{figure}[h!]
\begin{center}
\includegraphics[scale=0.45]{2dm}\,\,
\includegraphics[scale=0.43]{2dt}
\end{center}
\caption{ Numerical results for holographic mutual information (left plot) and tripartite information (right plot) as a function of the separation distance: for ${\ell}=1,\cdots, 5$. Note in all cases the mutual information is positive while for three regions the tripartite information remains negative.
}
\label{WDWP}
\end{figure}
\begin{figure}[h!]
\begin{center}
\includegraphics[scale=0.48]{3dt}\,\,\,\,\,\,\,\,\,\,
\includegraphics[scale=0.48]{3dm}
\end{center}
\caption{ The 3D plot of holographic mutual information and tripartite information as a function of the separation distance and the length of the strips.
}
\label{WDWP}
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%\section{Discussions}
In this paper, we have investigated some features of the scale-invariant gravity theory in four dimensions by computing the HEE for the logarithmic solution. Our numerical results indicate that in a such theory, the measure of entanglement for two systems is always positive. On the other hand for three regions the corresponding measure becomes negative. The latter quantity is named tripartite information, and the negativity of tripartite information imposes a condition on the mutual information. This statement can be understood from the relation \eqref{3par1}: when ${I^{[3]}}({A_1},{A_2},{A_3})$ is negative one receives $$I({A_1} ,{A_2} \cup {A_3})\ge I({A_1} , {A_2}) + I({A_1}, {A_3}). $$
It is worth mentioning that in a generic case in quantum field theory, the above inequality does not hold, and this inequality becomes very important in the context of quantum information theory. This indicates a well-known feature of quantum correlation, namely the monogamy of mutual information. This characteristic of measures of quantum entanglement completely comes from the nature of quantum mechanics. In other words, correlations in the holographic regimes arise from quantum entanglement rather than classical correlations.