A R\'enyi entropy perspective on topological order in classical toric code models

Concepts of information theory are increasingly used to characterize collective phenomena in condensed matter systems, such as the use of entanglement entropies to identify emergent topological order in interacting quantum many-body systems. Here we employ classical variants of these concepts, in particular R\'enyi entropies and their associated mutual information, to identify topological order in classical systems. Like for their quantum counterparts, the presence of topological order can be identified in such classical systems via a universal, subleading contribution to the prevalent volume and boundary laws of the classical R\'enyi entropies. We demonstrate that an additional subleading $O(1)$ contribution generically arises for all R\'enyi entropies $S^{(n)}$ with $n \geq 2$ when driving the system towards a phase transition, e.g. into a conventionally ordered phase. This additional subleading term, which we dub connectivity contribution, tracks back to partial subsystem ordering and is proportional to the number of connected parts in a given bipartition. Notably, the Levin-Wen summation scheme -- typically used to extract the topological contribution to the R\'enyi entropies -- does not fully eliminate this additional connectivity contribution in this classical context. This indicates that the distillation of topological order from R\'enyi entropies requires an additional level of scrutiny to distinguish topological from non-topological $O(1)$ contributions. This is also the case for quantum systems, for which we discuss which entropies are sensitive to these connectivity contributions. We showcase these findings by extensive numerical simulations of a classical variant of the toric code model, for which we study the stability of topological order in the presence of a magnetic field and at finite temperatures from a R\'enyi entropy perspective.Comment: 17 pages, 19 figure

Entanglement entropy scaling in the bilayer Heisenberg spin system

We examine the entanglement properties of the spin-half Heisenberg model on the two-dimensional square-lattice bilayer based on quantum Monte Carlo calculations of the second R\'enyi entanglement entropy. In particular, we extract the dominant area-law contribution to the bipartite entanglement entropy that shows a non-monotonous behavior upon increasing the inter-layer exchange interaction: a local maximum in the area-law coefficient is located at the quantum critical point separating the antiferromagnetically ordered region from the disordered dimer-singlet regime. Furthermore, we consider subleading logarithmic corrections to the R\'enyi entanglement entropy scaling. Employing different subregion shapes, we isolate the logarithmic corner term from the logarithmic contribution due to Goldstone modes that is found to be enhanced in the limit of decoupled layers. At the quantum critical point, we estimate a contribution of $0.016(1)$ due to each $90^{\circ}$ corner. This corner term at the SU(2) quantum critical point deviates from the Gaussian theory value, while it compares well with recent numerical linked cluster calculations on the bilayer model.Comment: 7 pages, 7 figure

Mott transition of fermionic atoms in a three-dimensional optical trap

We study theoretically the Mott metal-insulator transition for a system of fermionic atoms confined in a three-dimensional optical lattice and a harmonic trap. We describe an inhomogeneous system of several thousand sites using an adaptation of dynamical mean field theory solved efficiently with the numerical renormalization group method. Above a critical value of the on-site interaction, a Mott-insulating phase appears in the system. We investigate signatures of the Mott phase in the density profile and in time-of-flight experiments.Comment: 4 pages and 5 figure

Kondo proximity effect: How does a metal penetrate into a Mott insulator?

We consider a heterostructure of a metal and a paramagnetic Mott insulator using an adaptation of dynamical mean field theory to describe inhomogeneous systems. The metal can penetrate into the insulator via the Kondo effect. We investigate the scaling properties of the metal-insulator interface close to the critical point of the Mott insulator. At criticality, the quasiparticle weight decays as 1/x^2 with distance x from the metal within our mean field theory. Our numerical results (using the numerical renormalization group as an impurity solver) show that the prefactor of this power law is extremely small.Comment: 4 pages, 3 figure

Hadron-Hadron Interactions from Nf=2+1+1N_f=2+1+1 Lattice QCD: I=3/2I=3/2 πK\pi K Scattering Length

In this paper we report on results for the s-wave scattering length of the $\pi$-$K$ system in the $I=3/2$ channel from $N_f=2+1+1$ Lattice QCD. The calculation is based on gauge configurations generated by the European Twisted Mass Collaboration with pion masses ranging from about $230$ to $450\,\text{MeV}$ at three values of the lattice spacing. Our main result reads $M_{\pi}\,a_0^{3/2,\text{phys}} = -0.059(2)$. Using chiral perturbation theory we are also able to estimate $M_{\pi}\,a_0^{1/2,\text{phys}} = 0.163(3)$. The error includes statistical and systematic uncertainties, and for the latter in particular errors from the chiral and continuum extrapolations.Comment: 30 pages, 6 figures, 15 table