BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Stefan Wessel

Stefan Wessel appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.str-elcond-mat.quant-gascond-mat.mes-hallcond-mat.stat-mechcond-mat.mtrl-scicond-mat.dis-nncond-mat.othercond-mat.supr-conhep-thquant-ph

Catalog footprint

What is connected

45works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

45 published item(s)

preprint2022arXiv

Cluster quantum Monte Carlo study of two-dimensional weakly-coupled frustrated trimer antiferromagnets

We report results from spin trimer-based cluster quantum Monte Carlo simulations for the thermodynamic properties of two-dimensional frustrated quantum antiferromagnets that are composed of weakly-coupled three-spin (trimer) clusters. In particular, we consider the spin-1/2 kagome lattice with a strong breathing distortion, and the triangle-square lattice model proposed previously for the cuprate La${}_4$Cu${}_3$MoO${}_{12}$. For both cases, we demonstrate that an appropriately chosen trimer-based computational basis allows us to significantly reduce the quantum Monte Carlo sign problem down to the low-temperature regime. Besides exploring the thermodynamic behavior for the triangle-square lattice model we also assess a mean field theory-based prediction for the onset of chiral order. For the breathing distorted kagome lattice model, we observe a robust two-peak structure in the specific heat, both in the quantum spin liquid and the lattice-nematic regimes.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Evidence for pressure induced unconventional quantum criticality in the coupled spin ladder antiferromagnet C$_9$H$_{18}$N$_2$CuBr$_4$

Quantum phase transitions in quantum matter occur at zero temperature between distinct ground states by tuning a nonthermal control parameter. Often, they can be accurately described within the Landau theory of phase transitions, similarly to conventional thermal phase transitions. However, this picture can break down under certain circumstances. Here, we present a comprehensive study of the effect of hydrostatic pressure on the magnetic structure and spin dynamics of the spin-1/2 ladder compound C$_9$H$_{18}$N$_2$CuBr$_4$. Single-crystal heat capacity and neutron diffraction measurements reveal that the N$\rm \acute{e}$el-ordered phase breaks down beyond a critical pressure of $P_{\rm c}$$\sim$1.0 GPa through a continuous quantum phase transition. Estimates of the critical exponents suggest that this transition may fall outside the traditional Landau paradigm. The inelastic neutron scattering spectra at 1.3 GPa are characterized by two well-separated gapped modes, including one continuum-like and another resolution-limited excitation in distinct scattering channels, which further indicates an exotic quantum-disordered phase above $P_{\rm c}$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Hierarchical single-ion anisotropies in spin-1 Heisenberg antiferromagnets on the honeycomb lattice

We examine the thermal properties of the spin-1 Heisenberg antiferromagnet on the honeycomb lattice in the presence of an easy-plane single-ion anisotropy as well as the effects of an additional weak in-plane easy-axis anisotropy. In particular, using large-scale quantum Monte Carlo simulations, we analyze the scaling of the correlation length near the thermal phase transition into the ordered phase. This allows us to quantify the temperature regime above the critical point in which -- in spite of the additional in-plane easy-axis anisotropy -- characteristic easy-plane physics, such as near a Berezinskii-Kosterlitz-Thouless transition, can still be accessed. Our theoretical analysis is motivated by recent neutron scattering studies of the spin-1 compound BaNi${}_2$V${}_2$O${}_8$ in particular, and it addresses basic quantum spin models for generic spin-1 systems with weak anisotropies, which we probe over the full range of experimentally relevant correlation length scales.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Surrogate models for quantum spin systems based on reduced order modeling

We present a methodology to investigate phase-diagrams of quantum models based on the principle of the reduced basis method (RBM). The RBM is built from a few ground-state snapshots, i.e., lowest eigenvectors of the full system Hamiltonian computed at well-chosen points in the parameter space of interest. We put forward a greedy-strategy to assemble such small-dimensional basis, i.e., to select where to spend the numerical effort needed for the snapshots. Once the RBM is assembled, physical observables required for mapping out the phase-diagram (e.g., structure factors) can be computed for any parameter value with a modest computational complexity, considerably lower than the one associated to the underlying Hilbert space dimension. We benchmark the method in two test cases, a chain of excited Rydberg atoms and a geometrically frustrated antiferromagnetic two-dimensional lattice model, and illustrate the accuracy of the approach. In particular, we find that the ground-manifold can be approximated to sufficient accuracy with a moderate number of basis functions, which increases very mildly when the number of microscopic constituents grows - in stark contrast to the exponential growth of the Hilbert space needed to describe each of the few snapshots. A combination of the presented RBM approach with other numerical techniques circumventing even the latter big cost, e.g., Tensor Network methods, is a tantalising outlook of this work.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Exact Critical Casimir Amplitude of Anisotropic Systems from Conformal Field Theory and Self-Similarity of Finite-Size Scaling Functions in $d\geq 2$ Dimensions

