\section{Introduction}
\label{intro}
Entanglement entropy of a subregion in a conformal field theory is conjectured to be found holographically as the surface of the extension of subregion into the bulk (i.e. RT surface) provided that its variation vanishes \cite{Ryu:2006bv,Ryu:2006ef}. This quantity diverges near the boundary. Most of the time people are interested in the behavior of the divergent terms and their relationship with physical properties of the system. Especially the logarithmic singularity may be interpreted as some cut-off independent characteristics of the system. For instance, in three dimensions the coefficient of  logarithmic term has been shown to be proportional to the central charge of the  corresponding CFT at the boundary \cite{Casini.2007,Casini.2009,Takayanagi.2007, Moore.2006}. Such relationship are also aimed for higher dimensions in subsequent works \cite{Solodukhin:2008cdf}.

There is also some interest in the entanglement entropy of singular subregions  in recent years \cite{Bueno:2015rda,Bueno:2015xda,mollabash12,mollabash13}. The effect of singularities gained some attractions, first in three dimensions \cite{Takayanagi.2007} and then in higher dimensions \cite{Singh:2012ala}. The most interest is to find new singular terms in the expansion of the entanglement entropy due to singularity of the considered subregion \cite{Myers:2012vs}.

Another quantity of great interest is the complexity of a QFT living in the boundary. Besides the well-known approaches of $ complexity = volume $ \cite{Stanford:2014jda, Susskind:2014rva,Susskind:2014moa,Couch:2016exn} and $ comlexity =action $ \cite{Carmi:2016wjl, Lehner:2016vdi,Brown:2015bva,Brown:2015lvg}, there is also another conjecture to  consider the complexity as the volume enclosed by the Rio-Takayanagi surface \cite{Alishahiha:2015rta,Ben-Ami:2016qex}.This approach is known as subregion complexity.
In a recent paper \cite{Bakhshae:2017} we used this approach to study the subregion complexity of a number of singular surfaces, focusing on the singularities in terms of UV cutoff parameter.

Recently some interests have been arised to deform the CFT by imposing a relevant operator \cite{parvizi1,parvizi3,parvizi4,parvizi5,parvizi6,parvizi7}.
Using the standard AdS/CFT correspondence \cite{myers11}, this may be achieved, for example, by turning on a scalar field in the bulk. Hence, the geometry of the bulk is no longer pure AdS; however, it turns out to be asymptotically AdS near the boundary. It has been shown that this deformation may lead to appearing universal logarithmic correction in the entanglement entropy for definite values of the conformal dimension of the relevant operator \cite{myers}.
The effect of deformation of the CFT on the entanglement entropy is
investigated for regular subregions such as sphere in Refs. \cite{parvizi6,myers}. For singular subregions, this effect has  been recently studied only for kink in three dimensions in \cite{prvizi}. Some general aspects of the complexity of the deformed theories are studied in \cite{alishahiha11}.

In this paper we are mostly interested in calculating the entanglement entropy and complexity for a set of singular surfaces (including cones and creases) in a deformed conformal field theory.
After a brief review of a deformed conformal field theory in the next section, we calculate the entanglement entropy for singular subregions of a kink in $d=3$ and cones in $d=4, 5$ and $6$ for a deformed theory. We also consider creases in  $d=4, \ 5$ and $6$. This is done in section~\ref{Entanglement Entropies}. The essential calculations and technical points are given in more details for the case of kink and to some extent for cone $c_1$. For other singular surfaces we give only the important results.
In section~\ref{Complexity} we give our results for the subregion complexities of the same singular submanifolds. We discuss about our results in section~\ref{Discussions}.

It is important to distinguish among different types of singularities. The most familiar kind is that of ordinary UV divergences of the entanglement entropy, as well as the complexity, when we approach the boundary in the AdS/CFT framework. The next kind corresponds to geometrical singularities near the needle points or wedges of singular subregion. Finally we encounter new singular terms due to deforming a theory by a relevant operator.  Among different singular terms people are mostly interested in the logarithmic singular terms because of their universal characteristics due to independence of the regularization process.
\label{Introduction}

