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

Fakher F. Assaad

Fakher F. Assaad 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.mes-hallcond-mat.supr-conhep-latquant-phcond-mat.stat-mechcond-mat.mtrl-scicond-mat.quant-gasphysics.comp-ph

Catalog footprint

What is connected

49works
9topics
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

49 published item(s)

preprint2026arXiv

Revealing altermagnetic Fermi surfaces with two Kondo impurities

Motivated by recent advances in the study of altermagnetism, or unconventional magnetism, and in the realization and manipulation of two-impurity Kondo physics in real materials, we propose a phase-sensitive method to explore unconventional magnetic symmetries. Our method can be implemented with spin-resolved scanning tunneling microscopy to study two-impurity Kondo phenomena on altermagnetic metals by varying the distance and orientation between magnetic impurities. Using quantum Monte Carlo simulations, we analyze the spin splitting of the Kondo resonance, whose spatial distribution sensitively captures the symmetry of the underlying altermagnetic order. Furthermore, the impurity spin correlations reflects the anisotropy of the RKKY interaction due to the altermagnetic Fermi surface splitting. This work provides a framework for studying the competition between the Kondo effect, the RKKY interaction and altermagnetism, in the simplest possible system.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Dissipation-Induced Order: The $S=1/2$ Quantum Spin Chain Coupled to an Ohmic Bath

We consider an $S=1/2$ antiferromagnetic quantum Heisenberg chain where each site is coupled to an independent bosonic bath with ohmic dissipation. The coupling to the bath preserves the global SO(3) spin symmetry. Using large-scale, approximation-free quantum Monte Carlo simulations, we show that any finite coupling to the bath suffices to stabilize long-range antiferromagnetic order. This is in stark contrast to the isolated Heisenberg chain where spontaneous breaking of the SO(3) symmetry is forbidden by the Mermin-Wagner theorem. A linear spin-wave theory analysis confirms that the memory of the bath and the concomitant retarded interaction stabilize the order. For the Heisenberg chain, the ohmic bath is a marginal perturbation so that exponentially large system sizes are required to observe long-range order at small couplings. Below this length scale, our numerics is dominated by a crossover regime where spin correlations show different power-law behaviors in space and time. We discuss the experimental relevance of this crossover phenomena.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Metallic and Deconfined Quantum Criticality in Dirac Systems

Motivated by the physics of spin-orbital liquids, we study a model of interacting Dirac fermions on a bilayer honeycomb lattice at half filling, featuring an explicit global SO(3)$\times$U(1) symmetry. Using large-scale auxiliary-field quantum Monte Carlo (QMC) simulations, we locate two zero-temperature phase transitions as function of increasing interaction strength. First, we observe a continuous transition from the weakly-interacting semimetal to a different semimetallic phase in which the SO(3) symmetry is spontaneously broken and where two out of three Dirac cones acquire a mass gap. The associated quantum critical point can be understood in terms of a Gross-Neveu-SO(3) theory. Second, we subsequently observe a transition towards an insulating phase in which the SO(3) symmetry is restored and the U(1) symmetry is spontaneously broken. While strongly first order at the mean-field level, the QMC data is consistent with a direct and continuous transition. It is thus a candidate for a new type of deconfined quantum critical point that features gapless fermionic degrees of freedom.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Nematic Quantum Criticality in Dirac Systems

We investigate nematic quantum phase transitions in two different Dirac fermion models. The models feature twofold and fourfold, respectively, lattice rotational symmetries that are spontaneously broken in the ordered phase. Using negative-sign-free quantum Monte Carlo simulations and an $ε$-expansion renormalization group analysis, we show that both models exhibit continuous phase transitions. In contrast to generic Gross-Neveu dynamical mass generation, the quantum critical regime is characterized by large velocity anisotropies, with fixed-point values being approached very slowly. Hence both experimental and numerical investigations will not be representative of the infrared fixed point, but of a crossover regime characterized by drifting exponents.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Thermodynamic and Dynamical Signatures of a Quantum Spin-Hall Insulator to Superconductor Transition

Thermodynamic and dynamical properties of a model of Dirac fermions with a deconfined quantum critical point (DQCP) separating an interaction-generated quantum spin-Hall insulator from an s-wave superconductor [Nature Comm.~{\bf 10}, 2658 (2019)] are studied by quantum Monte Carlo simulations. Inside the deconfined quantum critical region bound by the single-particle gap, spinons and spinless charge-2e skyrmions emerge. Since the model conserves total spin and charge, and has a single length scale, these excitations lead to a characteristic linear temperature dependence of the uniform spin and charge susceptibilities. At the DQCP, the order parameter dynamic structure factors show remarkable similarities that support emergent Lorentz symmetry. Above a critical temperature, superconductivity is destroyed by the proliferation of spin-1/2 vortices.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Valence-bond solid to antiferromagnet transition in the two-dimensional Su-Schrieffer-Heeger model by Langevin dynamics