The exact critical Casimir amplitude is derived for anisotropic systems within the $d=2$ Ising universality class by combining conformal field theory (CFT) with anisotropic $φ^4$ theory. Explicit results are presented for the general anisotropic scalar $φ^4$ model and for the fully anisotropic triangular-lattice Ising model in finite rectangular and infinite strip geometries with periodic boundary conditions (PBC). These results demonstrate the validity of multiparameter universality for confined anisotropic systems and the nonuniversality of the critical Casimir amplitude. We find an unexpected complex form of self-similarity of the anisotropy effects near the instability where weak anisotropy breaks down. This can be traced back to the property of modular invariance of isotropic CFT for $d=2$. More generally, for $d>2$ we predict the existence of self-similar structures of the finite-size scaling functions of $O(n)$-symmetric systems with planar anisotropies and PBC both in the critical region for $ n \geq 1$ as well as in the Goldstone-dominated low-temperature region for $ n \geq 2$.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Quantum Monte Carlo simulations of highly frustrated magnets in a cluster basis: The two-dimensional Shastry-Sutherland model

Quantum Monte Carlo (QMC) simulations constitute nowadays one of the most powerful methods to study strongly correlated quantum systems, provided that no "sign problem" arises. However, many systems of interest, including highly frustrated magnets, suffer from an average sign that is close to zero in standard QMC simulations. Nevertheless, a possible sign problem depends on the simulation basis, and here we demonstrate how a suitable choice of cluster basis can be used to eliminate or at least reduce the sign problem in highly frustrated magnets that were so far inaccessible to efficient QMC simulations. We focus in particular on the application of a two-spin (dimer)-based QMC method to the thermodynamics of the spin-1/2 Shastry-Sutherland model for SrCu$_2$(BO$_3$)$_2$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Magnetism of Finite Graphene Samples: Mean-Field Theory compared with Exact Diagonalization and Quantum Monte Carlo Simulation

The magnetic properties of graphene on finite geometries are studied using a self-consistent mean-field theory of the Hubbard model. This approach is known to predict ferromagnetic edge states close to the zig-zag edges in single-layer graphene quantum dots and nanoribbons. In order to assess the accuracy of this method, we perform complementary exact diagonalization and quantum Monte Carlo simulations. We observe good quantitative agreement for all quantities investigated provided that the Coulomb interaction is not too strong.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Quantifying the fragility of unprotected quadratic band crossing points

We examine a basic lattice model of interacting fermions that exhibits quadratic band crossing points (QBCPs) in the non-interacting limit. In particular, we consider spinless fermions on the honeycomb lattice with nearest neighbor hopping $t$ and third-nearest neighbor hopping $t''$, which exhibits fine-tuned QBCPs at the corners of the Brillouin zone for ${t'' = t/2}$. In this situation, the density of states remains finite at the Fermi level of the half-filled band and repulsive nearest-neighbor interactions $V$ lead to a charge-density-wave (CDW) instability at infinitesimally small $V$ in the random-phase approximation or mean-field theory. We examine the fragility of the QBCPs against dispersion renormalizations in the ${t\mbox{-}t''\mbox{-}V}$ model using perturbation theory, and find that the $t''$-value needed for the QBCPs increases with $V$ due to the hopping renormalization. However, the instability toward CDW formation always requires a nonzero threshold interaction strength, i.e., one cannot fine-tune $t''$ to recover the QBCPs in the interacting system. These perturbative arguments are supported by quantum Monte Carlo simulations for which we carefully compare the corresponding threshold scales at and beyond the QBCP fine-tuning point. From this analysis, we thus gain a quantitative microscopic understanding of the fragility of the QBCPs in this basic interacting fermion system.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Reduction of the sign problem near $T=0$ in quantum Monte Carlo simulations

Building on a recent investigation of the Shastry-Sutherland model [S. Wessel et al., Phys. Rev. B 98, 174432 (2018)], we develop a general strategy to eliminate the Monte Carlo sign problem near the zero temperature limit in frustrated quantum spin models. If the Hamiltonian of interest and the sign-problem-free Hamiltonian---obtained by making all off-diagonal elements negative in a given basis---have the same ground state and this state is a member of the computational basis, then the average sign returns to one as the temperature goes to zero. We illustrate this technique by studying the triangular and kagome lattice Heisenberg antiferrromagnet in a magnetic field above saturation, as well as the Heisenberg antiferromagnet on a modified Husimi cactus in the dimer basis. We also provide detailed appendices on using linear programming techniques to automatically generate efficient directed loop updates in quantum Monte Carlo simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Spin vs. bond correlations along dangling edges of quantum critical magnets

