Researcher profile

Fakher F. Assaad

Fakher F. Assaad contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

14 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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.