The two-dimensional Su-Schrieffer-Heeger model of electrons coupled to quantum phonons is investigated using Langevin dynamics within the framework of auxiliary-field quantum Monte Carlo. Based on an explicit determination of the density of zeros of the fermion determinant, it is argued that the method is efficient in the challenging adiabatic limit. Large-scale simulations at the O(4)-symmetric point establish that the ground state of the 2D SSH model undergoes a transition from a $(π,π)$ valence bond solid to an antiferromagnet with increasing phonon frequency, yet still in the adiabatic regime. The single-particle spectrum illustrates the renormalization of the electronic band and suggests the existence of a gapped polaronic band, whereas the particle-hole channels show gapless modes associated with long-range bond and magnetic order, respectively. The simulations are supplemented with a mean-field analysis and a self-consistent Born proximation.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Automatic differentiation applied to excitations with Projected Entangled Pair States

The excitation ansatz for tensor networks is a powerful tool for simulating the low-lying quasiparticle excitations above ground states of strongly correlated quantum many-body systems. Recently, the two-dimensional tensor network class of infinite entangled pair states gained new ground state optimization methods based on automatic differentiation, which are at the same time highly accurate and simple to implement. Naturally, the question arises whether these new ideas can also be used to optimize the excitation ansatz, which has recently been implemented in two dimensions as well. In this paper, we describe a straightforward way to reimplement the framework for excitations using automatic differentiation, and demonstrate its performance for the Hubbard model at half filling.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Fermi-surface reconstruction without symmetry breaking

We present a sign-problem free quantum Monte Carlo study of a model that exhibits quantum phase transitions without symmetry breaking and associated changes in the size of the Fermi surface. The model is an Ising gauge theory on the square lattice coupled to an Ising matter field and spinful `orthogonal' fermions at half-filling, both carrying Ising gauge charges. In contrast to previous studies, our model hosts an electron-like, gauge-neutral fermion excitation providing access to Fermi liquid phases. One of the phases of the model is a previously studied orthogonal semi-metal, which has $\mathbb{Z}_2$ topological order, and Luttinger-volume violating Fermi points with gapless orthogonal fermion excitations. We elucidate the global phase diagram of the model: along with a conventional Fermi liquid phase with a large Luttinger-volume Fermi surface, we also find a `deconfined' Fermi liquid in which the large Fermi surface co-exists with fractionalized excitations. We present results for the electron spectral function, showing its evolution from the orthogonal semi-metal with a spectral weight near momenta $\{\pm π/2, \pm π/2\}$, to a large Fermi surface.

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

Phase diagram and dynamics of the SU($\boldsymbol N$) symmetric Kondo lattice model

In heavy-fermion systems, the competition between the local Kondo physics and intersite magnetic fluctuations results in unconventional quantum critical phenomena which are frequently addressed within the Kondo lattice model (KLM). Here we study this interplay in the SU($N$) symmetric generalization of the two-dimensional half-filled KLM by quantum Monte Carlo simulations with $N$ up to 8. While the long-range antiferromagnetic (AF) order in SU($N$) quantum spin systems typically gives way to spin-singlet ground states with spontaneously broken lattice symmetry, we find that the SU($N$) KLM is unique in that for each finite $N$ its ground-state phase diagram hosts only two phases -- AF order and the Kondo-screened phase. The absence of any intermediate phase between the $N=2$ and large-$N$ cases establishes adiabatic correspondence between both limits and confirms that the large-$N$ theory is a correct saddle point of the KLM fermionic path integral and a good starting point to include quantum fluctuations. In addition, we determine the evolution of the single-particle gap, quasiparticle residue of the doped hole at momentum $(π,π)$, and spin gap across the magnetic order-disorder transition. Our results indicate that increasing $N$ modifies the behavior of the coherence temperature: while it evolves smoothly across the magnetic transition at $N=2$ it develops an abrupt jump -- of up to an order of magnitude -- at larger but finite $N$. We discuss the magnetic order-disorder transition from a quantum-field-theoretic perspective and comment on implications of our findings for the interpretation of experiments on quantum critical heavy-fermion compounds.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Phases of the (2+1) dimensional SO(5) non-linear sigma model with topological term

We use the half-filled zeroth Landau level in graphene as a regularization scheme to study the physics of the SO(5) non-linear sigma model subject to a Wess-Zumino-Witten topological term in 2+1 dimensions. As shown by Ippoliti et al. [PRB 98, 235108 (2019)], this approach allows for negative sign free auxiliary field quantum Monte Carlo simulations. The model has a single free parameter, $U_0$, that monitors the stiffness. Within the parameter range accessible to negative sign free simulations, we observe an ordered phase in the large $U_0$ or stiff limit. Remarkably, upon reducing $U_0$ the magnetization drops substantially, and the correlation length exceeds our biggest system sizes, accommodating 100 flux quanta. The implications of our results for deconfined quantum phase transitions between valence bond solids and anti-ferromagnets are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Hubbard model on the honeycomb lattice: from static and dynamical mean-field theories to lattice quantum Monte Carlo simulations