Dangling edge spins of two-dimensional quantum critical antiferromagnets display strongly enhanced spin correlations with scaling dimensions that fall outside of the classical theory of surface critical phenomena. We provide large-scale quantum Monte Carlo results for both spin and bond correlation functions for the case of the columnar dimer model in particular. Unlike the spin correlations, we find the bond correlations to differ starkly between the spin-1/2 and spin-1 case. Furthermore, we compare the corresponding scaling dimensions to recent theoretical predictions. These predictions are, in part, supported by our numerical data, but cannot explain our findings completely. Our results thus put further constraints on completing the understanding of dangling edge correlations, as well as surface phenomena in strongly-correlated quantum systems in general.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Emergent symmetries and coexisting orders in Dirac fermion systems

The quantum phase diagram and critical behavior of two-dimensional Dirac fermions coupled to two compatible order-parameter fields with $O(N_1)\oplus O(N_2)$ symmetry is investigated. Recent numerical studies of such systems have reported evidence for non-Landau-Ginzburg-Wilson transitions and emergent $O(N_1+N_2)$ symmetry between the two ordered states, which has been interpreted within a scenario of deconfined quantum criticality in (2+1)-dimensional Dirac materials. Here, we provide two theoretical approaches to refine the phase diagrams of such systems. In the immediate vicinity of the multicritical point between the ordered phases and the semimetallic phase, we employ a non-perturbative field-theoretical analysis based on the functional renormalization group. For the particular case of $N_1=3$, $N_2=1$, we perform a large-scale quantum Monte Carlo analysis of the strong-coupling region, where both orders meet. Our findings support the robust emergence of enhanced symmetry at the multicritical point and suggest the transition between the two ordered phases to take place via a sequence of continuous transitions. In particular, we find that intermediate regimes of coexistence are present in the phase diagram for all values of $N_1$ and $N_2$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Competing pairing channels in the doped honeycomb lattice Hubbard model

Proposals for superconductivity emerging from correlated electrons in the doped Hubbard model on the honeycomb lattice range from chiral $d+id$ singlet to $p+ip$ triplet pairing, depending on the considered range of doping and interaction strength, as well as the approach used to analyze the pairing instabilities. Here, we consider these scenarios using large-scale dynamic cluster approximation (DCA) calculations to examine the evolution in the leading pairing symmetry from weak to intermediate coupling strength. These calculations focus on doping levels around the van Hove singularity (VHS) and are performed using DCA simulations with an interaction-expansion continuous-time quantum Monte Carlo cluster solver. We calculated explicitly the temperature dependence of different uniform superconducting pairing susceptibilities and found a consistent picture emerging upon gradually increasing the cluster size: while at weak coupling the $d+id$ singlet pairing dominates close to the VHS filling, an enhanced tendency towards $p$-wave triplet pairing upon further increasing the interaction strength is observed. The relevance of these systematic results for existing proposals and ongoing pursuits of odd-parity topological superconductivity are also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Thermal Ising transitions in the vicinity of two-dimensional quantum critical points

The scaling of the transition temperature into an ordered phase close to a quantum critical point as well as the order parameter fluctuations inside the quantum critical region provide valuable information about universal properties of the underlying quantum critical point. Here, we employ quantum Monte Carlo simulations to examine these relations in detail for two-dimensional quantum systems that exhibit a finite-temperature Ising-transition line in the vicinity of a quantum critical point that belongs to the universality class of either (i) the three-dimensional Ising model for the case of the quantum Ising model in a transverse magnetic field on the square lattice or (ii) the ${\it chiral}$ Ising transition for the case of a half-filled system of spinless fermions on the honeycomb lattice with nearest-neighbor repulsion. While the first case allows large-scale simulations to assess the scaling predictions to a high precision in terms of the known values for the critical exponents at the quantum critical point, for the later case we extract values of the critical exponents $ν$ and $η$, related to the order parameter fluctuations, which we discuss in relation to other recent estimates from ground state quantum Monte Carlo calculations as well as analytical approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Correlations and entanglement in quantum critical bilayer and necklace XY models

We analyze the critical properties and the entanglement scaling at the quantum critical points of the spin-half XY model on the two-dimensional square-lattice bilayer and necklace lattice, based on quantum Monte Carlo simulations on finite tori and for different subregion shapes. For both models, the finite-size scaling of the transverse staggered spin structure factor is found in accord with a quantum critical point described by the two-component, three-dimensional $ϕ^4$-theory. The second Rényi entanglement entropy in the absence of corners along the subsystem boundary exhibits area-law scaling in both models, with an area-law prefactor of $0.0674(7)$ [$0.0664(4)$] for the bilayer [necklace] model, respectively. Furthermore, the presence of $90^{\circ}$ corners leads to an additive logarithmic term in both models. We estimate a contribution of $-0.010(2)$ [$-0.009(2)$] due to each $90^{\circ}$ corner to the logarithmic correction for the bilayer [necklace] model, and compare our findings to recent numerical linked cluster calculations and series expansion results on related models.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Mott transition in the triangular lattice Hubbard model: a dynamical cluster approximation study

