Researcher profile

T. Schoof

T. Schoof contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

8 published item(s)

preprint2016arXiv

Ab Initio Quantum Monte Carlo Simulations of the Uniform Electron Gas without Fixed Nodes II: Unpolarized Case

In a recent publication [S. Groth \textit{et al.}, PRB (2016)], we have shown that the combination of two novel complementary quantum Monte Carlo approaches, namely configuration path integral Monte Carlo (CPIMC) [T. Schoof \textit{et al.}, PRL \textbf{115}, 130402 (2015)] and permutation blocking path integral Monte Carlo (PB-PIMC) [T. Dornheim \textit{et al.}, NJP \textbf{17}, 073017 (2015)], allows for the accurate computation of thermodynamic properties of the spin-polarized uniform electron gas (UEG) over a wide range of temperatures and densities without the fixed-node approximation. In the present work, we extend this concept to the unpolarized case, which requires non-trivial enhancements that we describe in detail. We compare our new simulation results with recent restricted path integral Monte Carlo data [E. Brown \textit{et al}., PRL \textbf{110}, 146405 (2013)] for different energy contributions and pair distribution functions and find, for the exchange correlation energy, overall better agreement than for the spin-polarized case, while the separate kinetic and potential contributions substantially deviate.

preprint2015arXiv

Ab Initio Quantum Monte Carlo Simulations of the Uniform Electron Gas without Fixed Nodes

The uniform electron gas (UEG) at finite temperature is of key relevance for many applications in the warm dense matter regime, e.g. dense plasmas and laser excited solids. Also, the quality of density functional theory calculations crucially relies on the availability of accurate data for the exchange-correlation energy. Recently, new benchmark results for the N = 33 spin-polarized electrons at high density, r_s = r/a_B <= 4 and low temperature, have been obtained with the configuration path integral Monte Carlo (CPIMC) method [T. Schoof et al., Phys. Rev. Lett. 115, 130402 (2015)]. To achieve these results, the original CPIMC algorithm [T. Schoof et al., Contrib. Plasma Phys. 51, 687 (2011)] had to be further optimized to cope with the fermion sign problem (FSP). It is the purpose of this paper to give detailed information on the manifestation of the FSP in CPIMC simulations of the UEG and to demonstrate how it can be turned into a controllable convergence problem. In addition, we present new thermodynamic results for higher temperatures. Finally, to overcome the limitations of CPIMC towards strong coupling, we invoke an independent method|the recently developed permutation blocking path integral Monte Carlo approach [T. Dornheim et al., accepted for publication in J. Chem Phys., arXiv:1508.03221]. The combination of both approaches is able to yield ab initio data for the UEG over the entire density range, above a temperature of about one half of the Fermi temperature. Comparison with restricted path integral Monte Carlo data [E. W. Brown et al., Phys. Rev. Lett. 110, 146405 (2013)] allows us to quantify the systematic error arising from the free particle nodes.

preprint2015arXiv

Ab initio thermodynamic results for the degenerate electron gas at finite temperature

The uniform electron gas (UEG) at finite temperature is of key relevance for many applications in dense plasmas, warm dense matter, laser excited solids and much more. Accurate thermodynamic data for the UEG are an essential ingredient for many-body theories, in particular, density functional theory. Recently, first-principle restricted path integral Monte Carlo results became available which, however, due to the fermion sign problem, had to be restricted to moderate degeneracy, i.e. low to moderate densities with $r_s={\bar r}/a_B \gtrsim 1$. Here we present novel first-principle configuration PIMC results for electrons for $r_s \leq 1$. We also present quantum statistical data within the $e^4$-approximation that are in good agreement with the simulations at small to moderate $r_s$.

preprint2014arXiv

Introduction to Configuration Path Integral Monte Carlo

