We calculate the entanglement and the universal boundary entropy (BE) in critical quantum spin chains, such as the transverse field Ising chain and the XXZ chain, with arbitrary direction of the boundary magnetic field (ADBMF). We determine the boundary universality class that an ADBMF induces. In particular, we show that the induced boundary conformal field theory depends on the point on the Bloch sphere where the boundary magnetic field directs. We show that the classification of the directions boils down to the simple fact that the boundary field breaks the bulk symmetry or does not. We present a procedure to estimate the universal BE, based on the finite-size corrections of the entanglement entropy, which apply to ADBMF. To calculate the universal BE in the XXZ chain, we use the density matrix renormalization group. The transverse field XY chain with ADBMF after Jordan-Wigner transformation is not a quadratic free fermion Hamiltonian. We map this model to a quadratic free fermion chain by introducing two extra ancillary spins coupled to the main chain at the boundaries, which makes the problem integrable. The eigenstates of the transverse field XY chain can be obtained by proper projection in the enlarged chain. Using this mapping, we are able to calculate the entanglement entropy of the transverse field XY chain using the usual correlation matrix technique up to relatively large sizes.
We investigate the entanglement entropy of a block of
L sites in quasifree translation-invariant spin chains concentrating on the effect of reflection-symmetry breaking. The Majorana two-point functions corresponding to the Jordan-Wigner transformed fermionic modes are determined in the most general case; from these, it follows that reflection symmetry in the ground state can only be broken if the model is quantum critical. The large
L asymptotics of the entropy are calculated analytically for general gauge-invariant models, which have, until now, been done only for the reflection-symmetric sector. Analytical results are also derived for certain nongauge-invariant models (e.g., for the Ising model with Dzyaloshinskii-Moriya interaction). We also study numerically finite chains of length
N with a nonreflection-symmetric Hamiltonian and report that the reflection symmetry of the entropy of the first
L spins is violated but the reflection-symmetric Calabrese-Cardy formula is recovered asymptotically. Furthermore, for noncritical reflection-symmetry-breaking Hamiltonians, we find an anomaly in the behavior of the saturation entropy as we approach the critical line. The paper also provides a concise but extensive review of the block-entropy asymptotics in translation-invariant quasifree spin chains with an analysis of the nearest-neighbor case and the enumeration of the yet unsolved parts of the quasifree landscape.
Block entanglement entropy in the ground state of a quantum spin chain is investigated. The spins have Kitaev-type nearest-neighbor interaction, of strength
Jx or
Jy, through either
x or
y components of the spins on alternating bonds, along with a transverse magnetic field
h. An exact solution is obtained through Jordan-Wigner fermionization, and it exhibits a macroscopically degenerate ground state for
h=0 and a nondegenerate ground state for $h\ensuremath{\ne}0$ for all interaction strengths. For a chain of
N spins, we study the block entropy of a partition of
L contiguous spins. The block entanglement entropy needs the eigenvalues of the
2L-dimensional reduced density matrix. We employ an efficient method that reduces this problem to evaluating eigenvalues of an
L-dimensional matrix, which enables us to calculate easily the block entanglement for large-
N chains numerically. The entanglement entropy grows as
logL at the degeneracy point
h=0 and only for
Jx=Jy. For $h\ensuremath{\ne}0$, the entropy becomes independent of the size, thus obeying the area law. For ${J}_{x}\ensuremath{\ne}{J}_{y}$, the block entropy shows a nonmonotonic behavior for
L<N/2.
We examine the entanglement entropy of the even half of a translationally invariant finite chain or lattice in its ground state. This entropy measures the entanglement between the even and odd halves (each forming a ``comb'' of
n/2 sites) and can be expected to be extensive for short-range couplings away from criticality. We first consider bosonic systems with quadratic couplings, where analytic expressions for arbitrary dimensions can be provided. The bosonic treatment is then applied to finite spin chains and arrays by means of the random-phase approximation. Results for first-neighbor anisotropic
XY couplings indicate that, while at strong magnetic fields this entropy is strictly extensive, at weak fields important deviations arise, stemming from parity-breaking effects and the presence of a factorizing field (in the vicinity of which it becomes size-independent and identical to the entropy of a contiguous half). Exact numerical results for small spin
s chains are shown to be in agreement with the bosonic random-phase approximation prediction.
We study general entanglement properties of the excited states of the one-dimensional translational invariant free fermions and coupled harmonic oscillators. In particular, using the integrals of motion, we prove that these Hamiltonians independent of the gap (mass) have infinite excited states that can be described by conformal field theories with integer or half-integer central charges. In the case of free fermions, we also show that because of the huge degeneracy in the spectrum, even a gapless Hamiltonian can have excited states with an area law. Finally, we study the universal average entanglement entropy over eigenstates of the Hamiltonian of the free fermions introduced recently in Vidmar et al. [Phys. Rev. Lett. 119, 020601 (2017)]. In particular, using a duality relation, we map the average bipartite entanglement entropy over all the exponential numbers of eigenstates of a generic free fermion Hamiltonian to the problem of the calculation of the average multibipartite entanglement entropy over a single eigenstate of the
XX chain. This relation can be useful for experimental measurement of the universal average entanglement entropy. Part of our conclusions can be extended to the quantum spin chains that are associated with the free fermions via Jordan-Wigner transformation.
The half-chain entanglement entropy (EE) of the one-dimensional
XX (isotropic
XY) chain and quantum Ising chain (QIC) at the critical point are calculated numerically for large lattices. The periodic, free-free, fixed-fixed (
++) and ($+\ensuremath{-}$), and fixed-free boundary conditions (BCs) are studied. The models are transformed into a fermion system with a Jordan-Wigner transformation. The correlation functions of these fermions are obtained numerically. With the correlation functions, the entanglement entropy is obtained. The numerical calculation is performed with quadruple precision. The numerical results agree with the exact relationship with an accuracy of ${10}^{\ensuremath{-}28}$ between the EEs of the
XX chain and QICs: (1) the EE of the
XX chain with periodic BCs is the sum of two QICs with periodic BCs; (2) the EE of the
XX chain with free-free BCs is the sum of two QICs with free-free BCs; (3) the EE of the
XX chain with fixed-free BCs is the sum of the EE of a QIC with fixed-free BCs and that with free-free BCs; (4) the EE of the
XX chain with fixed-fixed (
++) BCs is the sum of the EE of a QIC with fixed-fixed (
++) BCs and that with free-free BCs; and (5) the EE of the
XX chain with fixed-fixed ($+\ensuremath{-}$) BCs is the sum of the EE of a QIC with fixed-fixed (
++) BCs and that of the first excited state of a QIC with free-free BCs. We fit the half-chain EE with finite-size scaling expansion of terms ${L}^{\ensuremath{-}n}$ and ${L}^{\ensuremath{-}n}lnL$, where
L is the length of the spin chain. We find that the fitted parameters for the leading term
lnL and the constant term are consistent with conformal field theory and that the exact solutions are highly accurate. The fitting accuracy is so high that the fitted central charge of the Ising model is $0.5\ifmmode\pm\else\textpm\fi{}1\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}26}$ for the periodic boundary condition. From the higher-order corrections in the fitted expansion, we find that there exist ${L}^{\ensuremath{-}4}lnL$ and higher-order corrections for the fixed-fixed ($+\ensuremath{-}$) boundary conditions, whereas there exist ${L}^{\ensuremath{-}3}lnL$ and higher-order corrections for the fixed-free boundary conditions. These types of corrections mean corrections of the central charge.