Based on dynamical cluster approximation (DCA) quantum Monte Carlo simulations, we study the interaction-driven Mott metal-insulator transition (MIT) in the half-filled Hubbard model on the anisotropic two-dimensional triangular lattice, where the degree of frustration is varied between the unfrustrated case and the fully frustrated, isotropic triangular lattice. Upon increasing the DCA cluster size, we analyze the evolution of the MIT phase boundary as a function of frustration in the phase diagram spanned by the interaction strength and temperature, and provide a quantitative description of the MIT phase boundary in the triangular lattice Hubbard model. Qualitative differences in the phase boundary between the unfrustrated and fully frustrated cases are exhibited. In particular, a change in the sign of the phase boundary slope is observed, which via an impurity cluster eigenstate analysis, may be related to a change in the nature of the insulating state. We discuss our findings within the scenario that the triangular lattice electron system might exhibit a quantum critical Mott MIT with a possible quantum spin liquid insulating state, such as considered for the organic charge transfer salts $κ$-(BEDT-TTF)$_{2}$Cu$_{2}$(CN)$_{3}$ and EtMe$_{3}$Sb[Pd(dmit)$_{2}$]$_{2}$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Nonlocal density interactions in auxiliary-field quantum Monte Carlo simulations: application to the square lattice bilayer and honeycomb lattice

We consider an efficient scheme to simulate fermionic Hubbard models with nonlocal density-density interactions in two dimensions, based on bond-centered auxiliary-field quantum Monte Carlo. The simulations are shown to be sign-problem free within a finite, restricted parameter range. Using this approach, we first study the Hubbard model on the half-filled square lattice bilayer, including an interlayer repulsion term in addition to the local repulsion, and present the ground state phase diagram within the accessible parameter region. Starting from the antiferromagnetically ordered state in the absence of interlayer repulsion, the interlayer interactions are found to destabilize the antiferromagnetic order, leading to a band insulator state. Moreover, for sufficiently strong interlayer tunneling, we also observe the emergence of a direct dimer product state of mixed D-Mott and S-Mott character along the equal coupling line. We discuss the stability range of this state within strong-coupling perturbation theory. Furthermore, we consider the Hubbard model on the honeycomb lattice with next-nearest-neighbor interactions. Such an interaction is found to enhance both charge density and spin-current correlations within the semimetallic region. However, inside the accessible parameter region, they do not stabilize long-ranged charge density wave order nor a quantum spin Hall state, and the only insulating state that we observe exhibits long-range antiferromagnetism.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Corner contribution to the entanglement entropy of strongly-interacting O(2) quantum critical systems in 2+1 dimensions

In a D=2+1 quantum critical system, the entanglement entropy across a boundary with a corner contains a subleading logarithmic scaling term with a universal coefficient. It has been conjectured that this coefficient is, to leading order, proportional to the number of field components N in the associated O(N) continuum $ϕ^4$ field theory. Using density matrix renormalization group calculations combined with the powerful numerical linked cluster expansion technique, we confirm this scenario for the O(2) Wilson-Fisher fixed point in a striking way, through direct calculation at the quantum critical points of two very different microscopic models. The value of this corner coefficient is, to within our numerical precision, twice the coefficient of the Ising fixed point. Our results add to the growing body of evidence that this universal term in the Rényi entanglement entropy reflects the number of low-energy degrees of freedom in a system, even for strongly interacting theories.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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ényi 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ényi 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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Ground state phase diagram of the half-filled bilayer Hubbard model

Employing a combination of functional renormalization group calculations and projective determinantal quantum Monte Carlo simulations, we examine the Hubbard model on the square lattice bilayer at half filling. From this combined analysis, we obtain a comprehensive account on the ground state phase diagram with respect to the extent of the system's metallic and (antiferromagnetically ordered) Mott-insulating as well as band-insulating regions. By means of an unbiased functional renormalization group approach, we exhibit the antiferromagnetic Mott-insulating state as the relevant instability of the free metallic state, induced by any weak finite onsite repulsion. Upon performing a careful analysis of the quantum Monte Carlo data, we resolve the difficulty of identifying this antiferromagnetic ground state for finite interlayer hopping in the weak-coupling regime, where nonmonotonous finite-size corrections are shown to relate to the two-sheeted Fermi surface structure of the metallic phase. On the other hand, quantum Monte Carlo simulations are well suited to identify the transition between the Mott-insulating phase and the band insulator in the intermediate-to-strong coupling regime. Here, we compare our numerical findings to indications for the transition region obtained from the functional renormalization group procedure.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Quantum Nature of Edge Magnetism in Graphene