\section{Deformed CFT}
\label{Deformed CFT}
Consider a CFT living in the $d$-dimensional boundary of a $(d+1)$-dimensional space-time. As is well-known \cite{myers}, turning on a scalar field in the bulk, one can deform  the CFT by a relevant operator. To do this, the Hilbert-Einstein action with a negative cosmological constant $\Lambda=-\dfrac{d(d-1)}{L^2}$ is perturbed as follows
\begin{align}
I=\dfrac{1}{16\pi G_N}\int d^{d+1}x\sqrt{g}\left[R-\dfrac{d(d-1)}{L^2}\right]-\frac{1}{2}\int d^{d+1}x\sqrt{g}\left[(\partial \Phi)^2+M^2 \Phi ^2\right].
\label{a50}
\end{align}
The mass parameter $ M$ determines the conformal dimention of the boundary operator $O$ (dual to $\Phi$ ) as
\begin{align}
\Delta _\pm=\frac{d}{2}\pm \sqrt{(ML)^2+\frac{d^2}{4}},
\label{a51}
\end{align}
where $\Delta _+$ and $\Delta_-$ are valid for $ -\frac{d^2}{4}\leq (ML)^2 \leq-\frac{d^2}{4}+1 $;  while  for $ -\frac{d^2}{4}+1< (ML)^2 $ only $\Delta _+$ is valid \cite{nishi1}. Due to this deformation the AdS solution of the Einstein equation also deforms as follows
\begin{align}
ds^2=\dfrac{L^2}{z^2}\left[\dfrac{dz^2}{f(z)}-dt^2+d\rho ^2+ \rho^ 2(d\theta ^2+ \sin ^2\theta d\Omega _n ^2)\right],
\label{a52}
\end{align}
where $f(z)$ is the deformation function. The geometry is asymptotically AdS, i.e. $ f(z)\rightarrow1$ as $z\rightarrow0$. Considering the equations of motion for metric components as well as the scalar field $\Phi$, one can show directly that the expansion of $ f(z)$ near the boundary is as follows
\begin{equation}
f(z)=\left\{ \begin{array}{ll}
	1+ (\mu z)^{2\alpha}+..., \quad &\Delta \neq d/2\nonumber\\
( \mu z)^{d}(\log \mu z)^2+..., \quad &\Delta = d/2 \end{array} \right.
\end{equation}
where $\mu $ is a mass parameter determined from the parameter $\lambda$ of coupling  the relevant operator, and
$\alpha = d-\Delta_+=\Delta_-$ for  both $\Delta=\Delta_+$ and $ \Delta=\Delta_-$.

Deforming the bulk geometry leads to some changes in the entanglement entropy of the subregions as stated above. It is shown \cite{parvizi6} that to lowest order of $\mu$, the variation of the entanglement entropy for the sphere in a $d$-dimensional CFT  and $ \Delta \neq d/2$ reads
 \begin{align}
 \delta S=-\dfrac{\Gamma(\dfrac{d+1}{2})\Gamma(\dfrac{2-d+2 \alpha}{2})}{4\Gamma(\frac{3}{2}+\alpha)}K(\mu R)^{2\alpha}+\dfrac{K(\mu R)^{d-2}}{2(2-d+2 \alpha)}(\mu\delta)^{2-d+2\alpha}+\mathcal{O}(\delta ^{4-d+2\alpha}),
  \end{align}
 where $K=L^{d-1} \textstyle{Vol} (S^{d-2}) /4 G_N $, $\delta$ is the UV cutoff and $R$ is the radius of the sphere. The first term is the finite result while the subsequent terms are singular corrections. They are a limited number of terms which are present unless $2n-d+2\alpha \le 0$ where equality leads to logarithmic singularity. For $\Delta=\Delta_-$ where $d/2-1<\alpha<d/2$, we find $n<1$ which gives no contribution for $\alpha \ne d/2$.  For $\Delta=\Delta_+$ we have $d/2<\Delta_+<d$ which leads to $0<\alpha<d/2$. Hence, for $d=3$ we have a singular term as $\delta^{-1+2\alpha}$ for $\alpha<1/2$ a logarithmic singularity for $\alpha=1/2$ and no singular term for $\alpha>1/2$. For $d=4, 5$ and $6$  the results are consistent with what we will find in the next section (see the table~\ref{tab1} ).
 \section{Entanglement Entropies}
\label{Entanglement Entropies}
Let us rewrite the solution \eqref{a52} of the equations of motion of the action \eqref{a50} in the following form
\begin{align}
	ds^2=\frac{L^2}{z^2}\left[ - dt^2+\frac{dz^2}{f(z)}+d\rho ^2+\rho ^2 (d\theta ^2+ \sin^2 \theta d \Omega _n ^2) +\sum _{i=1}^{i=l}d x_i ^2\right],\label{a0}
\end{align}
In this form we have divided the flat $(d-1)$-dimensional manifold of the boundary at fixed time into a $(n+2)$ flat space described by spherical coordinates and a $l$-dimensional space described by Cartesian coordinates, such that $d=l+n+3$. This enables us to introduce the geometrical singularity by assuming $ \rho\rightarrow 0 $  when the angle $\theta$ is limited to the interval $[-\Omega,\Omega]$ for kink (i.e. $n=0$ and $l=0$) and the interval $[0,\Omega]$ for cones (i.e. $n \ge 1$ and $l=0$). Creases correspond to extensions $l \ge 1$. All of our singular subregions are restricted to $\rho \le H$ where $H$ is the IR cutoff.
Assuming $\rho=\rho(z, \theta)$ to describe the RT surface, the induced metric reads
\begin{equation}
h=\left(\begin{array}{c c c c c c}
\dfrac{L^2}{z^2}(\frac{1}{f(z)}+\rho '^2)&\dfrac{L^2}{z^2}\rho' \dot{\rho} &&&&\\
\dfrac{L^2}{z^2}\rho' \dot{\rho}&\dfrac{L^2}{z^2}(\rho ^2+\dot{\rho}^2)&&&&\\
&&\dfrac{L^2\rho ^2 \sin ^2(\theta)}{z^2 }g_{ab} (S^n)&&&\\
&&&\dfrac{L^2}{z^2}&&\\
&&&&\ddots &\\
&&&&&\dfrac{L^2}{z^2} \\
\end{array}\right).
\label{ab11}
\end{equation}
According to RT prescription, the entanglement entropy is proportional to the minimized  area of the RT surface as
\begin{equation}
S=\frac{2\pi}{l_p^{d-1}} L^{d-1} \Omega_n  \int dz  d\theta\,\dfrac{\rho^{n}}{z^{d-1}}\,\sin^n \theta\,\sqrt{{\rho ' }^{2}\rho ^2+\frac{\rho ^2 +\dot{\rho}^2}{f(z)}}, \label{ab10}
\end{equation}
where $ \Omega_n$ is the volume of the unit $n$-sphere,  $\dot{z}=\partial_ {\theta} z $ and $ z'=\partial _{\rho} z $. In the following subsections we  calculate the  entanglement entropy for different singular subregions.
\subsection{kink $ k $}
\label{kink }
 The entanglement entropy \eqref{ab10} for kink is given by the following integral
\begin{align}
S\vert_{k}=\frac{2\pi L^2}{l_p ^2}\int d z d\theta L(\rho, \dot{\rho}, \rho'), \label{a'4}
\end{align}
where
\begin{align}
L(\rho, \dot{\rho}, \rho')= \dfrac{\sqrt{\frac{\dot{\rho}^2+\rho ^2}{f(z)}+{\rho '}^2 \rho^2 }}
{z^2}.\label{a4}
\end{align}
 Minimizing the above integral gives the equation of motion of $ \rho(z,\theta) $ as follows
\begin{align}
 2 f z \rho(\rho ^2 +\dot{\rho}^2){\rho} '' +2 z \rho(1+f{\rho ' }^2)\ddot{\rho}-4 f z \rho \rho ' \dot{\rho}\dot{\rho}'+z \rho \rho'(\rho ^2 +\dot{\rho} ^2) f '  & \nonumber\\
-2z\left((1+f{\rho '}^2)\rho ^2 +2\dot{\rho}^2\right)-4 f \rho \rho '\left((1+f {\rho '}^2)\rho ^2 +\dot{\rho}^2\right) & =0. \label{ab2}
\end{align}
This partial differential equation should be solved according to the following boundary conditions
\begin{align}
& \frac{\partial\rho}{\partial \theta }(z,0)=0,  \label{d2}\\ & \rho(z, \Omega) =H. \label{d3}
\end{align}
Eq. \eqref{d2} shows that at arbitrary $z$ the coordinate $\rho$ on RT surface acquires its minimum value at $\theta=0$, while Eq. \eqref{d3} shows that at the limiting points $\theta= \pm \Omega$ the RT surface touches the boundary $\rho=H$. See figure~\ref{fig1}
 to get a geometrical feeling about the problem.
To solve Eq. \eqref{ab2} we consider a perturbative ensatz where in the limit $ \mu \rightarrow 0 $, we assume $\rho=z/h(\theta)$. This leads to the following equation for $h(\theta)$
\begin{align}
2+3 h^2+h^4+2 \dot{h}^2+h(1+h^2)\ddot{h}=0,\label{a6}
\end{align}
which give the following constant
\begin{align}
K_3=\dfrac{(1+h^2)}{h^2\sqrt{1+h^2+\dot{h}^2}}.\label{a10}
\end{align}

\begin{figure}[tbp]
\centering
\includegraphics[width=.45\textwidth]{photo2.pdf}
\caption{\label{fig1}Schematic figure of the RT extension of the kink (green region) into the bulk \cite{Myers:2012vs}. }
\end{figure}

The boundary condition \eqref{d2} gives $\dot{h}(0)=0$ at the turning point $\theta=0$. Hence, the constant $K_3$ can be written in terms of $h_0=h(0)$. In principal the Eq. \eqref{a6} can be solved to find the function $h(\theta)$. However, we use this equation to find $\dot{h}$ as a function of $h$. Inserting these results in Eq. \eqref{ab2}, one can find the entanglement entropy of kink \cite{Myers:2012vs}. Then we can complete our enzatz by considering the following expansion
\begin{align}
\rho (z,\theta)=\frac{z}{h(\theta)}+\mu^{2\alpha} z^{2\alpha+1} g_2(\theta)+\cdots .\label{a5}
\end{align}
Inserting the complete ansatz \eqref{a5} into the equation of motion\eqref{ab2} we find (in addition to Eq. \eqref{a6}) the following equations for $g_2$,
\begin{align}
& 2 h^3(1+h^2)(1+h^2 +\dot{h}^2)\ddot{g}_2+2 h^2 \dot{h}(2(5 +2 h^2 +2 \alpha)\dot{h}^2+(1+h^2)(10+5 h^2+4 \alpha+h \ddot{h}))\dot{g}_2\nonumber\\
&-2h (-2(-2+h^2 +\alpha+2 \alpha ^2)\dot{h}^4+\dot{h}^2(22 +h^4 +2 (7 -2\alpha)\alpha+ h^2(19 +4\alpha- 8 \alpha ^2)\nonumber\\
&+2 h(2 + h^2 +2 \alpha)\ddot{h})+(1 + h^2)(2 (9 + 8 \alpha) + h^2 (19 + 4 h^2 - 4 (-2 + \alpha) \alpha)  \nonumber\\
&+3 h (2 + h^2 + 2 \alpha) \ddot{h})) g_2+(2 (-1 + \alpha) \dot{h}^4 + (1 + h^2) (-8 + h^4 +
 2 h^2 (-2 + \alpha) \nonumber\\
&+ h (-2 + h^2)\ddot{h})+\dot{h}^2(h^4 + 2 (-5 + \alpha) + h^2 (-5 + 4 \alpha) + h (-1 + h^2)\ddot{h}))=0.\label{ab1}
\end{align}
As is expected, for the generic problem of the entanglement entropy, we will find UV divergences in the limit $z \rightarrow 0$ or equivalently as $\rho \rightarrow 0$ (for arbitrary $-\Omega <\theta <\Omega$) or $\theta \rightarrow \Omega$ (for finite $\rho$). In terms of the variable $h$ the latter limit is equivalent to $h \rightarrow 0$. The UV divergent terms of the entanglement entropy  originate from the divergences of the integrand as well as the limits of the integral  Eq. \eqref{a'4} as
\begin{align}
	S\vert_{k}=&-\frac{2\pi L^2}{l_p ^2}\int_{z_m} ^{\delta }dz \int_{-\Omega+\epsilon
	}^{\Omega -\epsilon} d\theta L(\rho, \dot{\rho}, \rho'),
\end{align}
which upon changing the integration limits of $\theta$ from $(-\Omega,\Omega)$ to $(0, \Omega)$ and replacing $d\theta$ in Eq.  \eqref{a'4}
by $dh/\dot{h}$,  reads
\begin{align}
S\vert_{k}=&-\frac{4\pi L^2}{l_p ^2} \int_{z_m} ^{\delta }\frac{dz}{z} \int_ {h_0}^{h_{1c}} \frac{dh}{\dot{h}} L(\rho, \dot{h}d\rho/dh, \rho') . \label{aa8}
\end{align}
Concerning the integral bounds in Eq. \eqref{aa8},  from  the boundary condition \eqref{d3}  we have $\rho (z, \Omega -\epsilon(z)) = H $; hence we can define $ h_{1c}(z) = h(\Omega - \epsilon(z))$. In the limit $z=\delta \rightarrow 0$ we have
$ h_{1c}(\delta) = h(\Omega - \epsilon(\delta))$. The limit $z_m$ is also achieved via the condition $\rho(z_m, 0) = H $. Inserting $\rho$ from Eq. \eqref{a5} into \eqref{aa8} gives the following expansion for the entanglement entropy with respect to $\mu$
\begin{align}
S\vert_{k}=-\frac{4\pi L^2}{l_p ^2}(I_0 +\mu ^{2\alpha} I_1+ \cdots ), \label{a9}
\end{align}
where
\begin{align}
&I_k=\int_{z_m} ^{\delta }dz z^{(2\alpha)k-1} \int_ {h_0}^{h_{1c}}dh G_k(h),\label{a11}
\end{align}
in which
\begin{align}
&G_0(h)=\dfrac{\sqrt{1+h^2 +\dot{h}^2}}{\dot{h} h^2},\label{aa10}\\
&G_1(h)=\dfrac{ -\dot{h}^2-h^2(1+2 \dot{h}\dot{g}_2)+4 (1+\alpha) h g_2+2 h^3 g_2}{2\dot{h} h^2 \sqrt{1+h^2+\dot{h}^2}}\, .\label{aa11}
\end{align}
Now to calculate the entanglement entropy one needs to solve Eq. \eqref{ab1}  for $ g_2 $. This is a difficult task. However, since we are mostly interested to find the UV divergent terms of the entanglement entropy, we just need to find the asymptotic behavior the corresponding function $g_2$ near the boundary, i.e. in the limit $h\rightarrow 0$. According to Eq. \eqref{a5} , in order to keep $\rho$ finite in the limit $ h \rightarrow 0 $ and $ \delta\rightarrow 0 $, the most singular term of $g_2$ should be of order $h^{-2\alpha-1}$. Assuming $g_2=a h^{-2\alpha-1}+a' h^{-2\alpha}+\cdots$ and inserting it in the corresponding equations \eqref{ab1} (as well as $\dot{g_2} =\dot{h}dg_2/dh$ etc.) and using Eq. \eqref{a10} for $\dot{h}$ we find the following results
\begin{align}
g_2=\dfrac{a_1}{h^{2\alpha +1}}+\dfrac{a_2}{h^{2\alpha-1}}+\dfrac{a_3}{h^{2\alpha-3}}+\frac{a_4}{ h}+a_5 h+\mathcal{O}(h^3),\label{a13}
\end{align}
where
\begin{align}
&a_2=\frac{a_1}{5} (2 + 5 \alpha - 2 \alpha^2),\\
&a_3=\frac{a_1}{70} (-12 + 10 K_3^2 - 7 \alpha+ 43 \alpha^2 - 28 \alpha ^3 + 4 \alpha^4),\\
&a_4=-\frac{1}{(6 + 4 \alpha)}\ ,\\
&a_5= \frac{2}{15 + 31 \alpha + 20 \alpha^2 + 4 \alpha^3}\ .
\end{align}
The constant $a_1$ would be fixed from the boundary conditions \eqref{d2} and \eqref{d3}.
Now we need to determine the limiting value $h_{1c}$ to find the divergencies of the integrals \eqref{a11}. This can be down from the same boundary condition $\rho(z, \Omega)=H$ in the limit $ z=\delta$, i.e. $\rho (\delta, \Omega-\epsilon)=H$ where $h_{1c} = h(\Omega-\epsilon)$. Using Eqs. \eqref{a10} and  \eqref{a13} we can find the expansion of  $ h_{1c} $ in terms of $ \delta $ as follows
\begin{align}
h_{1c}(\delta)=&\left( \frac{1}{H}+a_1 \mu ^{2\alpha} H^{2\alpha -1}\right)\delta +a_2\mu ^{2\alpha} H^{2\alpha -3}\delta ^3\nonumber\\ &+a_3 \mu ^{2\alpha} H^{2\alpha-5}\delta ^5+ \mu ^{2\alpha}\frac{a_4}{H}\delta ^{2\alpha+1} +\mu ^{2\alpha} \frac{a_5}{H^3}\delta ^{2\alpha+3}
+\mathcal{O}(\delta ^{2\alpha+5}).\label{a155}
\end{align}
From the Eqs. \eqref{a10} and \eqref{a13} we can find the following expansions for $G_0$ and $G_1$  in terms of the UV cut-off parameter $\delta$
{
\begin{align}
&G_0(h)\sim- \frac{1}{h^2}-\dfrac{k_3 ^2 h^2}{2}+\mathcal{O}(h^4)\label{m1},\\
&G_1(h)\sim -\frac{a_1(1+2\alpha)}{h^{2\alpha+2}}+\frac{2+\alpha}{(3 +2 \alpha) h^2}+\mathcal{O}(h^{-2\alpha}).\label{m2}
\end{align}
Separating the singular part of the integrals $ I_k $, we can divide them as $I_k=I_k'+I_k'' $, where the integrand of $I'_k$ is regular. Hence, we have
\begin{align}
I_0=&I_0'+I_0''\nonumber\\
=&\int_{z_m} ^{\delta}\frac{dz}{z}  \int _{h_0}^{h_{1c}} dh\left(G_0(h)+\frac{1}{h^2}\right)-\int_{z_m} ^{\delta}\frac{dz}{z}  \int _{h_0}^{h_{1c}} dh\frac{1}{h^2}\nonumber\\
=&\int_{z_m} ^{\delta}\frac{dz}{z}  \int _{h_0}^{h_{1c}} dh\left(G_0(h)+\frac{1}{h^2}\right)+\int_{z_m} ^{\delta}\frac{dz}{z}\left(\frac{1}{h_{1c}}-\frac{1}{h_0}\right),
\label{aa15}
\end{align}
\begin{align}
I_1&=I_1'+I_1''\nonumber\\
=&\int_{z_m} ^{\delta}z^{2\alpha-1} dz\int _{h_0}^{h_{1c}} dh\left(G_1(h)+\frac{a_1(1+2\alpha)}{h^{2\alpha+2}}-\frac{2+\alpha}{(3 +2 \alpha) h^2}\right)\nonumber\\
&-\int_{z_m} ^{\delta} z^{2\alpha-1}dz\int _{h_0}^{h_{1c}}dh\left(\frac{a_1(1+2\alpha)}{h^{2\alpha+2}}-\frac{2+\alpha}{(3 +2 \alpha) h^2}\right) .
\label{a16}
\end{align}
 Note that for instance in the second term of Eq. \eqref{a16} due to the term $  z^2 dz $, we just need to consider the first two terms of the expansion of $ G_1 $ in Eq. \eqref{m2}. The same task is done in the other calculations. In order to find the singular behavior in each case, let take the derivative of $I_k$'s with respect to $ \delta $, i.e.
\begin{align}
\frac{d I'_0}{d\delta}=&\frac{1}{\delta } \int _{h_0}^{h_{1c}} dh\left(G_0(h)+\frac{1}{h^2}\right)\nonumber\\
&=\frac{1}{\delta } \int _{h_0}^{0} dh\left(G_0(h)+\frac{1}{h^2}\right)+\frac{d h_{1c}}{d\delta}\left(G_0(h)+\frac{1}{h^2}\right)\vert _{h=h{1c}}+..\nonumber\\
&=\frac{1}{\delta } \int _{h_0}^{0} dh\left(G_0(h)+\frac{1}{h^2}\right)+\frac{d h_{1c}}{d\delta}\left(-\frac{k_3 ^2 h_{1c}^2}{2}\right)\nonumber\\
&=\frac{1}{\delta } \int _{h_0}^{0} dh\left(G_0(h)+\frac{1}{h^2}\right)+\mathcal{O}(\delta ^2),\label{a17}
\end{align}
\begin{align}
&\frac{d I_0''}{d\delta}=\frac{H}{\delta ^2}\left(1-a_1 \mu ^{2\alpha} H^{2\alpha }-\mu ^{2\alpha} a_4\delta ^{2\alpha}\right)-\frac{1}{h_0 \delta }+\mathcal{O}(\delta ^0),\label{a18}\\
&\frac{d I_1'}{d\delta}=\delta ^{2\alpha-1} \int _{h_0}^{0} dh\left(G_1(h)+\frac{a_1(1+2\alpha)}{h^{2\alpha+2}}-\frac{2+\alpha}{(3 +2 \alpha) h^2}\right)+\mathcal{O}(\delta ^0 ),\label{a19}\\
&\frac{d I_1''}{d\delta}=\dfrac{a_1 H^{1 + 2 \alpha}}{\delta^2} -\dfrac{ H (2 + \alpha)}{(3 + 2 \alpha) }\delta ^{2\alpha-2}+\mathcal{O}( \delta ^0).\label{a20}
\end{align}
Integrating with respect to $ \delta $ we can find the entanglement entropy for kink in a deformed CFT as
\begin{align}
S\vert _{k}=S_{(0,k)} +S_{(1,k)},
\end{align}
where
\begin{align}
S_{(0,k)}=-\frac{4\pi L^2}{l_p ^2}\left[\log\left(\frac{\delta}{H}\right)\left(\int _{h_0}^{0} dh\big(\dfrac{\sqrt{1+h^2+\dot{h}^2}}{ \dot{h}h^2}+\frac{1}{h^2})-\frac{1}{h_0}\right)-\dfrac{H}{\delta}\right]
\end{align}
is the he entanglement entropy of kink for the pure AdS case, consistent with the existing results \cite{Myers:2012vs}, and
\begin{equation}
S_{(1,k)}=\left\{ \begin{array}{ll}
	\frac{4\pi L^2}{l_p ^2}\dfrac{\mu^{2\alpha} H\delta ^{2\alpha-1}}{2(2 \alpha-1)} \quad & 0 <\alpha <\frac{1}{2}\nonumber\\
\frac{4\pi L^2}{l_p ^2}\dfrac{\mu^{2\alpha} H}{2} \log\left(\mu\delta\right) \quad &\alpha=\frac{1}{2}\nonumber\\
0\quad &  \alpha>\frac{1}{2},\alpha\neq\frac{3}{2} \end{array} \right.
\end{equation}
is the first order correction of the entanglement entropy due to deformation with a relevant operator. This result also is consistent with Ref. \cite{prvizi} (see table~\ref{tab1}). The important point is that here we have a new universal logarithmic correction for the case $\alpha=1/2$.  Considering the definition of $\alpha$ in the previous section, this shows that for a special tuning $(ML)^2=-5/2$ we have new logarithmic term in entanglement entropy. On the other hand we have no new UV correction for $\alpha >1/2$.

\subsection{cone $ c_n $}
\label{cone}
In this subsection, we give some details for the case $n=1$ corresponding to $d=4$, while for $n=2$ and $n=3$ we give only the results. The entanglement entropy of the cone $c_1$ is achieved
by minimizing the RT surface which may be formulated by  $ \rho =\rho(z, \theta) $ as follows
\begin{align}
S\vert_{c_1}=\frac{2\pi \Omega _1 L^3}{l_{p }^3}\int  d z d\theta\dfrac{\rho\sin(\theta)\sqrt{\frac{\dot{\rho}^2+\rho ^2}{f(z)}+{\rho '}^2 \rho^2 }}
{z^3},\label{b3}
\end{align}
where $ \dot{\rho}=\partial _{\theta} \rho $ and $ \rho '=\partial _{z} \rho $. For this reason we should solve the following equation of motion
\begin{align}
-6 z \sin(\theta) \rho \dot{\rho}^2
+ 2 z \cos(\theta) \dot{\rho}^3
 - 4 z \sin(\theta) \rho^3 (1 + f{ \rho '}^2)+  \sin(\theta) \rho^4 \big(z f' \rho'
- 6 f^2 {\rho '}^3 + f (-6 \rho '  &\nonumber\\ +2 z \rho '')\big)
+\rho ^2 \big(2 z \sin(\theta)\ddot{\rho} (1 + f {\rho ' }^2)+2 z \dot{\rho}
\left(\cos(\theta) + f \rho ' (\cos(\theta) \rho '
 - 2 \sin(\theta)\dot{\rho  }')\right) & \nonumber\\
 +\sin(\theta)\dot{\rho}^2\left(z f' \rho '
 +f(-6 \rho '+ 2 z\rho'')\right)\big)& =0.
\label{b4}
\end{align}
Using the general ensatz \eqref{a5}, we can find iterative equations for the unknown functions $h(\theta)$, $g_2(\theta)$. So the entropy functional $S\vert_{c_1}$ would be written peturbatively as
\begin{align}
S\vert_{c_1}=&-\frac{4\pi ^2 L^3}{l_p ^3}(I_0 +\mu^{2\alpha} I_1+\cdots),\label{b6}
\end{align}
where
\begin{align}
&I_k=\int_{z_m} ^{\delta }dz z^{2 \alpha k-1} \int_ {h_0}^{h_{1c}}dh G_k(h),\label{bb6}
\end{align}
in which
\begin{align}
&G_0(h)=\dfrac{\sin(\theta)\sqrt{1+h^2 +\dot{h}^2}}{\dot{h}h^3},\label{d4}\\
&G_1(h)=\dfrac{\sin(\theta)\left(-\dot{h}^2 -h^2(1+2 \dot{h}\dot{g}_2)+4 h^3 g_2 +2 h (3 +2 \alpha +\dot{h}^2) g_2\right)}{ 2 h^3 \dot{h}\sqrt{1+h^2 +\dot{h}^2}}, \label{d5}
\end{align}
where as before we have changed the variable $\theta$ to $h$.
Now let us insert the ensatz \eqref{a5} in the equation of motion \eqref{b4}. Assuming $ y =\sin(\theta) $, we have $ \dot{h}(y) = \sqrt{1-y^2}/y' $ and $ \ddot{h}(y) = -((1 - y^2)y''+
y{y'}^2)/{y'}^3 $ where $ y' = dy/dh $ and $y''=d^2y/dh^2$. Hence the zeroth order part of the equation of motion reads
\begin{align}
&(3+h^2) y (-1+y^2) y'+h(1+h^2)(-1+2 y^2){ y'}^2-(3 +5 h^2 +2 h^4) y {y'}^3 \nonumber\\
&-h (-1 +y^2)\left(-1 +y^2+(1+h^2) y y'' \right)=0.
\label{b7}
\end{align}
In the next order we can use $\dot{g_2}=g_2' \dot{h}$ and $\ddot{g_2}=g_2''\dot{h}^2 +g_2' \ddot{h}$, and relations of $\dot{h}$ and $\ddot{h}$ in terms of $y$, $y'$ and $y''$, to find the following equation
\begin{align}
&h g_2 \big(-2 h (-1 + y^2)^2 (2 + h^2 + 2 \alpha) +
    y (-1 + y^2) \left(5 + 3 h^4 + 12 \alpha+ 4 \alpha^2 +
       4 h^2 (1 + \alpha + \alpha^2)\right)\nonumber\\& \times y' +(1 + h^2)  y
      \left(9 + 2 h^4 + 6 \alpha +
       h^2 (11 + 4 \alpha - 4 \alpha^2)\right) {y'}^3\big) -
 h (-1 + y^2)^2 \left(1 +
    2 (h^2 + h^4) g_2'\right) \nonumber\\& + (1 +
    h^2)y
     \times {y'}^3\left(3 - h^2 (-3 + \alpha) +
    h^2 (3 + 5 h^2 + 2 h^4) g_2 '\right) +
 y (-1 + y^2) y'\big(-2 h^2 + \alpha + h^2 \nonumber\\& \times \alpha+
    h^2 (1 + h^2) (7 + 5 h^2 + 4 \alpha)
   g_2'+ h^3 (1 + h^2)^2 g_2''\big)=0.
\label{b8}
\end{align}
 The Eq. \eqref{b7} should be solved with the initial conditions  $ y = sin(\Omega) $ at $ h = 0 $ and $y'=0$ at $h=h_0$ where $h_0=h(\theta=0)$ is the value of $h$ at the turning point. This equation can not be solved exactly. However, we are only interested in the behavior of the quantities near the boundary $h=0$. Hence, we only need to know the expansion of $y(h)$ around $h=0$ which satisfy Eq. \eqref{b7} and the boundary conditions. The result turns out to be as follows
\begin{align}
y =\sin(\Omega) - \frac{1}{4} \cos(\Omega) \cot(\Omega) h^2 +\left(\frac{1}{64} (3 - \cos(2\Omega)) \cot(\Omega)^2 \csc(\Omega) \log(h)\right) h^4+\mathcal{O}(h^4).\label{b9}
\end{align}
 Now we want to find the expansion of $g_2(h)$ near the boundary. Since $\rho$ should be finite in the limit $h \rightarrow 0$ and $z \rightarrow 0$, the most sigular term in the expansion of $g_2$ should bo of order $h^{1-2\alpha}$. Inserting a power expansion for $g_2$ in terms of $h$ which begins from $h^{-5}$ into Eq. (\ref{b8}) gives the following result
\begin{align}
g_2=\dfrac{a_1}{h^{2\alpha -1}}+\dfrac{a_2}{h^{2\alpha -3}}+\frac{a_3}{ h}+a_4 h+\left(\dfrac{b_1}{h^{2\alpha -3}}+b_2 h\right)\log(h)+\mathcal{O}(h^3),
\label{b10}
\end{align}
where
\begin{align}
&a_2=\frac{a_1}{96} \csc(2\Omega)^2\left(39 + 48 \alpha- 16 \alpha ^2 + 76 \cos(2\Omega) + (37 - 48 \alpha + 16 \alpha ^2)\cos(4 \Omega) \right)\\
&a_3=-\dfrac{-2 + \alpha}{4 (-1 + \alpha^2)}\\
&b_1=-\frac{a_1}{8}\left(-3 + \cos(2 \Omega)\right) \csc(\Omega)^2\\
&b_2=\dfrac{3\left(5 + 4 \cos(2 \Omega) \right) \csc(2 \Omega)^2}{32 (2 + \alpha) (-1 + \alpha ^2)}
\end{align}
The undetermined constant ( $a_1$) in Eq. (\ref{b10}) should be fixed from the boundary conditions  \eqref{d2} and  \eqref{d3}.
Finally we should insert our results concerning $\sin \theta$, $\dot{h}$ and $g_2$  in terms of $h$ in Eqs. \eqref{d4} and  \eqref{d5} for $G_0$ and $G_1$  to find their needed singular terms as
\begin{align}
&G_0(h) \sim -\dfrac{\sin(\Omega)}{h^3}+\dfrac{\cos(\Omega)\cot(\Omega)}{8 h}+\mathcal{O}\left(h\right), \label{n3}\\
&G_1(h)\sim -\dfrac{2a_1\alpha \sin(\Omega)}{h^{2\alpha +1}}+\dfrac{(-3 + \alpha+ \alpha ^2) \sin(\Omega)}{2(-1 + \alpha^2)h^3}+\mathcal{O}\left(h^{-1}\right) .
 \label{n4}
\end{align}

Now, we can use Eqs.  \eqref{n3} and  \eqref{n4} in the integrands of Eq.  \eqref{bb6} to find their behavior near the boundary:
\begin{align}
I_0=&I_0'+I_0''\nonumber\\
&\int_{z_m} ^{\delta}\frac{dz}{z}  \int _{h_0}^{h_{1c}} dh\left(G_0(h)+\dfrac{\sin(\Omega)}{h^3}-\dfrac{\cos(\Omega)\cot(\Omega)}{8 h}\right)\nonumber\\
&+\int_{z_m} ^{\delta}\frac{dz}{z}  \int _{h_0}^{h_{1c}} dh\left(-\dfrac{\sin(\Omega)}{h^3}+\dfrac{\cos(\Omega)\cot(\Omega)}{8 h}\right), \label{n6}
\end{align}
\begin{align}
I_1=&I_1'+I_1''\nonumber\\
&\int_{z_m} ^{\delta}z^{2\alpha- 1} dz\int _{h_0}^{h_{1c}} dh\left(G_1(h)+\dfrac{2a_1\alpha \sin(\Omega)}{h^{2\alpha +1}}-\dfrac{(-3 + \alpha+ \alpha ^2) \sin(\Omega)}{2(-1 + \alpha^2)h^3}\right)\nonumber\\
&-\int_{z_m} ^{\delta} z^{2\alpha -1} dz\int _{h_0}^{h_{1c}}dh \left(\dfrac{2a_1\alpha \sin(\Omega)}{h^{2\alpha +1}}-\dfrac{(-3 + \alpha+ \alpha ^2) \sin(\Omega)}{2(-1 + \alpha^2)h^3}\right), \label{n7}
\end{align}
In order to find the singular terms of the integrals  \eqref{n6} and  \eqref{n7} we should find $h_{1c}$ similar to what did in Eq. \eqref{a155}. The result is
\begin{align}
h_{1c}(\delta)=&\frac{1}{H}\delta +a_1\mu^{2\alpha} H^{2\alpha -3}\delta ^3+a_2 \mu^{2\alpha} H^{2\alpha-5}\delta ^5+\mu^{2\alpha}\frac{a_3}{H}\delta ^{2\alpha+1}+\mu^{2\alpha} \frac{a_4}{H^3}\delta ^{2\alpha+3}\nonumber\\
&+\left(b_1 \mu^{2\alpha} H^{2\alpha-5}\delta ^5+\mu^{2\alpha} \frac{b_2}{H^3}\delta ^{2\alpha+3}\right)\log\delta+\mathcal{O}(\delta ^{2\alpha+5}).\label{b'11}
\end{align}
Following the same technique for calculating derivatives of the singular terms in terms of the cut-off parameter and integrating it again, we find the entanglement entropy of cone as follows

\begin{align}
S\vert_{c_1}=S\vert_{(0,c_1)}+S\vert_{(1,c_1)}
\end{align}
\begin{align}
&S\vert_{(0,c_1)}=-\frac{4\pi ^2 L^3}{l_p ^3}\left(-\dfrac{\sin(\Omega) H^2}{4\delta ^2}+\dfrac{\cos(\Omega)\cot(\Omega)}{16}{\log ^2\left(\frac{\delta}{H}\right)}+q_1\log\left(\frac{\delta}{H}\right)\right),\\
&S\vert_{(1,c_1)}=\left\{ \begin{array}{ll}
	-\frac{4\pi ^2 L^3}{l_p ^3}\left(\dfrac{H^2\mu^{2\alpha}  \sin(\Omega)}{8 (1 - \alpha)\delta ^{2 - 2 \alpha}}\right), \quad & 0 <\alpha <1\nonumber\\
-\frac{4\pi ^2 L^3}{l_p ^3}\left(-\dfrac{H^2 \mu^{2\alpha} \sin(\Omega)}{4} \log (\mu\delta) \right), \quad &\alpha=1 \nonumber\\
0\quad &  \alpha> 1 \end{array} \right.
\end{align}

\begin{table}[tbp]
\centering
\begin{tabular}{|c|c|c|c|c|c|}
			\hline
			$ d$ &Geometry  &Natural & Geometrical& $ \alpha$& New \\
			&   &  &  singulrities & & divergences for EE \\
			\hline
			&   &&&$ 0<\alpha<\frac{1}{2}  $ &  $1/\delta ^{1-2\alpha}$ \\
			3&$ k$ & $1/\delta $&$\log \delta$&$\alpha =\frac{1}{2}$& $\log\delta $\\
			&&&&$\alpha >\frac{1}{2}$& 0 \\
			\hline
			\hline
			& & &&$ 0<\alpha<1  $ &  $1/\delta ^{2-2\alpha}$ \\
			4&$ c_1$ &$1/\delta^2,\log\delta$&$\log^2\delta$ &$\alpha =1$& $\log\delta $\\
			&&&&$\alpha >1$& 0\\
			\hline
			\hline
			& &&& $ 0<\alpha<3/2,\alpha\neq1/2  $ &  $1/\delta ^{1-2\alpha}$,$1/\delta ^{3-2\alpha}$ \\
			5&$ c_2$&$1/\delta,1/\delta^3$&  $\log\delta $& $\alpha =1/2$ & $1/\delta^2$, $\log\delta $\\
			&& && $\alpha =3/2$ & $\log\delta $\\
			&&&&$\alpha >3/2$& 0\\
			\hline
			\hline
			& & &&$ 0<\alpha<2,\alpha\neq1  $ &  $1/\delta ^{2-2\alpha}$,$1/\delta ^{4-2\alpha}$ \\
			6&$ c_3$ &$1/\delta^2,1/\delta^4,\log\delta$&  $\log^2\delta $& $\alpha =1$ & $1/\delta^2$, $\log\delta $\\
			&& && $\alpha =2$ & $\log\delta $\\
			&&&&$\alpha >2$& 0\\
			\hline
\end{tabular}
\caption{\label{tab1}Entanglement entropy for kink-cones}
\end{table}
where
\begin{align}
&q_1=-\dfrac{\cos(\Omega)\cot(\Omega)}{8}\log(h_0)-\dfrac{\sin(\Omega)}{2 h_0^2}+ \int _{h_0}^{0} dh\bigg(\dfrac{\sin(\theta)\sqrt{1+h^2+\dot{h}^2}}{ \dot{h}h^3} +\dfrac{\sin(\Omega)}{h^3}\nonumber\\
&-\dfrac{\cos(\Omega)\cot(\Omega)}{8 h}\bigg) .
\end{align}

The above calculation can also be done for higher dimensional cones, i.e. for $d=5$ and $d=6$ corresponding to $n=2$ and $n=3$ in Eq. \eqref{a0}, respectively (see the appendix~\ref{Entropy1} for our final results).  In table~\ref{tab1} above, we have shown the nature of new singularities which emerge due to deformation of theory by a relevant operator. These singularities should be accompanied by the ordinary UV singularities, as well as singularities due to needle  points of the subregions considered. As is seen, for the case $c_1$ with the special tuning $M^2L^2=-3$ we have new logarithmic corrections in entanglement entropy. The same thing happens for example for $c_2$ with special tunings $M^2L^2=-9/4$ and $-21/4$.

\subsection{crease $ k\times R^l $}
\label{crease1}
Consider the case $n=0$ for arbitrary $l$  in the general metric (\ref{a0}),  where our singular subregion is limited to the region
 $ \theta \in [-\Omega, \Omega] $, $ \rho \in  [0,H] $  and $ x_i \in [-\frac{\tilde{H}}{2},\frac{\tilde{H}}{2}] $ where $H$ and $ \tilde{H} $ are IR cut-off. Assuming $ \rho =\rho (z ,\theta) $ as the RT surface, we need to optimize the following integral
\begin{align}
S\vert_{k \times R}=\frac{2\pi L^3\tilde{H}}{l_p ^3}\int d z d\theta \dfrac{\sqrt{\frac{\dot{\rho}^2+\rho ^2}{f(z)}+{\rho '}^2 \rho^2 }}
{z^3}.\label{c3}
\end{align}
Using the general ensatz \eqref{a5} and  inserting it in the equation  optimizing the integral (\ref{c3}), we can expand the expression of the  entanglement entropy as
\begin{align}
S\vert_{k\times R}=-\frac{4\pi L^3\tilde{H}}{l_p ^3}(I_1 +\mu^{2\alpha} I_2 + \cdots),\label{c6}
\end{align}
where
\begin{align}
I_1=\int_{z_m} ^{\delta }\frac{dz}{z^2} \int_ {h_0}^{h_{1c}} \frac{dh}{\dot{h}}\dfrac{\sqrt{1+h^2 +\dot{h}^2}}{ h^2},\label{s4}
\end{align}
\begin{align}
I_2=\int_{z_m} ^{\delta }dz z^{2\alpha -2} \int_ {h_0}^{h_{1c}} \frac{dh}{\dot{h}}\dfrac{ (-\dot{h}^2-h^2(1+2 \dot{h}\dot{g}_2)+4 (1+\alpha) h g_2+2 h^3 g_2)}{2 h^2 \sqrt{1+h^2+\dot{h}^2}}.\label{s5}
\end{align}

Similar to our treatment for kink in section 2, we find in zeroth order the following constant of motion (see Eq.  \eqref{a10})
\begin{align}
K_4=\dfrac{(1+h^2)^{\frac{3}{2}}}{h^3\sqrt{1+h^2+\dot{h}^2}}.\label{c5}
\end{align}
In the first order with respect to $\mu$ the singular behavior of $g_2$ near the boundary turns out to be given (similar to Eq. \eqref{a13}) by the following expansion with respect to the variable $h$
\begin{align}
g_2 =\dfrac{a_1}{h^{2\alpha +1}}+\frac{a_2}{ h^{2\alpha-1}}+\frac{a_3}{ h^{2\alpha-3}}+\frac{a_4}{ h}+a_5 h+\mathcal{O}(h^{3}).\label{s1}
\end{align}
where
\begin{align}
&a_2=-\frac{a_1}{6} (-3 - 7 \alpha + 2 \alpha^2),\label{h1}\\
&a_3=\frac{a_1}{96} (-18 - 3 \alpha+ 71 \alpha^2 - 36 \alpha^3 +
   4 \alpha^4),\label{h2}\\
&a_4=-\dfrac{1}{4(2 + \alpha)},\label{h3}\\
&a_5=\dfrac{9}{8 (2 + \alpha) (3 + 4 \alpha+ \alpha^2)}\label{h5}.
\end{align}
The quantity $h_{1c} = h(\Omega-\epsilon )$ turns out to have the following expansion (similar to Eq. \eqref{a155}) with respect to the UV cut-off
\begin{align}
h_{1c}(\delta)=&\left(\frac{1}{H}+a_1 \mu^{2\alpha} H^{2\alpha -1}\right)\delta +a_2 \mu^{2\alpha} H^{2\alpha -3}\delta ^3+a_3 \mu^{2\alpha} H^{2\alpha-5}\delta ^5\nonumber\\&+ \mu^{2\alpha}\frac{a_4}{H}\delta ^{2\alpha+1}+\mu^{2\alpha} \frac{a_5}{H^3}\delta ^{2\alpha+3}+
\mathcal{O}(\delta ^{2\alpha+5}).\label{ac10}
\end{align}
Using Eqs.  \eqref{c5},  \eqref{s1} and  \eqref{ac10} we can expand the integrand of Eqs. \eqref{s4} and \eqref{s5}. Then the singular terms with respect to the UV cut-off can be found similar to Eqs \eqref{aa15}--\eqref{a20}.
\begin{align}
S_{k\times R}=S_{(0,k\times R)}+S_{(1,k\times R)},
\end{align}

\begin{align}
S_{(0,k\times R)}=-\frac{4\pi L^3\tilde{H}}{l_p ^3}\left[-\frac{1}{\delta} \int_ {h_0}^{0} dh\left(\dfrac{\sqrt{1+h^2 +\dot{h}^2}}{\dot{h} h^2}+\frac{1}{h^2}\right)-\frac{H}{2\delta ^2}+\frac{1}{h_0 \delta}\right],
\label{c8}
\end{align}
\begin{align}
&S_{(1,k\times R),0<\alpha<1,\alpha\neq\frac{1}{2}}=\nonumber\\
&-\frac{4\pi L^3\tilde{H}}{l_p ^3}\mu^{2\alpha}\bigg[\frac{1}{(2\alpha-1)\delta ^{1-2\alpha}}\int_ {h_0}^{0}dh \bigg(\dfrac{ (-\dot{h}^2-h^2(1+2 \dot{h}\dot{g}_2)+4 (1+\alpha) h g_2+2 h^3 g_2)}{2 h^2\dot{h} \sqrt{1+h^2+\dot{h}^2}}\nonumber\\
&-\bigg(a_5 +  \dfrac{(1 - 2 a_4 )}{2 h^2} +
 h^{-2 \alpha} \big(a_2  (1 - 2 \alpha)  -\dfrac{a_1 (1 + 2 \alpha)}{h^2}\big)\bigg) \bigg)\nonumber\\
&-\dfrac{h_0^{-1-2\alpha}\left (2 a_2 h_0^2 + (-1 + 2 a_4) h_0^{2 \alpha} +   2 a_5 h_0^{2 + 2 \alpha} + 2 a_1\right) }{2(2\alpha-1)\delta ^{1-2\alpha }}- \dfrac{H}{4(\alpha -1)\delta ^{2-2\alpha}}\bigg],
\end{align}
\begin{align}
&S_{(1,k\times R),\alpha=\frac{1}{2}}=\nonumber\\&
\frac{4\pi L^3\tilde{H}\mu^{2\alpha}}{l_p ^3}\bigg[\dfrac{H}{2 \delta} -\dfrac{(-21 h_0 + 3 h_0^3 + 35 a_1 + 35 h_0^2 a_1}{35 h_0^2} \log(\mu\delta)+\nonumber\\
 &+\log (\mu\delta)\int_ {h_0}^{0}dh \bigg(\dfrac{ -\dot{h}^2-h^2(1+2 \dot{h}\dot{g}_2)+6 h g_2+2 h^3 g_2}{2 \dot{h}h^2 \sqrt{1+h^2+\dot{h}^2}}-\big(\dfrac{21 h + 3 h^3 - 70 a_1}{35 h^3}\big) \bigg)\bigg],\\
&S_{(1,k\times R),\alpha=1}=\frac{4\pi L^3\tilde{H}\mu^{2\alpha}}{l_p ^3}\dfrac{  H}{2}\log(\mu \delta), \\ & S_{(1,k\times R), \alpha>1}=0.
\end{align}
 The above calculation can also be done for higher dimensional cresaes, i.e. for $d=4$, $d=5$ and $d=6$ corresponding to $m=1$, $m=2$ and $m=3$ in Eq. \eqref{a0} respectively (see the appendix~\ref{Entropy1}  for our final results). In table~\ref{tab2} below, we have shown the nature of new singularities which emerge due to deformation of theory by a relevant operator.
\begin{table}[tbp]
\centering
\begin{tabular}{|c|c|c|c|c|c|}
\hline
$ d$ &Geometry & Natural & Geometrical& $ \alpha$ & New \\
 &   &   &singularities&&  divergences for EE \\
\hline
 &   &&& $ 0<\alpha<1,\alpha\neq\frac{1}{2} $&  $1/\delta ^{1-2\alpha}$, $1/\delta ^{2-2\alpha}$ \\
4&$ k\times R$ &$1/\delta ^2$&$1/\delta $& $\alpha =\frac{1}{2}$& $\log\delta $, $1/\delta$ \\
& &&& $\alpha =1$& $\log\delta $ \\
&&&&$\alpha >1$& 0 \\
\hline
 \hline
 &   & &&$ 0<\alpha<3/2,\alpha\neq1  $ &  $1/\delta ^{2-2\alpha}$, $1/\delta ^{3-2\alpha}$ \\
5&$ k\times R^2$&$1/\delta ^3$&$1/\delta ^2$ & $\alpha =1$& $\log\delta $, $1/\delta$ \\
& & &&$\alpha =\frac{3}{2}$& $\log\delta $ \\
&&&&$\alpha >\dfrac{3}{2}$& 0 \\
\hline
 \hline
 &   &&& $ 0<\alpha<2,\alpha\neq\frac{3}{2}  $ &  $1/\delta ^{3-2\alpha}$, $1/\delta ^{4-2\alpha}$ \\
6&$ k\times R^3$&$1/\delta ^4$&$1/\delta^3 $ & $\alpha =\frac{3}{2}$& $\log\delta $, $1/\delta$ \\
& &&& $\alpha =2$& $\log\delta $ \\
&&&&$\alpha >2$& 0 \\
\hline
\end{tabular}
\caption{\label{tab2}Entanglement entropy for $ k\times R^l $}
\end{table}

\subsection{crease $ c_n\times R^l $}
\label{crease2}
Now we consider some examples of the more general form of the bulk geometry introduced earlier in equation \eqref{a0}. For  $ l = 1 $ and $ n=1 $ the singular subregion in the boundary is defined as
$ \theta \in \left[ 0; \Omega \right] $, $ \phi \in \left[ 0; 2 \pi \right] $, $ \rho \in\left[ 0;  H \right] $ and $ x^l \in \left[ \frac{\tilde{H}}{2};- \frac{\tilde{H}}{2} \right] $. The entanglement entropy is obtained by minimizing the area of RT surface given by
\begin{align}
S\vert_{c_1\times R}=\frac{2\pi \Omega _1 L^4\tilde{H}}{l_{p }^4}\int  d z d\theta\dfrac{\rho\sin(\theta)\sqrt{\frac{\dot{\rho}^2+\rho ^2}{f(z)}+{\rho '}^2 \rho^2 }}
{z^4},\label{c10}
\end{align}
where $ \rho =\rho (z ,\theta) $ gives the extension of subregion to the bulk. As before this quantity can be expanded as given by the ensatz \eqref{a5}. Similar to case of $ c_n $ we find the expansions of the quantities $y=\sin \theta$ and $g_2$ near the boundary (i.e. in the limit $h\rightarrow 0$) as follows
\begin{align}
&y=\sin(\Omega) - \frac{1}{6}\cos(\Omega)\cot(\Omega) h^2 +  \frac{1}{432}\left(-19 + 5 \cos(2\Omega)\right) \cot^2(\Omega)\csc(\Omega) h^4+\mathcal{O}(h^5),\\
&g_2(h)=\frac{a_1}{h^{2 \alpha-2}}+\frac{a_2}{h^{2\alpha-4}}+\dfrac{a_3}{h} +a_4 h +\mathcal{O}(h^3).
\end{align}
where
\begin{align}
&a_2=\frac{1}{252} a_1 \left(-235 + 144 \alpha-
   36 \alpha^2 + (179 - 144 \alpha + 36 \alpha^2) \cos(2 \Omega) \right)\csc^2(\Omega),\\
&a_3=\dfrac{3-\alpha}{-6 - 2 \alpha+ 4 \alpha^2},\\
&a_4=\dfrac{\alpha \left(-14 + \alpha + (4 + \alpha) \cos(2 \Omega)\right) \csc^2(\Omega)}{9 (6 - 7 \alpha- 13 \alpha^2 + 4 \alpha^3 + 4 \alpha^4)}.
\end{align}
We also need to know the limiting value $h_{1c}$ in terms of the UV-cutoff. The result is
\begin{align}
h_{1c}(\delta)=\frac{\delta}{H} +a_1 \mu ^{2\alpha} H^{2\alpha -4}\delta ^4+a_2 \mu ^{2\alpha} H^{2\alpha-6}\delta ^6+ \mu ^{2\alpha}\frac{a_3}{H}\delta ^{2\alpha+1}+\mu ^{2\alpha} \frac{a_4}{H^3}\delta ^{2\alpha+3}+\mathcal{O}(\delta ^{2\alpha+5}).\label{a15}
\end{align}
Inserting all the data in the integral \eqref{c10} leads to the following result
\begin{align}
S_{c_1\times R}=S_{(0,c_1\times R)}+S_{(1,c_1\times R)},
\end{align}
where
\begin{align}
&S_{(1,c_1\times R),0<\alpha<\frac{3}{2},\alpha\neq\frac{1}{2}}=\nonumber \\&-\frac{4\pi ^2 L^4\tilde{H}\mu ^{2\alpha} }{l_p ^4}\bigg[\dfrac{H^2\sin(\Omega)}{(12 - 8 \alpha) \delta ^{3-2\alpha}}-\dfrac{h0^{-2 (1 + \alpha)}}{ \delta^{1-2\alpha}36 (1 -
     2 \alpha)^2 (1 + \alpha) (2 + \alpha) (-3 +
     2 \alpha)}\nonumber\\
 &\times\bigg(162 h_0^{2 \alpha} - 216 h_0^3 a_1 - 261 h_0^{2 \alpha} \alpha+
    16 h_0^{2 + 2 \alpha} \alpha+ 252 h_0^3 a_1 \alpha-
    171 h_0^{2 \alpha} \alpha^2 +
    4 h_0^{2 + 2 \alpha} \alpha^2 \nonumber\\
    &+ 468 h_0^3 a_1 \alpha^2 +
    72 h_0^{2 \alpha} \alpha^3 - 144 h_0^3 a_1\alpha^3 +
    36 h_0^{2 \alpha} \alpha^4 - 144 h_0^3 a_1 \alpha^4
     -  4 h_0^{2 + 2 \alpha} (-14 + \alpha) \alpha \times \nonumber\\
     &\csc^2(\Omega)+
    2 h_0^{2 + 2 \alpha}
     \cot^2(\Omega)\big(-24 + 4 \alpha +
       10 \alpha ^2
       + \left(18 - 37 \alpha - 7 \alpha^2 +
          16 \alpha^3 + 4 \alpha^4\right) \log(
         h_0) \nonumber\\&+ \left(-18 + 37 \alpha+ 7 \alpha^2 - 16 \alpha^3 -
          4 \alpha^4\right) \log(\delta/H)\big) \bigg)\sin(\Omega)\nonumber\\
       &+ \frac{1}{\delta ^{1-2\alpha}(-1 + 2 \alpha)}\int_{h_0} ^0 dh\dfrac{\sin(\theta) \left(-
\dot{h}^2 - h^2 (1 + 2 \dot{h}\dot{g_2}) +4 h^3 g_2 +  2 h (3 + 2 \alpha+ \dot{h}^2)g_2\right)}{2h^3\dot{h}\sqrt{1+h^2 +\dot{h}^2}}\nonumber\\
&-\big(-\dfrac{ (-9 + 5 \alpha+ 2 \alpha^2) \cos(\Omega) \cot(\Omega)}{18 h (-3 - \alpha + 2 \alpha^2)}+h^{-2 \alpha} a_1 (1 - 2 \alpha) \sin(\Omega)\nonumber\\
&+\dfrac{ (-9 + \alpha + 2 \alpha^2) \sin(\Omega)}{ h^3 (-6 - 2 \alpha + 4 \alpha^2)} \big)\bigg],
\end{align}

\begin{align}
&S_{(1,c_1\times R),\alpha=\frac{1}{2}}=\nonumber \\&-\frac{4\pi ^2 L^4\tilde{H}\mu ^{2\alpha} }{l_p ^4}\Bigg[\dfrac{H^2 \sin(\Omega)}{8\delta ^2}+  \frac{1}{72}\log(\mu\delta)\bigg(-\frac{10}{3} \cos(\Omega) \cot(\Omega) +    8 \cos(\Omega) \cot(\Omega) \nonumber\\ &\times\log(h_0)\-
    4 \cos(\Omega) \cot(\Omega)\log(\mu\delta)+
    9 \left(\frac{16}{3 h_0^2} - 8 a_1\right)  \sin(\Omega)\bigg)\nonumber\\
   +& \log(\mu \delta )\int_{h_0} ^0 dh\bigg(\dfrac{\sin(\theta) \left(-
\dot{h}^2 - h^2 (1 + 2 \dot{h}\dot{g_2}) +4 h^3 g_2 +  2 h (4+ \dot{h}^2)g_2\right)}{2h^3\dot{h}\sqrt{1+h^2 +\dot{h}^2}}\nonumber\\
&+\dfrac{h^2 \cos(\Omega) \cot(\Omega) - 12 \sin(\Omega)}{9 h^3}\bigg)\Bigg],\\
&S_{(1,c_1\times R),\alpha=\frac{3}{2}}=-\frac{4\pi ^2 L^4\tilde{H}\mu ^{2\alpha} }{l_p ^4}\left(-\frac{1}{4} H^2  \log(\mu\delta)\sin(\Omega)\right), \\& S_{(1,c_1\times R),\alpha>\frac{3}{2}}=0.
\end{align}
For crease $ c_1 \times R^2$ the final results are given in appendix~\ref{Entropy1}. Table \ref{tab3} below, shows the nature of new singular terms similar to previous cases.
\begin{table}[tbp]
\centering
\begin{tabular}{|c|c|c|c|c|c|}
\hline
$ d$ &Geometry &Natural & Geometrical& $ \alpha$ & New \\
 &   &  & singularities && divergences for EE \\
\hline
 &   & &&$ 0<\alpha<\frac{3}{2},\alpha\neq\frac{1}{2} $& $\log( \delta)/\delta ^{1-2\alpha}$ ,$1/\delta ^{1-2\alpha}$, $1/\delta ^{3-2\alpha}$ \\
5&$ c_1\times R$ &$1/\delta^3$,$1/\delta$&$\log(\delta)/\delta$& $\alpha =\frac{1}{2}$&$\log ^2(\delta) $ ,$\log(\delta) $, $1/\delta^2$ \\
& &&& $\alpha =\frac{3}{2}$& $\log(\delta) $ \\
&&&&$\alpha >\frac{3}{2}$& 0 \\
\hline
 &   &&& $ 0<\alpha<2,\alpha\neq1  $ &$\log\delta/\delta ^{2-2\alpha}$ , $1/\delta ^{2-2\alpha}$, $1/\delta ^{4-2\alpha}$ \\
6&$ c_1\times R^2$& $1/\delta^4$,$1/\delta^2,\log(\delta )$&$\log(\delta)/\delta$& $\alpha =1$& $\log ^2(\delta) $,$\log(\delta) $, $1/\delta^2$ \\
& & &&$\alpha =2$& $\log\delta $ \\
&&&&$\alpha >2$& 0 \\
\hline
 \hline
 &   &&& $ 0<\alpha<2,\alpha\neq\frac{1}{2},1  $ &  $1/\delta ^{1-2\alpha}$,$1/\delta ^{2-2\alpha}$,$1/\delta ^{4-2\alpha}$ \\
6&$ c_2\times R$ &$1/\delta^4$,$1/\delta^2,\log(\delta )$&$1/\delta$& $\alpha =\frac{1}{2}$& $\log(\delta) $, $1/\delta$,$1/\delta^3$ \\
& & &&$\alpha =1$& $\log(\delta) $, $1/\delta^2$\\
& & &&$\alpha =2$& $\log(\delta) $ \\
&&&&$\alpha >2$& 0 \\
\hline
  &   & &&$ 0<\alpha<\frac{5}{2},\alpha\neq1,\frac{3}{2}  $ &$1/\delta ^{2-2\alpha}$  $1/\delta ^{3-2\alpha}$, $1/\delta ^{5-2\alpha}$ \\
7&$ c_2\times R^2$ &$1/\delta^5$,$1/\delta^3,1/\delta$&$1/\delta^2$& $\alpha =1$& $\log(\delta) $, $1/\delta$ ,$1/\delta^3$\\
& & &&$\alpha =\frac{3}{2}$& $\log(\delta) $ ,$1/\delta^2$\\
& & &&$\alpha =\frac{5}{2}$& $\log(\delta) $ \\
&&&&$\alpha >\frac{5}{2}$& 0 \\
\hline
 \hline
\end{tabular}
\caption{\label{tab3}Entanglement entropy for $ c_n\times R^l $}
\end{table}

\section{Subregion Complexity}
\label{Complexity}
{As stated in the introduction, the subregion approach for complexity of a given static state \cite{Alishahiha:2015rta,Ben-Ami:2016qex} in a conformal theory concerns the volume enclosed by the Ryu-Takayanagi surface, i.e.
\begin{align}
\mathcal{C}=\dfrac{V(\gamma)}{8 \pi l G_N},\label{a22}
\end{align}
where $l$ is a characteristic length scale of the bulk geometry, and $ \gamma $ is the RT surface corresponding to the subregion specified in the boundary. Now let us calculate this quantity for the kink given by the metric \eqref{a0}. The volume surrounded by the surface with the induced metric \eqref{ab11} can be given as
\begin{align}
V(\gamma)=&L^3 \int dz d\theta d\rho \dfrac{\rho}{\sqrt{f(z)}z^3}\nonumber\\
&=L^3 \int dz d\theta \dfrac{1}{\sqrt{f(z)}z^3}\int _{\rho(z,\theta)} ^H d\rho \rho\nonumber\\
&=\frac{L^3}{2}\left(-\int dz d\theta \dfrac{1}{\sqrt{f(z)}z^3} \rho ^2+H^2\int dz d\theta \dfrac{1}{\sqrt{f(z)}z^3}\right)\nonumber\\
&\equiv V_1 +V_2,\label{a23}
\end{align}
Using the ansatz (\ref{a5}) in the integrand we find
\begin{align}
V_2=\dfrac{L^3 H^2 \Omega}{2 \delta ^2}+\dfrac{H^2 L^3 \mu ^{2\alpha} \Omega \delta ^ {2\alpha-2}}{4(\alpha-1)}+\mathcal{O}(\delta ^0),
\end{align}
\begin{align}
V_1=&L^3\int_{z_m}^\delta \frac{dz}{z} \int_{h_0}^{h_{1c}} dh \dfrac{1}{h^2 \dot{h}}+L^3\mu ^{2\alpha} \int_{z_m}^\delta dz z^{2\alpha -1}\int _{h_0}^{h_{1c}} dh\dfrac{-1+4 h g_2}{2 h^2\dot{h}}+ \cdots \\. \label{a24}
\end{align}
Let us search for the divergent terms of $V_1$ using
\begin{align}
\dfrac{1}{h^2 \dot{h}}\sim -k_3 +k_3 h^2+\mathcal{O}(h^3),
\end{align}
\begin{align}
\dfrac{-1+4 h g_2}{2 h^2\dot{h}}\sim -\dfrac{2a_1 k_3}{h^{2 \alpha}} +\frac{5+2\alpha}{2(3+2 \alpha)}k_3+\mathcal{O}(h^{2 -2\alpha}),
\end{align}
where $ a_1 $ is the  coefficient  defined in Eq. \eqref{a13}. In terms of different powers of $\mu$ we have
\begin{align}
V_1=L^3( I_1 +\mu ^{2\alpha} I_2+\cdots),
\end{align}
where
\begin{align}
I_1=\int_{z_m}^\delta \frac{dz}{z} \int_{h_0}^{h_{1c}} dh \dfrac{1}{h^2 \dot{h}},
\end{align}
\begin{align}
I_2=&\int_{z_m}^\delta dz z^{2\alpha -1}\int _{h_0}^{h_{1c}} dh\left(\dfrac{-1+4 h g_2}{2 h^2\dot{h}}+\dfrac{2a_1 k_3}{h^{2 \alpha}} -\frac{5+2\alpha}{2(3+2 \alpha)}k_3\right)\nonumber\\
&-\int_{z_m}^\delta dz z^{2\alpha-1}\int _{h_0}^{h_{1c}} dh\left(+\dfrac{2a_1 k_3}{h^{2 \alpha}} -\frac{5+2\alpha}{2(3+2 \alpha)}k_3\right)\nonumber\\
&=I_1'+I_2'.
\end{align}
In the expansions of the above integrands we encounter terms of the form $\int dz z^n \int dh h^{-m}$ for positive $n$ and $m$. To find the singularities we need just to keep terms with $m \ge n+2$. In the first integral $V_1$, which is just the complexity for the undeformed CFT, the integral over $h$ is finite, so we have a logarithmic divergent term for integration over $z$. However, in the expansions of the terms in the integrals $V_1$ over $h$, we do not find any term for which $m \ge n+2$. The final result reads

\begin{align}
\mathcal{C}_{k}=\mathcal{C}_{(0,k)}+\mathcal{C}{(1,k)}
\end{align}
\begin{align}
\mathcal{C}_{(0,k)}=\frac{L^3}{8 \pi l G}\left(\dfrac{ H^2 \Omega}{2 \delta ^2}+\log(\delta/H)\int_{h_0}^{0} dh \dfrac{1}{h^2 \dot{h}}\right)+finite,
\end{align}
\begin{align}
\mathcal{C}_{(1,k)}=\left\{ \begin{array}{ll}
	\frac{L^3 \mu ^{2\alpha}}{8 \pi l G}\left(\dfrac{H^2  \Omega}{4(-1+\alpha)\delta ^{2- 2 \alpha}}\right), \quad & 0 <\alpha <1\nonumber\\
\frac{L^3 \mu ^{2\alpha}}{8 \pi l G}\left(\dfrac{H^2  \Omega}{2}\log(\mu \delta) \right), \quad &\alpha=1 \nonumber\\
0 \quad &  \alpha>1\end{array} \right.
\end{align}

We can also find the subregion complexity for all of the cones discussed in the previous section. The final result for $c_1$ is
\begin{align}
\mathcal{C}_{c_1}=\mathcal{C}_{(0,c_1)}+\mathcal{C}{(1,c_1)}
\end{align}
where
\begin{align}
\mathcal{C}_{(0,c_1)}=\dfrac{L^3}{8 G_N}\left[\dfrac{2\left(1-\cos(\Omega)\right)}{9}\dfrac{H^3}{\delta ^3}-\frac{\cos(\Omega)}{3}\dfrac{ H}{\delta}+\dfrac{\beta(h_0)}{3}\log\left(\dfrac{\delta}{H}\right) \right],
\end{align}
\begin{table}
\centering
\begin{tabular}{|c|c|c|c|c|c|}
\hline
$ d$ &Geometry  &  Natural&Geometrical& $ \alpha$ & New \\
 &   &   & singularities& & divergences for $\mathcal{C}$ \\
\hline
 &    &&&$0<\alpha<1$ & $1/\delta ^{2-2\alpha}$\\
3&$ k$ &$1/\delta^2$&$\log(\delta) $&$\alpha=1$ &$\log(\delta) $\\
&&&&$\alpha>1$ &0 \\
\hline
 & &&&$ 0<\alpha<\frac{3}{2}  $,$\alpha\neq \frac{1}{2},\frac{3}{2} $ & $1/\delta ^{1-2\alpha}$,$1/\delta ^{3-2\alpha}$\\
4&$ c_1$ &$1/\delta^3$,$1/\delta$&$\log(\delta )$&$\alpha=\frac{1}{2}$ &$\log(\delta) $,$1/\delta^2$\\
&&&&$\alpha=\frac{3}{2}$ &$\log(\delta) $\\
&&&&$\alpha>3/2$ &0 \\
\hline
 & &&&$ 0<\alpha<2  $,$\alpha\neq 1,2 $& $1/\delta ^{2-2\alpha}$,$1/\delta ^{4-2\alpha}$\\
5&$ c_2$ &$1/\delta^4$,$1/\delta^2$,$\log(\delta) $&$\log^2(\delta) $&$\alpha=1$ &$\log(\delta) $,$1/\delta^2$\\
&&&&$\alpha=2$ &$\log(\delta) $\\
&&&&$\alpha>2$ &0 \\
\hline
 & &&&$ 0<\alpha<\frac{5}{2}  $,$\alpha\neq \frac{1}{2},\frac{3}{2},\frac{5}{2} $& $1/\delta ^{1-2\alpha}$,$1/\delta ^{3-2\alpha}$,$1/\delta ^{5-2\alpha}$\\
6&$ c_3$ &$1/\delta^5$,$1/\delta^3$,$1/\delta $&$\log(\delta )$&$\alpha=1/2$ &$\log(\delta) $,$1/\delta^2$,$1/\delta^4$\\
&&&&$\alpha=3/2$ &$\log(\delta) $,$1/\delta^2$\\
&&&&$\alpha=5/2$ &$\log(\delta) $\\
&&&&$\alpha>5/2$ &0 \\
\hline
\end{tabular}
\caption{\label{tab4}Complexity for $k-c_n$}
\end{table}
 \begin{table}
\centering
 \begin{tabular}{|c|c|c|c|c|c|}
\hline
$ d$ &Geometry & Natural&Geometrical& $ \alpha$ & New \\
 &   &  & singularities&& divergences for $\mathcal{C}$\\
\hline
 &   &&& $ 0<\alpha<\frac{3}{2} ,\alpha\neq \frac{1}{2}$&  $1/\delta ^{1-2\alpha}$, $1/\delta ^{3-2\alpha}$ \\
4&$ k\times R$ &$1/\delta$ ,$1/\delta^3$ &-& $\alpha =\frac{1}{2}$& $\log\delta $, $1/\delta^2$ \\
& & &&$\alpha =\frac{3}{2}$& $\log\delta $ \\
&&&&$\alpha >\frac{3}{2}$& 0 \\
\hline
 &   &&& $ 0<\alpha<2,\alpha\neq1  $ &  $1/\delta ^{2-2\alpha}$, $1/\delta ^{4-2\alpha}$ \\
5&$ k\times R^2$&$1/\delta^2$ ,$1/\delta^4$&- & $\alpha =1$& $\log\delta $, $1/\delta^2$ \\
& & &&$\alpha =2$& $\log\delta $ \\
&&&&$\alpha >2$& 0 \\
\hline
 &   &&& $ 0<\alpha<\frac{5}{2},\alpha\neq\frac{3}{2}  $ &  $1/\delta ^{3-2\alpha}$, $1/\delta ^{5-2\alpha}$ \\
6&$ k\times R^3$ &$1/\delta^3$ ,$1/\delta^5$&-& $\alpha =\frac{3}{2}$& $\log\delta $, $1/\delta^2$ \\
& & &&$\alpha =\frac{5}{2}$& $\log\delta $ \\
&&&&$\alpha >\dfrac{5}{2}$& 0 \\
\hline
\end{tabular}
\caption{\label{tab5}Complexity for $ k\times R^l $}
\end{table}
\begin{table}[tbp]
\centering
\begin{tabular}{ |c|c|c|c|c|}
\hline
Geometry & Natural&Geometrical &$ \alpha$ & New \\
   &  & singularities&&  divergences for $\mathcal{C}$ \\
\hline
    &&& $ 0<\alpha<2,\alpha\neq\frac{1}{2},1 $& $1/\delta ^{i-2\alpha}$ ,$i=1,2,4$ \\
$ c_1\times R$ &$\log(\delta) $,$1/\delta^i$,$ i=2,4$& $1/\delta$& $\alpha =\frac{1}{2}$&$\log(\delta) $,$1/\delta $ , $1/\delta^3$ \\
 & &&$\alpha =1$& $\log(\delta) $ , $1/\delta^2$\\
 &&& $\alpha =2$& $\log(\delta) $ \\
&&&$\alpha >2$& 0 \\
\hline
\hline
 &&& $ 0<\alpha<\frac{5}{2},\alpha\neq\frac{1}{2},\frac{3}{2}  $ &  $\log(\delta)/\delta ^{1-2\alpha}$,$1/\delta ^{i-2\alpha}i=1,3,5$ \\
$ c_2\times R$ &$1/\delta^i$,$i=1,3,5$&$\log(\delta)/\delta$& $\alpha =\frac{1}{2}$&  $\log^2(\delta )$,$\log(\delta )$, $1/\delta^2$, $1/\delta^4$ \\
 &&& $\alpha =\frac{3}{2}$& $\log(\delta) $,$1/\delta^2$ \\
 &&& $\alpha =\frac{5}{2}$& $\log(\delta) $\\
&&&$\alpha >\frac{5}{2}$& 0 \\
\hline
 &&  & $ 0<\alpha<3,\alpha\neq1,2  $ &   $\log(\delta)/\delta ^{2-2\alpha}$,$1/\delta ^{i-2\alpha}$, $i=2,4,6$ \\
$ c_2\times R^2$ &$\log(\delta)$,$1/\delta^i$,$ i=2,4,6$&$\log(\delta)/\delta^2$& $\alpha =1$& $\log^2(\delta) $,$\log(\delta) $, $1/\delta^2$, $1/\delta^4$ \\
 &&& $\alpha =2$& $\log(\delta) $, $1/\delta^2$\\
 &&& $\alpha =3$& $\log(\delta) $\\
&&&$\alpha >3$& 0 \\
\hline
\end{tabular}
\caption{\label{tab6} Complexity for $ c_n\times R^l $}
\end{table}
and
\begin{align}
\mathcal{C}_{(1,c_1)}=\left\{ \begin{array}{ll}
	\dfrac{L^3\mu ^{2\alpha}}{8 G_N}\left[-\dfrac{H (-3 + \alpha + \alpha ^2) \cos(\Omega) }{6 (-1 + 2 \alpha) (-1 + \alpha ^2)\delta ^{1-2\alpha}}+\dfrac{ H^3  (-1 + \cos(\Omega))}{\delta^{3 - 2 \alpha}(9-6 \alpha)}\right], \quad & 0 <\alpha <\frac{3}{2},\alpha \neq\frac{1}{2}\nonumber\\
\dfrac{L^3\mu ^{2\alpha}}{8 G_N}\left[\dfrac{ H^3(\cos(\Omega)-1)}{6 \delta^{2}}-\frac{1}{2}H  \cos(\Omega) \log(\mu\delta)\right], \quad &\alpha=\frac{1}{2} \nonumber\\
\dfrac{L^3\mu ^{2\alpha}}{8 G_N}\left[-\frac{1}{3} H^3 (-1 +\cos(\Omega))\log(\mu\delta)\right],\quad &  \alpha= \dfrac{3}{2}\nonumber\\
0, \quad &  \alpha> \dfrac{3}{2}\end{array} \right.
\end{align}

in which
\begin{align}
\beta(h_0)=&2 \int_ {h_0}^{0} dh\bigg(\dfrac{\sin(\theta)}{h^3\dot{h}}+\dfrac{1}{2} \dfrac{\cos( \Omega)}{h^2}-\frac{1}{8}\cot ^2 (\Omega ) \sin (\Omega ) \csc (2\Omega)(3 - \cos (2\Omega)) \log (h) \nonumber\\ &- \dfrac{1}{8}\cos(\Omega)\cot ^2(\Omega)\bigg)
-\dfrac{\cos(\Omega)}{ h_0} -\dfrac{h_0}{4}\cos(\Omega) \cot ^2(\Omega)+\frac{1}{4} h_0( 1-\log(h_0))\cot ^2 (\Omega )  \nonumber\\ &\times\sin (\Omega )
\csc (2\Omega)(3 - \cos (2\Omega)).
\end{align}
In the tables~\ref{tab4},\ref{tab5} and \ref{tab6} we give the nature of UV singular terms for different cases.

\section{Discussions}
\label{Discussions}
In this paper we concentrated on the effect of deformation of the boundary CFT due to a relevant  operator achieved by turning on a scalar field in the bulk. The entanglement entropy, as well as the subregion complexity, are important quantities  which demonstrate physical properties of the deformed CFT. The singular behavior of these quantities near the boundary or near the singular point (or wedge) of an assumed singular subregion may be viewed as key points toward investigating the properties of the corresponding CFT.

Considering the asymptotic  AdS-theory, we introduced a set of subregions characterized by singular points or wedges. For the pure AdS space-time in the bulk (i.e. undeformed theory) two kinds of singularity appears in entanglement entropy and complexity. The first category is the natural  UV singularity of these quantities. However, it is showed that for singular subregions, such as kink, cone or cresae, new kinds of singularity may appear as one  approaches the singular point or wedge. Deformation of the CFT by means of a relevant operator, imposes a third kind of singularity in the final expressions.

Our  results are summarized in tables  tables~\ref{tab1}--\ref{tab6} in the text. For completeness we have shown the known results for natural as well as geometrical singularities in independent columns (forth and fifth columns). The new kind of singularities due to deformation are $ \mu$-dependent terms which are shown in the sixth column of the corresponding tables. Tables \ref{tab1},\ref{tab2} and \ref{tab3} give the entanglement entropy for kink, cones, and creases $k\times R^l$ and  $C_n\times R^l$, respectively. The tables \ref{tab4},\ref{tab5} and \ref{tab6} give the subregion complexity for the same singular subregions. In the tables we just have indicated the kind of UV divergences. However, in appendices \ref{Entropy1} and \ref{Entropy2} we have given the complete singular terms for entanglement entropy and complexity respectively.

Among different results we point on the following important ones.\\
i) In all cases we have new logarithmic divergent terms in entanglement entropy for particular values of conformal dimension of the relevant operator which corresponds to some special tunings of the parameter $M$ of the scalar field in the bulk.\\
ii) The same thing happens for complexity but for different values of conformal dimension. In fact for deformed theory we may have logarithmic corrections in all dimensions for particular values of the conformal dimension.\\
ii) For crease $C_1\times R^l$ we found for entanglement entropy divergent terms of the form $\delta^{-k} \log \delta$ or $(\log \delta)^2$ depending on the values of $l$ and the conformal dimension. The same thing happens for complexity for crease $C_2\times R^l$.

\acknowledgments

 We would like to thank Mohsen Alishahiha, Amin Faraji-Astaneh  and Ali Naseh for fruitful discussions and M. Reza Mohammadi-Mozaffar and Ali Mollabashi for careful reading of the manuscript.


\section{Introduction}

The ability to perceive complex visual scenes in a structured and compositional way is crucial for humans to understand the seemingly chaotic world \cite{lake2017building}. Finding the underlying mechanisms of grouping basic visual elements is usually termed as the binding \cite{treisman1996binding,wolfe1999psychophysical}, or more specifically, the perceptual grouping \cite{grossberg1997visual} problem, and has been studied extensively in the fields of neuroscience, cognitive science and psychology. Designing more human-like AI systems which learn compositional and disentangled representations \cite{bengio2013representation} is desirable, because of their superior expressiveness and generalization ability compared to those systems learning a single complex representation. Complex visual scenes are composed of numerous relatively simpler primitive entities. By hierarchically decomposing the observed visual scenes \cite{bienenstock1997compositionality}, the raw observations can be summarized into organized and compact knowledge which is generalizable to an infinite number of novel scenes constructed by the combinations of these primitive entities \cite{biederman1987recognition,hummel1992dynamic,van2015part}.

Inspired by the synchronization theory \cite{milner1974model}, various approaches have been proposed to tackle the perceptual grouping problem by using neuronal synchrony as the grouping mechanism \cite{wang1995locally,rao2008unsupervised,reichert2013neuronal}. Despite of the  satisfactory grouping results already obtained, these methods do not learn separate representations for individual objects in the visual scene, which limits the expressiveness and generalization ability of the learned models. In order to learn disentangled and compositional representations of visual scenes, several methods which integrate spatial mixture models with neural networks have been proposed in recent years \cite{greff2015binding,greff2016tagger,premont2017recurrent,greff2017neural}. In these approaches, different objects are model by different mixture components, and each component is encouraged to focus on an individual object. The underlying intuition is that the image of visual scene can be decomposed into the pixel-wise weighted sum (modeled by mixture weights) of reconstructed images of individual objects (modeled by mixture components).

The perceptual grouping results of these recently proposed spatial mixture models are gratifying. Images of individual objects can be reconstructed accurately if objects do not overlap each other. However, there exist two potential downsides: 1) different attributes (e.g. shape, appearance) of objects cannot be well disentangled, and 2) the background is either not explicitly modeled by a mixture component or considered as an ordinary component, which in turn either complicates the determination of background region or increases the modeling complexity because the background region is much harder to represent due to the diverse combinations of foreground objects.

In this paper, we propose a new type of spatial mixture models which places a learnable deep prior on the model parameters and considers modeling the background as a special component complementary to the foreground objects. Parameters of the learnable deep prior consist of both mixture weight parameters and mixture component parameters, and are computed based on the outputs of two \emph{learnable} neural networks with latent representations of objects as inputs. The basic idea is that, in the proposed model, the mixture weights should focus on modeling the spatial dependencies (i.e., shapes) of the objects while the mixture components deal with variations (i.e., appearance) within objects. The inherent structure and learned parameters of the network regularize the model to correctly disentangle and estimate the ``shapes'' and ``appearances'' of objects in the visual scenes. By constructing the mixture weights in an innovative manner inspired by the truncated stick-breaking process, foreground objects and the background which is harder to represent are modeled differently. Because the background is considered as a special component whose mixture weights are implicitly determined by estimations of less varying foreground objects, the proposed method can thus achieve better compositionality and generalizes better to novel visual scenes.

We evaluate our method on two perceptual grouping datasets, in which images are composed of simple shapes or handwritten images, under different experimental configurations. Extensive empirical results suggest that representing the complex regions of background pixels in a compositional manner is crucial to high-quality grouping results. Compared with two state-of-the-art perceptual grouping algorithms, the proposed method not only achieves comparable or higher grouping accuracies, but also better estimates individual objects from the observed visual scene.
The learned conceptual entities are generalizable to novel visual scenes and insensitive to the diversity of objects.

\section{Preliminaries}

\subsection{Perceptual Grouping}

Perceptual grouping is a type of binding problem that studies how brains arrange elements of visual scenes into conceptual entities like objects or patterns. Based on the inductive biases that are either innate or learned from experiences, humans tend to segregate, for example, the upper-left corner image in Figure \ref{fig:shapes_tagger} into three white hollow shapes and one black background although these shapes are connected and partially overlapped. The underlying mechanisms have been studied extensively for years. Gestalt psychologists suggested that humans perceive meaningful contents from the sensory information based on inborn mental laws, and summarized their theories into several rule-based principles of perception (e.g. proximity, similarity, continuity, closure, and common fate) which is often known as Gestalt laws of grouping \cite{goldstein2016sensation}. The principle of common fate differs from others in that it requires temporal information and is only applicable when relative motions exist. In this work, we only consider utilization of spatial information and focus on building models which follow the principles of proximity, similarity, continuity, and closure by learning prior knowledge of spatial dependencies from data in an unsupervised setting.

\subsection{Spatial Mixture Model}

In spatial mixture models, pixels are assumed to be generated independently from the mixture. Correlations between pixels can be achieved by assuming mixture components, or mixture weights, or both of them to be spatially dependent. Let $X_m$ be the $m$th pixel of image $\boldsymbol{X}$, and $Z_m$ be the latent variable specifying the mixture component of $X_m$. $\boldsymbol{\theta}_{m,k} = \{\boldsymbol{\theta}_{m,k}^c, \boldsymbol{\theta}_{m,k}^w\}$ represents parameters for generating $X_m$ in the $k$th component. For the sake of notational simplicity, we let $p_{m,k}=P(X_m|Z_m\!\!=\!k; \boldsymbol{\theta}_{m}^c)$ and $\pi_{m,k}=P(Z_m\!\!=\!k; \boldsymbol{\theta}_{m}^w)=\theta_{m,k}^w$. The log probability of observing image $\boldsymbol{X}$ is given by
\begin{equation}
	\label{equ:spatial_mixture_model}
	\log{P(\boldsymbol{X}; \boldsymbol{\theta})} = \sum_{m}{\log{\sum_{k}{p_{m,k} \pi_{m,k}}}}
\end{equation}

Computing the maximum likelihood estimate (MLE) of $\boldsymbol{\theta}$ directly is difficult due to the existence of unobserved component assignments. Approximated estimates can be obtained iteratively by applying the Expectation-Maximization (EM) algorithm. Let $\gamma_{m,k}^{(t)}$ denote the posterior probability at the $t$th iteration $P(Z_m\!\!=\!k|X_m; \boldsymbol{\theta}_{m}^{(t)})$, which is computed by
\begin{equation}
	\label{equ:posterior}
	\gamma_{m,k}^{(t)} = \frac{p_{m,k}^{(t)} \pi_{m,k}^{(t)}}{\sum_{k'}{p_{m,k'}^{(t)} \pi_{m,k'}^{(t)}}}
\end{equation}
Parameters are updated using $\boldsymbol{\theta}^{(t+1)} \!=\! \argmax_{\boldsymbol{\theta}}{Q(\boldsymbol{\theta}; \boldsymbol{\theta}^{(t)})}$, where $Q(\boldsymbol{\theta}; \boldsymbol{\theta}^{(t)})$ is the Q-function defined as
\begin{equation}
	\label{equ:q_function}
	Q(\boldsymbol{\theta}; \boldsymbol{\theta}^{(t)}) = \sum_{m}{\sum_{k}{\gamma_{m,k}^{(t)}(\log{p_{m,k}} + \log{\pi_{m,k}})}}
\end{equation}

\subsection{Neural Expectation Maximization}

Neural Expectation Maximization (N-EM) \cite{greff2017neural} is a state-of-the-art spatial mixture model for solving the perceptual grouping problem. In N-EM, mixture components and mixture weights are modeled as spatially dependent and spatially independent, respectively. Parameters $\boldsymbol{\theta}_{m,k}^c \!\!=\!\! f_{\phi}(\boldsymbol{s}_k)_m$ are generated by a neural network $f_{\phi}$ with latent representations of objects $\boldsymbol{s}_k$ as inputs, and mixture weights $\pi_{m,k}$ which are identical at all pixels are either computed based on the EM update rule or fixed to $1 / K$ ($K$ is the number of mixture components) for simplicity \cite{greff2015binding}.

Because of the non-linearity of neural networks, the optimal latent representation $\boldsymbol{s}_k^*$ that maximizes $Q(\boldsymbol{\theta}; \boldsymbol{\theta}^{(t)})$ does not have a closed-form solution. Two types of generalized EM algorithms were proposed to solve this problem. One is to utilize the gradient descent to optimize $\boldsymbol{s}_k$ iteratively
\begin{equation}
	\boldsymbol{s}_k^{(t+1)} = \boldsymbol{s}_k^{(t)} + \eta \sum_{m}{\gamma_{m,k}^{(t)} \frac{\partial \log{p_{m,k}}}{\partial \boldsymbol{\theta}_{m,k}^c} \frac{\partial \boldsymbol{\theta}_{m,k}^c}{\partial \boldsymbol{s}_k}}\bigg|_{\boldsymbol{s}_k = \boldsymbol{s}_k^{(t)}}
\end{equation}
The other is to mimic the process of gradient descent by a recurrent neural network (RNN). The neural network $f_{\phi}$ is trained by the backpropagation through time (BPTT) algorithm \cite{robinson1987utility,werbos1988generalization}.

The loss function of each image is given by
\begin{align}
	\label{equ:loss_nem}
	L = & -\sum_{m}{\sum_{k}{\gamma_{m,k} (\log{p_{m,k}} + \log{\pi_{m,k}})}} \nonumber \\
	& + \lambda \sum_{m}{\sum_{k}{(1 - \gamma_{m,k}) D_{\text{KL}}(p_{\text{prior}}||p_{m,k})}}
\end{align}
This loss function consists of two parts. The first part is the negative Q-function which measures the reconstruction errors. The second part is the weighted sum of the Kullback-Leibler (KL) divergences between the mixture components $p_{m,k}$ and a distribution $p_{\text{prior}}$ predefined by the prior knowledge of the image (e.g. the background intensities). $p_{m,k}$ is encouraged to be close to the distribution of background when the $k$th component is unconfident that the $m$th pixel should be assigned to it ($\gamma_{m,k}$ is small).

\section{Spatial Mixture Model with\\Learnable Deep Priors}

In N-EM, only spatially dependent mixture components are transformed from latent representations of objects. Spatially independent mixture weights are modeled as coefficients which measure the importances of mixture components in the mixture model.
As a result, different attributes are entangled in mixture components and the region occupied by each object cannot be determined solely based on the latent representation. Moreover, the background is not explicitly modeled by a mixture component, and cannot be determined via the Maximum \emph{a posteriori} (MAP) estimation.

To better represent individual objects and model the background pixels, we integrate neural networks and spatial mixture models in a different way.
Both mixture weights and mixture components are transformed from latent representations that fully characterize attributes of objects.
Attributes are disentangled into ``shapes'' and ``appearances'', and modeled separately by mixture weights and mixture components. To follow the principle of similarity in the Gestalt laws of grouping, all pixels share the same mixture component parameters. In the following, we first introduce our approach in a general framework named Learnable Deep Priors and then make a concrete example for clarity.

\subsection{Learnable Deep Priors}

Let $P(\boldsymbol{\theta}_m'; \boldsymbol{\theta}_m) \!=\! \sum_k{\pi_{m,k} \delta(\boldsymbol{\theta}_m' - \boldsymbol{\theta}_{m,k}^c)}$, where $\delta(\cdot)$ is the Dirac delta function. The log probability \eqref{equ:spatial_mixture_model} of the observed image $\boldsymbol{X}$ can be rewritten as
\begin{equation*}
	\log{P(\boldsymbol{X}; \boldsymbol{\theta})} = \sum_{m}{\log{\int{P(X_m|\boldsymbol{\theta}_m') P(\boldsymbol{\theta}_m'; \boldsymbol{\theta}_m) d\boldsymbol{\theta}_m'}}}
\end{equation*}
where the prior $P(\boldsymbol{\theta}_m'; \boldsymbol{\theta}_m)$ of model parameters $\boldsymbol{\theta}_m'$ is defined with both \emph{mixture weights} $\pi_{m,k}$ and \emph{mixture component parameters} $\boldsymbol{\theta}_{m,k}^c$. Because $\pi_{m,k}$ and $\boldsymbol{\theta}_{m,k}^c$ are transformed from latent representations of objects by \emph{learnable} neural networks $f_{\phi}$ and $g_{\psi}$, $P(\boldsymbol{\theta}_m'; \boldsymbol{\theta}_m)$ is referred to as the \emph{learnable deep prior} in our method.

\textbf{Truncated Stick-Breaking Mixture Weights}.
The spatially dependent mixture weights $\pi_{m,k}$ are computed based on outputs of a neural network $f_{\phi}$ and are constrained by the structure and parameters of the network. Let $K$ be the number of components, $\boldsymbol{s}_k$ the latent representation of the $k$th component, and $\sigma(\cdot)$ the sigmoid function. $c_{k,m}$ is defined as $\sigma(f_{\phi}(\boldsymbol{s}_k)_m)$. The computation of $\pi_{m,k}$ which shares a similar form to the truncated stick-breaking process gives
\begin{equation}
	\label{equ:prior}
	\pi_{m,k} =
	\begin{cases}
		c_{k,m} \prod_{k'<k}{(1 - c_{k',m})}, & k < K \\
		\prod_{k'<K}{(1 - c_{k',m})}, & k = K
	\end{cases}
\end{equation}
\eqref{equ:prior} differs from the truncated stick-breaking process in that the length of stick $c_{k,m}$ is computed by a neural network instead of sampled from the beta distribution. Only the first $K - 1$ components which model foreground objects are assigned with latent representations. The complex-shaped background region is estimated based on the predicted shapes of these foreground objects.

Applying the softmax function to the outputs of neural networks is another viable way of modeling mixture weights. The main consideration of choosing \eqref{equ:prior} over the softmax function is \emph{compositionality}. The purpose of applying the mixture model is to decompose the complex visual scene into simpler components. Generally speaking, objects in the visual scene are relatively easy to represent. The background pixels, however, are almost as hard to model as the visual scene due to the diverse combinations of foreground objects. Because the softmax function treats each component equally, the complex-shaped background region need to be represented directly and the modeling complexity is increased significantly.

\textbf{Mixture Component Parameters}.
To follow the principle of similarity in the Gestalt laws of grouping, mixture component parameters $\boldsymbol{\theta}_{m,k}^c$ are modeled as spatially independent. If no additional constraint is added to $\boldsymbol{\theta}_{m,k}^c$, the optimal solution that maximizes the Q-function can be obtained by setting $\sum_{m}{\partial Q(\boldsymbol{\theta}; \boldsymbol{\theta}^{(t)}) / \partial \boldsymbol{\theta}_{m,k}^c}$ to zero. To fully characterize attributes of objects with latent representations, all $\boldsymbol{\theta}_{m,k}^c$ with $k < K$ are outputs of a neural network $g_{\psi}$ with latent representations $\boldsymbol{s}_k$ as inputs and are optimized by updating $\boldsymbol{s}_k$. Parameters $\boldsymbol{\theta}_{m,K}^c$ for the background component are modeled as learnable variables that are optimized by gradient descent in the M-step of the EM algorithm.

\subsection{Optimization of Mixture Model Parameters}

Latent representations $\boldsymbol{s}_{k}$ can be updated by either gradient descent with a learnable learning rate or an RNN which imitate the behavior of gradient descent. For the sake of notational simplicity, outputs of the neural networks $f_{\phi}(\boldsymbol{s}_k)_m$ and $g_{\psi}(\boldsymbol{s}_k)$ are denoted by $f_{k,m}$ and $g_k$, respectively. If using gradient descent as the update rule, latent representations $\boldsymbol{s}_k$ with $1 \leq k < K$ are updated by
\begin{multline}
	\label{equ:update_proposed}
	\boldsymbol{s}_k^{(t+1)} = \boldsymbol{s}_k^{(t)} + \eta_s \sum_{m}\Big(\gamma_{m,k}^{(t)} \frac{\partial \log{p_{m,k}}}{\partial g_k} \frac{\partial g_k}{\partial \boldsymbol{s}_k} \\
	+ \sum_{k' \geq k}{\gamma_{m,k'}^{(t)} \frac{\partial \log{\pi_{m,k'}}}{\partial f_{k,m}} \frac{\partial f_{k,m}}{\partial \boldsymbol{s}_k}}\Big)\bigg|_{\boldsymbol{s}_k = \boldsymbol{s}_k^{(t)}}
\end{multline}

If applying an RNN to learn the procedure of gradient descent, inputs to the RNN are features extracted from the concatenation of $\sum_{m}{\gamma_{m,k}^{(t)} (\partial \log{p_{m,k}} / \partial g_k)}$ and all $\sum_{k' \geq k}{\gamma_{m,k'}^{(t)} (\partial \log{\pi_{m,k'}} / \partial f_{k,m})}$ with $1 \leq m \leq M$ by an encoder network.

The update rule for the parameters $\boldsymbol{\theta}_{m,K}^c$ is given by
\begin{equation}
	\label{equ:update_last}
	{\boldsymbol{\theta}_{m,K}^{c^{(t+1)}}} \!\!= {\boldsymbol{\theta}_{m,K}^{c^{(t)}}} + \eta_{\theta} \sum_{m'}{\gamma_{m'\!,K}^{(t)} \frac{\partial \log{p_{m'\!,K}}}{\partial \boldsymbol{\theta}_{m'\!,K}^c}}\bigg|_{\boldsymbol{\theta}_{m'\!\!,K}^c = {\boldsymbol{\theta}_{m'\!\!,K}^{c^{(t)}}}}
\end{equation}

\subsection{Learning of Neural Networks}

Neural networks $f_{\phi}$ and $g_{\psi}$ are trained by BPTT algorithm. The loss function is chosen as
\begin{align}
	\label{equ:loss_proposed}
	L = & -\sum_{m}{\sum_{k}{\gamma_{m,k} (\log{p_{m,k}} + \log{\pi_{m,k}})}} \nonumber \\
	& + \lambda \sum_{m}{D_{\text{KL}}(p_{\text{prior}}||p_{m,K})}
\end{align}
As with the loss function \eqref{equ:loss_nem} used in N-EM, \eqref{equ:loss_proposed} also consists of the Q-function term and a regularization term. This regularization acts like a simplified version of the inductive biases required for human to solve the figure-ground organization problem (distinguishing a figure from the background). In N-EM, prior information of the background $p_{\text{prior}}$ is used to regularize all mixture components. In our method, the last mixture component $p_{m,K}$ which models the background is encouraged to be similar to $p_{\text{prior}}$, so that the network can better focus on modeling the foreground objects.

\subsection{Example: Gaussian Mixture Components}

The conditional probability distribution of a pixel given its component assignment is usually assumed to be Gaussian for real-valued images, and Bernoulli for binary images. We consider only Gaussian distribution in this section (detailed procedures for each image in an epoch are described in Algorithm \ref{alg:proposed}). Formulas for other types of probability distributions can be derived with little modification. Only the means of the Gaussian distributions are learnable parameters, and the variances are assumed to be constant and set to $\alpha^{-1}$ with $\alpha$ being a hyperparameter tuned by cross-validation. Except for the last component, the means of the Gaussian distributions are computed by the neural network $g_{\psi}$. For the $K$th component, it is modeled as a parameter $\mu_K$ that is learned by gradient descent. Let $x_m$ be the intensity of the $m$th pixel. The mixture component $p_{m,k}$ is computed by
\begin{equation}
	\label{equ:conditional_gaussian}
	p_{m,k} =
	\begin{cases}
		\mathcal{N}(x_m; g_{\psi}(\boldsymbol{s}_k), \alpha^{-1}), & k < K \\
		\mathcal{N}(x_m; \mu_K, \alpha^{-1}), & k = K
	\end{cases}
\end{equation}

\begin{algorithm}[t]
	\caption{Proposed Method (Gaussian Distribution)}
	\label{alg:proposed}
	\begin{algorithmic}
		\State // Initialization:
		\State Draw $\boldsymbol{s}_k^{(0)}$ randomly from $\mathcal{N}(\boldsymbol{0}, \sigma_{\text{init}}^2 \boldsymbol{I})$; \quad $\mu_K^{(0)} \gets \mu_{\text{prior}}$;
		\For{$t \gets 0$ to $T-1$}
		\State // E-Step:
		\State $f_{k,m}^{(t)} \gets f_{\phi}(\boldsymbol{s}_k^{(t)})_m$; \quad $g_k^{(t)} \gets g_{\psi}(\boldsymbol{s}_k^{(t)})$;
		\State Compute $p_{m,k}^{(t)}$, $\pi_{m,k}^{(t)}$ and $\gamma_{m,k}^{(t)}$ using \eqref{equ:conditional_gaussian}, \eqref{equ:prior} and \eqref{equ:posterior};
		\State // M-Step:
		\If{using gradient descent}
		\State Update $\boldsymbol{s}_k^{(t+1)}$ using \eqref{equ:update_proposed_gaussian};
		\Else
		\State $u \gets \alpha \sum_{m}{\gamma_{m,k}^{(t)} (x_m - g_k^{(t)})}$;
		\State $v_m \gets \gamma_{m,k}^{(t)} - \sum_{k' \geq k}{\gamma_{m,k'}^{(t)}} \sigma(f_{k,m}^{(t)})$;
		\State $r_{\text{in}} \gets enc([u, v_1, v_2, \dots, v_M])$;
		\State $\boldsymbol{s}_k^{(t+1)} \gets \text{RNN}(r_{\text{in}}, \boldsymbol{s}_k^{(t)})$;
		\EndIf
		\State Update $\mu_K^{(t+1)}$ using
		\eqref{equ:update_last_gaussian};
		\State // Loss:
		\State Compute $L^{(t)}$ using \eqref{equ:loss_proposed};
		\EndFor
		\State // Neural Networks and Learning Rates:
		\State Compute the weighted sum of losses $L_{\text{sum}} \gets \sum_{t}{w_t L^{(t)}}$;
		\State Update $\phi$, $\psi$, $\eta_s$ and $\eta_{\theta}$ based on $L_{\text{sum}}$;
	\end{algorithmic}
\end{algorithm}

Substituting \eqref{equ:conditional_gaussian} and \eqref{equ:prior} into \eqref{equ:update_proposed}, the gradient descent update rule for latent representations becomes
\begin{multline}
	\label{equ:update_proposed_gaussian}
	\boldsymbol{s}_k^{(t+1)} = \boldsymbol{s}_k^{(t)} + \eta_s \sum_{m}\Big(\alpha \gamma_{m,k}^{(t)} (x_m - g_k^{(t)}) \frac{\partial g_k}{\partial \boldsymbol{s}_k} \\
	+ \big(\gamma_{m,k}^{(t)} - \sum_{k' \geq k}{\gamma_{m,k'}^{(t)}} \sigma(f_{k,m}^{(t)})\big) \frac{\partial f_{k,m}}{\partial \boldsymbol{s}_k}\Big)\bigg|_{\boldsymbol{s}_k = \boldsymbol{s}_k^{(t)}}
\end{multline}

If updating $\boldsymbol{s}_k$ with an RNN, the encoder network $enc$ transforms the concatenation of $\alpha \sum_{m}{\gamma_{m,k}^{(t)} (x_m - g_k^{(t)})}$ and $\gamma_{m,k}^{(t)} - \sum_{k' \geq k}{\gamma_{m,k'}^{(t)}} \sigma(f_{k,m}^{(t)}), 1 \leq m \leq M$ to features $r_{\text{in}}$ as the inputs to the RNN.

Substituting \eqref{equ:conditional_gaussian} into \eqref{equ:update_last} and replacing $\boldsymbol{\theta}_{m,K}^c$ with $\mu_K$, the learnable parameter for the last component is updated by
\begin{equation}
	\label{equ:update_last_gaussian}
	\mu_K^{(t+1)} = \mu_K^{(t)} + \eta_{\theta} \alpha \sum_{m}{\gamma_{m,K}^{(t)} (x_m - \mu_K^{(t)})}
\end{equation}

\section{Related Work}

Several approaches have been proposed to solve the perceptual grouping problem and related tasks in recent years. Tagger \cite{greff2016tagger} combines the iterative amortized grouping (TAG) mechanism and the Ladder Network \cite{rasmus2015semi} to learn perceptual grouping in an unsupervised manner. It utilizes multiple copies of the same neural network to model different groups in the visual scene and iteratively refine the reconstruction result. RTagger \cite{premont2017recurrent} replaces the Ladder Network with the Recurrent Ladder Network and extends Tagger to sequential data. Neural Expectation Maximization (N-EM) \cite{greff2017neural} tackles the problem based on the Expectation-Maximization (EM) framework \cite{dempster1977maximum}, and achieves comparable performance to Tagger with much fewer parameters. Relational Neural Expectation Maximization (R-NEM) \cite{van2018relational} integrates N-EM with a type of Message Passing Neural Network \cite{gilmer2017neural} to learn common-sense physical reasoning based on the compositional object-representations extracted by N-EM.

Our method is mostly related to N-EM. Same as N-EM, the update rules for the distributed representations of objects in the visual scene are derived from the EM algorithm. Different from N-EM, the foreground objects and background are modeled separately, and attributes of objects are disentangled into ``shapes'' and ``appearances'' and modeled by mixture weights and mixture components, respectively.

\section{Experiments}

\subsubsection{Datasets}
Our method is evaluated on two datasets derived from a set of publicly released perceptual grouping datasets provided by \cite{greff2015binding,greff2016tagger,greff2017neural}. We refer to these two datasets as the \textit{Multi-Shapes} dataset and the \textit{Multi-MNIST} dataset. The Multi-Shapes dataset consists of $4$ subsets, which differs from each other in either the image size ($20 \!\times\! 20$ or $28 \!\times\! 28$) or the number of objects in each image ($2$, $3$ or $4$). Each subset contains $\num[group-separator={,}]{70000}$ binary images with simple shapes located at random positions. The Multi-MNIST dataset is composed of $3$ subsets with different degrees of object variations. In each subset, there are $\num[group-separator={,}]{70000}$ grayscale $48 \!\times\! 48$ images which contain $2$ handwritten digits. Images in both Multi-Shapes and Multi-MNIST datasets may contain overlapped objects.

\subsubsection{Comparison Methods}
The proposed method is compared with two state-of-the-art perceptual grouping algorithms, namely, \textit{N-EM} \cite{greff2017neural} and \textit{Tagger} \cite{greff2016tagger}. To assess the effectiveness of computing mixture weights using \eqref{equ:prior}, we substitute it with the softmax function and evaluate this modification (refered as \textit{softmax method}) on the Multi-Shapes dataset. Neural networks are trained with the Adam optimization algorithm \cite{kingma2014adam} for all approaches. To make a fair comparison, we do not specially design the neural networks $f_{\phi}$ and $g_{\psi}$ for the proposed method. The structure of $f_{\phi}$ which generates mixture weights in our method is identical to the network $f_{\phi}$ in N-EM. The extra $g_{\psi}$ for computing mixture component parameters is simply a one-layer fully-connected neural network, in which the number of parameters equals the dimension of latent representations and negligible compared to the total number of network parameters.

\subsubsection{Evaluation Metrics}
We use two metrics to evaluate the grouping results. The first is \textit{Adjusted Mutual Information} (AMI), which measures the accuracy of mixture assignments. To be consistent with previous work \cite{greff2015binding,greff2016tagger,greff2017neural}, AMI scores are not computed in the background and overlapping regions. The second is \textit{Mean Squared Error} (MSE), which computes the smallest mean squared differences between the estimated images of individual objects and the ground truth images over all possible permutations (formulated as the assignment problem and solved efficiently by the Hungarian method). MSE scores are evaluated at all pixels.

\subsection{Perceptual Grouping Performance Comparisons}

Performances of Tagger, N-EM and the proposed method are compared on two subsets of the Multi-Shapes dataset which differ from each other in image size ($20 \!\times\! 20$ and $28 \!\times\! 28$). Images in both subsets contain $3$ objects. Because objects in both subsets are of the same size, $20 \!\times\! 20$ images which consist of less pixels are more likely to contain overlapped objects, and performing perceptual grouping on them is not necessarily easier. For Tagger, parameters of mixture models are updated via a $3$-layer Ladder Network. For N-EM and the proposed method, the encoder and decoder networks are convolutional neural networks (CNNs) with $2$ convolutional, $2$ fully-connected and $3$ layer normalization \cite{ba2016layer} layers. The qualitative results are presented in Figure \ref{fig:shapes_compare_methods}, and the AMI/MSE scores are presented as follows.

\begin{center}
	\begin{tabular}{c C{1in} C{1in}}
		\hline
		& Size $20 \times 20$ & Size $28 \times 28$ \\ \hline
		Tagger & \textbf{0.933} / 0.71e-2 & 0.820 / 1.49e-2 \\
		N-EM & 0.824 / 1.59e-2 & 0.897 / 0.82e-2 \\
		Proposed & 0.920 / \textbf{0.40e-2} & \textbf{0.941} / \textbf{0.28e-2} \\
		\hline
	\end{tabular}
\end{center}

On the subset consisting of $20 \times 20$ images, Tagger achieves the highest grouping accuracy (AMI) and the proposed method best estimates the images of individual objects (MSE). The AMI score of the proposed method is slightly lower than Tagger. On the other subset, The proposed method achieves both the best AMI and MSE scores.

According to Figure \ref{fig:shapes_tagger} (\textit{Tagger}), the learned spatial dependencies of pixels are entangled. Both mixture components (row 2; cols 2--5, 7--10) and mixture weights (row 3; cols 2--5, 7--10) contain partial and imperfect information of shape and appearance. It is caused by modeling both mixture components and mixture weights to be spatially dependent without adding additional regularizations. Another observation is that objects are mainly distinguished based on mixture components, and the mixture weights assist the segregation when ambiguity exists. Although the mixture component assignments are predicted accurately for foreground objects (row3; cols 1, 6), the background pixels are not well inferred based on posterior probabilities $\gamma$ (row2; cols 1, 6). The possible reason is that the diverse combinations of foreground objects prevent background pixels from being modeled well using non-combinational representations.

In Figure \ref{fig:shapes_nem} (\textit{N-EM}), the colors of background pixels in the $2$nd row are gray (red+green+blue) because the posterior probabilities $\gamma$ are almost identical for all mixture components. Similar to Tagger, background pixels cannot be well determined via the Maximum \emph{a posteriori} (MAP) estimation in N-EM. This is because N-EM does not use a mixture component to model the background.

As demonstrated in Figure \ref{fig:shapes_proposed} (\textit{proposed method}), objects are well segregated based on the learned knowledge of spatial dependencies even if objects are heavily overlapped (cols 2, 6). Although the reconstructed images of individual objects are not completely correct, by modeling the background complementary to the foreground objects and disentangling different attributes of objects, the predicted results follow the principles of proximity, similarity, continuity, and closure in the Gestalt laws of grouping.

\begin{figure}[!ht]
	\centering
	\subfloat[Tagger (Size $20 \times 20$)]{\includegraphics[width=\linewidth]{fig_shapes_tagger
		\label{fig:shapes_tagger}}
	\\
	\subfloat[N-EM (Size $28 \times 28$)]{\includegraphics[width=\linewidth]{fig_shapes_nem
		\label{fig:shapes_nem}}
	\\
	\subfloat[Proposed (Size $28 \times 28$)]{\includegraphics[width=\linewidth]{fig_shapes_proposed
		\label{fig:shapes_proposed}}
	\\
	\subfloat[Softmax (Size $28 \times 28$)]{\includegraphics[width=\linewidth]{fig_shapes_softmax
		\label{fig:shapes_softmax}}
	\caption{Component assignments and reconstructed images of individual objects evaluated on the Multi-Shapes dataset.}
	\label{fig:shapes_compare_methods}
\end{figure}

\subsection{Effectiveness of Modeling Background Separately}

Besides the observation from Figure \ref{fig:shapes_tagger} that modeling the background as an ordinary component results in imperfect separation of foreground and background based on posterior probabilities, we also assess the effectiveness of modeling the background separately by evaluating the performances of the softmax method on the second subset (size $28 \!\times\! 28$) used in the previous experiment. The best AMI and MSE scores achieved are $0.875$ and $1.10 \!\times\! 10^{-2}$, which are significantly worse than the results of the proposed method ($0.941$ and $0.28 \!\times\! 10^{-2}$). The burden of directly representing the complex-shaped background region prevents the model to focus on spatial dependencies of foreground pixels. The qualitative results of the softmax method is shown in Figure \ref{fig:shapes_softmax}. The estimated images of individual objects are not visually plausible although the results of posterior inferences (row 2) are satisfactory.

\subsection{Generalizability and Sensitivity Analyses}

\subsubsection{Generalizability of Compositional Representations}

To illustrate the generalizability of the learned compositional representations, the proposed method are evaluated on three subsets of the Multi-Shapes dataset that differ in the number of objects in images. Results of N-EM are also presented. Models of both methods are trained \textit{only} on the second subset. The AMI/MSE scores are shown below.
\begin{center}
	\begin{tabular}{c C{0.9in} C{0.9in}}
		\hline
		& N-EM & Proposed \\ \hline
		2 Objects & \textit{0.906} / 0.98e-2 & \textit{0.961} / \textit{0.11e-2} \\
		3 Objects (Train) & 0.897 / \textit{0.82e-2} & 0.941 / 0.28e-2 \\
		4 Objects & 0.838 / 1.44e-2 & 0.901 / 0.69e-2 \\
		\hline
	\end{tabular}
\end{center}
For both approaches, the learned representations of conceptual entities generalize well to visual scenes containing more or less objects. The accuracies of estimated mixture component assignment are slightly higher when the task is simpler, and moderately lower when dealing with more challenging visual scenes. The MSE score of N-EM becomes slightly worse on the simpler task.

\subsubsection{Sensitivity to Choice of Conditional Distributions}

\begin{table}
	\centering
	\begin{tabular}{|c|c|c||C{0.37in}|C{0.42in}||C{0.37in}|C{0.42in}|}
		\hline
		& \multirow{2}{*}{$F$} & \multirow{2}{*}{$N$} & \multicolumn{2}{c||}{N-EM} & \multicolumn{2}{c|}{Proposed} \\ \cline{4-7}
		& & & AMI & MSE & AMI & MSE \\ \hline\hline
		\multirow{5}{*}{\STAB{\rotatebox[origin=c]{90}{Bernoulli}}}
		& 16 & 16 & 0.876 & 0.95e-2 & 0.856 & 1.28e-2 \\ \cline{2-7}
		& 32 & 16 & \textit{0.897} & \textit{0.82e-2} & 0.868 & 1.23e-2 \\ \cline{2-7}
		& 32 & 32 & 0.816 & 1.32e-2 & 0.926 & 0.43e-2 \\ \cline{2-7}
		& 64 & 32 & 0.814 & 1.34e-2 & 0.920 & 0.40e-2 \\ \cline{2-7}
		& 64 & 64 & 0.749 & 1.70e-2 & \textbf{0.941} & \textbf{0.28e-2} \\
		\hline\hline
		\multirow{5}{*}{\STAB{\rotatebox[origin=c]{90}{Gaussian}}}
		& 16 & 16 & \textit{0.687} & 2.35e-2 & 0.915 & 0.68e-2 \\ \cline{2-7}
		& 32 & 16 & 0.676 & \textit{2.34e-2} & 0.911 & 0.63e-2 \\ \cline{2-7}
		& 32 & 32 & 0.566 & 2.87e-2 & 0.914 & 0.51e-2 \\ \cline{2-7}
		& 64 & 32 & 0.503 & 3.06e-2 & 0.935 & 0.30e-2 \\ \cline{2-7}
		& 64 & 64 & 0.353 & 3.77e-2 & \textbf{0.939} & \textbf{0.29e-2} \\
		\hline
	\end{tabular}
	\caption{Comparisons of AMI and MSE scores with different mixture component distributions. $F$ and $N$ are respective dimensions of RNN inputs and latent representations.}
	\label{tab:shapes}
\end{table}

We assess the influences of different choices of mixture component distributions on the second subset ($3$ objects) used in the previous experiment. Detailed performances of N-EM and the proposed method are shown in Table \ref{tab:shapes}. If intensities of binary images are modeled by Gaussian distributions, the grouping results of N-EM deteriorate drastically. The proposed method, on the other hand, are less sensitive to the form of mixture component distributions.

\subsubsection{Sensitivity to Diversity of Objects}

\begin{table}
	\centering
	\begin{tabular}{|c|c|c||C{0.37in}|C{0.42in}||C{0.37in}|C{0.42in}|}
		\hline
		& \multirow{2}{*}{$F$} & \multirow{2}{*}{$N$} & \multicolumn{2}{c||}{N-EM} & \multicolumn{2}{c|}{Proposed} \\ \cline{4-7}
		& & & AMI & MSE & AMI & MSE \\ \hline\hline
		\multirow{4}{*}{\#1}
		& 32 & 16 & \textit{0.780} & \textbf{1.17e-2} & 0.748 & 1.89e-2 \\ \cline{2-7}
		& 32 & 32 & 0.573 & 1.54e-2 & 0.765 & 1.64e-2 \\ \cline{2-7}
		& 64 & 32 & 0.521 & 1.54e-2 & 0.787 & 1.35e-2 \\ \cline{2-7}
		& 64 & 64 & 0.391 & 1.79e-2 & \textbf{0.790} & \textit{1.30e-2} \\
		\hline\hline
		\multirow{4}{*}{\#2}
		& 32 & 16 & \textit{0.317} & 2.86e-2 & \textbf{0.719} & 2.07e-2 \\ \cline{2-7}
		& 32 & 32 & 0.254 & 2.75e-2 & 0.704 & 2.04e-2 \\ \cline{2-7}
		& 64 & 32 & 0.194 & \textit{2.71e-2} & 0.719 & \textbf{1.61e-2} \\ \cline{2-7}
		& 64 & 64 & 0.152 & 2.73e-2 & 0.676 & 1.63e-2 \\
		\hline\hline
		\multirow{4}{*}{\#3}
		& 32 & 16 & \textit{0.335} & 2.93e-2 & \textbf{0.715} & 2.16e-2 \\ \cline{2-7}
		& 32 & 32 & 0.223 & 3.15e-2 & 0.702 & 2.07e-2 \\ \cline{2-7}
		& 64 & 32 & 0.185 & 2.94e-2 & 0.705 & 1.69e-2 \\ \cline{2-7}
		& 64 & 64 & 0.158 & \textit{2.84e-2} & 0.683 & \textbf{1.62e-2} \\
		\hline
	\end{tabular}
	\caption{AMI and MSE scores evaluated on different subsets of the Multi-MNIST dataset. $F$ and $N$ are respective dimensions of RNN inputs and latent representations.}
	\label{tab:static_mnist}
\end{table}

The Multi-MNIST dataset is much more challenging than the Multi-Shapes dataset for the larger image size ($48 \!\times\! 48$ versus $20 \!\times\! 20$ or $28 \!\times\! 28$), as well as higher degrees of both inter-object variations (the number of unique objects) and intra-object variations (intensities of pixels belonging to the same object). There are $3$ subsets in this dataset. Subsets 1, 2 and 3 are constructed using $20$, $500$ and $\num[group-separator={,}]{70000}$ unique digits, respectively. Different from the first two subsets, the digits composing images in the test set of subset 3 do not exist in the training and validation set. We evaluate the performances of N-EM and the proposed method on these subsets, and compare the sensitivities to object diversities of the two approaches. Experimental results are presented in Table \ref{tab:static_mnist}.

When the diversity of objects is small, N-EM and the proposed method perform similarly well. N-EM achieves slightly higher AMI score and the proposed method attains lower MSE score. As the variations of objects increase, the performance of the proposed method drops moderately, and the grouping results of N-EM decline significantly. On subsets 2 and 3, both methods achieve the highest AMI scores when dimensions of RNN inputs and latent representations are small, and the best MSE scores when these dimensions are slighly larger. The possible reason is that higher dimensions allow models to focus on more details and help reconstruct individual digits. The large diversities of objects, on the other hand, misguide models to represent parts of different objects in the same mixture component and lower the grouping accuracies as the capacities of models increase.

Figure \ref{fig:mnist} demonstrates the predicted component assignments and reconstructed images of individual objects produced by both methods. On subset 1, the results generated by N-EM are slightly more visually plausible. Except for the occluded regions, details of each object are accurately reconstructed. This is because N-EM allows the mixture components to vary among pixels, which helps better reconstruct the intra-object variations of each digit. Both methods can accurately segregate objects from some complex visual scenes like the last two images in Figure \ref{fig:mnist_20_nem} and \ref{fig:mnist_20_proposed}. On subsets 2 and 3, the estimated images of individual objects generated by N-EM are also sharper than the proposed method. However, each component contains pixels of multiple objects. The proposed method models the spatial dependencies of pixels with mixture weights, and utilizes the principle of similarity implied by the spatially independent mixture components to assist the learning of correlations between pixels. Empirical results suggest that this scheme is less sensitive to the diversities of objects in the visual scenes.

\begin{figure}[t]
	\centering
	\subfloat[N-EM (20 variants)]{\includegraphics[width=0.5\linewidth]{fig_mnist_20_nem
		\label{fig:mnist_20_nem}}
	\hfil
	\subfloat[Proposed (20 variants)]{\includegraphics[width=0.5\linewidth]{fig_mnist_20_proposed
		\label{fig:mnist_20_proposed}}
	\\
	\subfloat[N-EM (500 variants)]{\includegraphics[width=0.5\linewidth]{fig_mnist_500_nem
		\label{fig:mnist_500_nem}}
	\hfil
	\subfloat[Proposed (500 variants)]{\includegraphics[width=0.5\linewidth]{fig_mnist_500_proposed
		\label{fig:mnist_500_proposed}}
	\\
	\subfloat[N-EM (70,000 variants)]{\includegraphics[width=0.5\linewidth]{fig_mnist_all_nem
		\label{fig:mnist_all_nem}}
	\hfil
	\subfloat[Proposed (70,000 variants)]{\includegraphics[width=0.5\linewidth]{fig_mnist_all_proposed
		\label{fig:mnist_all_proposed}}
	\caption{Component assignments and reconstructed images of individual objects evaluated on the Multi-MNIST dataset.}
	\label{fig:mnist}
\end{figure}

\section{Conclusion}

In this paper, we have proposed a framework called Learnable Deep Priors (LDP) which defines priors of model parameters with compositional latent representations and learnable neural networks for tackling the perceptual grouping problem. By modeling the background as a special component in a compositional manner and decomposing the attributes of conceptual entities into shapes and appearances that are separately represented, the proposed method better learns spatial dependencies than existing methods for perceptual grouping. We have demonstrated that the proposed method is insensitive to diversities of objects, and the learned compositional representations of individual entities generalize well to visual scenes constructed by novel combinations of these entities. Further research in this direction could be considering learnable deep priors with latent representations of conceptual entities organized in a structured form, integrating hierarchical spatial mixture models with neural networks, and solving higher-level tasks like relational inferences based on the learned compositional representations.

\section{Acknowledgments}

This work was supported in part by National Key R\&D Program of China (No.2017YFC0803700), NSFC under Grant No.61572138 \& No.U1611461, and STCSM Projects under Grant No.16JC1420400 \& No.18511103104.