Source author record

Maksim Ulybyshev

Maksim Ulybyshev 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-el hep-lat hep-th

Catalog footprint

What is connected

3works

3topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

Mitigating spikes in fermion Monte Carlo methods by reshuffling measurements

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.

preprint2020arXiv

Lefschetz thimbles decomposition for the Hubbard model on the hexagonal lattice

We propose a framework to study the properties of the Lefschetz thimbles decomposition for lattice fermion models approaching the thermodynamic limit. The proposed set of algorithms includes the Schur complement solver and the exact computation of the derivatives of the fermion determinant. It allows us to solve the gradient flow (GF) equations taking into account the fermion determinant exactly, with high performance. We can find both real and complex saddle points and describe the structure of the Lefschetz thimbles decomposition for large enough lattices to extrapolate the results to the thermodynamic limit. The algorithms are described for a general lattice fermion model, with emphasis on two types of lattice discretizations for relativistic fermions (staggered and Wilson), as well as on interacting tight-binding models for condensed matter systems. We apply these algorithms to the Hubbard model on a hexagonal lattice, dealing with lattice volumes as large as 12x12 with 256 steps in Euclidean time, in order to capture the properties of the thimbles decomposition as the thermodynamic, low-temperature, and continuum limits are approached. The complexity of the thimbles decomposition appears to be dependent on the form of the Hubbard-Stratonovich (HS) transformation. In particular, the evidence is provided for the existence of an optimal regime for the hexagonal lattice Hubbard model, with a reduced number of thimbles being important in the overall sum. We have performed quantum Monte Carlo (QMC) simulations using GF to deform the integration contour into the complex plane and demonstrated the agreement with exact diagonalization on small volumes (8 sites in space) if optimal setup for HS transformation is used. The residual sign problem is compared with the state-of-the-art BSS-QMC. We show that the average sign can be kept substantially higher using the Lefschetz thimbles approach.

preprint2016arXiv

Competing order in the fermionic Hubbard model on the hexagonal graphene lattice

We study the phase diagram of the fermionic Hubbard model on the hexagonal lattice in the space of on-site and nearest neighbor couplings with Hybrid-Monte-Carlo simulations. With pure on-site repulsion this allows to determine the critical coupling strength for spin-density wave formation with the standard approach of introducing a small mass term, explicitly breaking the sublattice symmetry. The analogous mass term for charge-density wave formation above a critical nearest-neighbor repulsion, on the other hand, would introduce a fermion sign problem. The competition between the two and the phase diagram in the space of the two coouplings can however be studied in simulations without explicit sublattice symmetry breaking. Our results compare qualitatively well with the Hartree-Fock phase diagram. We furthermore demonstrate how spin-symmetry breaking by the Euclidean time discretization can be avoided also, when using an improved fermion action based on an exponetial transfer matrix with exact sublattice symmetry.