BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Dario Alfè

Dario Alfè appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.mtrl-sciphysics.chem-phphysics.comp-phcond-mat.othercond-mat.str-elphysics.geo-phcond-mat.mes-hallcond-mat.softquant-ph

Catalog footprint

What is connected

14works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

14 published item(s)

preprint2016arXiv

Boosting the accuracy and speed of quantum Monte Carlo: size-consistency and time-step

Diffusion Monte Carlo (DMC) simulations for fermions are becoming the standard to provide high quality reference data in systems that are too large to be investigated via quantum chemical approaches. DMC with the fixed-node approximation relies on modifications of the Green function to avoid singularities near the nodal surface of the trial wavefunction. We show that these modifications affect the DMC energies in a way that is not size-consistent, resulting in large time-step errors. Building on the modifications of Umrigar {\em et al.} and of DePasquale {\em et al.} we propose a simple Green function modification that restores size-consistency to large values of time-step; substantially reducing the time-step errors. The new algorithm also yields remarkable speedups of up to two orders of magnitude in the calculation of molecule-molecule binding energies and crystal cohesive energies, thus extending the horizons of what is possible with DMC.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Evidence for Stable Square Ice from Quantum Monte Carlo

Recent experiments on ice formed by water under nanoconfinement provide evidence for a two-dimensional (2D) `square ice' phase. However, the interpretation of the experiments has been questioned and the stability of square ice has become a matter of debate. Partially this is because the simulation approaches employed so far (force fields and density functional theory) struggle to accurately describe the very small energy differences between the relevant phases. Here we report a study of 2D ice using an accurate wave-function based electronic structure approach, namely Diffusion Monte Carlo (DMC). We find that at relatively high pressure square ice is indeed the lowest enthalpy phase examined, supporting the initial experimental claim. Moreover, at lower pressures a `pentagonal ice' phase (not yet observed experimentally) has the lowest enthalpy, and at ambient pressure the `pentagonal ice' phase is degenerate with a `hexagonal ice' phase. Our DMC results also allow us to evaluate the accuracy of various density functional theory exchange correlation functionals and force field models, and in doing so we extend the understanding of how such methodologies perform to challenging 2D structures presenting dangling hydrogen bonds.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Toward Accurate Adsorption Energetics on Clay Surfaces

Clay minerals are ubiquitous in nature, and the manner in which they interact with their surroundings has important industrial and environmental implications. Consequently, a molecular-level understanding of the adsorption of molecules on clay surfaces is crucial. In this regard computer simulations play an important role, yet the accuracy of widely used empirical force fields (FF) and density functional theory (DFT) exchange-correlation functionals is often unclear in adsorption systems dominated by weak interactions. Herein we present results from quantum Monte Carlo (QMC) for water and methanol adsorption on the prototypical clay kaolinite. To the best of our knowledge, this is the first time QMC has been used to investigate adsorption at a complex, natural surface such as a clay. As well as being valuable in their own right, the QMC benchmarks obtained provide reference data against which the performance of cheaper DFT methods can be tested. Indeed using various DFT exchange-correlation functionals yields a very broad range of adsorption energies, and it is unclear a priori which evaluation is better. QMC reveals that in the systems considered here it is essential to account for van der Waals (vdW) dispersion forces since this alters both the absolute and relative adsorption energies of water and methanol. We show, via FF simulations, that incorrect relative energies can lead to significant changes in the interfacial densities of water and methanol solutions at the kaolinite interface. Despite the clear improvements offered by the vdW-corrected and the vdW-inclusive functionals, absolute adsorption energies are often overestimated, suggesting that the treatment of vdW forces in DFT is not yet a solved problem.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Tuning dissociation using isoelectronically doped graphene and hexagonal boron nitride: water and other small molecules

Novel uses for 2-dimensional materials like graphene and hexagonal boron nitride (h-BN) are being frequently discovered especially for membrane and catalysis applications. Still however, a great deal remains to be understood about the interaction of environmentally and industrially elevant molecules such as water with these materials. Taking inspiration from advances in hybridising graphene and h-BN, we explore using density functional theory, the dissociation of water, hydrogen, methane, and methanol on graphene, h-BN, and their isoelectronic doped counterparts: BN doped graphene and C doped h-BN. We find that doped surfaces are considerably more reactive than their pristine counterparts and by comparing the reactivity of several small molecules we develop a general framework for dissociative adsorption. From this a particularly attractive consequence of isoelectronic doping emerges: substrates can be doped to enhance their reactivity specifically towards either polar or non-polar adsorbates. As such, these substrates are potentially viable candidates for selective catalysts and membranes, with the implication that a range of tuneable materials can be designed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Water on BN doped benzene: A hard test for exchange-correlation functionals and the impact of exact exchange on weak binding

Density functional theory (DFT) studies of weakly interacting complexes have recently focused on the importance of van der Waals dispersion forces whereas, the role of exchange has received far less attention. Here, by exploiting the subtle binding between water and a boron and nitrogen doped benzene derivative (1,2-azaborine) we show how exact exchange can alter the binding conformation within a complex. Benchmark values have been calculated for three orientations of the water monomer on 1,2-azaborine from explicitly correlated quantum chemical methods, and we have also used diffusion quantum Monte Carlo. For a host of popular DFT exchange-correlation functionals we show that the lack of exact exchange leads to the wrong lowest energy orientation of water on 1,2-azaborine. As such, we suggest that a high proportion of exact exchange and the associated improvement in the electronic structure could be needed for the accurate prediction of physisorption sites on doped surfaces and in complex organic molecules. Meanwhile to predict correct absolute interaction energies an accurate description of exchange needs to be augmented by dispersion inclusive functionals, and certain non-local van der Waals functionals (optB88- and optB86b-vdW) perform very well for absolute interaction energies. Through a comparison with water on benzene and borazine (B$_3$N$_3$H$_6$) we show that these results could have implications for the interaction of water with doped graphene surfaces, and suggest a possible way of tuning the interaction energy.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Water on hexagonal boron nitride from diffusion Monte Carlo

Despite a recent flurry of experimental and simulation studies, an accurate estimate of the interaction strength of water molecules with hexagonal boron nitride is lacking. Here we report quantum Monte Carlo results for the adsorption of a water monomer on a periodic hexagonal boron nitride sheet, which yield a water monomer interaction energy of -84 +/- 5 meV. We use the results to evaluate the performance of several widely used density functional theory (DFT) exchange correlation functionals, and find that they all deviate substantially. Differences in interaction energies between different adsorption sites are however better reproduced by DFT.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Benchmarking the performance of Density Functional Theory and Point Charge Force Fields in their Description of sI Methane Hydrate against Diffusion Monte Carlo

High quality reference data from diffusion Monte Carlo calculations are presented for bulk sI methane hydrate, a complex crystal exhibiting both hydrogen-bond and dispersion dominated interactions. The performance of some commonly used exchange-correlation functionals and all-atom point charge force fields is evaluated. Our results show that none of the exchange-correlation functionals tested are sufficient to describe both the energetics and the structure of methane hydrate accurately, whilst the point charge force fields perform badly in their description of the cohesive energy but fair well for the dissociation energetics. By comparing to ice Ih, we show that a good prediction of the volume and cohesive energies for the hydrate relies primarily on an accurate description of the hydrogen bonded water framework, but that to correctly predict stability of the hydrate with respect to dissociation to ice Ih and methane gas, accuracy in the water-methane interaction is also required. Our results highlight the difficulty that density functional theory faces in describing both the hydrogen bonded water framework and the dispersion bound methane.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Comparison of polynomial approximations to speed up planewave-based quantum Monte Carlo calculations

The computational cost of quantum Monte Carlo (QMC) calculations of realistic periodic systems depends strongly on the method of storing and evaluating the many-particle wave function. Previous work [A. J. Williamson et al., Phys. Rev. Lett. 87, 246406 (2001); D. Alfè and M. J. Gillan, Phys. Rev. B 70, 161101 (2004)] has demonstrated the reduction of the O(N^3) cost of evaluating the Slater determinant with planewaves to O(N^2) using localized basis functions. We compare four polynomial approximations as basis functions -- interpolating Lagrange polynomials, interpolating piecewise-polynomial-form (pp-) splines, and basis-form (B-) splines (interpolating and smoothing). All these basis functions provide a similar speedup relative to the planewave basis. The pp-splines have eight times the memory requirement of the other methods. To test the accuracy of the basis functions, we apply them to the ground state structures of Si, Al, and MgO. The polynomial approximations differ in accuracy most strongly for MgO and smoothing B-splines most closely reproduce the planewave value for of the variational Monte Carlo energy. Using separate approximations for the Laplacian of the orbitals increases the accuracy sufficiently to justify the increased memory requirement, making smoothing B-splines, with separate approximation for the Laplacian, the preferred choice for approximating planewave-represented orbitals in QMC calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the Accuracy of van der Waals Inclusive Density-Functional Theory Exchange-Correlation Functionals for Ice at Ambient and High Pressures

Density-functional theory (DFT) has been widely used to study water and ice for at least 20 years. However, the reliability of different DFT exchange-correlation (xc) functionals for water remains a matter of considerable debate. This is particularly true in light of the recent development of DFT based methods that account for van der Waals (vdW) dispersion forces. Here, we report a detailed study with several xc functionals (semi-local, hybrid, and vdW inclusive approaches) on ice Ih and six proton ordered phases of ice. Consistent with our previous study [Phys. Rev. Lett. 107, 185701 (2011)] which showed that vdW forces become increasingly important at high pressures, we find here that all vdW inclusive methods considered improve the relative energies and transition pressures of the high-pressure ice phases compared to those obtained with semi-local or hybrid xc functionals. However, we also find that significant discrepancies between experiment and the vdW inclusive approaches remain in the cohesive properties of the various phases, causing certain phases to be absent from the phase diagram. Therefore, room for improvement in the description of water at ambient and high pressures remains and we suggest that because of the stern test the high pressure ice phases pose they should be used in future benchmark studies of simulation methods for water.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Melting curve of face-centred-cubic nickel from first principles calculations

The melting curve of Ni up to 100 GPa has been calculated using first principles methods based on density functional theory (DFT). We used two complementary approaches: i) coexistence simulations with a reference system and then free energy corrections between DFT and the reference system, and ii) direct DFT coexistence using simulation cells including 1000 atoms. The calculated zero pressure melting temperature is slightly underestimated at $1637 \pm 10$ K (experimental value is 1728 K), and at high pressure is significantly higher than recent measurements in diamond anvil cell experiments [Phys. Rev. B {\bf 87}, 054108 (2013)]. The zero pressure DFT melting slope is calculated to be $30 \pm 2$ K, in good agreement with the experimental value of 28 K.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

High-temperature behaviour of supported graphene: electron-phonon coupling and substrate-induced doping

One of the salient features of graphene is the very high carrier mobility that implies tremendous potential for use in electronic devices. Unfortunately, transport measurements find the expected high mobility only in freely suspended graphen. When supported on a surface, graphene shows a strongly reduced mobility, and an especially severe reduction for temperatures above 200 K. A temperature-dependent mobility reduction could be explained by scattering of carriers with phonons, but this is expected to be weak for pristine, weakly-doped graphene. The mobility reduction has therefore been ascribed to the interaction with confined ripples or substrate phonons. Here we study the temperature-dependent electronic structure of supported graphene by angle-resolved photoemission spectroscopy, a technique that can reveal the origin of the phenomena observed in transport measurements. We show that the electron-phonon coupling for weakly-doped, supported graphene on a metal surface is indeed extremely weak, reaching the lowest value ever reported for any material. However, the temperature-dependent dynamic interaction with the substrate leads to a complex and dramatic change in the carrier type and density that is relevant for transport. Using ab initio molecular dynamics simulations, we show that these changes in the electronic structure are mainly caused by fluctuations in the graphene-substrate distance.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Thermal and electrical conductivity of iron at Earth's core conditions

The Earth acts as a gigantic heat engine driven by decay of radiogenic isotopes and slow cooling, which gives rise to plate tectonics, volcanoes, and mountain building. Another key product is the geomagnetic field, generated in the liquid iron core by a dynamo running on heat released by cooling and freezing to grow the solid inner core, and on chemical convection due to light elements expelled from the liquid on freezing. The power supplied to the geodynamo, measured by the heat-flux across the core-mantle boundary (CMB), places constraints on Earth's evolution. Estimates of CMB heat-flux depend on properties of iron mixtures under the extreme pressure and temperature conditions in the core, most critically on the thermal and electrical conductivities. These quantities remain poorly known because of inherent difficulties in experimentation and theory. Here we use density functional theory to compute these conductivities in liquid iron mixtures at core conditions from first principles- the first directly computed values that do not rely on estimates based on extrapolations. The mixtures of Fe, O, S, and Si are taken from earlier work and fit the seismologically-determined core density and inner-core boundary density jump. We find both conductivities to be 2-3 times higher than estimates in current use. The changes are so large that core thermal histories and power requirements must be reassessed. New estimates of adiabatic heat-flux give 15-16 TW at the CMB, higher than present estimates of CMB heat-flux based on mantle convection; the top of the core must be thermally stratified and any convection in the upper core driven by chemical convection against the adverse thermal buoyancy or lateral variations in CMB heat flow. Power for the geodynamo is greatly restricted and future models of mantle evolution must incorporate a high CMB heat-flux and explain recent formation of the inner core.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Transport properties for liquid silicon-oxygen-iron mixtures at Earth's core conditions

We report on the thermal and electrical conductivities of two liquid silicon-oxygen-iron mixtures (Fe$_{0.82}$Si$_{0.10}$O$_{0.08}$ and Fe$_{0.79}$Si$_{0.08}$O$_{0.13}$), representative of the composition of the Earth's outer core at the relevant pressure-temperature conditions, obtained from density functional theory calculations with the Kubo-Greenwood formulation. We find thermal conductivities $k$ =100 (160) W m$^{-1}$ K$^{-1}$, and electrical conductivities $σ= 1.1 (1.3) \times 10^6 Ω^{-1}$ m$^{-1}$ at the top (bottom) of the outer core. These new values are between 2 and 3 times higher than previous estimates, and have profound implications for our understanding of the Earth's thermal history and the functioning of the Earth's magnetic field, including rapid cooling rate for the whole core or high level of radiogenic elements in the core. We also show results for a number of structural and dynamic properties of the mixtures, including the partial radial distribution functions, mean square displacements, viscosities and speeds of sound.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Hydrogen bonds and van der Waals forces in ice at ambient and high pressures

The first principles approaches, density functional theory (DFT) and quantum Monte Carlo, have been used to examine the balance between van der Waals (vdW) forces and hydrogen (H) bonding in ambient and high pressure phases of ice. At higher pressure, the contribution to the lattice energy from vdW increases and that from H bonding decreases, leading vdW to have a substantial effect on the transition pressures between the crystalline ice phases. An important consequence, likely to be of relevance to molecular crystals in general, is that transition pressures obtained from DFT functionals which neglect vdW forces are greatly overestimated.

0 citations0 reviewsNo code linkedOpen access