arXiv:cond-mat/0101452 · Bulletin of National University of Uzbekistan Mathematics and Natural Sciences

Complex diffusion Monte-Carlo method: test by the simulation of the 2D fermions

Bakhodir Abdullaev, Mirzayousuf Musakhanov, Atsushi Nakamura

Published Sep 15, 2020• 1 citation

Abstract

On the base of the diffusion Monte-Carlo method we develop the method allowing to simulate the quantum systems with complex wave function. The method is exact and there are no approximations on the simulations of the module and the phase of the system's wave function. In our method averaged value of any quantity have no direct contribution from the phase of distribution function but only from the phase of the Green function of diffusion equation. We test the method by the simulations of the ground state of fermions in two-dimensional parabolic well. Anyons are used for the representation of the two-dimensional (2D) fermions. We vary the number of fermions from two to ten and find a good agreement of the numerical results with analytical ones for the numbers of the particles N>4.

DOI: 10.56017/2181-1318.1122

Related papers

Complex Diffusion Monte-Carlo method for the systems with complex wave function: test by the simulation of 2D electron in uniform magnetic field

B. Abdullaev, M.M. Musakhanov, Atsushi Nakamura

On the base of Diffusion Monte-Carlo method it is developed a new Complex Diffusion Monte-Carlo (CDMC) method allowing to simulate the quantum systems with complex wave function. There are no approximations on the calculation of modulus and phase of wave function in contrast to other methods. We find that the averaged value of any quantity in CDMC will have no direct contribution from the phase of the distribution function but only from the phase of the Green function of the diffusion equation. This is most important and crucial point of CDMC. We are testing CDMC by the calculation of the wave function and the ground state energy of two-dimensional electron placed into the external uniform magnetic field. There is an excellent agreement between simulations results and an analytical ones.

arXiv:cond-mat/0101330Jan 22, 20011 citation

Complex Diffusion Monte-Carlo method: tests by the simulations of 2D electron in magnetic field and 2D fermions-anyons in parabolic well

B. Abdullaev, M.M. Musakhanov, Atsushi Nakamura

We propose a new Complex Diffusion Monte Carlo (CDMC) method for the simulation of quantum systems with complex wave function. In CDMC the modulus and phase of wave function are simulated both in contrast to other methods. We successfully test CDMC by the simulation of the ground state for 2D electron in magnetic field and 2D fermions-anyons in parabolic well.

arXiv:cond-mat/0202256Feb 15, 20020 citations

Real-Time Diffusion Monte Carlo Method

Ilkka Ruokosenmäki, Tapio T. Rantala

Direct sampling of multi-dimensional systems with quantum Monte Carlo methods allows exact account of many-body effects or particle correlations. The most straightforward approach to solve the Schrödinger equation, Diffusion Monte Carlo, has been used in several benchmark cases for other methods to pursue. Its robustness is based on direct sampling of a positive probability density for diffusion in imaginary time. It has been argued that the corresponding real time diffusion can not be realised, because the corresponding oscillating complex valued distribution can not be used to drive diffusion. Here, we demonstrate that this can be done by turning the distribution piecewise positive and normalizable, and also, by using four types of walkers. This study is a proof of concept demonstration using the well-known and transparent case: one-dimensional harmonic oscillator. Furthermore, we show that our novel method can be used to find not only the ground state but also excited states and even the time evolution of a given wave function. Considering fermionic systems, this method may turn out to be feasible for finding the wave function nodes for other approaches.

arXiv:1703.05558Communications in Computational PhysicsOct 11, 20181 citation

Grid-based diffusion Monte Carlo for fermions without the fixed-node approximation

Alexander A. Kunitsa, So Hirata

A diffusion Monte Carlo algorithm is introduced that can determine the correct nodal structure of the wave function of a few-fermion system and its ground-state energy without an uncontrolled bias. This is achieved by confining signed random walkers to the points of a uniform infinite spatial grid, allowing them to meet and annihilate one another to establish the nodal structure without the fixed-node approximation. An imaginary-time propagator is derived rigorously from a discretized Hamiltonian, governing a non-Gaussian, sign-flipping, branching, and mutually annihilating random walk of particles. The accuracy of the resulting stochastic representations of a fermion wave function is limited only by the grid and imaginary-time resolutions and can be improved in a controlled manner. The method is tested for a series of model problems including fermions in a harmonic trap as well as the He atom in its singlet or triplet ground state. For the latter case, the energies approach from above with increasing grid resolution and converge within 0.015E_{h} of the exact basis-set-limit value for the grid spacing of 0.08 a.u. with a statistical uncertainty of 10^{-5}E_{h} without an importance sampling or Jastrow factor.

arXiv:1811.03000Physical review. EJan 27, 202012 citations

Diffusion Monte Carlo for fermions with replica reduction

Myung-Hoon Chung, D. P. Landau

We propose an improved framework of diffusion Monte Carlo method for quantum many-fermion systems. This approach enables us to find the ground state of the system without any a priori information on the nodal structure of the wave function. The main idea in this approach is to use sign-carrying replicas which are the basis vectors in a restricted Hilbert space. For efficiency, we propose to reduce the number of replicas in a proper manner. We test our method on a spinless fermion system, and compare the results with the exact values. We also discuss this method in relation to exact diagonalization.

Physical Review BMar 15, 20123 citations

Monte Carlo simulations of one-dimensional fermion systems

J. E. Hirsch, Robert L. Sugar, Douglas J. Scalapino +1 more

We discuss a new method to perform numerical simulations of one-dimensional systems with fermion and boson degrees of freedom. The method is based on a direct-space, imaginary-time representation of the fermion field. It is fast so that systems having up to 100 sites can easily be simulated. In addition, the method provides an intuitive physical "picture" of the ground state of a one-dimensional many-body system. We discuss in detail how to implement the method and how to compute various physical quantities. In particular, we show how to extend the method to study averages of off-diagonal quantities in an occupation-number representation. To assess the accuracy of our procedure, we apply it to free fermions in one dimension and compare with exact results. We then study a model of spinless interacting fermions and obtain the expected phase structure and behavior of correlation functions. We also consider the extended Hubbard model at various points in its phase diagram and study the behavior of spin-density, charge-density, and pairing correlation functions. We then study the Gross-Neveu model and show how the behavior depends on the number of fermion flavors. Finally, we consider an electron-phonon model and study its behavior both in the one-particle polaron sector and in the half-filled-band case. Along the way we show pictures of the ground-state configurations that give physical insight into the properties of the systems, like charge-density-wave, spin-density-wave, and superconducting states, "fractional charges," and solitons. We conclude by comparing our method with other methods and discuss the possibility of extending it to higher dimensions.

Physical review. B, Condensed matterNov 1, 1982487 citations