We study the one-band Hubbard model on the honeycomb lattice using a combination of quantum Monte Carlo (QMC) simulations and static as well as dynamical mean-field theory (DMFT). This model is known to show a quantum phase transition between a Dirac semi-metal and the antiferromagnetic insulator. The aim of this article is to provide a detailed comparison between these approaches by computing static properties, notably ground-state energy, single-particle gap, double occupancy, and staggered magnetization, as well as dynamical quantities such as the single-particle spectral function. At the static mean-field level local moments cannot be generated without breaking the SU(2) spin symmetry. The DMFT approximation accounts for temporal fluctuations, thus captures both the evolution of the double occupancy and the resulting local moment formation in the paramagnetic phase. As a consequence, the DMFT approximation is found to be very accurate in the Dirac semi-metallic phase where local moment formation is present and the spin correlation length small. However, in the vicinity of the fermion quantum critical point the spin correlation length diverges and the spontaneous SU(2) symmetry breaking leads to low-lying Goldstone modes in the magnetically ordered phase. The impact of these spin fluctuations on the single-particle spectral function -- \textit{waterfall} features and narrow spin-polaron bands -- is only visible in the lattice QMC approach.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A search for correlation-induced adiabatic paths between distinct topological insulators

Correlations in topological states of matter provide a rich phenomenology, including a reduction in the topological classification of the interacting system compared to its non-interacting counterpart. This happens when two phases that are topologically distinct on the non-interacting level become adiabatically connected once interactions are included. We use a quantum Monte Carlo method to study such a reduction. We consider a 2D charge-conserving analog of the Levin-Gu superconductor whose classification is reduced from $\mathbb{Z}$ to $\mathbb{Z}_4$. We may expect any symmetry-preserving interaction that leads to a symmetric gapped ground state at strong coupling, and consequently a gapped symmetric surface, to be sufficient for such reduction. Here, we provide a counter example by considering an interaction which (i) leads to a symmetric gapped ground state at sufficient strength and (ii) does not allow for any adiabatic path connecting the trivial phase to the topological phase with $w=4$. The latter is established by numerically mapping the phase diagram as a function of the interaction strength and a parameter tuning the topological invariant. Instead of the adiabatic connection, the system exhibits an extended region of spontaneous symmetry breaking separating the topological sectors. Frustration reduces the size of this long-range ordered region until it gives way to a first order phase transition. Within the investigated range of parameters, there is no adiabatic path deforming the formerly distinct free fermion states into each other. We conclude that an interaction which trivializes the surface of a gapped topological phase is necessary but not sufficient to establish an adiabatic path within the reduced classification. In other words, the class of interactions which trivializes the surface is different from the class which establishes an adiabatic connection in the bulk.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Universal Fermi-surface anisotropy renormalization for interacting Dirac fermions with long-range interactions

Recent experimental and numerical evidence suggest an intriguing universal relationship between the Fermi surface anisotropy of the non-interacting parent two-dimensional electron gas and the strongly correlated composite Fermi liquid formed in a strong magnetic field close to half-filing. Inspired by these observations, we explore more generally the question of anisotropy renormalization in interacting 2D Fermi systems. Using a recently developed non-perturbative and numerically-exact projective quantum Monte Carlo simulation as well as other numerical and analytic techniques, only for Dirac fermions with long-range Coulomb interactions do we find a universal square-root decrease of the Fermi-surface anisotropy. For the half-filled composite Fermi liquid, this result is surprising since a Dirac fermion ground state was only recently proposed as an alternative to the usual HLR state. The importance of the long-range interaction, expected for Dirac systems, is also consistent with recent transport measurements. Our proposed universality can be tested in several anisotropic Dirac materials including graphene, topological insulators organic conductors, and magic-angle twisted bilayer graphene.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Continuous-time quantum Monte Carlo for fermion-boson lattice models: Improved bosonic estimators and application to the Holstein model

We extend the continuous-time interaction-expansion quantum Monte Carlo method with respect to measuring observables for fermion-boson lattice models. Using generating functionals, we express expectation values involving boson operators, which are not directly accessible because simulations are done in terms of a purely fermionic action, as integrals over fermionic correlation functions. We also demonstrate that certain observables can be inferred directly from the vertex distribution, and present efficient estimators for the total energy and the phonon propagator of the Holstein model. Furthermore, we generalize the covariance estimator of the fidelity susceptibility, an unbiased diagnostic for phase transitions, to the case of retarded interactions. The new estimators are applied to half-filled spinless and spinful Holstein models in one dimension. The observed renormalization of the phonon mode across the Peierls transition in the spinless model suggests a soft-mode transition in the adiabatic regime. The critical point is associated with a minimum in the phonon kinetic energy and a maximum in the fidelity susceptibility.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Edge instabilities of topological superconductors

