Journal of Physics Conference Series

Path integral Monte Carlo simulations of nanowires and quantum point contacts

John B. Shumway, Matthew J. Gilbert

Published Apr 1, 2006• 4 citations

Abstract

We study quantum wires and point contacts with a fixed-node path integral Monte Carlo technique. The fixed-node technique uses a variational principle to map fermionic problems into effective bosonic problems, which are then evaluated with standard quantum Monte Carlo techniques. While fixed-node is an approximation, it has the useful properties of being variational and being able to recover the exact answer when the exact nodes or phases of the density matrix are known. In these finite-temperature simulations we use the free particle density matrix as a fixed-node constraint to efficiently simulate hundreds of interacting electrons. We have calculated charge densities, pair correlation functions, and the currentcurrent Matsubura Green's functions for quantum wires and quantum point contacts.

DOI: 10.1088/1742-6596/35/1/017

Related papers

Path Integral Monte Carlo Methods for Fermions

Ethan Brown

In general, Quantum Monte Carlo methods suffer from a sign problem when simulating fermionic systems. This causes the efficiency of a simulation to decrease exponentially with the number of particles and inverse temperature. To circumvent this issue, a nodal constraint is often implemented, restricting the Monte Carlo procedure from sampling paths that cause the many-body density matrix to change sign. Unfortunately, this high-dimensional nodal surface is not a priori known unless the system is exactly solvable, resulting in uncontrolled errors. We will discuss two possible routes to extend the applicability of finite-temperatue path integral Monte Carlo. First we extend the regime where signful simulations are possible through a novel permutation sampling scheme. Afterwards, we discuss a method to variationally improve the nodal surface by minimizing a free energy during simulation. Applications of these methods will include both free and interacting electron gases, concluding with discussion concerning extension to inhomogeneous systems.

Jan 1, 20140 citations

Electronic structure quantum Monte Carlo

Michal Bajdich, Luboš Mitáš

Electronic structure quantum Monte Carlo Quantum Monte Carlo (QMC) is an advanced simulation methodology for studies of manybody quantum systems. The QMC approaches combine analytical insights with stochastic computational techniques for efficient solution of several classes of important many-body problems such as the stationary Schrödinger equation. QMC methods of various flavors have been applied to a great variety of systems spanning continuous and lattice quantum models, molecular and condensed systems, BEC-BCS ultracold condensates, nuclei, etc. In this review, we focus on the electronic structure QMC, i.e., methods relevant for systems described by the electron-ion Hamiltonians. Some of the key QMC achievements include direct treatment of electron correlation, accuracy in predicting energy differences and favorable scaling in the system size. Calculations of atoms, molecules, clusters and solids have demonstrated QMC applicability to real systems with hundreds of electrons while providing 90-95% of the correlation energy and energy differences typically within a few percent of experiments. Advances in accuracy beyond these limits are hampered by the so-called fixed-node approximation which is used to circumvent the notorious fermion sign problem. Many-body nodes of fermion states and their properties have therefore become one of the important topics for further progress in predictive power and efficiency of QMC calculations. Some of our recent results on the wave function nodes and related nodal domain topologies will be briefly reviewed. This includes analysis of few-electron systems and descriptions of exact and approximate nodes using transformations and projections of the highly-dimensional nodal hypersurfaces into the 3D space. Studies of fermion nodes offer new insights into topological properties of eigenstates such as explicit demonstrations that generic fermionic ground states exhibit the minimal number of two nodal domains. Recently proposed trial wave functions based on Pfaffians with pairing orbitals are presented and their nodal properties are tested in calculations of first row atoms and molecules. Finally, backflow "dressed" coordinates are introduced as another possibility for capturing correlation effects and for decreasing the fixed-node bias.

arXiv:1008.2369Acta Physica Slovaca Reviews and TutorialsApr 1, 200933 citations

Fixed-node Monte Carlo study of the two-dimensional Hubbard model

Guozhong An, J. M. J. van Leeuwen