In low-temperature high-density plasmas quantum effects of the electrons are becoming increasingly important. This requires the development of new theoretical and computational tools. Quantum Monte Carlo methods are among the most successful approaches to first-principle simulations of many-body quantum systems. In this chapter we present a recently developed method---the configuration path integral Monte Carlo (CPIMC) method for moderately coupled, highly degenerate fermions at finite temperatures. It is based on the second quantization representation of the $N$-particle density operator in a basis of (anti-)symmetrized $N$-particle states (configurations of occupation numbers) and allows to tread arbitrary pair interactions in a continuous space. We give a detailed description of the method and discuss the application to electrons or, more generally, Coulomb-interacting fermions. As a test case we consider a few quantum particles in a one-dimensional harmonic trap. Depending on the coupling parameter (ratio of the interaction energy to kinetic energy), the method strongly reduces the sign problem as compared to direct path integral Monte Carlo (DPIMC) simulations in the regime of strong degeneracy which is of particular importance for dense matter in laser plasmas or compact stars. In order to provide a self-contained introduction, the chapter includes a short introduction to Metropolis Monte Carlo methods and the second quantization of quantum mechanics.

preprint2014arXiv

Towards ab initio thermodynamics of the electron gas at strong degeneracy

Recently a number of theoretical studies of the uniform electron gas (UEG) at finite temperature have appeared that are of relevance for dense plasmas, warm dense matter and laser excited solids and thermodynamic density functional theory simulations. In particular, restricted path integral Monte Carlo (RPIMC) results became available which, however, due to the Fermion sign problem, are confined to moderate quantum degeneracy, i.e. low to moderate densities. We have recently developed an alternative approach---configuration PIMC [T. Schoof {\em et al.}, Contrib. Plasma Phys. {\bf 51}, 687 (2011)] that allows one to study the so far not accessible high degeneracy regime. Here we present the first step towards UEG simulations using CPIMC by studying implementation and performance of the method for the model case of $N=4$ particles. We also provide benchmark data for the total energy.

preprint2013arXiv

Comment on &#34;Discussion on `Novel attractive force between ions in quantum plasmas -- failure of simulations based on a density functional approach&#34;

In a recent article [P.K. Shukla, B. Eliasson and M. Akbari-Moghanjoughi, Physica Scripta {\bf 87}, 018202 (2013)] the authors criticized our analysis of the screened proton potential in dense hydrogen that was based on {\em ab initio} density functional theory (DFT) simulations [M. Bonitz, E. Pehlke, and T. Schoof, Phys. Rev. E {\bf 87}, 037102 (2013)]. In particular, they attributed the absence of the Shukla-Eliasson attractive force between protons in the DFT simulations to a failure of DFT. Here we discuss in detail their arguments and show that their conclusions are incorrect.

preprint2013arXiv

On &#34;Novel attractive forces&#34; between ions in quantum plasmas -- failure of linearized quantum hydrodynamics

In a recent letter [P.K. Shukla and B. Eliasson, Phys. Rev. Lett. 108, 165007 (2012)] the discovery of a new attractive force between protons in a hydrogen plasma was reported that would be responsible for the formation of molecules and of a proton lattice. Here we show, based on ab initio density functional theory calculations, that these predictions are wrong and caused by using linearized quantum hydrodynamics beyond the limits of its applicability.

preprint2013arXiv

Reply to comment of Shukla et al. on &#39;On &#34;Novel attractive forces&#34; between ions in quantum plasmas - failure of linearized quantum hydrodynamics&#39;

In an earlier paper we demonstrated that the &#34;novel attractive force&#34; between protons in dense hydrogen reported by Shukla and Eliasson based on linearized quantum hydrodynamics (LQHD) is wrong. As benchmark results we used state of the art density functional theory (DFT) simulations. In their response to that paper Shukla et al. (arXiv:1112.5556) claimed that the disagreement is due to a failure of DFT. Here we show - by comparing the properties and limitations of DFT and QHD - that their statements have no scientific basis and that DFT is a suitable benchmark for more approximate approaches such as QHD and LQHD.