Shannon and Renyi mutual information in quantum critical spin chains
Version 1 Released on 04 April 2015 under Creative Commons Attribution 4.0 International LicenseAuthors' affiliations
- Topology and Correlations in Condensed Matter Research Group, Department of Condensed Matter, Max Planck Institute for the Physics of Complex Systems
Keywords
- Conformal field theory
- Critical point phenomena
- Spins
Abstract
We study the Shannon mutual information in one-dimensional critical spin chains, following a recent conjecture (Phys. Rev. Lett. 111, 017201 (2013)), as well as Rényi generalizations of it. We combine conformal field theory arguments with numerical computations in lattice discretizations with central charge $c=1$ and $c=1/2$. For a periodic system of length $L$ cut into two parts of length $\ell$ and $L-\ell$, all our results agree with the general shape-dependence $I_n(\ell,L)=(b_n/4)\ln \left(\frac{L}{\pi}\sin \frac{\pi \ell}{L}\right)$, where $b_n$ is a universal coefficient. For the free boson CFT we show from general arguments that $b_n=c=1$. At $c=1/2$ we conjecture a result for $n>1$. We perform extensive numerical computations in Ising chains to confirm this, and also find $b_1\simeq 0.4801629(2)$, a nontrivial number which we do not understand analytically. Open chains at $c=1/2$ and $n=1$ are even more intriguing, with a shape-dependent logarithmic divergence of the Shannon mutual information.
Introduction
The entanglement entropy (EE) has emerged as a particularly convenient tool in the study of quantum many-body systems. For example in one-dimensional critical systems it is known to be universal [26,60,8], with a slow logarithmic divergence proportional to the central charge of the underlying conformal field theory [14,4](CFT). It is not the only information-theoretic quantity that exhibits universal behavior: various types of other entanglement[9] or fidelity[16] measures have been shown to be universal. However one particularly attractive feature of the EE is that it only depends to leading order on the central charge, not on more refined properties of the CFT. The Shannon entropy, a measure of disorder in a certain basis, is also an interesting quantity in its own right. It is defined as \begin{equation} S=-\sum_{\sigma} p_\sigma \ln p_\sigma\qquad,\qquad p_\sigma=\left|\left\langle \sigma|\psi \right\rangle\right|^2, \end{equation} where $\left|\psi\right\rangle$ will be for us the ground state wavefunction. The sum runs over all configurations in the chosen basis of the Hilbert space. This entropy has been studied in a variety of contexts, from the Anderson localization [18,19], multifractality[3] and quantum chaos [30,49], to quantum quenches[47] and critical phenomena. It is is also related [53] to the entanglement entropy of certain 2d Rokhsar-Kivelson states [21,27,53,56,45,62,57,31,55], as well as the classical mutual information in 2d systems [54]. In the following we will only consider a basis obtained from tensor products of local degrees of freedom. Typically in such cases the Shannon entropy obeys a “volume” law1 : it is proportional to the volume $L^d$ of the $d-$dimensional quantum system. Subleading terms are also interesting, as they have been shown to be universal. For example in periodic 1d chains whose low-energy properties are governed by a Luttinger liquid theory, a subleading term gives access to the Luttinger parameter [53]. In other systems or geometries it can be used to extract universal properties [56,62,57] of the underlying CFT, and thus identify it. We refer to Ref. [37] for a review. Universal terms can also be accessed in higher dimensions [36,35,37]. Inspired by all the results obtained for the one-dimensional entanglement entropy, it is natural to study the Shannon entropy of a subsystem [59]. Let us cut our chain in two parts $A$ and $B$. Subsystem $A$ is described by the reduced density matrix $\rho_A=\textrm{Tr}_B \left|\psi\right\rangle\left\langle \psi\right|$, and the subsystem Shannon entropy is \begin{equation} S(A)=-\sum_{\mu} p_\mu \ln p_\mu\qquad,\qquad p_\mu=\left\langle \mu|\rho_A|\mu \right\rangle. \end{equation} Alternatively, $p_\mu$ can be seen as the marginal probability obtained from summing over all spin configurations outside of $A$. Then, a most natural object is the mutual information between them, defined as \begin{equation}\label{eq:smidef} I(A,B)=S(A)+S(B)-S(A\cup B). \end{equation} $S(A\cup B)$ is the Shannon entropy of the total system. Note that the mutual information is symmetric with respect to $A$ and $B$, and that the leading contributions proportional to $L$ cancel in the definition (\ref{eq:smidef}). It has been first studied for the Ising universality class [33,59], and the scaling argued to be universal. More precisely, it was found that the Shannon mutual information (SMI) is well described by the following formula: \begin{equation}\label{eq:smiscaling} I(\ell,L)=\frac{b}{4}\ln \left(\frac{L}{\pi}\sin \frac{\pi \ell}{L}\right)+O(1), \end{equation} for a periodic system of length $L$ cut into two parts of respective lengths $\ell$ and $L-\ell$. The $O(1)$ term includes a non-universal constant as well as subleading corrections. This scaling form is almost identical to the celebrated EE result [26,8], upon substituting $c/3$ with $b/4$. In the following, we will refer to such a scaling as the conformal scaling. In an interesting recent development, Alcaraz and Rajabpour [2] further conjectured that the coefficient $b$ appearing in (\ref{eq:smiscaling}) is nothing but the central charge $c$ of the underlying CFT in general. This conjecture was supported by exact diagonalizations in a variety of spin chains corresponding to different universality classes, as well as on an exact correspondence with the second Rényi entanglement entropy for a particular model of harmonic oscillators with $c=1$. This conjecture, while very reasonable, is intriguing for two reasons. First, there is a small mismatch (of the order of $2$ percent) between the numerical estimate of $b$ in Ref. [2] compared to Refs. [33,52,59], for different models and limits belonging to the Ising universality class. On the conceptual level also, the result appears much simpler than previous findings in the Shannon entropy of the full chain [53,56,62]. Indeed, highly nontrivial transitions in the Rényi entropy, \begin{equation} S_n=\frac{1}{1-n}\ln \left(\sum_\sigma \left[p_{\sigma}\right]^n\right), \end{equation} were observed as a function of the Rényi index $n$. For Ising, this transition occurs precisely at $n=1$, where $S_n$ reduces to the Shannon entropy $\lim_{n\to 1}S_n=S_1=S$. The above considerations motivate us to clarify this issue and revisit this problem, generalizing the study to that of the Rényi mutual information (RMI) \begin{equation}\label{eq:rmidef} I_n(A,B)=S_n(A)+S_n(B)-S_n(A\cup B) \end{equation} for general $n$. A unifying conclusion from our study will be that the RMI of periodic systems generically obeys the conformal scaling \begin{equation}\label{eq:general_rmi} I_n(\ell,L)=\frac{b_n}{4}\ln\left(\frac{L}{\pi}\sin \frac{\pi \ell}{L}\right)+O(1), \end{equation} where $b_n$ is a universal – and in general nontrivial – coefficient. We also explain why such a result should be, typically, independent of the choice of local basis. We derive Eq. (\ref{eq:general_rmi}) for the free boson CFT, exploiting the gaussian form of the action, and combining this with boundary CFT arguments. We show that $b_n$ is proportional to the central charge in this case. At $n=1$ we recover the numerical result of Ref. [2]. Crucial to our derivation is the gaussian nature of the action. Already for the Ising CFT such arguments do not apply anymore. To investigate the scaling of the RMI we performed extensive numerical simulations on the example of the XY chain in transverse field. We used the free fermion structure to access larger system sizes than considered in previous works, and confirm (\ref{eq:general_rmi}). We then argue for a formula proportional to the central charge when $n>1$. Using the very accurate free fermion data and combining with simple extrapolation techniques, we also managed to obtain \begin{equation} b=b_1=0.4801629(2) \end{equation} at the Shannon point. This differs from the central charge ($c=1/2$), and convincingly disproves the general conjecture of Ref. [2]. However the conformal scaling still holds. We do not know how to derive this from first principles, but conjecture that all minimal models give similar results to that of Ising, with (Shannon) prefactor close but not quite identical to the central charge. The paper is organized as follows. In Sec. 2. we study the path-integral representation of the Rényi mutual information. We derive the conformal scaling for the free boson CFT ($c=1$), and show that (\ref{eq:general_rmi}) holds (with $b_n=1$) for all Rényi entropies, provided $n$ is not too large. We also show how the phase transition scenario of Ref. [56,62] allows us to derive the scaling of the mutual information for the Ising universality class provided $n>1$. These predictions are checked numerically in Sec. 3. for two spin chains, the XXZ and XY chain in transverse field, where very good agreement is found. Sec. 4.2. focuses on a numerical extraction of the universal coefficient $b_n$ in free fermionic systems, as well as on the “transition” points where our analytical arguments do not apply. Some additional details are gathered in three appendices. The first shows how the conformal scaling can be derived from standard boundary CFT techniques (App. 7.). The second (App. 8.) shows some numerical computations of correlation functions, in support of our arguments for Ising. Finally, we perform in App. 9. some exact computations of Rényi entropies for integer $n$ in the full XX chain, one of the simplest models described by a free boson CFT.
Mutual information in CFT
In this section we study the Rényi mutual information for all Rényi indices $n$ in a CFT setup.
Infinite Renyi limit
Let us start by considering the case $n\to \infty$. This limit turns out to be simplest, and will be useful later on. The entropy is given by \begin{equation} S_\infty(A)=-\ln \left(\left\langle \textrm{max}|\rho_A|\textrm{max} \right\rangle\right). \end{equation} A similar expression holds for $S_\infty(B)$ and $S_{\infty}(A\cup B)$. Here $\left|\textrm{max}\right\rangle$ denotes the spin configuration(s) with the biggest probability(ies). In most spin chains these are usually attained by homogeneous states. For example in the (ferromagnetic) Ising chain in transverse field $\left|\textrm{max}\right\rangle=\left|\uparrow\ldots\uparrow\right\rangle$ or $\left|\textrm{max}\right\rangle=\left|\downarrow\ldots\downarrow\right\rangle$, and in the antiferromagnetic XXZ chain $\left|\textrm{ max}\right\rangle=\left|\uparrow\downarrow\ldots\uparrow\downarrow\right\rangle$ or $\left|\textrm{max}\right\rangle=\left|\downarrow\uparrow\ldots\downarrow\uparrow\right\rangle$. In a euclidean picture this probability is given by \begin{equation} \left\langle \textrm{max}|\rho_A|\textrm{max} \right\rangle=\lim_{\tau\to \infty}\frac{\left\langle a|e^{-\tau H} \delta (\sigma-\sigma_\textrm{max})e^{-\tau H}|a \right\rangle}{\left\langle a|e^{-2\tau H}|a \right\rangle} \end{equation} $H$ is the Hamiltonian of the system, $\left|a\right\rangle$ is a state at infinity that has non-zero overlap with the ground state, and the $\tau \to \infty$ limit ensures the projection onto the ground-state of the Hamiltonian $H$. $\delta(\sigma-\sigma_\textrm{max})$ selects the spin configuration(s) with maximum probability(ies) on the segment $\tau=0,0\leq x\leq \ell$. In a transfer matrix picture, this can be represented as the following ratio of partition functions \begin{equation} \left\langle \textrm{max}|\rho_A|\textrm{max} \right\rangle=\frac{\mathcal{Z}_{\textrm{slit}}}{\mathcal{Z}}. \end{equation} For a periodic (open) system $\mathcal{Z}_{\textrm{slit}}$ is the partition function of an infinite cylinder (strip) with a slit, as is shown in Fig. 1(a,b).
Let us also mention that subleading corrections can be computed in these two simple cases. The leading corrections [50,51] to (\ref{eq:efp_periodic},\ref{eq:efp_open}) take the form of a series in $\beta_k(\ln L)L^{-k}$, where $\beta_k(x)$ is a polynomial of degree at most $k$ in $x$.
The free boson case
This section is devoted to the free compact field with action \begin{equation}\label{eq:action} \mathcal{A}=\frac{g}{4\pi}\int_{0}^{L} dx \int_{-\infty}^{\infty} d\tau \left(\nabla \varphi\right)^2\quad,\quad \varphi\equiv \varphi+2\pi r. \end{equation} The field $\varphi$ is compactified on a circle of radius $r$. This is the euclidean time version of the Luttinger liquid CFT; the Luttinger parameter, that controls the decay of correlation functions, is $K=(2gr^2)^{-1}$. For a periodic (open) chain the field lives on a infinite cylinder (strip) of circumference (width) $L$. In the continuum limit, each spin configuration $\mu$ in subsystem $A$ is replaced by a corresponding field configuration \begin{equation} \phi(x)=\varphi(x,\tau=0)\qquad,\qquad 0 \leq x\leq \ell. \end{equation} The Renyi entropy is \begin{equation} \mathcal{S}_n=\frac{\ln \mathcal{Z}^{(n)}}{1-n}\qquad,\qquad \mathcal{Z}^{(n)}=\int \left[\mathcal{D}\phi\right] \left(p_g(\phi)\right)^n. \end{equation} We kept track of the stiffness $g$, for reasons that will become apparent shortly. The probability $p_g(\phi)$ can be evaluated following Ref. [57]. We decompose the field $\varphi$ into the sum of two terms. The first is a harmonic function $\varphi_\phi$ that satisfies the boundary condition \begin{equation} \varphi_\phi(x,\tau=0)=\phi(x)\qquad ,\forall x\in [0;\ell], \end{equation} and the second an oscillator part $\varphi_0$ that satisfies a Dirichlet boundary condition. We now exploit the gaussian nature of the action, as well as the fact that $\varphi_\phi$ has a vanishing laplacian, to get $\mathcal{A}[\varphi_\phi+\varphi_0]=\mathcal{A}[\varphi_\phi]+\mathcal{A}[\varphi_0]$. Hence \begin{equation}\label{eq:proba} p_g(\phi)=\exp(-S_g[\varphi_\phi])\frac{\mathcal{Z}_g^D}{\mathcal{Z}_g}. \end{equation} $\mathcal{Z}_g$ is the partition function of an infinite cylinder (resp. strip), and $\mathcal{Z}^D$ is the partition function of an infinite cylinder (resp. strip) with a Dirichlet defect line at $\tau=0$, $0\leq x\leq \ell$, just as in Fig. 1.
Now we raise $p_g(\phi)$ to the power $n$, to get \begin{equation}\label{eq:proban} \left[p_g(\phi)\right]^n=\exp(-n S_g[\varphi^\phi])\left(\frac{\mathcal{Z}_g^D}{\mathcal{Z}_g}\right)^n. \end{equation} The crucial point is that the exponential prefactor in (\ref{eq:proban}) can be interpreted as a Boltzmann factor in a system with stiffness $g^\prime=ng$. Hence we get \begin{equation}\label{eq:changedstiffness} [p_g(\varphi^\phi)]^n\propto p_{ng}(\varphi^\phi). \end{equation} The Rényi index $n$ therefore changes the stiffness to $g^\prime=ng$ (or alternatively, the Luttinger parameter to $K^\prime=K/n$) near the slit. Keeping track of the proportionality coefficients and using the normalization of the $p_{ng}$, we finally arrive at \begin{equation}\label{eq:smp} \mathcal{Z}^{(n)}=\frac{\mathcal{Z}_{ng}}{\mathcal{Z}_{ng}^D}\left(\frac{\mathcal{Z}_g^D}{\mathcal{Z}_g}\right)^n. \end{equation} This result generalizes Ref. [57] to an arbitrary subsytem of length $\ell$. A crucial property of the two cut-cylinder (a) and cut-strip (b) geometries we focus on is that they can be conformally mapped to the upper half-plane (see. 7.). The leading universal shape-dependent piece is exactly that given in 2.1., and the result only depends on the central charge $c$, not the stiffness $g$ 2 . Hence the stiffness dependence in the partition functions can be discarded, and we get \begin{equation}\label{eq:freebosonZreplicas} \mathcal{Z}^{(n)}=\left(\frac{\mathcal{Z}^D}{\mathcal{Z}}\right)^{n-1}, \end{equation} so that \begin{equation} \mathcal{S}_n=-\ln \left(\mathcal{Z}^D/\mathcal{Z}\right). \end{equation} This ratio is formally identical to the one studied in the previous section. Therefore we obtain \begin{equation}\label{eq:freeboson} I_n(\ell,L)=\frac{c}{4}\ln \left[\frac{L}{\pi}\sin \frac{\pi \ell}{L}\right]+O(1)\quad,\quad c=1 \end{equation} for a periodic system, and half that in a open system. The result (\ref{eq:freeboson}) should also apply to orbifolds of the free boson theory, as well as in the non compact limit $r\to\infty$. When $n=1$ our derivation recovers the numerical results of Ref. [2] for the $XXZ$ chain and the $Q=4$ state quantum Potts model. We emphasize that even though the central charge appears in Eq. (\ref{eq:freeboson}), the derivation relies crucially on Eq. (\ref{eq:smp}), which is specific to the free field ($c=1$). The above derivation implicitly assumes that the boundary is still critical at stiffness $g^\prime=ng$, and this may not necessarily be so. For example lattice effects may change the result if $n$ gets too large in a compact theory. These effects can be tackled by adding vertex operators to the action (\ref{eq:action}). The least irrelevant is $V_d=\cos( \frac{d}{r}\varphi)$, where $d$ is the smallest integer allowed by the lattice symmetries. It is irrelevant provided $d^2>2gr^2=K^{-1}$, which is the case since we are studying a critical system. However, in presence of a stiffness $ng$ this condition becomes $d^2>2ngr^2=nK^{-1}$, which can be rewritten as [57] \begin{equation}\label{eq:nc} n<n_c\qquad,\qquad n_c=K d^{\,2}. \end{equation} Therefore, when $n>n_c$ a phase transition takes place at the boundary, and the field gets locked into one of the $d$ minima of the cosine potential $V_{d}$. In this boundary-locked phase the universal contributions to the entropy are given by \begin{equation} S_n\sim \frac{1}{1-n}\ln \left[d (p_{max})^n\right], \end{equation} where $p_{\textrm{max}}=\left\langle \textrm{max}|\rho_A|\textrm{max} \right\rangle$. The universal shape-dependence becomes \begin{equation} \mathcal{S}_{n>n_c}=\frac{n}{n-1}\ln \left(\frac{\mathcal{Z}^D}{\mathcal{Z}}\right). \end{equation} Using the result established in Sec. 2.1., we get \begin{equation}\label{eq:freeboson_afternc} I_{n>n_c}(\ell,L)=\frac{n}{n-1}\frac{c}{4}\ln \left[\frac{L}{\pi}\sin \frac{\pi \ell}{L}\right]+O(1) \end{equation} for a periodic system, and once again half that in a open system. We recover the result Eq. (\ref{eq:Iinfty}) in the limit $n\to \infty$. Interestingly the conformal scaling holds in both phases ($n<n_c$) and ($n>n_c$), but with different prefactors proportional to $c$. From our arguments however, it is not obvious how the mutual information behaves at the transition point $n=n_c$, where the vertex operator is marginal. We will come back to this point in Sec. 4.2..
Path-integral and replicas
For more general models or CFTs it is not so easy to use the explicit form of the action as in the Luttinger liquid. The standard method is to introduce $n$ replicas, with $n$ an integer greater than one. We have \begin{equation}\label{eq:replicas} \mathcal{Z}^{(n)}=\frac{\mathcal{Z}_{\textrm{rep}}}{\mathcal{Z}^n}, \end{equation} where $\mathcal{Z}$ is the partition function of the infinite cylinder (strip). $\mathcal{Z}_{\textrm{rep}}$ is the partition function of a replicated system: we have $n$ independent copies of the cut-cylinder (cut-strip) geometries shown in Fig. 1, but where the configurations of all copies are identified along the cut. Note that for the Shannon entropy of the full chain this geometry can be folded, and seen as a system of $2n$ semi-infinite cylinders (strips) glued along their common boundary. It is also instructive to derive the free boson result before the transition ($n<n_c$), using replicas as in Ref. [21]. We introduce $n$ copies $\varphi_i$, for $i=1,\ldots,n$ of the boson field, with total action \begin{equation} \mathbf{A}=\sum_{i=1}^n \mathcal{A}_i\quad,\quad \mathcal{A}_i=\frac{g}{4\pi}\int_{0}^{L} dx_i \int_{-\infty}^{\infty} d\tau_i \left(\nabla \varphi_i\right)^2. \end{equation} All these fields are independent in the bulk, but they are stitched together at the slit. Neglecting the ambiguities in the compactification that do not matter at the leading order in this geometry (see the discussion after (\ref{eq:smp})), one can introduce the orthogonal transformation [21] $\tilde{\varphi}_1=\frac{1}{\sqrt{n}}\sum_i \varphi_i$ as well as $\tilde{\varphi}_i=\frac{1}{\sqrt{2}}\left(\varphi_{i}-\varphi_{i-1}\right)$ for $i=2,\ldots,n$. In terms of the new fields the action decouples in the bulk, $\mathbf{A}=\sum_i \tilde{\mathcal{A}}_i$. Since all the fields $\varphi_{i\geq 2}$ have to match on the slit, the new fields $\tilde{\varphi}_i$ vanish on it, hence they obey a Dirichlet boundary condition. The remaining field $\tilde{\varphi}_1$ fluctuates freely at the slit and cancels with one of the $n$ normalization partition functions $\mathcal{Z}$. In the end we recover Eq. (\ref{eq:freebosonZreplicas}), and the result (\ref{eq:freeboson}) follows. For the Ising CFT, such a gluing of CFT already becomes nontrivial. As can be seen in a Landau-Ginzburg (“$\phi^4$”) point of view (see e.g. Ref. [45]), the total bulk action does not decouple in terms of similar orthogonal transformation of the fields. Any unitary minimal model has a similar representation, and would suffer from this problem. However, previous numerical simulations of the Rényi entropy of the full chain [53,56,62] have found results consistent with a boundary transition already at $n=n_c=1$. This implies that the entropy should be dominated by ordered configurations for any $n>1$, similar to what happens in the free boson case for $n>n_c$. Hence the copies decouple, and we expect the RMI to be given by the simple formula \begin{equation}\label{eq:rmi_conj} I_{n>1}(\ell,L)=\frac{c}{4}\frac{n}{n-1}\ln \left[\frac{L}{\pi}\sin \frac{\pi \ell}{L}\right]. \end{equation} Although it can be justified for large $n$ [10], it is not obvious from perturbative RG arguments why this should happen already for $n>1$ 3 . The couplings near the slit are enhanced in the replica picture, or equivalently the temperature is lowered. One would naively expect these to be responsible for the ordering of the boundary. However, a typical scenario for a boundary in 2d classical critical systems with positive local degrees of freedom is that of the ordinary transition, where a change of couplings near a boundary does not order it [15]. Here the gluing of $n$ copies appears to escape this scenario, and the slit on which the copies are stiched does order. In surface critical phenomena language this is an extraordinary transition [15]. We have at present no analytical understanding of this observation. Similar gluing of 2d Ising models have already been studied in the literature (see e.g. [29,58]), but in slightly different limits. Let us finally mention that universal scaling forms proportional to $n/(n-1)$ have been found in the 2d classical Rényi mutual information [54] for Ising, as well as in the Rényi entropy of the 2d quantum transverse field Ising model [36]. In both cases the underlying ordering assumption is easier to justify: for the former the critical system is coupled to a bulk in the ordered phase [28], while for the latter the higher dimensionality makes an extraordinary transition more likely. To confirm our ordering assumption, we have computed numerically the spin-spin correlation function along the slit in the replicated geometry corresponding to $\mathcal{Z}^{(n)}$ ($n=1$ is the usual bulk spin-spin correlation function). For $n>1$ the data shows that this correlation becomes ordered, and the copies effectively decouple. We refer to Appendix 8. for the details. Such a behavior is also strikingly different from the free boson case, where the boundary is still critical with a modified stiffness $g^\prime=ng$, see Eq. (\ref{eq:changedstiffness}). We therefore still expect critical correlations for $n<n_c$ for the free boson, and ordered ones for $n>n_c$. This has already been shown [32] numerically in lattice discretizations such as the XXZ spin chain. We present numerical results for the mutual information that support our conjecture (\ref{eq:rmi_conj}) in Sec. 3.. As in the free boson case, our arguments do not predict the behavior of the mutual information exactly at the transition point ($n=1$ here). This study is deferred to Sec. 4., where we rely on numerics.
Numerical checks of the CFT
In this section we perform various numerical checks of the conformal scaling in lattice discretizations of CFTs with central charge $c=1$ and $c=1/2$.
Lattice computations
We focus here on two spin chains. The first is the antiferromagnetic XXZ chain \begin{equation} H=\sum_{i=1}^L \left(\sigma_i^x\sigma_{i+1}^x+\sigma_i^y\sigma_{i+1}^y+\Delta \sigma_i^z \sigma_{i+1}^z\right), \end{equation} which is (for $-1<\Delta\leq 1$) a Luttinger liquid CFT ($c=1$). Here the basis generated by the eigenstates of the $\sigma_j^z$ is most natural, as the bosonic field is diagonal in such a basis [23]. Even though the chain is integrable, it is in practice not easy to exploit the additional structure to compute the Shannon entropy (see however Refs. [7,46] for an explicit computation of $S_{\infty}(L,L)=-\ln p_{\uparrow\downarrow\ldots\uparrow\downarrow}$ in the periodic case). Here we generate the ground-state wave function using the Lanczos algorithm. For a periodic system one can use translational invariance, the reflection symmetry, and combine these with the $U(1)$ and particle-hole symmetry. Doing so system sizes up to $L=36$ can be reached with a reasonable amount of computer effort.
The other is the quantum XY chain in transverse field: \begin{equation} H=-\sum_i\left[ \left(\frac{1+\gamma}{2}\right)\sigma_i^x\sigma_{i+1}^x+\left(\frac{1-\gamma}{2}\right)\sigma_i^y\sigma_{i+1}^y+h\sigma_i^z\right] \end{equation} For $h=1$ and $\gamma>0$, its long distance behavior is described by the Ising $c=1/2$ CFT ($h=1$ and $\gamma=1$ is the Ising chain in transverse field –ICTF). Note that for $\gamma=0,h=0$ we recover the XXZ model at $\Delta=0$. The main advantage of this chain is that it can be written in terms of free fermions, and numerical computations are somewhat simplified. Indeed, performing a Jordan-Wigner transformation \begin{eqnarray}\label{eq:jw1} \sigma_j^z&=&2c_j^\dagger c_j-1,\\ \label{eq:jw2} \frac{\sigma_j^x+i\sigma_j^y}{2}&=& \exp\left(i\pi \sum_{l=1}^{j-1}c_l^\dagger c_l\right)c_j^\dagger, \end{eqnarray} allows to express $H$ as a quadratic form in the $c_i,c_i^\dagger$, which may be diagonalized by a Bogoliubov transformation
. Each probability in subsystem $A$ can be obtained exactly, following e.g. Ref. [56]. We have \begin{equation}\label{eq:probadet} p_{\mu}=\det_{1\leq i,j\leq \ell}\left(\,\left\langle f_i^\dagger \left(f_j^\dagger+f_j\right)\right\rangle_{\!L}\,\right), \end{equation} where $\langle \ldots\rangle_L$ denotes the average in the ground-state of the chain of size $L$, and $f_i^\dagger=c_i^\dagger$ (resp. $f_i^\dagger=c_i$) if $\sigma_i^z=\uparrow$ (resp. $\sigma_i^z=\downarrow$). There are analogous formulae for systems $B$ and $A\cup B$, respectively as $(L-\ell)\times(L-\ell)$ and $L\times L$ determinants. Using this, the entropy in the $\sigma^z$ basis can be obtained by brute-force summation over the $2^\ell$ configurations, with total complexity $\ell^3 2^\ell$. The number of determinants to compute can be slightly reduced by making use of the lattice symmetries. Since each of the probabilities can be computed independently, parallelization is also trivial. In practice we can push the numerics up to $\ell \simeq 40$ independent on $L$, which significantly improves on Ref. [2], another advantage being that the data is exact up to machine-precision. For the Ising chain the most natural local basis is generated by the $\sigma_j^x$, as these correspond to the actual spins in the classical two-dimensional model. In a periodic chain, it is also possible[53] to obtain the entropy in this basis, using the Kramers-Wannier duality $\sigma_j^z\to\tilde{\sigma}_{j-1}^x\tilde{\sigma}_j^x$, $\sigma_j^x\sigma_{j+1}^x\to \tilde{\sigma}_j^z$, which maps $H$ onto itself at the critical point. We get \begin{equation} p(\sigma_1^x,\ldots,\sigma_L^x)=\frac{1}{2}p(\tilde{\sigma}_1^z,\ldots,\tilde{\sigma}_L^z), \end{equation} and evaluate $p(\tilde{\sigma}_1^z,\ldots,\tilde{\sigma}_L^z)$ using Eq. (\ref{eq:probadet}). The factor $1/2$ is due to the fact that the mapping is two to one. The complexity for the subsystem entropy is slightly greater ($2^L L^3$) than in the $z$ basis, as one needs to generate all the probabilities of the full chain to access the subsystem (marginal) probabilities.
For open chains the Kramers-Wannier duality does not map $H$ onto itself, and this trick does not work anymore. The free fermions method does not outperform the Lanczos algorithm in this case.
Results for the free boson
We first focus our attention on the periodic XXZ chain in the $z$ basis. In this chain the least irrelevant vertex operator allowed by the lattice symmetries is $V_2=\cos \left(\frac{2}{r}\varphi\right)$, so that (see Eq. (\ref{eq:nc})) \begin{equation} n_c=4K. \end{equation} The Luttinger parameter is known to be [23] \begin{equation} K=\left(2-\frac{2}{\pi}\arccos \Delta\right)^{-1}. \end{equation} We test Eq. (\ref{eq:freeboson}) for different values of $n$ and $\Delta$ before the transition ($n<n_c$), and plot the ratio \begin{equation}\label{eq:cestimate} \frac{I_n(\ell,L)-I_n(L/2,L)}{\frac{1}{4}\ln \sin \frac{\pi \ell}{L}} \end{equation} in Fig. 2. Provided (\ref{eq:freeboson}) is correct, this should give $c$ for all $\ell/L$. The results are, within at worst a few percents, compatible with this. Note that in principle the conformal limit is reached only when $\ell\gg 1$, so that we expect the result to deteriorate at small aspect ratios $\ell/L$. Finite-size effects are bigger when $\Delta<0$ and $n<1$, or $\Delta>0$ and $n>1$. However, as illustrated in the inset, the accuracy improves as we increase the system size, so we expect a slow convergence towards $c=1$.
Results for the Ising CFT
As an example of $n>1$ simulation in the Ising universality class, we computed $I_2(\ell,L)$ in the Ising chain in transverse field (ICTF). The simulations were performed both in the $x-$ and $z-$ basis. Assuming the transition argument of Sec 2.3., the CFT predicts a scaling \begin{equation}\label{eq:I2Ising} I_2(\ell,L)=\frac{c}{2}\ln \left(\frac{L}{\pi}\sin \frac{\pi \ell}{L}\right)+O(1). \end{equation} In Fig. 3 we show $I_2(\ell,L)-I_2(L/2,L)$ for several system sizes, and compare it to the CFT result. As can be seen the agreement is excellent, and improves as $L$ gets larger. In the inset we also perform a central charge extraction similar to that done in Fig. 2. We plot the ratio \begin{equation} \frac{I_2(\ell,L)-I_2(L/2,L)}{\frac{1}{2}\ln \sin \frac{\pi \ell}{L}} \end{equation} as a function of $\ell/L$. The central charge extracted from this method gives an excellent agreement with $c=1/2$ in the $z-$basis. In the $x-$ basis there are slightly bigger finite-size effects, and the numerical data for our largest system size agrees up to $1.5$ percents. However the trend towards $c=1/2$ is clear. All these numerical results nicely agree with our predictions, whichever the local basis ($x$ or $z$) we choose. This basis independence can be justified in the following way, assuming the phase transition scenario. For $n>1$ (universal terms in) the entropy will be dominated by the ordered configurations. In the $x$-basis these are $\left|\uparrow\ldots \uparrow\right\rangle_x$ and $\left|\downarrow\ldots\downarrow\right\rangle_x$. Since these correspond to the actual classical spin configuration in the 2d Ising model, the expected CIBC is the fixed boundary condition. In the $z-$ basis only one configuration dominates, namely $\left|\uparrow\ldots\uparrow\right\rangle_z$. In the classical model this corresponds to a superposition of all possible spin configurations. Hence the appropriate CIBC is the free boundary condition. However we have seen in Sec. 2.1. that the conformal scaling is insensitive to the precise nature of the (conformal) boundary conditions, and this justfies why (\ref{eq:I2Ising}) holds in both basis. We have also checked that the shape function is half the periodic one in open chains, as predicted by CFT.
Transitions and transition points
The aim of this section is to perform a detailed study of the transitions happening as a function of the Rényi parameter $n$ for the Ising and free boson universality class. In particular, we try and extract the universal prefactor $b_n$ with the highest accuracy possible. Recall we expect the mutual information to behave as \begin{equation}\label{eq:RMIscaling2} I_n(\ell,L)=\frac{b_n}{4}\ln\left[\frac{L}{\pi}\sin\left(\frac{\pi \ell}{L}\right)\right]. \end{equation} There are several different methods to extract $b_n$ from the numerical data. The main two we use are the following. Since the results of Figs. 2,3 appear most accurate near $\ell=L/2$, we first look at a discrete derivative of the mutual information around this point \begin{equation}\label{eq:deltadef} \delta_n(L)=\frac{2L^2}{\pi^2}\left[I_n(L/2,L)-I_n(L/2+2,L)\right] \end{equation} for large $L$. We consider only total system sizes $L$ multiple of four, so that all sizes are even, and potential parity effects are avoided. Assuming (\ref{eq:RMIscaling2}), $\delta_n(L)$ should scale as \begin{equation}\label{eq:derivative} \delta_n(L)=b_n+o(L^0), \end{equation} and so can be used to extract $b_n$. In case the structure of the subleading corrections to (\ref{eq:derivative}) can be determined, extremely accurate values of $b_n$ can be extracted. For example including corrections of the form $\sum_{p=1}^3\alpha_p L^{-p}$ in (\ref{eq:derivative}) yields $b_{1/2}=c=1.001(2)$ for the worst data set ($n=1.5,\Delta=0.4$) in Fig. 2. We call this method the “discrete derivative method”. An alternative procedure is to look at a finite subsystem of size $\ell$ in an infinite system ($L\to \infty$) [52]. The Rényi entropy is given in this case by \begin{equation}\label{eq:infscaling} S_n(\ell,L\to \infty)=a_n\ell+\frac{b_n}{8}\ln \ell+\mathcal{O}(1), \end{equation} where $a_n$ is a line free energy. Numerical computations for free fermionic systems simulated in $z$-basis are possible, as is explained in Sec 3.1.. The logarithm is subleading, however the linear term can be substracted off by studying $S_n(\ell,\infty)-S_n(\ell,\ell)$, because $a_n$ does not depend on the precise geometry. As in the previous method, the results can be improved by extrapolation. In the following we will use both methods to extract $b_n$; their respective merits depend on the physical systems considered, and on our ability to guess the correct structure of the subleading corrections.
Transitions
Let us now demonstrate our results for the universal coefficient $b_n$. From the arguments presented above, we expect a phase transition in both the XX (4.1.1.) and Ising (4.1.3.) chains.
XX chain
We start with the XX chain, as an example of the free boson described in Sec. 2.2.. Our arguments predict a transition at $n=n_c=4$. The expected result from CFT is \begin{equation}\label{eq:xx_bn} b_n=\left\{ \begin{array}{ccc} c&,&0<n<4\\ \\ \frac{n}{n-1}c&,&n>4 \end{array} \right. \end{equation} with a central charge $c=1$.
Ising chain
We now turn our attention to the ICTF, and perform the same extraction of $b_n$. As in the previous section, the entropy is studied both in the $z$-basis and $x$ basis (where the infinite method is not available). The numerical results are shown in Fig. 5.
The transition points
The exact transition points $n=n_c$ are potentially the most interesting, and seem nontrivial from a boundary CFT perspective. It is not even guarantied that the conformal scaling survives. As we shall see the numerical results suggest that it does in periodic chains, albeit with a nontrivial value of $b_{n_c}$.
XX chain
We first start with the XX chain. We study the Rényi mutual information for $n=n_c=4$. The numerical results for several system sizes are shown in Fig. 6.
Ising chain
We now finally arrive at the Shannon $n=1$ point in the ICTF, where $b_1=c=1/2$ has been conjectured [2]. The numerical results in the $x$ and $z$ basis are shown in Fig. 7.
The peculiar case of open chains
It is also tempting to look at the shape dependence of the mutual information with open boundary conditions. In systems or Rényi indices where the CFT arguments of Sec 2. do apply, the mutual information is expected to be half that of a periodic system. Our numerical results show that this is not true for Ising at $n=1$. The numerical simulations clearly suggest the following form for the Shannon entropy \begin{equation}\label{eq:shannonopen_scaling} S_1(\ell,L)=a\ell+b^\prime\left(\ln \ell\right)^2+f(\ell/L)\ln \ell+O(1), \end{equation} which implies \begin{equation}\label{eq:smi_open} I_1(\ell,L)=\left[f(\ell/L)+f(1-\ell/L)\right]\ln \ell+O(1). \end{equation} This unusual scaling, with a shape-dependent logarithmic term, is very different from what we have seen before, and defies our simple CFT arguments. In practice, it is difficult to extract the shape-dependent $\ln \ell$ term in the mutual information. This is because an extrapolation with fixed aspect ratio $\ell/L$ is required, and due to commensurability effects this may require unreachable system sizes. However, we are able to extract $f(\ell/L)$ from the Shannon subsystem entropy in the $z$ basis, where slightly bigger system sizes are available. In practice we extrapolate $S_1(\ell,L)$ to $a\ell+b^\prime (\ln \ell)^2+\beta_0\ln \ell+\alpha_0+\beta_1\ell^{-1}\ln \ell+\alpha_1\ell^{-1}$ for fixed aspect ratio $\ell/L$. The $\ell^{-1}\ln \ell$ is a natural correction here, as can be seen by making the substitution $\ell\to \ell+\epsilon$ for some ultraviolet cutoff $\epsilon$, and expanding again in $\ell$. Using this method, we found very stable values of $f(\ell/L)$, which is good evidence for the correctness of our scaling ansatz. The results are shown in Fig. 9 for three different values of $\gamma$ in the XY chain, and strongly suggest that the shape function $f(\ell/L)-f(\ell/L=1/2)$ is universal.
Conclusion
We have studied in this paper the Rényi and Shannon mutual information of one-dimensional quantum critical systems, using a combination of CFT and numerical techniques. Our main result is that the RMI for $n\geq 1$ follows the simple conformal scaling formula \begin{equation} I_n(\ell,L)=\frac{b_n}{4}\ln \left[\frac{L}{\pi}\sin \left(\frac{\pi \ell}{L}\right)\right], \end{equation} for a periodic system cut into two subsystems of sizes $\ell$ and $L-\ell$. We have derived this result from CFT in the case of the free boson, and shown that $b_n$ is proportional to the central charge, except at a special point $n=n_c$ where it is unknown. At $n=1$ our analytical formula recovers the numerical result of Ref. [2]. The Ising universality class proved more tricky, as we were unable to formally derive the conformal scaling. However, inspired by the phase transitions observed in [56,62], we showed numerically that the system orders at the slit in a replica picture for $n>1$. This implies the conformal scaling, and the result $b_n=c\frac{n}{n-1}$, independent of the local ($x$ or $z$) spin basis. We checked this prediction numerically, and found very good agreement. By construction this argument breaks down when $n\to1$. However the conformal scaling survives, with a value for $b_1$ that we computed numerically (see Eq. (\ref{eq:b_estimate})). The Shannon point proves even more mysterious in open chains, where the full Shannon entropy has a $\left(\ln L\right)^2$ contribution, and the SMI diverges logarithmically with a shape-dependent prefactor (\ref{eq:smi_open}). We have currently no analytical understanding of these observations, and it would be highly desirable to make progress on this vexing problem. The effect of having a Rényi index $n\neq 1$ is intuitively similar to a change of temperature in the 2d classical model, as can be seen from the difference between the XXZ/six-vertex and the Ising results. Indeed, the former models possess a line of critical points, while the latter have a critical point that separates two gapped phases. All the competing orders are present in the (critical) ground-state wavefunction, and changing $n$ allows to distinguich between them. Such an interpretation should not be taken too litteraly, as is discussed in Sec. 2.3.. It would even give a wrong $n_c$ for the Luttinger liquid, by neglecting the fact that the transition is a boundary transition, not a bulk transition. Note however that $n$ is a true inverse temperature if we interpret the $p_i$ as Boltzmann weights for a 1d classical chain. In this case transitions at finite temperature are possible because the interactions are long-range (see Appendix. 9.2. for an illustration in the XX chain). Finally and due to the correspondence with the Rényi entanglement entropy of certain Rokhsar-Kivelson states [53], these transitions also give us a lot of information about the entanglement spectrum in such states [13,56]. From our results for the Ising universality class it is tempting to conjecture a similar behavior for all minimal models (or any subset that closes under fusion rules). For example, Ref. [2] found $b_1 \simeq 0.79$ in the $3$-state quantum Potts model ($c=4/5=0.8$), and our findings indicate that this probably differs slightly from the central charge. Interestingly all these numbers should be universal: even though their exact value is not known exactly, they can still be used to identify the universality class. The fact that they are so close to the actual central charge is also very intriguing. Another interesting project would be to study the possible transitions as a function of $n$. For more complicated models the set of allowed conformal boundary conditions becomes larger, and one could imagine several transitions in the $b_n$ exponent or the Rényi entropy, with the conformal scaling still holding. Another interesting direction would be to study the Rényi entropy when $n<1$ for non free bosonic theories. For such values the replica approach brings some results for the entropy of the full chain. For example $Z_{1/2}$ becomes the partition function of a single half-sheet and boundary CFT applies [56]. This is not the case anymore for the subsystem entropy, and one needs to rely on other methods. Even the numerical results show big finite-size effects, and it is difficult to distinguish between a conformal scaling with strong subleading corrections and a slightly enhanced logarithmic scaling.
Acknowledgments
I wish to thank Fabien Alet, Bertrand Berche, Jérôme Dubail, Paul Fendley, Grégoire Misguich and Vincent Pasquier for several stimulating discussions. I especially thank John Cardy for sharing some of his notes on the subject with me. The simulations were performed on the UVACSE ITS cluster. This work was supported by the US National Science Foundation under the grant DMR/MPS1006549.
Appendix
Conformal mappings
In this appendix we recall the derivation of the conformal scaling of the emptiness formation probability in the periodic and open geometries (see Fig. 1 (a) and (b)). We follow the method of Ref. [12].
Cylinder with a slit
We start with the cylinder with a slit, relevant to periodic systems.
Strip with a slit
Correlation functions for Ising in the replica systems
The aim of this appendix is to show some numerical results for spin-spin correlations at the slit in a replica picture. We focus on the Ising case, and show that the correlations are ordered when $n>1$. For convenience we consider the periodic case with $\ell=L$, relevant to the total Rényi entropy $S_n(L,L)$. However the ordering can also be observed for any $\ell/L$. Since $\ell=L$ we can fold the system, so that we are studying correlations at the binding of a $2n$-sheeted “book”, where each of the sheets are independent except at the binding. To compute them we use the hamiltonian limit. The ground-state can be written, in the spin basis, as \begin{equation} \left|\psi\right\rangle=\sum_{\sigma} \psi_\sigma \left|\sigma\right\rangle \end{equation} The correlation we look at is \begin{equation} C_n(\ell,L)=\left\langle \sigma_0^x\sigma_\ell^x \right\rangle_{(n)}, \end{equation} where $\left\langle \ldots \right\rangle_{(n)}$ denotes the average in the “rényified” ground state \begin{equation} \left|\psi\right\rangle_{(n)}=\frac{1}{\sqrt{Z_n}}\sum_{\sigma}\left(\psi_\sigma\right)^n \left|\sigma\right\rangle. \end{equation} Recall that the $\sigma^x$ correspond to the classical Ising spins in the 2d model, so we use this basis for the numerical computations. For critical correlations conformal invariance implies \begin{equation}\label{eq:corr_conf} C_n(\ell,L)\sim \left|\frac{L}{\pi}\sin \frac{\pi \ell}{L}\right|^{-\alpha_n}, \end{equation} where $\alpha_n$ is the critical exponent. There are three values of $n$ for which the exponent can easily be determined:
- The book with one sheet ($n=1/2$). The correlation is that of spins living at the boundary of a semi-infinite cylinder with free boundary condition. In this case the exponent is $\alpha_{1/2}=1$.
- The book with two sheets ($n=1$). This corresponds to the usual ground-state correlations, with the well known Onsager exponent $\alpha_1=1/4$.
- The book with an infinite number of sheets ($n\to\infty$). Here the rényified ground-state is dominated by the two ordered states, and $\alpha_\infty=0$.
Some exact results for the Rényi entropy of the full $XX$ chain
This last appendix is devoted to the study of the Rényi entropy $S_n(L,L)$ of the full XX chain. In this simple model and geometry the Rényi entropy is related to a partition function for a 2d gas of particles on a ring, a discrete analog of the Dyson gas [17]. We use this connection to derive several exact results for integer $n$ Rényi indices. These are in agreement with the arguments presented in the main text, and offer supplementary evidence to the phase transition scenario for the free boson CFT. The section is organized as follows. We start in 9.1. by presenting some determinant and combinatorial identities that will be useful throughout. The connection between the entropy and the discrete gas is explained in 9.2., and a summary of the results in given in 9.3.. Some of the technical details of the computations are finally gathered in 9.4..
Determinants and constant term identities
Vandermonde determinants
Here we start by presenting some useful determinant identities. The first is the Vandermonde determinant, given by \begin{eqnarray}\label{eq:vand} V(u_1,u_2,\ldots,u_N)&=& \det_{1\leq j,k\leq N} \left(u_j^{k-1}\right)\\ &=&\prod_{1\leq j<k\leq N} (u_j-u_k). \end{eqnarray} A similar determinant, sometimes called symplectic Vandermonde, describes the ground-state of the open XX chain. \begin{eqnarray}\label{eq:svand} W(u_1,\ldots ,u_N)&=& \det_{1\leq j,k\leq N}\left(u_j^k-u_j^{-k}\right)\\ \nonumber &=&\prod_{j=1}^N u_j^{-N}(1-u_j^2)\prod_{k>j}(u_j-u_k)(1-u_ju_k) \end{eqnarray} Certain powers of $V$ and $W$ can also be expressed as determinants. For example the fourth power of $V$ is related to a $2N\times 2N$ determinant: \begin{equation}\label{eq:cvand} V(x_1,\ldots,x_N)^4=\det_{\begin{array}{c} \scriptstyle 1\leq j\leq N \\ \scriptstyle 1\leq k\leq 2N\end{array}} \left(\begin{array}{c} u_j^{k-1}\\ (k-1)u_j^{k-2} \end{array} \right). \end{equation} (\ref{eq:cvand}) can obtained by considering the limit \begin{equation} \lim_{v_1\to u_1,\ldots,v_N\to u_N} \;\frac{V(u_1,v_1,\ldots,u_N,v_N)}{(u_1-v_1)\ldots(u_N-v_N)}, \end{equation} and performing elementary row manipulation that leave the determinant invariant. A similar trick can be used on $W$. We get \begin{equation}\label{eq:csvand} \det_{\begin{array}{c} \scriptstyle 1\leq j\leq N \\ \scriptstyle 1\leq k\leq 2N\end{array}} \left(\begin{array}{c} u_j^{k}-u_j^{-k}\\ k u_j^{k-1}+k u_j^{-k-1} \end{array} \right)=\frac{W(u_1,\ldots,u_N)^4}{\prod_{j=1}^{N}(u_j^2-1)}. \end{equation}
Dyson and Macdonald constant term identities
Many of the partition functions shown in this appendix will follow from a constant term identity due to MacDonald [38]: \begin{equation}\label{eq:macdonald} \textrm{CT} \left[\prod_{\alpha\in R}\left(1-e^\alpha\right)^k\right]=\prod_{i=1}^{N} \left(\begin{array}{c}kd_i \\ k\end{array}\right) \end{equation} where the product runs over vectors of $\mathbb{R}^N$, elements of a crystallographic root system. $\textrm{CT}$ stands for the constant term in the expansion of the product. The $d_i$ are a set of integers characteristic of the root system[38]. The simplest (infinite) family of root systems is given by the series \begin{equation} A_{N-1}=\{\pm\left(t_i-t_j\right)\;,1\leq i< j\leq N\}, \end{equation} for which $d_i=i$ and (\ref{eq:macdonald}) reduces to the better known Dyson constant term identity [17,24] \begin{equation} \textrm{CT} \left[\prod_{j\neq i}^{N} \left(1-e^{t_i-t_j}\right)^{n}. \right]=\frac{(nN)!}{n!^{N}}. \end{equation} In the following, we will also get the $C_N$ root systems, given by \begin{equation} C_N=\{\pm \left(t_i\pm t_j\right),1\leq i<j\leq N\}\cup\{\pm 2t_i,1\leq i\leq N\}. \end{equation} In this case, (\ref{eq:macdonald}) holds with \begin{equation}\label{eq:di} d_i=2i. \end{equation}
Rényi-Shannon entropy and Dyson-Gaudin gas
We now establish the relation between the entropy and the thermodynamics of a discrete log-gas. The XX Hamiltonian we consider is given by \begin{equation} H=\sum_{j=1}^{L-1}\left(\sigma_j^x\sigma_{j+1}^x+\sigma_j^y\sigma_{j+1}^y\right)+H_{bound}, \end{equation} where $H_{bound}$ encodes the boundary conditions. In the following we wish to study both periodic boundary conditions ($H_{bound}=\sigma_L^x\sigma_1^x+\sigma_L^y\sigma_1^y$) and open boundary conditions ($H_{bound}=0$). In the latter case the relevant external conformal boundary condition is Dirichlet. To illustrate the effect of boundary changing operators, we also considered a case where the external boundary condition is Neuman. This can be achieved by applying any finite magnetic field along $x$ at the two boundary spins[1]. Here we focus on $H_{bound}=\sqrt{2}(\sigma_1^x+\sigma_L^x)$: the amplitude $\sqrt{2}$ (instead of a general $h$) is chosen because the diagonalization of $H$ simplifies considerably at this special point[6,61]. An important difference with the first two cases is that the total magnetization $M=\sum_{j=1}^L \sigma_j^z$ is not conserved anymore.
Summary of the results
We have seen (see above) that the Rényi entropy for three different boundary conditions was given by \begin{equation} S_n=\frac{\ln Z_n}{1-n}, \end{equation} where $Z_n$ was the partition function of a Dyson-Gaudin gas with certain symmetries (see Fig. ref.{C\tmspace +\thinmuskip {.1667em}2}(a,b,c)). This correspondence can be exploited to derive exact results in the XX chain. We present them for an XX chain of length $L$, assuming –except for the open Neuman case– a fixed filling fraction $\rho=N/L$, where $N$ is the number of particles. We also implicitly assume $\rho \leq 1/2$ throughout. The results for $\rho>1/2$ can simply be deduced from the particle-hole symmetry $N\to L-N$.
Periodic
For any integer Rényi index $n\leq \rho^{-1}$ we have \begin{equation}\label{eq:gaudin} Z_n^{per}(N,L)=\frac{(Nn)!}{N! L^{N(n-1)}n!^N}. \end{equation} This result is known [22] and was exploited in Ref. [53]. Note that all these values belong to the “replica” region $n<n_c=\rho^{-2}$. When $n=2$ we also obtained the spin-spin correlation functions. Using the notations of App. 8., we have for $N\leq L/2$ \begin{eqnarray}\nonumber \left\langle \sigma_j^z\sigma_{j+l}^z \right\rangle_{\!(2)}&=&\frac{2N\cos \frac{2N\pi l}{L}-\cot \frac{\pi l}{L}\sin \frac{2N\pi l}{L}}{L^2\sin \frac{\pi l}{L}}g_N\!\left(\textstyle{\frac{\pi l}{L}}\right) \\ &-&\left[\frac{\sin \frac{2N\pi l}{L}}{L\sin \frac{\pi l}{L}}\right]^2, \label{eq:corrxx} \end{eqnarray} where \begin{equation} g_N(\theta)=\sum_{k=1}^{N} \frac{\sin\left[ \left(2k-1\right)\theta\right]}{k-1/2}. \end{equation} When $N=L/2$ the result is known [39], as $\left|\psi\right\rangle_{(2)}$ coincides at half-filling with the ground-state of the Haldane-Shastry (HS) chain[25,48]. Hence $Z_2$ is also the norm the (unnormalized) HS ground-state. In the thermodynamic limit the sum converges to $\pi/2$ for any $l/L$ and $N/L$, so that $\left\langle \sigma_j^z\sigma_{j+l}^z \right\rangle\approx (L\sin\frac{\pi l}{L})^{-1}$. We recover the conformal scaling with a nontrivial exponent $\alpha_2=1$. For a Luttinger liquid the usual one is $\alpha=2K$, and $K=1$ in the XX chain. The exponent $\alpha_2$ we find here is consistent with the image of a modified Luttinger parameter $K^\prime= K/n=K/2=1/2$ near the binding in a replica picture (see Eq. (\ref{eq:changedstiffness}) and Ref. [32]). We were also able to derive a new result for $n=4$ at half-filling: \begin{equation}\label{eq:periodicspecial} Z_4^{per}(L/2,L)=\frac{2^{L/2}}{L^L}\left(\frac{\Gamma\left(\frac{L}{2}+\frac{1}{2}\right)}{\Gamma\left(\frac{1}{2}\right)}\right)^2. \end{equation} This is interesting because it corresponds to the critical value of the Rényi index. Unfortunately, we did not manage to compute correlation functions at this point. Let us also mention \begin{equation}\label{eq:sinfty} S_{\infty}=(\rho \ln 2) L \end{equation} for $\rho\leq 1/2$. All subleading (constant) terms in (\ref{eq:gaudin}, \ref{eq:periodicspecial}, \ref{eq:sinfty}) can be computed, and agree with the predictions of Refs. [53,56].
Open geometry
Similar results can be derived in the open geometry. When $n\leq \rho^{-1}$ we have \begin{equation}\label{eq:openmacdonald} Z_n^{open}(N,L)=\left(\frac{n^{n}}{n!(L+1)^{n-1}}\right)^N\times \prod_{k=1}^{n-1}\frac{\Gamma\left(N+\frac{1}{2}+\frac{k}{2n}\right)}{\Gamma\left(\frac{1}{2}+\frac{k}{2n}\right)} \end{equation} which, after asymptotic expansion, gives a subleading power-law term $L^{(n-1)/4}$ in $Z_n$, and therefore a $-\frac{1}{4}\ln L$ contribution to the Rényi entropy. This is the expected result from CFT [21,62]. Here also, $Z_2$ coincides with the norm of the Haldane-Shastry chain, with open boundary conditions[5]. Just as in the periodic case, we were also able to access $Z_4$ at half-filling. It is given by \begin{equation}\label{eq:openspecial} Z_4^{open}(L/2,L)=\frac{2^{L/2}}{(L+1)^L}\left(\frac{\Gamma \left(\frac{L}{2}+\frac{3}{4}\right)}{\Gamma\left(\frac{3}{4}\right)}\right)^2. \end{equation} An asymptotic expansion yields a power-law term $L^{1/2}$, and therefore a contribution $-\frac{1}{6}\ln L$ to the Rényi entropy. This is discussed in Sec. 4.2.3..
Open Neuman geometry
The open Neuman chain is slightly more complicated, as the number $N$ of particles is not conserved anymore. We managed to compute $Z_2$, which is given by \begin{equation}\label{eq:openneuman} Z_2=\frac{2^{L/2}}{(L+1)^{L/2+1}}\times \frac{\Gamma\left(\frac{L}{2}+\frac{5}{4}\right)}{\Gamma\left(\frac{5}{4}\right)}. \end{equation} An asymptotic expansion gives a power-law term $L^{-1/4}$, and therefore a $+\frac{1}{4}\ln L$ contribution to the second Rényi entropy. This is again consistent with CFT [62]. Also $S_\infty=-\ln p_{\textrm{max}}$ can be computed exactly, with \begin{eqnarray}\nonumber \\ p_{\textrm{max}}&=&\frac{2^{L/2(L/2-1)}}{(L+1)^{L/2}}\left(\cos \frac{\pi}{2L+2}\right)^{-L/2}\\ &\times&\prod_{k=1}^{L/4} \left(\sin \frac{2k\pi}{L+1}\right)^{L-2k}\left(\sin \frac{[2k+1]\pi}{L+1}\right)^{2k} \end{eqnarray} (we have assumed $L/2$ even). The asymptotic expansion once again yields a $\frac{1}{4}\ln L$ term, consistent with the Cardy-Peschel formula[12], with $h_{ND}=1/16$, where $h_{ND}$ is the dimension of the operator that changes the boundary condition from Dirichlet to Neuman[1].
Derivations
Let us now finally explain how the results of the previous subsection can be derived. The techniques are standard in random matrix theory [20], but there are here a few subtleties [22] due to the discrete nature of the problem. We focus on the most complicated case, the open geometry, where the results are to our knowledge not known. We have seen that the partition function of the gas is given by \begin{equation} Z_n=\frac{(2L+2)^{-Nn}}{N!}\sum_{\{x\}}\left[ \det_{1\leq j,k\leq N}\left(\omega^{x_j k}-\omega^{-x_j k}\right) \right]^{2n} \end{equation} The determinant inside the sum can be evaluated using (\ref{eq:svand}), but we leave it in this form for now. It is easy to check that $Z_1=1$. Indeed for $n=1$ the square can be written as a product of two determinants: explicitly writing the two sum over permutations, and then exchanging the order of the sums yields the result. The case $n=2$ can be obtained for any number of particles. To perform this calculation we rewrite the fourth power of the $\det$ using (\ref{eq:csvand}). We have
Now the sum over all the positions of the particles can be performed, and we obtain \begin{equation}\label{eq:pfaffian} Z_2=\frac{1}{2^NN!}\sum_P(-1)^P M_{p_1 p_2}M_{p_3p_4}\ldots M_{p_{2N-1}p_{2N}}, \end{equation} where the matrix elements are given by \begin{equation}\label{eq:truc} M_{pq}=\sum_{x=1}^{L}A(x,p,q). \end{equation} Eq. (\ref{eq:pfaffian}) is formally the Pfaffian of the antisymmetric matrix given in (\ref{eq:truc}). It may be calculated by first computing the $2\times 2$ determinant of Eq. (\ref{eq:2by2det}), and noticing that the matrix elements $M_{pq}$ are non zero only when $p=q\pm 1$. In the end we get, after some further algebra, \begin{equation} Z_2(N,L)=\prod_{k=1}^N \frac{4k-1}{2L+2} =\left(\frac{2}{L+1}\right)^N\frac{\Gamma\left(N+\frac{3}{4}\right)}{\Gamma\left(\frac{3}{4}\right)} \label{eq:open2result} \end{equation} when $N\leq L/2$, which is compatible with (\ref{eq:openmacdonald}). When $N>L/2$ we get \begin{equation} Z_2(N,L)=Z_2(L-N,L). \end{equation} In the periodic case Eqs. (\ref{eq:pfaffian}) and (\ref{eq:truc}) also hold, with [22] \begin{equation} A(x,p,q)=\frac{q-p}{L^2}z^{x(p+q-2N-1)}. \end{equation} The matrix elements are non zero only when $p+q=2N+1$, and using this we reproduce Eq. (\ref{eq:gaudin}) for $n=2$. When $n$ is greater than two we use a different method, and rely on the result of Macdonald explained in Sec. 9.1.2.. In an open system the partition function can be written, using (\ref{eq:svand}), as
Now if we expand the $n$-th power of the product in (\ref{eq:part}), we get a sum of terms of the form $\omega^{\sum_i m_i x_i}$ for some integers $m_i$. The crucial point is that each $m_i$ has its module bounded by $|m_i|\leq 2nN$. Since $\sum_{x=1}^{L} \omega^{mx}=0$ unless $m$ is zero or a multiple of $2L+2$, only the constant term of the expanded product has a non vanishing contribution in the sum, provided $n<(L+1)/N$. This constant term is exactly given by the Macdonald result for the $C_N$ root systems, as is explained in section 9.1.2.. Plugging (\ref{eq:di}) and (\ref{eq:macdonald}) in (\ref{eq:part}) gives our result Eq. (\ref{eq:openmacdonald}). The periodic case follows from a similar argument [22] with the constant term identity corresponding to the $A_{N-1}$ root system, also known as the Dyson conjecture. We recover Eq. (\ref{eq:gaudin}), provided $n<L/(N-1)$. The restriction on the allowed values of $n$ is an important difference with the continuous Dyson gas. At filling half or less, $n=2$ is in the allowed values of $n$, and we also reproduce Eq. (\ref{eq:open2result}) which we derived first using the determinant approach. The latter method has an interesting advantage, in that it can easily be extended to compute correlation functions at $n=2$. Indeed, let us for example look at the two point correlation function of two particles at position $x_1$ and $x_2$: \begin{equation} C_{2}(x_1,x_2)=\frac{\sum_{\{x\}}^\prime \psi(x_1,\ldots,x_N)^4}{\sum_{\{x\}}\psi(x_1,\ldots,x_N)^4} \end{equation} where in $\sum^\prime_{\{x\}}$ we sum over all the position $x_3,x_4,\ldots,x_N$ of the remaining particles. The calculation is easiest in the periodic case. Using the approach leading to (\ref{eq:someeq}), only few permutations give a nonzero contribution. For that to happen $(p_1,p_2,p_3,p_4)$ has to be a permutation of $(k,2N+1-k,k^\prime,2N+1-k^\prime)$ for some choice of $1\leq k,k^\prime\leq N$. We get \begin{equation}\label{eq:twopoints1} C_{2}(x_1,x_2)=\frac{L^2}{4}\sum_{1\leq k<k^\prime\leq N}\frac{\Upsilon_{kk^\prime}(x_1,x_2)}{(2N+1-2k)(2N+1-2k^\prime)}, \end{equation} where \begin{equation}\label{eq:twopoints2} \Upsilon_{kk^\prime}(x_1,x_2)=\sum_{Q\in S_4} (-1)^Q A(x_1,q_1,q_2)A(x_2,q_3,q_4). \end{equation} Here the sum runs over all $4!=24$ permutations of $(k,2N+1-k,k^\prime,2N+1-k^\prime)$, and can be computed explicitly. In the end we recover (\ref{eq:corrxx}), using $\left\langle \sigma_{x_1}^z\sigma_{x_2}^z \right\rangle_{\!(2)}=4C_{2}(x_1,x_2)-4N^2/L^2$, and further simplifications.
The determinant method also allows to access $Z_4$, using the confluent versions of Vandermonde (\ref{eq:cvand},\ref{eq:csvand}) and expanding the product of the two determinants similar to the calculation of $Z_1$. The computations are somewhat cumbersome, but we can nevertheless recover the results obtained through the Macdonald conjecture (when $\rho \leq 1/4$). Outside of this region we could not obtain a simple closed-form formula, except in the special case of half-filling. This is interesting, as the Macdonald conjecture does not give the result. The method works for both periodic and open (Dirichlet) systems. The results are given by Eq. (\ref{eq:periodicspecial}) and Eq. (\ref{eq:openspecial}), and correspond exactly to the transition point $n=n_c$ discussed in the text. In the open Neuman case this method is the only one available, as the grand-canonical nature of the partition function necessarily breaks the condition $n\leq \rho^{-1}$, and the constant term results cannot be used anymore. The result for $Z_2$ is given by (\ref{eq:openneuman}). Let us finally mention that $S_\infty$ can be obtained by noticing that the most likely configurations are attained for particles as far apart from each other as possible. For example at half-filling this corresponds to the configuration $(1,3,5,\ldots ,L-1)$ and $(2,4,6,\ldots ,L)$ for the positions of the particles. In spin language, these are the homogeneous states $\left|\uparrow\downarrow\ldots\uparrow\downarrow\right\rangle$ and $\left|\downarrow\uparrow\ldots\downarrow\uparrow\right\rangle$ discussed in the main text.
References
- Ian Affleck. J. Phys. A: Math. Gen., 31:2761, 1998.
- F. C. Alcaraz and M. A. Rajabpour. Universal behavior of the shannon mutual information of critical quantum chains. Phys. Rev. Lett., 111:017201, 2013.
- Y. Y Atas and E Bogomolny. Multifractality of eigenfunctions in spin chains. Phys. Rev. E, 86:021104, 2012.
- A. A Belavin, A. M Polyakov, and A. B Zamolodchikov. Infinite conformal symmetry in two-dimensional quantum field theory. Nucl. Phys. B, 241:333, 1984.
- Denis Bernard, Vincent Pasquier, and Didina Serban. Exact solution of long-range interacting spin chains with boundaries. Europhys. Lett., 300:301, 1995.
- Ulrich Bilstein and Birgit Wehefritz. The xx–model with boundaries. part i: Diagonalization of the finite chain. J. Phys. A: Math. Gen., 324:191, 1999.
- Michael Brockmann, Jacopo De Nardis, Bram Wouters, and Jean-Sébastien Caux. A generalized gaudin determinant for overlaps of néel and xxz bethe states. J. Phys. A: Math. Theor., 47:145003, 2014.
- Pasquale Calabrese and John Cardy. Entanglement entropy and quantum field theory. J.Stat.Mech. P06002, 2004.
- Pasquale Calabrese, John Cardy, and Erik Tonni. Entanglement negativity in quantum field theory. Phys. Rev. Lett., 109:130502, October 2012.
- J. L Cardy. The ising model in a random boundary field. J. Phys. A: Math. Gen., 24:L1315, 1991.
- John L. Cardy. Boundary conditions, fusion rules and the verlinde formula. Nucl. Phys. B, 324:581, 1989.
- John L. Cardy and Ingo Peschel. Finite-size dependence of the free energy in two-dimensional critical systems. Nuclear Physics B, 300(0):377–392, 1988.
- Anushya Chandran, Vedika Khemani, and S. L. Sondhi. How universal is the entanglement spectrum? 2014.
- Philippe {Di Franceso}, Pierre Mathieu, and David Sénéchal. Conformal Field Theory. Springer-verlag, 1997.
- H. Diehl. Field-theoretic approach to critical behaviour at surfaces, volume Phase transitions and critical phenomena. Academic Press, New-York, domb and lebowitz edition, 1986.
- Jérome Dubail and Jean-Marie Stéphan. Universal behavior of a bipartite fidelity at quantum criticality. J. Stat. Mech. L03002, 2011.
- F. J. Dyson. J. Math. Phys., 3:140, 1962.
- F. Evers and A. D. Mirlin. Fluctuations of the inverse participation ratio at the anderson transition. Phys. Rev. Lett., 84:3690, 2000.
- F. Evers and A. D. Mirlin. Anderson transitions. Rev. Mod. Phys. 1355, 80:1355, 2008.
- P. J Forrester. Log-gases and Random Matrices. Princeton University Press, 2010.
- Eduardo Fradkin and Joel E. Moore. Entanglement entropy of 2d conformal quantum critical points: hearing the shape of a quantum drum. Phys. Rev. Lett., 97:050404, 2006.
- M Gaudin. Gaz coulombien discret a une dimension. J. Phys, 34:511, 1973.
- T Giamarchi. Quantum Physics in One Dimension. Oxford University Press, New-York, 2004.
- I.J. Good. Short proof of a conjecture by dyson. J. Math. Phys., 11:1884, 1972.
- F. D. M Haldane. Exact jastrow-gutzwiller resonating-valence-bond ground state of the spin-(1/2 antiferromagnetic heisenberg chain with 1/r2 exchange. Phys. Rev. Lett., 60:635, 1988.
- C. Holzhey, F. Larsen, and F. Wilczek. Geometric and renormalized entropy in conformal field theory. Nucl.Phys. B, 424:443, 1994.
- Benjamin Hsu, Michael Mulligan, Eduardo Fradkin, and Eun-Ah Kim. Universal entanglement entropy in 2d conformal quantum critical points. Phys. Rev. B, 79:115421, 2009.
- Jason Iaconis, Stephen Inglis, Ann B. Kallin, and Roger G. Melko. Detecting classical phase transitions with renyi mutual information. Phys. Rev. B, 87:195134, 2013.
- F Igloi, L Turban, and B Berche. Critical behavior near multiple junctions and dirty surfaces in the two-dimensional ising model. J. Phys. A: Math. Gen., 24:L1031, 1991.
- F. M. Israilev. Simple models of quantum chaos: Spectrum and eigenfunctions. Physics Reports, 196:299, 1990.
- Hyejin Ju, Ann B. Kallin, Paul Fendley, Matthew B. Hastings, and Roger G. Melko. Entanglement scaling in two-dimensional gapless systems. Phys. Rev. B 85, 85:165121, 2012.
- Yuta Kumano. unpublished, 2013.
- Hon Wai Lau and Peter Grassberger. Information theoretic aspects of the two-dimensional ising model. Phys. Rev. E, 87:022128, 2012.
- A. Lazarescu. PhD thesis, 2013.
- David J. Luitz, Fabien Alet, and Nicolas Laflorencie. Shannon-rényi entropies and participation spectra across 3d $o(3)$ criticality. Phys. Rev. B, 89:165106, 2014.
- David J. Luitz, Fabien Alet, and Nicolas Laflorencie. Universal behavior beyond multifractality in quantum many-body systems. Phys. Rev. Lett., 112:057203, 2014.
- David J. Luitz, Nicolas Laflorencie, and Fabien Alet. Participation spectroscopy and entanglement hamiltonian of quantum spin models. 2014.
- I. G. Macdonald. SIAM J. Math. Anal., 13:998, 1982.
- Anne E. B. Nielsen, J. Ignacio Cirac, and German Sierra. Quantum spin hamiltonians for the su(2)_k wzw model. J. Stat. Mech. P11014, 2011.
- The volume law may be violated if the locality requirement is dropped. For example the Shannon entropy evaluated in the basis of the eigenstates of the Hamiltonian –which is usually highly nonlocal– is $0$ if the ground-state is nondegenerate.
- Note that this is not true for the full periodic chain, where there are no logarithms, and the leading universal piece depends on the stiffness [53,57].
- I thank John Cardy for pointing that out to me.
- In this precise case the method can be slightly improved using the exact result for $S_4(L,L)$ derived in the appendix 9., from which the leading term $a_4$ in Eq. (\ref{eq:infscaling}) can be extracted.
- Indeed, for a fixed aspect ratio $\ell /L=p/q$, we need to study even system sizes also multiple of $p$.
- Masaki Oshikawa. Boundary conformal field theory and entanglement entropy in two-dimensional quantum lifshitz critical point. preprint, 2010.
- B. Pozsgay. Overlaps between eigenstates of the xxz spin-1/2 chain and a class of simple product states. preprint, 2013.
- Lea F. Santos and Marcos Rigol. Onset of quantum chaos in one-dimensional bosonic and fermionic systems and its relation to thermalization. Phys. Rev. E, 81:036206, 2010.
- B. S. Shastry. Exact solution of an s=1/2 heisenberg antiferromagnetic chain with long-ranged interactions. Phys. Rev. Lett., 60:639, 1988.
- Y. G. Sinai. Funct. Anal. Appl., 24:61, 1968.
- J.-M. Stéphan and J. Dubail. Logarithmic corrections to the free energy from sharp corners with angle 2{$\pi$}. J. Stat. Mech P09002, 2013.
- Jean-Marie Stéphan. Emptiness formation probability, toeplitz determinants, and conformal field theory. J. Stat. Mech P05010, 2014.
- Jean-Marie Stéphan. Intrication dans des systèmes quantiques à basse dimension. PhD thesis, p87-89, in French, 2011.
- Jean-Marie Stéphan, Shunsuke Furukawa, Grégoire Misguich, and Vincent Pasquier. Shannon and entanglement entropies of one- and two-dimensional critical wave functions. Phys. Rev. B, 80:184421, 2009.
- Jean-Marie Stéphan, Stephen Inglis, Paul Fendley, and Roger G. Melko. Geometric mutual information at classical critical points. Phys. Rev. Lett., 112:127204, 2014.
- Jean-Marie Stéphan, Hyejin Ju, Paul Fendley, and Roger G. Melko. Entanglement in gapless resonating valence bond states. New. J. Phys., 15:015004, 2013.
- Jean-Marie Stéphan, Grégoire Misguich, and Vincent Pasquier. Rényi entropy of a line in two-dimensional ising models. Phys. Rev. B, 82:125455, 2010.
- Jean-Marie Stéphan, Grégoire Misguich, and Vincent Pasquier. Phase transition in the rényi-shannon entropy of luttinger liquids. Phys. Rev. B, 84:195128, 2011.
- A. M. Tsvelik. Phys. Rev. Lett., 110:147202, 2013.
- Jaegon Um, Hyunggyu Park, and Haye Hinrichsen. Entanglement versus mutual information in quantum spin chains. J. Stat. Mech. P10026, 2012.
- G. Vidal, J. I. Latorre, E. Rico, and A. Kitaev. Entanglement in quantum critical phenomena. Phys.Rev.Lett., 90:227902, 2003.
- Birgit Wehefritz-Kaufman. J. Phys. A: Math. Theor., 40:217, 2007.
- Michael P. Zaletel, Jens H. Bardarson, and Joel E. Moore. Logarithmic terms in entanglement entropies of 2d quantum critical points and shannon entropies of spin chains. Phys. Rev. Lett., 107:020402, 2011.
Footnotes
1. The volume law may be violated if the locality requirement is dropped. For example the Shannon entropy evaluated in the basis of the eigenstates of the Hamiltonian –which is usually highly nonlocal– is $0$ if the ground-state is nondegenerate.
2. Note that this is not true for the full periodic chain, where there are no logarithms, and the leading universal piece depends on the stiffness [53,57].
3. I thank John Cardy for pointing that out to me.
4. In this precise case the method can be slightly improved using the exact result for $S_4(L,L)$ derived in the appendix 9., from which the leading term $a_4$ in Eq. (\ref{eq:infscaling}) can be extracted.
5. Indeed, for a fixed aspect ratio $\ell/L=p/q$, we need to study even system sizes also multiple of $p$.