Abstract This study extended the fixed-node QMC method to a lattice model for hopping fermions using Green’s function techniques. In general, the modifications of the procedures for continuous systems to adapt them to a lattice system were minor, and there appeared to be no surprises. Applications were made to the two-dimensional Hubbard model with different on-site interactions, lattice dimensions, boundary conditions, node specifications, and trial wavefunctions for importance sampling. Results were compared with those of Hartree-Fock and variational QMC calculations. For the HF calculations with zero interaction, the exact solutions are known, and the Green’s function QMC results agreed with those within statistical error. Trial wavefunctions specifying the nodes were obtained directly as those of HF determinantal solutions for zero interaction, or as variationally optimized determinantal functions with correlation terms. The procedure for fixed nodes was to eliminate walkers crossing the boundaries of a nodal region, a procedure not quite correct, but probably a reasonable approximation. The results in all cases were consistent with the HF and variational QMC calculations, giving lower energies than either for nonzero interaction energies. The existence of a domain-wall phase suggested by the HF and variational QMC results was confirmed by the fixed-node Green’s function results. It was noted that the fixed-node method would be applicable to systems other than the Hubbard model.

Jun 18, 20070 citations

Quantum Monte Carlo simulations of solids

W. M. C. Foulkes, Luboš Mitáš, Richard J. Needs +1 more

This article describes the variational and fixed-node diffusion quantum Monte Carlo methods and how they may be used to calculate the properties of many-electron systems. These stochastic wave-function-based approaches provide a very direct treatment of quantum many-body effects and serve as benchmarks against which other techniques may be compared. They complement the less demanding density-functional approach by providing more accurate results and a deeper understanding of the physics of electronic correlation in real materials. The algorithms are intrinsically parallel, and currently available high-performance computers allow applications to systems containing a thousand or more electrons. With these tools one can study complicated problems such as the properties of surfaces and defects, while including electron correlation effects with high precision. The authors provide a pedagogical overview of the techniques and describe a selection of applications to ground and excited states of solids and clusters.

Reviews of Modern PhysicsJan 5, 20012,344 citations

Generalized Pairing Wave Functions and Nodal Properties for Electronic Structure Quantum Monte Carlo

Michal Bajdich

The quantum Monte Carlo (QMC) is one of the most promising many-body electronic structure approaches. It employs stochastic techniques for solving the stationary Schr\" odinger equation and for evaluation of expectation values. The key advantage of QMC is its capability to use the explicitly correlated wave functions, which allow the study of many-body effects beyond the reach of mean-field methods. The most important limit on QMC accuracy is the fixed-node approximation, which comes from necessity to circumvent the fermion sign problem. The size of resulting fixed-node errors depends on the quality of the nodes (the subset of position space where the wave function vanishes) of a used wave function. In this dissertation, we analyze the nodal properties of the existing fermionic wave functions and offer new types of variational wave functions with improved nodal structure.

arXiv:0712.3066Dec 18, 20070 citations

Monte Carlo studies in weakly interacting quantum systems

Ryan Springall

This thesis is composed of two parts. In part one we report results of two quantum Monte Carlo methods - variational Monte Carlo and diffusion Monte Carlo - on the potential energy curve of the helium dimer, the prototypical van der Waals system. In contrast to previous quantum Monte Carlo calculations on this system, we have employed trial wave functions of the Slater-Jastrow form and used the fixed node approximation for the fermion nodal surface. We find both methods to be in excellent agreement with the best theoretical results at short range. In addition, the diffusion Monte Carlo results give very good agreement across the whole potential energy curve, while the Slater-Jastrow wave function fails to bind the dimer at all. In part two we switch to investigations of many-body systems at finite temperature. We use the path integral representation of statistical mechanics to investigate the symmetry properties of the canonical ensemble partition function and find a representation in terms of irreducible representations of the symmetric group. We use this as a foundation to propose a novel technique for the stochastic sampling of the bosonic partition function for quantum gases. It is shown that in principle we are able to use two operators which enable us to construct a Markov chain through a graph of the irreducible representation of the symmetric group. As an illustration of this method, a test calculation of four particles in a harmonic trap is performed.

Jun 27, 20120 citations