BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

M. J. Gillan

M. J. Gillan contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.mtrl-sciphysics.chem-phcond-mat.stat-mechBiological PhysicsBiomoleculescond-mat.othercond-mat.softmtrl-thphysics.comp-phphysics.geo-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
17works
0followers
10topics
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

17 published item(s)

preprint2016arXiv

Perspective: How good is DFT for water?

Kohn-Sham density functional theory (DFT) has become established as an indispensable tool for investigating aqueous systems of all kinds, including those important in chemistry, surface science, biology and the earth sciences. Nevertheless, many widely used approximations for the exchange-correlation (XC) functional describe the properties of pure water systems with an accuracy that is not fully satisfactory. The explicit inclusion of dispersion interactions generally improves the description, but there remain large disagreements between the predictions of different dispersion-inclusive methods. We present here a review of DFT work on water clusters, ice structures and liquid water, with the aim of elucidating how the strengths and weaknesses of different XC approximations manifest themselves across this variety of water systems. Our review highlights the crucial role of dispersion in describing the delicate balance between compact and extended structures of many different water systems, including the liquid. By referring to a wide range of published work, we argue that the correct description of exchange-overlap interactions is also extremely important, so that the choice of semi-local or hybrid functional employed in dispersion-inclusive methods is crucial. The origins and consequences of beyond-2-body errors of approximate XC functionals are noted, and we also discuss the substantial differences between different representations of dispersion. We propose a simple numerical scoring system that rates the performance of different XC functionals in describing water systems, and we suggest possible future developments.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Density-functional theory study of gramicidin A ion channel geometry and electronic properties

Understanding the mechanisms underlying ion channel function from the atomic-scale requires accurate ab initio modelling as well as careful experiments. Here, we present a density functional theory (DFT) study of the ion channel gramicidin A, whose inner pore conducts only monovalent cations and whose conductance has been shown to depend on the side chains of the amino acids in the channel. We investigate the ground-state geometry and electronic properties of the channel in vacuum, focusing on their dependence on the side chains of the amino acids. We find that the side chains affect the ground state geometry, while the electrostatic potential of the pore is independent of the side chains. This study is also in preparation for a full, linear scaling DFT study of gramicidin A in a lipid bilayer with surrounding water. We demonstrate that linear scaling DFT methods can accurately model the system with reasonable computational cost. Linear scaling DFT allows ab initio calculations with 10,000 to 100,000 atoms and beyond, and will be an important new tool for biomolecular simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Energy benchmarks for water clusters and ice structures from an embedded many-body expansion

We show how an embedded many-body expansion (EMBE) can be used to calculate accurate \emph{ab initio} energies of water clusters and ice structures using wavefunction-based methods. We use the EMBE described recently by Bygrave \emph{et al.} (J. Chem. Phys. \textbf{137}, 164102 (2012)), in which the terms in the expansion are obtained from calculations on monomers, dimers, etc. acted on by an approximate representation of the embedding field due to all other molecules in the system, this field being a sum of Coulomb and exchange-repulsion fields. Our strategy is to separate the total energy of the system into Hartree-Fock and correlation parts, using the EMBE only for the correlation energy, with the Hartree-Fock energy calculated using standard molecular quantum chemistry for clusters and plane-wave methods for crystals. Our tests on a range of different water clusters up to the 16-mer show that for the second-order Møller-Plesset (MP2) method the EMBE truncated at 2-body level reproduces to better than 0.1 m$E_{\rm h}$/monomer the correlation energy from standard methods. The use of EMBE for computing coupled-cluster energies of clusters is also discussed. For the ice structures Ih, II and VIII, we find that MP2 energies near the complete basis-set limit reproduce very well the experimental values of the absolute and relative binding energies, but that the use of coupled-cluster methods for many-body correlation (non-additive dispersion) is essential for a full description. Possible future applications of the EMBE approach are suggested.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

First-principles energetics of water: a many-body analysis

Standard forms of density-functional theory (DFT) have good predictive power for many materials, but are not yet fully satisfactory for solid, liquid and cluster forms of water. We use a many-body separation of the total energy into its 1-body, 2-body (2B) and beyond-2-body (B2B) components to analyze the deficiencies of two popular DFT approximations. We show how machine-learning methods make this analysis possible for ice structures as well as for water clusters. We find that the crucial energy balance between compact and extended geometries can be distorted by 2B and B2B errors, and that both types of first-principles error are important.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The kinetics of homogeneous melting beyond the limit of superheating

