Source author record

M. D. Towler

M. D. Towler appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

6works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

6 published item(s)

preprint2012arXiv

Assessing the accuracy of quantum Monte Carlo and density functional theory for energetics of small water clusters

We present a detailed study of the energetics of water clusters (H$_2$O)$_n$ with $n \le 6$, comparing diffusion Monte Carlo (DMC) and approximate density functional theory (DFT) with well converged coupled-cluster benchmarks. We use the many-body decomposition of the total energy to classify the errors of DMC and DFT into 1-body, 2-body and beyond-2-body components. Using both equilibrium cluster configurations and thermal ensembles of configurations, we find DMC to be uniformly much more accurate than DFT, partly because some of the approximate functionals give poor 1-body distortion energies. Even when these are corrected, DFT remains considerably less accurate than DMC. When both 1- and 2-body errors of DFT are corrected, some functionals compete in accuracy with DMC; however, other functionals remain worse, showing that they suffer from significant beyond-2-body errors. Combining the evidence presented here with the recently demonstrated high accuracy of DMC for ice structures, we suggest how DMC can now be used to provide benchmarks for larger clusters and for bulk liquid water.

preprint2011arXiv

Timescales for dynamical relaxation to the Born rule

We illustrate through explicit numerical calculations how the Born-rule probability densities of non-relativistic quantum mechanics emerge naturally from the particle dynamics of de Broglie-Bohm pilot-wave theory. The time evolution of a particle distribution initially not equal to the absolute square of the wave function is calculated for a particle in a two-dimensional infinite potential square well. Under the de Broglie-Bohm ontology, the box contains an objectively-existing 'pilot wave' which guides the electron trajectory, and this is represented mathematically by a Schroedinger wave function composed of a finite out-of-phase superposition of M energy eigenstates (with M ranging from 4 to 64). The electron density distributions are found to evolve naturally into the Born-rule ones and stay there; in analogy with the classical case this represents a decay to 'quantum equilibrium'. The proximity to equilibrium is characterized by the coarse-grained subquantum H-function which is found to decrease roughly exponentially towards zero over the course of time. The timescale tau for this relaxation is calculated for various values of M and the coarse-graining length epsilon. Its dependence on M is found to disagree with an earlier theoretical prediction. A power law - tau inversely proportional to M - is found to be fairly robust for all coarse-graining lengths and, although a weak dependence of tau on epsilon is observed, it does not appear to follow any straightforward scaling. A theoretical analysis is presented to explain these results. This improvement in our understanding of timescales for relaxation to quantum equilibrium is likely to be of use in the development of models of relaxation in the early universe, with a view to constraining possible violations of the Born rule in inflationary cosmology.

preprint2010arXiv

Continuum variational and diffusion quantum Monte Carlo calculations

This topical review describes the methodology of continuum variational and diffusion quantum Monte Carlo calculations. These stochastic methods are based on many-body wave functions and are capable of achieving very high accuracy. The algorithms are intrinsically parallel and well-suited to petascale computers, and the computational cost scales as a polynomial of the number of particles. A guide to the systems and topics which have been investigated using these methods is given. The bulk of the article is devoted to an overview of the basic quantum Monte Carlo methods, the forms and optimisation of wave functions, performing calculations within periodic boundary conditions, using pseudopotentials, excited-state calculations, sources of calculational inaccuracy, and calculating energy differences and forces.

preprint2010arXiv

Fundamental high pressure calibration from all-electron quantum Monte Carlo calculations

We develop an all-electron quantum Monte Carlo (QMC) method for solids that does not rely on pseudopotentials, and use it to construct a primary ultra-high pressure calibration based the equation of state of cubic boron nitride(c-BN). We compute the static contribution to the free energy with QMC, and obtain the phonon contribution from density functional theory, yielding a high-accuracy calibration up to 900 GPa usable directly in experiment. Furthermore, we compute the anharmonic Raman frequency shift with QMC as a function of pressure and temperature, allowing optical pressure calibration in table-top experiments. In contrast to present experimental approaches, small systematic errors in the theoretical EOS do not increase with pressure, and no extrapolation is needed. This all-electron methodology is generally applicable to first-row solids, and can be used to provide a new reference for ab initio calculations of solids and to benchmark pseudopotential accuracy.

preprint2010arXiv

Quantum Monte Carlo for minerals at high pressure: Phase stability, equations of state, and elasticity of silica

Silica is an abundant component of the Earth whose crystalline polymorphs play key roles in its structure and dynamics. As the simplest silicates, understanding pure silica is a prerequisite to understanding the rocky part of the Earth, its majority. First principle density functional theory (DFT) methods have often been used to accurately predict properties of silicates. Here, we study silica with quantum Monte Carlo (QMC), which until now was not computationally possible for such complex materials, and find that QMC overcomes the failures of DFT. QMC is a benchmark method that does not rely on density functionals, but rather explicitly treats the electrons and their interactions via a stochastic solution of Schrodinger's equation. Using ground state QMC plus phonons within the quasiharmonic approximation from density functional perturbation theory, we obtain the thermal pressure and equations of state of silica phases up to Earth's core-mantle boundary. Our results provide the most well-constrained equations of state and phase boundaries available for silica. QMC indicates a transition to the most dense alpha-PbO2 structure above the core-insulating D" layer, suggesting the absence of significant free silica in the bulk lower mantle, which has been assumed but never proven. We also find an accurate shear elastic constant and its geophysically important softening with pressure.

preprint2004arXiv

Diamond and $β$-tin structures of Si studied with quantum Monte Carlo calculations

We have used diffusion quantum Monte Carlo (DMC) calculations to study the pressure-induced phase transition from the diamond to $β$-tin structure in silicon. The calculations employ the pseudopotential technique and systematically improvable B-spline basis sets. We show that in order to achieve a precision of 1 GPa in the transition pressure the non-cancelling errors in the energies of the two structures must be reduced to 30 meV/atom. Extensive tests on system size errors, non-local pseudopotential errors, basis-set incompleteness errors, and other sources of error, performed on periodically repeated systems of up to 432 atoms, show that all these errors together can be reduced to well below 30 meV/atom. The calculated DMC transition pressure is about 3-4 GPa higher than the accepted experimental range of values, and we argue that the discrepancy may be due to the fixed-node error inherent in DMC techniques.