Source author record

J. R. Trail

J. R. Trail 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

7works
3topics
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

7 published item(s)

preprint2016arXiv

Trail-Needs pseudopotentials in quantum Monte Carlo calculations with plane-wave/blip basis sets

We report a systematic analysis of the performance of a widely used set of Dirac-Fock pseudopotentials for quantum Monte Carlo (QMC) calculations. We study each atom in the periodic table from hydrogen (Z=1) to mercury (Z=80), with the exception of the 4f elements (57 <= Z <= 70). We demonstrate that ghost states are a potentially serious problem when plane-wave basis sets are used in density functional theory (DFT) orbital-generation calculations, but that this problem can be almost entirely eliminated by choosing the s channel to be local in the DFT calculation; the d channel can then be chosen to be local in subsequent QMC calculations, which generally leads to more accurate results. We investigate the achievable energy variance per electron with different levels of trial wave function and we determine appropriate plane-wave cutoff energies for DFT calculations for each pseudopotential. We demonstrate that the so-called "T-move" scheme in diffusion Monte Carlo is essential for many elements. We investigate the optimal choice of spherical integration rule for pseudopotential projectors in QMC calculations. The information reported here will prove crucial in the planning and execution of QMC projects involving beyond-first-row elements.

preprint2013arXiv

Pseudopotentials for correlated electron systems

A scheme is developed for creating pseudopotentials for use in correlated-electron calculations. Pseudopotentials for the light elements H, Li, Be, B, C, N, O, and F, are reported, based on data from high-level quantum chemical calculations. Results obtained with these correlated electron pseudopotentials (CEPPs) are compared with data for atomic energy levels and the dissociation energies, molecular geometries and zero-point vibrational energies of small molecules obtained from coupled cluster single double triple (CCSD(T)) calculations with large basis sets. The CEPPs give better results in correlated-electron calculations than Hartree-Fock-based pseudopotentials available in the literature.

preprint2010arXiv

Methods for calculating forces within quantum Monte Carlo simulations

Atomic force calculations within the variational and diffusion quantum Monte Carlo (VMC and DMC) methods are described. The advantages of calculating DMC forces with the "pure" rather than the "mixed" probability distribution are discussed. An accurate and practical method for calculating forces using the pure distribution is presented and tested for the SiH molecule. The statistics of force estimators are explored and violations of the Central Limit Theorem are found in some cases.

preprint2009arXiv

Accurate structure factors from pseudopotential methods

Highly accurate experimental structure factors of silicon are available in the literature, and these provide the ideal test for any \emph{ab initio} method for the construction of the all-electron charge density. In a recent paper [J. R. Trail and D. M. Bird, Phys. Rev. B {\bf 60}, 7863 (1999)] a method has been developed for obtaining an accurate all-electron charge density from a first principles pseudopotential calculation by reconstructing the core region of an atom of choice. Here this method is applied to bulk silicon, and structure factors are derived and compared with experimental and Full-potential Linear Augmented Plane Wave results (FLAPW). We also compare with the result of assuming the core region is spherically symmetric, and with the result of constructing a charge density from the pseudo-valence density + frozen core electrons. Neither of these approximations provide accurate charge densities. The aspherical reconstruction is found to be as accurate as FLAPW results, and reproduces the residual error between the FLAPW and experimental results.

preprint2009arXiv

Core reconstruction in pseudopotential calculations

A new method is presented for obtaining all-electron results from a pseudopotential calculation. This is achieved by carrying out a localised calculation in the region of an atomic nucleus using the embedding potential method of Inglesfield [J.Phys. C {\bf 14}, 3795 (1981)]. In this method the core region is \emph{reconstructed}, and none of the simplifying approximations (such as spherical symmetry of the charge density/potential or frozen core electrons) that previous solutions to this problem have required are made. The embedding method requires an accurate real space Green function, and an analysis of the errors introduced in constructing this from a set of numerical eigenstates is given. Results are presented for an all-electron reconstruction of bulk aluminium, for both the charge density and the density of states.

preprint2009arXiv

Density-functional embedding using a plane-wave basis

The constrained electron density method of embedding a Kohn-Sham system in a substrate system (first described by P. Cortona, Phys. Rev. B {\bf 44}, 8454 (1991) and T.A. Wesolowski and A. Warshel, J. Phys. Chem {\bf 97}, 8050 (1993)) is applied with a plane-wave basis and both local and non-local pseudopotentials. This method divides the electron density of the system into substrate and embedded electron densities, the sum of which is the electron density of the system of interest. Coupling between the substrate and embedded systems is achieved via approximate kinetic energy functionals. Bulk aluminium is examined as a test case for which there is a strong interaction between the substrate and embedded systems. A number of approximations to the kinetic-energy functional, both semi-local and non-local, are investigated. It is found that Kohn-Sham results can be well reproduced using a non-local kinetic energy functional, with the total energy accurate to better than 0.1 eV per atom and good agreement between the electron densities.

preprint2009arXiv

Electronic damping of molecular motion at metal surfaces

A method for the calculation of the damping rate due to electron-hole pair excitation for atomic and molecular motion at metal surfaces is presented. The theoretical basis is provided by Time Dependent Density Functional Theory (TDDFT) in the quasi-static limit and calculations are performed within a standard plane-wave, pseudopotential framework. The artificial periodicity introduced by using a super-cell geometry is removed to derive results for the motion of an isolated atom or molecule, rather than for the coherent motion of an ordered over-layer. The algorithm is implemented in parallel, distributed across both ${\bf k}$ and ${\bf g}$ space, and in a form compatible with the CASTEP code. Test results for the damping of the motion of hydrogen atoms above the Cu(111) surface are presented.