Molecular dynamics simulation is used to study the time-scales involved in the homogeneous melting of a superheated crystal. The interaction model used is an embedded-atom model for Fe developed in previous work, and the melting process is simulated in the microcanonical $(N, V, E)$ ensemble. We study periodically repeated systems containing from 96 to 7776 atoms, and the initial system is always the perfect crystal without free surfaces or other defects. For each chosen total energy $E$ and number of atoms $N$, we perform several hundred statistically independent simulations, with each simulation lasting for between 500 ps and 10 ns, in order to gather statistics for the waiting time $τ_{\rm w}$ before melting occurs. We find that the probability distribution of $τ_{\rm w}$ is roughly exponential, and that the mean value $<τ_{\rm w} >$ depends strongly on the excess of the initial steady temperature of the crystal above the superheating limit identified by other researchers. The mean $<τ_{\rm w}>$ also depends strongly on system size in a way that we have quantified. For very small systems of $\sim 100$ atoms, we observe a persistent alternation between the solid and liquid states, and we explain why this happens. Our results allow us to draw conclusions about the reliability of the recently proposed Z method for determining the melting properties of simulated materials, and to suggest ways of correcting for the errors of the method.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Ab initio statistical mechanics of surface adsorption and desorption: II. Nuclear quantum effects

We show how the path-integral formulation of quantum statistical mechanics can be used to construct practical {\em ab initio} techniques for computing the chemical potential of molecules adsorbed on surfaces, with full inclusion of quantum nuclear effects. The techniques we describe are based on the computation of the potential of mean force on a chosen molecule, and generalise the techniques developed recently for classical nuclei. We present practical calculations based on density functional theory with a generalised-gradient exchange-correlation functional for the case of H$_2$O on the MgO~(001) surface at low coverage. We note that the very high vibrational frequencies of the H$_2$O molecule would normally require very large numbers of time slices (beads) in path-integral calculations, but we show that this requirement can be dramatically reduced by employing the idea of thermodynamic integration with respect to the number of beads. The validity and correctness of our path-integral calculations on the H$_2$O/MgO~(001) system are demonstrated by supporting calculations on a set of simple model systems for which quantum contributions to the free energy are known exactly from analytic arguments.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Bulk and surface energetics of lithium hydride crystal: benchmarks from quantum Monte Carlo and quantum chemistry

We show how accurate benchmark values of the surface formation energy of crystalline lithium hydride can be computed by the complementary techniques of quantum Monte Carlo (QMC) and wavefunction-based molecular quantum chemistry. To demonstrate the high accuracy of the QMC techniques, we present a detailed study of the energetics of the bulk LiH crystal, using both pseudopotential and all-electron approaches. We show that the equilibrium lattice parameter agrees with experiment to within 0.03 %, which is around the experimental uncertainty, and the cohesive energy agrees to within around 10 meV per formula unit. QMC in periodic slab geometry is used to compute the formation energy of the LiH (001) surface, and we show that the value can be accurately converged with respect to slab thickness and other technical parameters. The quantum chemistry calculations build on the recently developed hierarchical scheme for computing the correlation energy of a crystal to high precision. We show that the hierarchical scheme allows the accurate calculation of the surface formation energy, and we present results that are well converged with respect to basis set and with respect to the level of correlation treatment. The QMC and hierarchical results for the surface formation energy agree to within about 1 %.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Melting properties of a simple tight-binding model of transition metals: I.The region of half-filled d-band

We present calculations of the free energy, and hence the melting properties, of a simple tight-binding model for transition metals in the region of d-band filling near the middle of a d-series, the parameters of the model being designed to mimic molybdenum. The melting properties are calculated for pressures ranging from ambient to several Mbar. The model is intended to be the simplest possible tight-binding representation of the two basic parts of the energy: first, the pairwise repulsion due to Fermi exclusion; and second, the d-band bonding energy described in terms of an electronic density of states that depends on structure. In addition to the number of d-electrons, the model contains four parameters, which are adjusted to fit the pressure dependent d-band width and the zero-temperature pressure-volume relation of Mo. We show that the resulting model reproduces well the phonon dispersion relations of Mo in the body-centred-cubic structure, as well as the radial distribution function of the high-temperature solid and liquid given by earlier first-principles simulations. Our free-energy calculations start from the free energy of the liquid and solid phases of the purely repulsive pair-potential model, without d-band bonding. The free energy of the full tight-binding model is obtained from this by thermodynamic integration. The resulting melting properties of the model are quite close to those given by earlier first-principles work on Mo. An interpretation of these melting properties is provided by showing how they are related to those of the purely repulsive model.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

The Schottky defect formation energy in MgO calculated by diffusion Monte Carlo

