Fakher F. Assaad

The fermion disorder operator has been shown to reveal the entanglement information in 1D Luttinger liquids and 2D free and interacting Fermi and non-Fermi liquids emerging at quantum critical points(QCP). Here we study, by means of large-scale quantum Monte Carlo simulation, the scaling behavior of disorder operator in correlated Dirac systems. We first demonstrate the logarithmic scaling behavior of the disorder operator at the Gross-Neveu (GN) chiral Ising and Heisenberg QCPs, where consistent conformal field theory (CFT) content of the GN-QCP in its coefficient is found. Then we study a 2D monopole free deconfined quantum critical point (DQCP) realized between a quantum-spin Hall insulator and a superconductor. Our data point to negative values of the logarithmic coefficients such that the DQCP does not correspond to a unitary CFT. Density matrix renormalization group calculations of the disorder operator on a 1D DQCP model also detect emergent continuous symmetries.

The fermion disorder operator has been shown to reveal the entanglement information in 1D Luttinger liquids and 2D free and interacting Fermi and non-Fermi liquids emerging at quantum critical points (QCPs) [W. Jiang et al., arXiv:2209.07103]. Here we study, by means of large-scale quantum Monte Carlo simulation, the scaling behavior of the disorder operator in correlated Dirac systems. We first demonstrate the logarithmic scaling behavior of the disorder operator at the Gross-Neveu (GN) chiral Ising and Heisenberg QCPs, where consistent conformal field theory (CFT) content of the GN-QCP in its coefficient is found. Then we study a 2D monopole-free deconfined quantum critical point (DQCP) realized between a quantum-spin Hall insulator and a superconductor. Our data point to negative values of the logarithmic coefficients such that the DQCP does not correspond to a unitary CFT. Density matrix renormalization group calculations of the disorder operator on a 1D DQCP model also detect emergent continuous symmetries.

In this paper, we consider a path integral formulation of the Hubbard model based on a Hubbard-Stratonovich transformation that couples the auxiliary field to the local electronic density. This decoupling is known to have a saddle-point structure that shows a remarkable regularity: The field configuration at each saddle point can be understood in terms of a set of elementary field configurations localized in space and imaginary time which we coin instantons. The interaction between instantons is short ranged. Here, we formulate a classical partition function for the instanton gas that has predictive power. For a given set of physical parameters, we can predict the distribution of instantons and show that the instanton number is sharply defined in the thermodynamic limit, thereby defining a unique dominant saddle point. Decoupling in the charge channel conserves SU(2) spin symmetry for each field configurations. Hence, the instanton approach provides an SU(2) spin-symmetric approximation to the Hubbard model. It fails, however, to capture the magnetic transition inherent to the Hubbard model on the honeycomb lattice despite being able to describe local moment formation. In fact, the instanton itself corresponds to local moment formation and concomitant short-ranged antiferromagnetic correlations. This aspect is also seen in the single particle spectral function that shows clear signs of the upper and lower Hubbard bands. Our instanton approach bears remarkable similarities to local dynamical approaches, such as dynamical mean-field theory, in the sense that it has the unique property of allowing for local moment formation without breaking the SU(2) spin symmetry. In contrast to local approaches, it captures short-ranged magnetic fluctuations. Furthermore, it also offers possibilities for systematic improvements by taking into account fluctuations around the dominant saddle point. Finally, we show that the saddle point structure depends upon the choice of lattice geometry. For the square lattice at half filling, the saddle-point structure reflects the itinerant to localized nature of the magnetism as a function of the coupling strength. The implications of our results for Lefschetz thimble approaches to alleviate the sign problem are also discussed.

We propose a method to mitigate heavy-tailed distributions in fermion quantum Monte Carlo simulations originating from zeros of the fermion determinant. In this case, the second moment of the observables might not be well defined, and we show that by merely changing the synchronization between local updates and computation of observables, one can reduce the prefactor of the heavy-tailed distribution, thus substantially suppressing statistical fluctuations of observables. We also show that the average, or the first moment, is well defined and hence is independent on our measuring scheme. The method is especially suitable for local observables similar to, e.g., double occupancy, where the resulting speedup can reach two orders of magnitude. For observables containing spatial correlators, the speedup is more moderate but still ranges between five and ten. Our results are independent on the nature of the auxiliary field, discrete or continuous, and pave the way to improve measurement strategies for hybrid Monte Carlo simulations.

We propose a method to mitigate heavy-tailed distributions in fermion Quantum Monte Carlo simulations originating from zeros of the fermion determinant. In this case the second moment of the observables might be not well defined, and we show that by merely changing the synchronization between local updates and computation of observables, one can reduce the prefactor of the heavy-tailed distribution, thus substantially suppressing statistical fluctuations of observables. We also show that the average, or the first moment, is well defined and hence is independent on our measuring scheme. The method is especially suitable for local observables similar to e.g. double occupancy, where the resulting speedup can reach two orders of magnitude. For observables, containing spatial correlators, the speedup is more moderate, but still ranges between five and ten. Our results are independent on the nature of the auxiliary field, discrete or continuous, and pave the way to improve measurement strategies for Hybrid Monte Carlo simulations.

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.