Nodal topological superconductors display zero-energy Majorana flat bands at generic edges. The flatness of these edge bands, which is protected by time-reversal and translation symmetry, gives rise to an extensive ground-state degeneracy. Therefore, even arbitrarily weak interactions lead to an instability of the flat-band edge states towards time-reversal and translation-symmetry-broken phases, which lift the ground-state degeneracy. We examine the instabilities of the flat-band edge states of d_{xy}-wave superconductors by performing a mean-field analysis in the Majorana basis of the edge states. The leading instabilities are Majorana mass terms, which correspond to coherent superpositions of particle-particle and particle-hole channels in the fermionic language. We find that attractive interactions induce three different mass terms. One is a coherent superposition of imaginary s-wave pairing and current order, and another combines a charge-density-wave and finite-momentum singlet pairing. Repulsive interactions, on the other hand, lead to ferromagnetism together with spin-triplet pairing at the edge. Our quantum Monte Carlo simulations confirm these findings and demonstrate that these instabilities occur even in the presence of strong quantum fluctuations. We discuss the implications of our results for experiments on cuprate high-temperature superconductors.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Mott Quantum Criticality in the Anisotropic 2D Hubbard Model

We present evidence for Mott quantum criticality in an anisotropic two-dimensional system of coupled Hubbard chains at half-filling. In this scenario emerging from variational cluster approximation and cluster dynamical mean-field theory, the interchain hopping $t_{\perp}$ acts as a control parameter driving the second-order critical end point $T_c$ of the metal-insulator transition down to zero at $t_{\perp}^{c}/t\simeq 0.2$. Below $t_{\perp}^{c}$, the volume of the hole and electron Fermi pockets of a compensated metal vanishes continuously at the Mott transition. Above $t_{\perp}^{c}$, the volume reduction of the pockets is cut off by a first-order transition. We discuss the relevance of our findings to a putative quantum critical point in layered organic conductors, whose location remains elusive so far.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Thermodynamic and spectral properties of adiabatic Peierls chains

We present exact numerical results for the effects of thermal fluctuations on the experimentally relevant thermodynamic and spectral properties of Peierls chains. To this end, a combination of classical Monte Carlo sampling and exact diagonalization is used to study adiabatic half-filled Holstein and Su-Schrieffer-Heeger models. The classical nature of the lattice displacements in combination with parallel tempering permit simulations on large system sizes and a direct calculation of spectral functions in the frequency domain. Most notably, the long-range order and the associated Peierls gap give rise to a distinct low-temperature peak in the specific heat. The closing of the gap and suppression of order by thermal fluctuations involves in-gap excitations in the form of soliton-antisoliton pairs, and is also reflected in the dynamic density and bond structure factors as well as in the optical conductivity. We compare our data to the widely used mean-field approximation, and highlight relations to symmetry-protected topological phases and disorder problems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Dynamical mean-field study of partial Kondo screening in the periodic Anderson model on the triangular lattice

The competition between Kondo screening and indirect magnetic exchange is studied for a system with geometrical frustration using dynamical mean-field theory (DMFT). We systematically scan the weak- to strong-coupling regime of the periodic Anderson model on the triangular lattice for a wide range of fillings n. The magnetic phase diagram is derived using a site-dependent DMFT approach by self-consistent mapping onto three independent single-impurity models corresponding to the three correlated f orbitals in the unit cell. At half-filling, the system is a non-magnetic Kondo insulator for all considered interaction strengths U>0 which immediately develops into a non-magnetic metallic Kondo-singlet phase for fillings slightly below half-filling. On the other hand, indirect magnetic exchange between the f moments results in antiferromagnetic order at lower fillings. The antiferromagnetic and the Kondo-singlet phases are separated in the U-n phase diagram by an extended region of partial Kondo screening, i.e., a phase where the magnetic moment at one site in the unit cell is Kondo screened while the remaining two are coupled antiferromagnetically. At even lower fillings the system crosses over from a local-moment to a mixed-valence regime where the minimization of the kinetic energy in a strongly correlated system gives rise to a metallic and partially polarized ferromagnetic state.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Excitation spectra and correlation functions of quantum Su-Schrieffer-Heeger models

We study one-dimensional Su-Schrieffer-Heeger (SSH) models with quantum phonons using a continuous-time quantum Monte Carlo method. Within statistical errors, we obtain identical results for the SSH model with acoustic phonons, and a related model with a coupling to an optical bond phonon mode. Based on this agreement, we first study the Peierls metal-insulator transition of the spinless SSH model, and relate it to the Kosterlitz-Thouless transition of a spinless Luttinger liquid. In the Peierls phase, the spectral functions reveal the single-particle and charge gap, and a central peak related to long-range order. For the spinful SSH model, which has a dimerized ground state for any nonzero coupling, we reveal a symmetry-related degeneracy of spin and charge excitations, and the expected spin and charge gaps as well as a central peak. Finally, we study the SSH-$UV$ model with electron-phonon and electron-electron interaction. We observe a Mott phase with critical spin and bond correlations at weak electron-phonon coupling, and a Peierls phase with gapped spin excitations at strong coupling. We relate our findings to the extended Hubbard model, and discuss the physical origin of the agreement between optical and acoustic phonons.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fermionic quantum criticality in honeycomb and $π$-flux Hubbard models: Finite-size scaling of renormalization-group-invariant observables from quantum Monte Carlo