It is argued that the subtle crossover from decoherence-dominated classical magnetism to fluctuation-dominated quantum magnetism is experimentally accessible in graphene nanoribbons. We show that the width of a nanoribbon determines whether the edge magnetism is on the classical side, on the quantum side, or in between. In the classical regime, decoherence is dominant and leads to static spin polarizations at the ribbon edges, which are well described by mean-field theories. The quantum Zeno effect is identified as the basic mechanism which is responsible for the spin polarization and thereby enables the application of graphene in spintronics. On the quantum side, however, the spin polarization is destroyed by dynamical processes. The great tunability of graphene magnetism thus offers a viable route for the study of the quantum-classical crossover.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Spin dynamics of the bilinear-biquadratic $S=1$ Heisenberg model on the triangular lattice: a quantum Monte Carlo study

We study thermodynamic properties as well as the dynamical spin and quadrupolar structure factors of the O(3)-symmetric spin-1 Heisenberg model with bilinear-biquadratic exchange interactions on the triangular lattice. Based on a sign-problem-free quantum Monte Carlo approach, we access both the ferromagnetic and the ferroquadrupolar ordered, spin nematic phase as well as the SU(3)-symmetric point which separates these phases. Signatures of Goldstone soft-modes in the dynamical spin and the quadrupolar structure factors are identified, and the properties of the low-energy excitations are compared to the thermodynamic behavior observed at finite temperatures as well as to Schwinger-boson flavor-wave theory.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Magnetic Correlations in Short and Narrow Graphene Armchair Nanoribbons

Electronic states at the ends of a narrow armchair nanoribbon give rise to a pair of non-locally entangled spins. We propose two experiments to probe these magnetic states, based on magnetometry and tunneling spectroscopy, in which correlation effects lead to a striking, nonlinear response to external magnetic fields. On the basis of low-energy theories that we derive here, it is remarkably simple to assess these nonlinear signatures for magnetic edge states. The effective theories are especially suitable in parameter regimes where other methods such as quantum Monte-Carlo simulations are exceedingly difficult due to exponentially small energy scales. The armchair ribbon setup discussed here provides a promisingly well-controlled (both experimentally and theoretically) environment for studying the principles behind edge magnetism in graphene-based nano-structures.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Quantum Monte Carlo studies of edge magnetism in chiral graphene nanoribbons

We investigate chiral graphene nanoribbons using projective quantum Monte Carlo simulations within the local Hubbard model description and study the effects of electron-electron interactions on the electronic and magnetic properties at the ribbon edges. Static and dynamical properties are analyzed for nanoribbons of varying width and edge chirality, and compared to a self-consistent Hartee-Fock mean-field approximation. Our results show that for chiral ribbons of sufficient width, the spin correlations exhibit exceedingly long correlation lengths, even between zigzag segments that are well separated by periodic armchair regions. Characteristic enhancements in the magnetic correlations for distinct ribbon widths and chiralities are associated with energy gaps in the tight-binding limit of such ribbons. We identify specific signatures in the local density of states and low- energy modes in the local spectral function which directly relate to enhanced electronic correlations along graphene nanoribbons and which can be accessed scanning tunneling spectroscopy.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Thermal phase transitions in a honeycomb lattice gas with three-body interactions

We study the thermal phase transitions in a classical (hard-core) lattice gas model with nearest-neighbor three-body interactions on the honeycomb lattice, based on parallel tempering Monte Carlo simulations. This system realizes incompressible low-temperature phases at fractional fillings of 9/16, 5/8 and 3/4 that were identified in a previous study of a related quantum model. In particular, both the 9/16 and the 5/8 phase exhibit an extensive ground state degeneracy reflecting the frustrated nature of the three-body interactions on the honeycomb lattice. The thermal melting of the 9/16 phase is found to be a first-order, discontinuous phase transition. On the other hand, from the thermodynamic behavior we obtain indications for a four-states Potts-model thermal transition out of the 5/8 phase. Employing an exact mapping to a hard-core dimer model on an embedded honeycomb super-lattice, we find that this thermal Potts-model transition relates to the selection of one out of four extensive sectors within the low-energy manifold of the 5/8 phase.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Adiabatic loading of one-dimensional SU(N) alkaline earth fermions in optical lattices

Ultracold fermionic alkaline earth atoms confined in optical lattices realize Hubbard models with internal SU(N) symmetries, where N can be as large as ten. Such systems are expected to harbor exotic magnetic physics at temperatures below the superexchange energy scale. Employing quantum Monte Carlo simulations to access the low-temperature regime, we show that after adiabatically loading a weakly interacting gas into the strongly interacting regime of an optical lattice, the final temperature decreases with increasing N. Furthermore, we estimate the temperature scale required to probe correlations associated with low-temperature SU(N) magnetism. Our findings are encouraging for the exploration of exotic large-N magnetic states in ongoing experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Anomalous Behavior of Spin Systems with Dipolar Interactions

