Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userIntroduction to Quantum Monte Carlo simulations for fermionic systems
77
0
0.0
(
0
)
Added by
Raimundo R. dos Santos
Publication date
2003
fields
Physics
and research's language is
English
Authors
Raimundo R. dos Santos
Ask ChatGPT about the research
No Arabic abstract
We tutorially review the determinantal Quantum Monte Carlo method for fermionic systems, using the Hubbard model as a case study. Starting with the basic ingredients of Monte Carlo simulations for classical systems, we introduce aspects such as importance sampling, sources of errors, and finite-size scaling analyses. We then set up the preliminary steps to prepare for the simulations, showing that they are actually carried out by sampling discrete Hubbard-Stratonovich auxiliary fields. In this process the Greens function emerges as a fundamental tool, since it is used in the updating process, and, at the same time, it is directly related to the quantities probing magnetic, charge, metallic, and superconducting behaviours. We also discuss the as yet unresolved minus-sign problem, and two ways to stabilize the algorithm at low temperatures.
rate research
Read More
We introduce and compare three different Monte Carlo determinantal algorithms that allow one to compute dynamical quantities, such as the self-energy, of fermionic systems in their thermodynamic limit. We show that the most efficient approach expresses the sum of a factorial number of one-particle-irreducible diagrams as a recursive sum of determinants with exponential complexity. By comparing results for the two-dimensional Hubbard model with those obtained from state-of-the-art diagrammatic Monte Carlo, we show that we can reach higher perturbation orders and greater accuracy for the same computational effort.
Quantum Monte Carlo simulations, while being efficient for bosons, suffer from the negative sign problem when applied to fermions - causing an exponential increase of the computing time with the number of particles. A polynomial time solution to the sign problem is highly desired since it would provide an unbiased and numerically exact method to simulate correlated quantum systems. Here we show, that such a solution is almost certainly unattainable by proving that the sign problem is NP-hard, implying that a generic solution of the sign problem would also solve all problems in the complexity class NP (nondeterministic polynomial) in polynomial time.
This review summarizes recent developments in the study of fermionic quantum criticality, focusing on new progress in numerical methodologies, especially quantum Monte Carlo methods, and insights that emerged from recently large-scale numerical simulations. Quantum critical phenomena in fermionic systems have attracted decades of extensive research efforts, partially lured by their exotic properties and potential technology applications and partially awaked by the profound and universal fundamental principles that govern these quantum critical systems. Due to the complex and non-perturbative nature, these systems belong to the most difficult and challenging problems in the study of modern condensed matter physics, and many important fundamental problems remain open. Recently, new developments in model design and algorithm improvements enabled unbiased large-scale numerical solutions to be achieved in the close vicinity of these quantum critical points, which paves a new pathway towards achieving controlled conclusions through combined efforts of theoretical and numerical studies, as well as possible theoretical guidance for experiments in heavy-fermion compounds, Cu-based and Fe-based superconductors, ultra-cold fermionic atomic gas, twisted graphene layers, etc., where signatures of fermionic quantum criticality exist.
We discuss a projector Monte Carlo method for quantum spin models formulated in the valence bond basis, using the S=1/2 Heisenberg antiferromagnet as an example. Its singlet ground state can be projected out of an arbitrary basis state as the trial state, but a more rapid convergence can be obtained using a good variational state. As an alternative to first carrying out a time consuming variational Monte Carlo calculation, we show that a very good trial state can be generated in an iterative fashion in the course of the simulation itself. We also show how the properties of the valence bond basis enable calculations of quantities that are difficult to obtain with the standard basis of Sz eigenstates. In particular, we discuss quantities involving finite-momentum states in the triplet sector, such as the dispersion relation and the spectral weight of the lowest triplet.
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.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