The energetics of point defects in oxide materials plays a major role in determining their high-temperature properties, but experimental measurements are difficult, and calculations based on density functional theory (DFT) are not necessarily reliable. We report quantum Monte Carlo (QMC) calculations of the formation energy $E_{\rm S}$ of Schottky defects in MgO, which demonstrate the feasibility of using this approach to overcome the deficiencies of DFT. In order to investigate system-size errors, we also report DFT calculations of $E_{\rm S}$ on repeating cells of up to $\sim 1000$ atoms, which indicate that QMC calculations on systems of only 54 atoms should yield high precision. The DFT calculations also provide the relaxed structures used in the variational and diffusion Monte Carlo calculations. For MgO, we find $E_{\rm S}$ to be in close agreement with results from DFT and from model interaction potentials, and consistent with the scattered experimental values. The prospects for applying the same approach to transition metal oxides such as FeO are indicated.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

An efficient localized basis set for quantum Monte Carlo calculations on condensed matter

We present an efficient scheme for representing many-body wavefunctions in quantum Monte Carlo (QMC) calculations. The scheme is based on B-splines (blip functions), which consist of localized cubic splines centred on the points of a regular grid. We show that blip functions are unbiased, systematically improvable, and conveniently obtained from any standard plane-waves density functional theory (PW-DFT) code, and therefore provide a convenient and natural interface between PW-DFT and QMC calculations. We present tests on a 16-atom system of Si in the $β$-tin structure, and on 2- and 8- atoms systems of MgO in the NaCl structure. We show that already with such small systems the speed-up of blip functions with respect to plane-waves is between one and two order of magnitudes, without compromising the accuracy.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Linear-scaling quantum Monte Carlo with non-orthogonal localized orbitals

We have reformulated the quantum Monte Carlo (QMC) technique so that a large part of the calculation scales linearly with the number of atoms. The reformulation is related to a recent alternative proposal for achieving linear-scaling QMC, based on maximally localized Wannier orbitals (MLWO), but has the advantage of greater simplicity. The technique we propose draws on methods recently developed for linear-scaling density functional theory. We report tests of the new technique on the insulator MgO, and show that its linear-scaling performance is somewhat better than that achieved by the MLWO approach. Implications for the application of QMC to large complex systems are pointed out.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

The particle-in-cell model for ab initio thermodynamics: implications for the elastic anisotropy of the Earth&#39;s inner core

We assess the quantitative accuracy of the particle-in-cell (PIC) approximation used in recent ab initio predictions of the thermodynamic properties of hexagonal-close-packed iron at the conditions of the Earth&#39;s inner core. The assessment is made by comparing PIC predictions for a range of thermodynamic properties with the results of more exact calculations that avoid the PIC approximation. It is shown that PIC gives very accurate results for some properties, but that it gives an incorrect treatment of anharmonic lattice vibrations. In addition, our assessment does not support recent PIC-based predictions that the hexagonal c/a ratio increases strongly with increasing temperature, and we point out that this casts doubt on a proposed re-interpretation of the elastic anisotropy of the inner core.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Electron correlation and the phase diagram of Si

Previous first-principles calculations of the melting properties of Si, based on the local-density approximation (LDA) for electronic exchange-correlation energy, under-predict the melting temperature by ~ 20%. We present new first-principles results demonstrating that this problem is due to non-cancellation of exchange-correlation errors between the semiconducting solid and the metallic liquid. It is shown that other sources of error, particularly those due to system size and Brillouin-zone sampling, can be made negligible. The same LDA errors cause an underprediction of the pressure of the diamond-Si --> beta-tin-Si transition. The generalized-gradient approximation largely corrects both features of the Si phase diagram.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Thermodynamic stability of Fe/O solid solution at inner-core conditions

We present a new technique which allows the fully {\em ab initio} calculation of the chemical potential of a substitutional impurity in a high-temperature crystal, including harmonic and anharmonic lattice vibrations. The technique uses the combination of thermodynamic integration and reference models developed recently for the {\em ab initio} calculation of the free energy of liquids and anharmonic solids. We apply the technique to the case of the substitutional oxygen impurity in h.c.p. iron under Earth&#39;s core conditions, which earlier static {\em ab initio} calculations indicated to be thermodynamically very unstable. Our results show that entropic effects arising from the large vibrational amplitude of the oxygen impurity give a major reduction of the oxygen chemical potential, so that oxygen dissolved in h.c.p. iron may be stabilised at concentrations up a few mol % under core conditions.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

A Self-Consistent First-Principles Technique Having Linear Scaling

An algorithm for first-principles electronic structure calculations having a computational cost which scales linearly with the system size is presented. Our method exploits the real-space localization of the density matrix, and in this respect it is related to the technique of Li, Nunes and Vanderbilt. The density matrix is expressed in terms of localized support functions, and a matrix of variational parameters, L, having a finite spatial range. The total energy is minimized with respect to both the support functions and the elements of the L matrix. The method is variational, and becomes exact as the ranges of the support functions and the L matrix are increased. We have tested the method on crystalline silicon systems containing up to 216 atoms, and we discuss some of these results.

0 citations0 reviewsNo code linkedOpen access