We numerically investigate the critical behavior of the Hubbard model on the honeycomb and the $π$-flux lattice, which exhibits a direct transition from a Dirac semimetal to an antiferromagnetically ordered Mott insulator. We use projective auxiliary-field quantum Monte Carlo simulations and a careful finite-size scaling analysis that exploits approximately improved renormalization-group-invariant observables. This approach, which is successfully verified for the three-dimensional XY transition of the Kane-Mele-Hubbard model, allows us to extract estimates for the critical couplings and the critical exponents. The results confirm that the critical behavior for the semimetal to Mott insulator transition in the Hubbard model belongs to the Gross-Neveu-Heisenberg universality class on both lattices.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Phonon spectral function of the one-dimensional Holstein-Hubbard model

We use the continuous-time interaction expansion (CT-INT) quantum Monte Carlo method to calculate the phonon spectral function of the one-dimensional Holstein-Hubbard model at half-filling. Our results are consistent with a soft-mode Peierls transition in the adiabatic regime, and the existence of a central peak related to long-range order in the Peierls phase. We explain a previously observed feature at small momenta in terms of a hybridization of charge and phonon excitations. Tuning the system from a Peierls to a metallic phase with a nonzero Hubbard interaction suppresses the central peak, but a significant renormalization of the phonon dispersion remains. In contrast, the dispersion is only weakly modified in the Mott phase. We discuss finite-size effects, the relation to the dynamic charge structure factor, as well as additional sum rules and their implications. Finally, we reveal the existence of a discrete symmetry in a continuum field theory of the Holstein model, which is spontaneously broken in the Peierls phase.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Ring-exchange periodic Anderson model for 3He bilayers

We present numerical results of a model calculation for the 3He bilayer system, which captures the interplay between fast and slow dynamics of the different layers and incorporates an independent scale for the three body ring exchange. By means of Cluster-DMFT+CTQMC, we find remarkable similarities with the experiments: a suppression of the coherence temperature upon approaching the solidification point of the first layer accompanied by the onset of strong ferromagnetic correlations within the first layer. Based on the consistent results for different cluster sizes we conjecture a first-order transition cutting short the Kondo breakdown, which allows us to interpret the experimental observation of an intervening phase preempting the quantum critical point in terms of a phase separation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Spin and charge dynamics of a quasi-one-dimensional antiferromagnetic metal

We use quantum Monte Carlo simulations to study a finite-temperature dimensional-crossover-driven evolution of spin and charge dynamics in weakly coupled Hubbard chains with a half-filled band. The low-temperature behavior of the charge gap indicates a crossover between two distinct energy scales: a high-energy one-dimensional (1D) Mott gap due to the umklapp process and a low-energy gap which stems from long-range antiferromagnetic (AF) fluctuations. Away from the 1D regime and at temperature scales above the charge gap, the emergence of a zero-frequency Drude-like feature in the interchain optical conductivity $σ_{\perp}(ω)$ implies the onset of a higher-dimensional metal. In this metallic phase, enhanced quasiparticle scattering off finite-range AF fluctuations results in incoherent single-particle dynamics. The coupling between spin and charge fluctuations is also seen in the spin dynamical structure factor $S({\pmb q},ω)$ displaying damped spin excitations (paramagnons) close to the AF wave-vector ${\pmb q}=(π,π)$ and particle-hole continua near 1D momentum transfers spanning quasiparticles at the Fermi surface. We relate our results to the charge deconfinement in quasi-1D organic Bechgaard-Fabre salts.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Stable Quantum Monte Carlo Simulations for Entanglement Spectra of Interacting Fermions

We show that the two recently proposed methods to compute Renyi entanglement entropies in the realm of determinant quantum Monte Carlo methods for fermions are in principle equivalent, but differ in sampling strategies. The analogy allows to formulate a numerically stable calculation of the entanglement spectrum at strong coupling. We demonstrate the approach by studying static and dynamical properties of the entanglement hamiltonian across the interaction driven quantum phase transition between a topological insulator and quantum antiferromagnet in the Kane-Mele Hubbard model. The formulation is not limited to fermion systems and can readily be adapted to world-line based simulations of bosonic systems.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dynamically generated edge states in topological Kondo insulators

Kondo insulators combine strong electronic correlations with spin orbit coupling and thereby provide a potential realization of correlated topological insulators. We present model calculations which allow us to study the onset of bulk coherence and concomitant topological edge states from the mixed valence to local moment regimes. Our real-space dynamical mean-field results include the detailed temperature dependence of the single particle spectral function on slab geometries as well as the temperature dependence of the topological invariant. The relevance of our calculations for candidate materials like SmB6 is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Equivalence of Rashba-Hubbard and Hubbard chains

We review the fact that U(1) gauge symmetry enables the mapping of one-dimensional Hubbard chains with Rashba-type spin orbit coupling to renormalized Hubbard Hamiltonians. The existence of the mapping has important consequences for the interpretation of ARPES experiments on one-dimensional chains subject to Rashba spin orbit coupling and can be exploited to check for the applicability of the mapping. We show numerical applications of the mapping and consider the implications for bosonization as well as for the Heisenberg-limit of the Hubbard model. In addition we point out the consequences for various generalized Hubbard models.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Kondo screening of spin-charge separated fluxons by a helical liquid

