Dirac zero-modes in compact U(1) gauge theory

We study properties of the zero and near-zero eigenmodes of the overlap Dirac operator in compact U(1) gauge theory. In the confinement phase the exact zero-modes are localized as found by studying the values of the inverse participation ratio and other features. Non-zero-eigenmodes are less localized in the confinement phase. In the Coulomb phase no zero-modes are observed and the eigenmodes show no localization at all.Comment: Minor corrections, 15 pages, 5 figures, LaTeX styl

Wilderness Solitude in the 21st Century

Recent advances in mobile communication technology have led to a decrease in opportunities for individuals to experience alone-time within daily life. As a result, the solitude offered by wilderness landscapes has become all the more valuable. Past research on wilderness solitude has been divided into two distinct frameworks: the Social-Spatial Perspective and the Humanistic Perspective. This distinction has severely limited the development of a comprehensive research model that incorporates all the possible conditions relating to wilderness solitude. This study synthesized past research and theory to create a quantitative model of wilderness solitude which includes elements from both research perspectives, while incorporating novel conditions that relate to digital connectivity. Study participants were wilderness visitors to Montana’s Bob Marshall Wilderness Complex during the summer and fall of 2017. Exploratory factor analysis revealed four components of wilderness solitude. These components suggest that our interpretation of the “opportunities for solitude” clause within the Wilderness Act of 1964 ought to consider the themes of Societal Release, Introspection, Physical Separation and De-tethering from Digital Connectivity

Entanglement Spectra of Interacting Fermions in Quantum Monte Carlo Simulations

In a recent article T. Grover [Phys. Rev. Lett. 111, 130402 (2013)] introduced a simple method to compute Renyi entanglement entropies in the realm of the auxiliary field quantum Monte Carlo algorithm. Here, we further develop this approach and provide a stabilization scheme to compute higher order Renyi entropies and an extension to access the entanglement spectrum. The method is tested on systems of correlated topological insulators.Comment: 7+ pages, 5 figure

Interaction induced Dirac fermions from quadratic band touching in bilayer graphene

We revisit the effect of local interactions on the quadratic band touching (QBT) of Bernal stacked bilayer graphene models using renormalization group (RG) arguments and quantum Monte Carlo simulations of the Hubbard model. We present an RG argument which predicts, contrary to previous studies, that weak interactions do not flow to strong coupling even if the free dispersion has a QBT. Instead they generate a linear term in the dispersion, which causes the interactions to flow back to weak coupling. Consistent with this RG scenario, in unbiased quantum Monte Carlo simulations of the Hubbard model we find compelling evidence that antiferromagnetism turns on at a finite $U/t$, despite the $U=0$ hopping problem having a QBT. The onset of antiferromagnetism takes place at a continuous transition which is consistent with a dynamical critical exponent $z=1$ as expected for 2+1 d Gross-Neveu criticality. We conclude that generically in models of bilayer graphene, even if the free dispersion has a QBT, small local interactions generate a Dirac phase with no symmetry breaking and that there is a finite-coupling transition out of this phase to a symmetry-broken state

Effective models for strong electronic correlations at graphene edges

We describe a method for deriving effective low-energy theories of electronic interactions at graphene edges. Our method is applicable to general edges of honeycomb lattices (zigzag, chiral, and even disordered) as long as localized low-energy states (edge states) are present. The central characteristic of the effective theories is a dramatically reduced number of degrees of freedom. As a consequence, the solution of the effective theory by exact diagonalization is feasible for reasonably large ribbon sizes. The quality of the involved approximations is critically assessed by comparing the correlation functions obtained from the effective theory with numerically exact quantum Monte-Carlo calculations. We discuss effective theories of two levels: a relatively complicated fermionic edge state theory and a further reduced Heisenberg spin model. The latter theory paves the way to an efficient description of the magnetic features in long and structurally disordered graphene edges beyond the mean-field approximation.Comment: 13 pages, 9 figure

Z2 topological invariants in two dimensions from quantum Monte Carlo

We employ quantum Monte Carlo techniques to calculate the $Z_2$ topological invariant in a two-dimensional model of interacting electrons that exhibits a quantum spin Hall topological insulator phase. In particular, we consider the parity invariant for inversion-symmetric systems, which can be obtained from the bulk's imaginary-time Green's function after an appropriate continuation to zero frequency. This topological invariant is used here in order to study the trivial-band to topological-insulator transitions in an interacting system with spin-orbit coupling and an explicit bond dimerization. We discuss the accessibility and behavior of this topological invariant within quantum Monte Carlo simulations.Comment: 7 pages, 6 figure

Spontaneous particle-hole symmetry breaking of correlated fermions on the Lieb lattice

We study spinless fermions with nearest-neighbor repulsive interactions ($t$-$V$ model) on the two-dimensional three-band Lieb lattice. At half-filling, the free electronic band structure consists of a flat band at zero energy and a single cone with linear dispersion. The flat band is expected to be unstable upon inclusion of electronic correlations, and a natural channel is charge order. However, due to the three-orbital unit cell, commensurate charge order implies an imbalance of electron and hole densities and therefore doping away from half-filling. Our numerical results show that below a finite-temperature Ising transition a charge density wave with one electron and two holes per unit cell and its partner under particle-hole transformation are spontaneously generated. Our calculations are based on recent advances in auxiliary-field and continuous-time quantum Monte Carlo simulations that allow sign-free simulations of spinless fermions at half-filling. It is argued that particle-hole symmetry breaking provides a route to access levels of finite doping, without introducing a sign problem.Comment: 9 pages, 6 figures, added data for strong Coulomb repulsion and classical Ising-limi

Dimerized Solids and Resonating Plaquette Order in SU(N)-Dirac Fermions

We study the quantum phases of fermions with an explicit SU(N)-symmetric, Heisenberg-like nearest-neighbor flavor exchange interaction on the honeycomb lattice at half-filling. Employing projective (zero temperature) quantum Monte Carlo simulations for even values of N, we explore the evolution from a weak-coupling semimetal into the strong-coupling, insulating regime. Furthermore, we compare our numerical results to a saddle-point approximation in the large-N limit. From the large-N regime down to the SU(6) case, the insulating state is found to be a columnar valence bond crystal, with a direct transition to the semimetal at weak, finite coupling, in agreement with the mean-field result in the large-N limit. At SU(4) however, the insulator exhibits a subtly different valence bond crystal structure, stabilized by resonating valence bond plaquettes. In the SU(2) limit, our results support a direct transition between the semimetal and an antiferromagnetic insulator.Comment: 5 pages, 6 figure