We study the properties of spin systems realized by cold polar molecules interacting via dipole-dipole interactions in two dimensions. Using a spin wave theory, that allows for the full treatment of the characteristic long-distance tail of the dipolar interaction, we find several anomalous features in the ground state correlations and the spin wave excitation spectrum, which are absent in their counterparts with short range interaction. The most striking consequence is the existence of true long-range order at finite temperature for a two-dimensional phase with a broken U(1) symmetry.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Antiferromagnetism in the Hubbard Model on the Bernal-stacked Honeycomb Bilayer

Using a combination of quantum Monte Carlo simulations, functional renormalization group calculations and mean-field theory, we study the Hubbard model on the Bernal-stacked honeycomb bilayer at half-filling as a model system for bilayer graphene. The free bands consisting of two Fermi points with quadratic dispersions lead to a finite density of states at the Fermi level, which triggers an antiferromagnetic instability that spontaneously breaks sublattice and spin rotational symmetry once local Coulomb repulsions are introduced. Our results reveal an inhomogeneous participation of the spin moments in the ordered ground state, with enhanced moments at the three-fold coordinated sites. Furthermore, we find the antiferromagnetic ground state to be robust with respect to enhanced interlayer couplings and extended Coulomb interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Critical scales in anisotropic spin systems from functional renormalization

We apply a recently developed functional renormalization group (fRG) scheme for quantum spin systems to the spin-1/2 antiferromagnetic XXZ model on a two-dimensional square lattice. Based on an auxiliary fermion representation we derive flow equations which allow a resummation of the perturbation series in the spin-spin interactions. Spin susceptibilities are calculated for different values of the anisotropy parameter. The phase transition between planar and axial ordering at the isotropic point is reproduced correctly. The results for the critical scales from the fRG as quantitative measures for the ordering temperatures are in good agreement with the exact solution in the Ising limit. On the easy-plane side, the deviations from critical temperatures obtained with quantum Monte Carlo are larger but still acceptable. However, at the isotropic point the Mermin-Wagner theorem is violated such that a precise description of the behavior in the vicinity of the phase transition is not possible. We discuss possible reasons for these discrepanies.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Geometry fluctuations in a two-dimensional quantum antiferromagnet

The paper considers the effects of random fluctuations of the local spin connectivities (fluctuations of the geometry) on ground state properties of a two-dimensional quantum antiferromagnet. We analyse the behavior of spins described by the Heisenberg model as a function of what we call phason flip disorder, following a terminology used for aperiodic systems. The calculations were carried out both within linear spin wave theory and using quantum Monte Carlo simulations. An "order by disorder" phenomenon is observed in this model, wherein antiferromagnetism is found to be enhanced by phason disorder. The value of the staggered order parameter increases with the number of defects, accompanied by an increase in the ground state energy of the system.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Half-Vortex Unbinding and Ising Transition in Constrained Superfluids

We analyze the thermodynamics of the atomic and (nematic) pair superfluids appearing in the attractive two-dimensional Bose-Hubbard model with a three-body hard-core constraint that has been derived as an effective model for cold atoms subject to strong three-body losses in optical lattices. We show that the thermal disintegration of the pair superfluidity is governed by the proliferation of fractional half-vortices leading to a Berezinskii-Kosterlitz-Thousless transition with unusual jump in the helicity modulus. In addition to the (conventional) Berezinskii-Kosterlitz-Thousless transition out of the atomic superfluid, we furthermore identify a direct thermal phase transition separating the pair and the atomic superfluid phases, and show that this transition is continuous with critical scaling exponents consistent with those of the two-dimensional Ising universality class. We exhibit a direct connection between the partial loss of quasi long-range order at the Ising transition between the two superfluids and the parity selection in the atomic winding number fluctuations that distinguish the atomic from the pair superfluid.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Interaction-range effects for fermions in one dimension

Experiments on quasi-one-dimensional systems such as quantum wires and metallic chains on surfaces suggest the existence of electron-electron interactions of substantial range and hence physics beyond the Hubbard model. We therefore investigate one-dimensional, quarter-filled chains with a Coulomb potential with variable screening length by quantum Monte Carlo methods and exact diagonalization. The Luttinger liquid interaction parameter K_rho decreases with increasing interaction strength and range. Experimentally observed values close to 1/4 require strong interactions and/or large screening lengths. As predicted by bosonization, we find a metal-insulator transition at K_rho=1/4. Upon increasing the screening length, the charge and spin correlation functions reveal the crossover from dominant 2k_F spin correlations to dominant 4k_F charge correlations, and a strong enhancement of the charge velocity. In the metallic phase, the signatures of spin-charge separation in the single-particle spectrum, spinon and holon bands, remain robust even for rather long-ranged interactions. The charge-density-wave state exhibits backfolded shadow bands.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dynamical Signatures of Edge-State Magnetism on Graphene Nanoribbons