The insertion of a magnetic $π$ flux into a quantum spin Hall insulator creates four localized, spin-charge separated states: the charge and spin fluxons with either charge $Q=\pm1$ or spin $S_z=\pm1/2$, respectively. In the presence of repulsive Coulomb interactions, the charged states are gapped out and a local moment is formed. We consider the Kane-Mele-Hubbard model on a ribbon with zigzag edges to construct an impurity model where the spin fluxon is screened by the helical edge liquid. In the noninteracting model, the hybridization between fluxon and edge states is dominated by the extent of the latter. It becomes larger with increasing spin-orbit coupling $λ$ but only has nonzero values for even distances between the $π$ flux and the edge. For the interacting system, we use the continuous-time quantum Monte Carlo method, which we have extended by global susceptibility measurements to reproduce the characteristic Curie law of the spin fluxon. However, due to the finite extent of the fluxons, the local moment is formed at rather low energies. The screening of the spin fluxon leads to deviations from the Curie law that follow the universal behavior obtained from a data collapse. Additionally, the Kondo resonance arises in the local spectral function between the two low-lying Hubbard peaks.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Rashba coupling and magnetic order in correlated helical liquids

We study strongly correlated helical liquids with and without Rashba coupling using quantum Monte Carlo simulations of the Kane-Mele model with a Hubbard interaction at the edge. Independent of the Rashba coupling, we find that interactions enhance spin correlations and suppress the spectral weight at the Fermi level. For sufficiently strong interactions, a gap can be observed in the single-particle spectral function. However, based on a finite-size scaling analysis and theoretical arguments, we argue that this gap is closed by order parameter fluctuations in the Luttinger liquid phase even at zero temperature, and filled in by thermally induced kinks in the order parameter in the Mott phase at finite temperatures. While the bosonization suggests an umklapp-driven Mott transition only in the presence of Rashba coupling and hence a significant impact of the latter, our numerical results are almost unaffected by Rashba coupling even at low temperatures.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Zero-modes and global antiferromagnetism in strained graphene

A novel magnetic ground state is reported for the Hubbard Hamiltonian in strained graphene. When the chemical potential lies close to the Dirac point, the ground state exhibits locally both the Néel and ferromagnetic orders, even for weak Hubbard interaction. Whereas the Néel order parameter remains of the same sign in the entire system, the magnetization at the boundary takes the opposite sign from the bulk. The total magnetization this way vanishes, and the magnetic ground state is globally only an antiferromagnet. This peculiar ordering stems from the nature of the strain-induced single particle zero-energy states, which have support on one sublattice of the honeycomb lattice in the bulk, and on the other sublattice near the boundary of a finite system. We support our claim with the self-consistent numerical calculation of the order parameters, as well as by the Monte Carlo simulations of the Hubbard model in both uniformly and non-uniformly strained honeycomb lattice. The present result is contrasted with the magnetic ground state of the same Hubbard model in the presence of a true magnetic field (and for vanishing Zeeman coupling), which is exclusively Néel ordered, with zero local magnetization everywhere in the system.

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

Excitation spectra and spin gap of the half-filled Holstein-Hubbard model

Single-- and two-particle excitation spectra of the one-dimensional, half-filled Holstein-Hubbard model are calculated using the continuous-time quantum Monte Carlo method. In the metallic phase, the results are consistent with a Luther-Emery liquid that has gapped spin and single-particle excitations but a gapless charge mode. However, given the initially exponential dependence of the spin gap on the backscattering matrix element, the numerical excitation spectra appear gapless in the weak-coupling regime, and therefore resemble those of a Luttinger liquid. The Mott phase has the expected charge gap and gapless spin excitations. The Peierls state shows a charge, spin and single-particle gap, a soft phonon mode, backfolded shadow bands and soliton excitations. Arguments and numerical evidence for the existence of a nonzero spin gap throughout the metallic phase are provided in terms of equal-time spin and charge correlation functions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Interaction-driven transition between topological states in a Kondo insulator

Heavy fermion materials naturally combine strong spin-orbit interactions and electronic correlations. When there is precisely one conduction electron per impurity spin, the coherent heavy fermion state is insulating. This Kondo insulating state has recently been argued to belong to the class of quantum spin Hall states. Motivated by this conjecture and a very recent experimental realization of this state, we investigate a model for Kondo insulators with spin-orbit coupling. Using DMFT, we observe an interaction-driven transition between two distinct topological states, indicated by a closing of the bulk gap and a simultaneous change of the Z2 topological invariant. At large interaction strength we find a topological heavy fermion state, characterized by strongly renormalized heavy bulk bands, hosting a pair of zero-energy edge modes. The model allows a detailed understanding of the temperature dependence of the single particle spectral function and in particular the energy scales at which one observes the appearance of edge states within the bulk gap.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Pinning the order: the nature of quantum criticality in the Hubbard model on honeycomb lattice

In numerical simulations, spontaneously broken symmetry is often detected by computing two-point correlation functions of the appropriate local order parameter. This approach, however, computes the square of the local order parameter, and so when it is {\it small}, very large system sizes at high precisions are required to obtain reliable results. Alternatively, one can pin the order by introducing a local symmetry breaking field, and then measure the induced local order parameter infinitely far from the pinning center. The method is tested here at length for the Hubbard model on honeycomb lattice, within the realm of the projective auxiliary field quantum Monte Carlo algorithm. With our enhanced resolution we find a direct and continuous quantum phase transition between the semi-metallic and the insulating antiferromagnetic states with increase of the interaction. The single particle gap in units of the Hubbard $U$ tracks the staggered magnetization. An excellent data collapse is obtained by finite size scaling, with the values of the critical exponents in accord with the Gross-Neveu universality class of the transition.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Single magnetic impurities in the Kane-Mele model

The realization of the spin-Hall effect in quantum wells has led to a plethora of studies regarding the properties of the edge states of a 2D topological insulator. These edge states constitute a class of one-dimensional liquids, called the helical liquid, where an electron's spin quantization axis is tied to its momentum. In contrast to one dimensional conductors, magnetic impurities - below the Kondo temperature - cannot block transport and one expects the current to circumvent the impurity. To study this phenomenon, we consider the single impurity Anderson model embedded into an edge of a Kane-Mele ribbon with up to 512x80 sites and use the numerically exact continuous time QMC method to study the Kondo effect. We present results on the temperature dependence of the spectral properties of the impurity and the bulk system that show the behaviour of the system in the various regimes of the Anderson model. A view complementary to the single particle spectral functions can be obtained using the spatial behaviour of the spin-spin correlation functions. Here we show the characteristic, algebraic decay in the edge channel near the impurity.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Spinon confinement: dynamics of weakly coupled Hubbard chains

Using large-scale determinant quantum Monte Carlo simulations in combination with the stochastic analytical continuation, we study two-particle dynamical correlation functions in the anisotropic square lattice of weakly coupled one-dimensional (1D) Hubbard chains at half-filling and in the presence of weak frustration. The evolution of the static spin structure factor upon increasing the interchain coupling is suggestive of the transition from the power-law decay of spin-spin correlations in the 1D limit to long-range antiferromagnetic order in the quasi-1D regime and at $T=0$. In the numerically accessible regime of interchain couplings, the charge sector remains gapped. The low-energy momentum dependence of the spin excitations is well described by the linear spin-wave theory with the largest intensity located around the antiferromagnetic wave vector. This magnon mode corresponds to a bound state of two spinons. At higher energies the spinons deconfine and we observe signatures of the two-spinon continuum which progressively fade away as a function of interchain hopping.

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

Dimensional-crossover-driven Mott transition in the frustrated Hubbard model

We study the Mott transition in a frustrated Hubbard model with next-nearest neighbor hopping at half-filling. The interplay between interaction, dimensionality and geometric frustration closes the one-dimensional Mott gap and gives rise to a metallic phase with Fermi surface pockets. We argue that they emerge as a consequence of remnant one-dimensional Umklapp scattering at the momenta with vanishing interchain hopping matrix elements. In this pseudogap phase, enhanced d-wave pairing correlations are driven by antiferromagnetic fluctuations. Within the adopted cluster dynamical mean-field theory on the $8\times 2$ cluster and down to our lowest temperatures the transition from one to two dimensions is continuous.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Effect of Electron-Phonon Interaction Range for a Half-Filled Band in One Dimension

We demonstrate that fermion-boson models with nonlocal interactions can be simulated at finite band filling with the continuous-time quantum Monte Carlo method. We apply this method to explore the influence of the electron-phonon interaction range for a half-filled band in one dimension, covering the full range from the Holstein to the Fröhlich regime. The phase diagram contains metallic, Peierls, and phase-separated regions, which we characterize in terms of static and dynamical correlation functions. In particular, our results reveal a suppression of $2k_F$ charge correlations with increasing interaction range, allowing for a power-law decay comparable to the pairing correlations.

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
preprint2012arXiv

Peierls to superfluid crossover in the one-dimensional, quarter-filled Holstein model

We use continuous-time quantum Monte Carlo simulations to study retardation effects in the metallic, quarter-filled Holstein model in one dimension. Based on results which include the one- and two-particle spectral functions as well as the optical conductivity, we conclude that with increasing phonon frequency the ground state evolves from one with dominant diagonal order---2k_F charge correlations---to one with dominant off-diagonal fluctuations, namely s-wave pairing correlations. In the parameter range of this crossover, our numerical results support the existence of a spin gap for all phonon frequencies. The crossover can hence be interpreted in terms of preformed pairs corresponding to bipolarons, which are essentially localised in the Peierls phase, and "condense" with increasing phonon frequency to generate dominant pairing correlations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Understanding the Josephson current through a Kondo-correlated quantum dot

We study the Josephson current 0-$π$ transition of a quantum dot tuned to the Kondo regime. The physics can be quantitatively captured by the numerically exact continuous time quantum Monte Carlo method applied to the single-impurity Anderson model with BCS superconducting leads. For a comparison to an experiment the tunnel couplings are determined by fitting the normal-state linear conductance. Excellent agreement for the dependence of the critical Josephson current on the level energy is achieved. For increased tunnel couplings the Kondo scale becomes comparable to the superconducting gap and the regime of the strongest competition between superconductivity and Kondo correlations is reached; we predict the gate voltage dependence of the critical current in this regime.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A time and spatially resolved quench of the fermionic Hubbard model showing restricted equilibration