We investigate the edge-state magnetism of graphene nanoribbons using projective quantum Monte Carlo simulations and a self-consistent mean-field approximation of the Hubbard model. The static magnetic correlations are found to be short ranged. Nevertheless, the correlation length increases with the width of the ribbon such that already for ribbons of moderate widths we observe a strong trend towards mean-field-type ferromagnetic correlations at a zigzag edge. These correlations are accompanied by a dominant low-energy peak in the local spectral function and we propose that this can be used to detect edge-state magnetism by scanning tunneling microscopy. The dynamic spin structure factor at the edge of a ribbon exhibits an approximately linearly dispersing collective magnonlike mode at low energies that decays into Stoner modes beyond the energy scale where it merges into the particle-hole continuum.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Generic First-Order vs. Continuous Quantum Nucleation of Supersolidity

We analyze the nucleation of supersolid order out of the superfluid ground state of bosons on the triangular lattice. While the stability of supersolidity against phase separation in this system is by now well established for nearest-neighbor and long-range dipolar interactions, relevant for two-dimensional arrays of ultra-cold polar molecules, here we address directly the nature of the superfluid-to-supersolid transition. Based on symmetry arguments and quantum Monte Carlo simulations, we conclude that this quantum phase transition is driven first-order beyond the line of particle-hole symmetry. Along this line, the transition is continuous and its scaling behavior consistent with the three-dimensional (3D) XY universality class. We relate this finding to a 3D Z6 clock model description of the enlarged symmetry of the solid order parameter field. In the generic case however, the symmetry reduces to that of a 3D Z3 clock model, which reflects the first-order nature of the generic superfluid-to-supersolid quantum phase transition on the triangular lattice.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Pair Superfluidity of Three-Body Constrained Bosons in Two Dimensions

We examine the equilibrium properties of lattice bosons with attractive on-site interactions in the presence of a three-body hard-core constraint that stabilizes the system against collapse and gives rise to a dimer superfluid phase formed by virtual hopping processes of boson pairs. Employing quantum Monte Carlo simulations, the ground state phase diagram of this system on the square lattice is analyzed. In particular, we study the quantum phase transition between the atomic and dimer superfluid regime and analyze the nature of the superfluid-insulator transitions. Evidence is provided for the existence of a tricritical point along the saturation transition line, where the transition changes from being first-order to a continuous transition of the dilute bose gas of holes. The Berzinskii-Kosterlitz-Thouless transition from the dimer superfluid to the normal fluid is found to be consistent with an anomalous stiffness jump, as expected from the unbinding of half-vortices.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Supersolid Polar Molecules beyond Pairwise Interactions

We explore the phase diagram of ultracold bosonic polar molecules confined to a planar optical lattice of triangular geometry. External static electric and microwave fields can be employed to tune the effective interactions between the polar molecules into a regime of extended two- and three-body repulsions of comparable strength, leading to a rich quantum phase diagram. In addition to various solid phases, an extended supersolid phase is found to persist deep into the three-body dominated regime. While three-body interactions break particle-hole symmetry explicitly, a characteristic supersolid-supersolid quantum phase transition is observed, which indicates the restoration of particle-hole symmetry at half-filling. We revisit the spatial structure of the supersolid at this filling, regarding the existence of a further supersolid phase with three inequivalent sublattices, and provide evidence that this state is excluded also at finite temperatures.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Critical Entropy of Quantum Heisenberg Magnets on Simple-Cubic Lattices

We analyze the temperature dependence of the entropy of the spin-1/2 Heisenberg model on the three-dimensional simple-cubic lattice, for both the case of antiferromagnetic and ferromagnetic nearest neighbor exchange interactions. Using optimized extended ensemble quantum Monte Carlo simulations, we extract the entropy at the critical temperature for magnetic order from a finite-size scaling analysis. For the antiferromagnetic case, the critical entropy density equals 0.341(5)$k_B$, whereas for the ferromagnet, a larger value of 0.401(5) $k_B$ is obtained. We compare our simulation results to estimates put forward recently in studies assessing means of realizing the antiferromagnetic Néel state in ultra-cold fermion gases in optical lattices.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Finite temperature ordering of dilute graphene antiferromagnets

We employ large-scale quantum Monte Carlo simulations to study the magnetic ordering transition among dilute magnetic moments randomly localized on the graphene honeycomb lattice, induced by long-ranged RKKY interactions at low charge carrier concentration. In this regime the effective exchange interactions are ferromagnetic within each sublattice, and antiferromagnetic between opposite sublattices, with an overall cubic decay of the interaction strength with the separation between the moments. We verify explicitly, that this commensurability leads to antiferromagnetic order among the magnetic moments below a finite transition temperature in this two-dimensional system. Furthermore, the ordering temperature shows a crossover in its power-law scaling with the moments' dilution from a low- to a high-concentration regime.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A phason disordered two dimensional quantum antiferromagnet