We investigate the quench of half-filled 1D and 2D fermionic Hubbard models to models without Coulomb interaction. Since the time propagation is gaussian we can use a variety of time-dependent quantum Monte Carlo methods to tackle this problem without generating a dynamical sign problem. Using a continuous time quantum Monte Carlo method (CTQMC) we achieve a system size of 128 sites in 1D, and using a Blankenbecler-Scalapino-Sugar (BSS) type algorithm we were able to simulate 20 x 20 square lattices. Applying these methods to study the dynamics after the quench, we observe that the final state of the system can be reasonably well described by a thermal single-particle density matrix that takes the initial single particle conservation laws into account. The characteristic decay towards this limit is found to be oscillatory with an additional power law decay that depends on the dimensionality. This numerically exact result is shown to compare favorable to mean-field approximations as well as to perturbation theory. Furthermore we observe the information propagation in the 1D-case in the charge charge and spin spin correlations and find that it is linear with a velocity of roughly v = 4 in units of the hopping amplitude.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dynamic charge correlations near the Peierls transition

The quantum phase transition between a repulsive Luttinger liquid and an insulating Peierls state is studied in the framework of the one-dimensional spinless Holstein model. We focus on the adiabatic regime but include the full quantum dynamics of the phonons. Using continuous-time quantum Monte Carlo simulations, we track in particular the dynamic charge structure factor and the single-particle spectrum across the transition. With increasing electron-phonon coupling, the dynamic charge structure factor reveals the emergence of a charge gap, and a clear signature of phonon softening at the zone boundary. The single-particle spectral function evolves continuously across the transition. Hybridization of the charge and phonon modes of the Luttinger liquid description leads to two modes, one of which corresponds to the coherent polaron band. This band acquires a gap upon entering the Peierls phase, whereas the other mode constitutes the incoherent, high-energy spectrum with backfolded shadow bands. Coherent polaronic motion is a direct consequence of quantum lattice fluctuations. In the strong-coupling regime, the spectrum is described by the static, mean-field limit. Importantly, whereas finite electron density in general leads to screening of polaron effects, the latter reappear at half filling due to charge ordering and lattice dimerization.

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

Exact diagonalization study of the tunable edge magnetism in graphene

The tunable magnetism at graphene edges with lengths of up to 48 unit cells is analyzed by an exact diagonalization technique. For this we use a generalized interacting one-dimensional model which can be tuned continuously from a limit describing graphene zigzag edge states with a ferromagnetic phase, to a limit equivalent to a Hubbard chain, which does not allow ferromagnetism. This analysis sheds light onto the question why the edge states have a ferromagnetic ground state, while a usual one-dimensional metal does not. Essentially we find that there are two important features of edge states: (a) umklapp processes are completely forbidden for edge states; this allows a spin-polarized ground state. (b) the strong momentum dependence of the effective interaction vertex for edge states gives rise to a regime of partial spin-polarization and a second order phase transition between a standard paramagnetic Luttinger liquid and ferromagnetic Luttinger liquid.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Melting of stripe phases and its signature in the single-particle spectral function

Motivated by the recent experimental data [Phys. Rev. B 79, 100502 (2009)] indicating the existence of a pure stripe charge order over unprecedently wide temperature range in La_{1.8-x}Eu_{0.2}Sr_xCuO_4, we investigate the temperature-induced melting of the metallic stripe phase. In spite of taking into account local dynamic correlations within a real-space dynamical mean-field theory of the Hubbard model, we observe a mean-field like melting of the stripe order irrespective of the choice of the next-nearest neighbor hopping. The temperature dependence of the single-particle spectral function shows the stripe induced formation of a flat band around the antinodal points accompanied by the opening a gap in the nodal direction.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A weak coupling CTQMC study of the single impurity and periodic Anderson models with s-wave superconducting baths

We apply the unbiased weak-coupling continuous time quantum Monte Carlo (CTQMC) method to review the physics of a single magnetic impurity coupled to s-wave superconducting leads described by the BCS reduced Hamiltonian. As a function of the superconducting gap $Δ$, we study the signature of the first order transition between the singlet and doublet (local moment) states on various quantities. In particular we concentrate on the Josephson current with 0 to $π$ phase shift, the crossing of the Andreev bound states in the single particle spectral function, as well as the local dynamical spin structure factor. Within DMFT, this impurity problem provides a link to the periodic Anderson model with superconducting conduction electrons (BCS-PAM). The first order transition observed in the impurity model is reproduced in the BCS-PAM and is signalized by the crossing of the low energy excitations in the local density of states. The momentum resolved single particle spectral function in the singlet state reveals the coherent, Bloch-like, superposition of Andreev bound states. In the doublet or local moment phase the single particle spectral function is characterized by incoherent quasiparticle excitations.

0 citations0 reviewsNo code linkedOpen access