We examine a novel type of disorder in quantum antiferromagnets. Our model consists of localized spins with antiferromagnetic exchanges on a bipartite quasiperiodic structure, which is geometrically disordered in such a way that no frustration is introduced. In the limit of zero disorder, the structure is the perfect Penrose rhombus tiling. This tiling is progressively disordered by augmenting the number of random "phason flips" or local tile-reshuffling operations. The ground state remains Néel ordered, and we have studied its properties as a function of increasing disorder using linear spin wave theory and quantum Monte Carlo. We find that the ground state energy decreases, indicating enhanced quantum fluctuations with increasing disorder. The magnon spectrum is progressively smoothed, and the effective spin wave velocity of low energy magnons increases with disorder. For large disorder, the ground state energy as well as the average staggered magnetization tend towards limiting values characteristic of this type of randomized tilings.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Polar Molecules with Three-Body Interactions on the Honeycomb Lattice

We study the phase diagram of ultra-cold bosonic polar molecules loaded on a two-dimensional optical lattice of hexagonal symmetry controlled by external electric and microwave fields. Following a recent proposal in Nature Physics \textbf{3}, 726 (2007), such a system is described by an extended Bose-Hubbard model of hard-core bosons, that includes both extended two- and three-body repulsions. Using quantum Monte-Carlo simulations, exact finite cluster calculations and the tensor network renormalization group, we explore the rich phase diagram of this system, resulting from the strongly competing nature of the three-body repulsions on the honeycomb lattice. Already in the classical limit, they induce complex solid states with large unit cells and macroscopic ground state degeneracies at different fractional lattice fillings. For the quantum regime, we obtain effective descriptions of the various phases in terms of emerging valence bond crystal states and quantum dimer models. Furthermore, we access the experimentally relevant parameter regime, and determine the stability of the crystalline phases towards strong two-body interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Phases and magnetization process of an anisotropic Shastry-Sutherland model

We examine ground state properties of the spin-1/2 easy-axis Heisenberg model on the Shastry-Sutherland lattice with ferromagnetic transverse spin exchange using quantum Monte Carlo and degenerate perturbation theory. For vanishing transverse exchange, the model reduces to an antiferromagnetic Ising model that besides Neel order harbors regions of extensive ground state degeneracy. In the quantum regime, we find a dimerized phase of triplet states, separated from the Neel ordered phase by a superfluid. The quantum phase transitions between these phases are characterized. The magnetization process shows a magnetization plateau at 1/3 of the saturation value, that persists down to the Ising limit, and a further plateau at 1/2 only in the quantum regime. For both plateaus, we determine the crystalline patterns of the localized triplet excitations. No further plateaus nor supersolid phases are found in this model.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Penrose Quantum Antiferromagnet

The Penrose tiling is a perfectly ordered two dimensional structure with fivefold symmetry and scale invariance under site decimation. Quantum spin models on such a system can be expected to differ significantly from more conventional structures as a result of its special symmetries. In one dimension, for example, aperiodicity can result in distinctive quantum entanglement properties. In this work, we study ground state properties of the spin-1/2 Heisenberg antiferromagnet on the Penrose tiling, a model that could also be pertinent for certain three dimensional antiferromagnetic quasicrystals. We show, using spin wave theory and quantum Monte Carlo simulation, that the local staggered magnetizations strongly depend on the local coordination number z and are minimized on some sites of five-fold symmetry. We present a simple explanation for this behavior in terms of Heisenberg stars. Finally we show how best to represent this complex inhomogeneous ground state, using the "perpendicular space" representation of the tiling.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Quantum Antiferromagnetism in Quasicrystals

The antiferromagnetic Heisenberg model is studied on a two-dimensional bipartite quasiperiodic lattice. The distribution of local staggered magnetic moments is determined on finite square approximants with up to 1393 sites, using the Stochastic Series Expansion Quantum Monte Carlo method. A non-trivial inhomogeneous ground state is found. For a given local coordination number, the values of the magnetic moments are spread out, reflecting the fact that no two sites in a quasicrystal are identical. A hierarchical structure in the values of the moments is observed which arises from the self-similarity of the quasiperiodic lattice. Furthermore, the computed spin structure factor shows antiferromagnetic modulations that can be measured in neutron scattering and nuclear magnetic resonance experiments. This generic model is a first step towards understanding magnetic quasicrystals such as the recently discovered Zn-Mg-Ho icosahedral structure.

0 citations0 reviewsNo code linkedOpen access