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

D. Alfè

D. Alfè contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.mtrl-sciphysics.chem-phcond-mat.str-elcond-mat.softcond-mat.stat-mechphysics.geo-ph

Trust snapshot

Quick read

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

13 published item(s)

preprint2020arXiv

Electronic correlations and transport in iron at Earth's core conditions

The transport properties of iron under Earth's inner core conditions are essential input for the geophysical modelling but are poorly constrained experimentally. Here we show that the thermal and electrical conductivity of iron at those conditions remains high even if the electron-electron-scattering (EES) is properly taken into account. This result is obtained by ab initio simulations taking into account consistently both thermal disorder and electronic correlations. Thermal disorder suppresses the non-Fermi-liquid behavior of the body-centered cubic iron phase, hence, reducing the EES; the total calculated thermal conductivity of this phase is 220 Wm$^{-1}$K$^{-1}$ with the EES reduction not exceeding 20%. The EES and electron-lattice scattering are intertwined resulting in breaking of the Matthiessen's rule with increasing EES. In the hexagonal close-packed iron the EES is also not increased by thermal disorder and remains weak. Our main finding thus holds for the both likely iron phases in the inner core.

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

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
preprint2012arXiv

Assessment of density functional theory for iron(II) molecules across the spin-crossover transition

Octahedral Fe$^{2+}$ molecules are particularly interesting as they often exhibit a spin-crossover transition. In spite of the many efforts aimed at assessing the performances of density functional theory for such systems, an exchange-correlation functional able to account accurately for the energetic of the various possible spin-states has not been identified yet. Here we critically discuss the issues related to the theoretical description of this class of molecules from first principles. In particular we present a comparison between different density functionals for four ions, namely [Fe(H$_2$O)$_6$]$^{2+}$, [Fe(NH$_3$)$_6$]$^{2+}$, [Fe(NCH)$_6$]$^{2+}$ and [Fe(CO)$_6$]$^{2+}$. These are characterized by different ligand-field splittings and ground state spin multiplicities. Since no experimental data are available for the gas phase, the density functional theory results are benchmarked against those obtained with diffusion Monte Carlo, one of the most accurate methods available to compute ground state total energies of quantum systems. On the one hand, we show that most of the functionals considered provide a good description of the geometry and of the shape of the potential energy surfaces. On the other hand, the same functionals fail badly in predicting the energy differences between the various spin states. In the case of [Fe(H$_2$O)$_6$]$^{2+}$, [Fe(NH$_3$)$_6$]$^{2+}$, [Fe(NCH)$_6$]$^{2+}$, this failure is related to the drastic underestimation of the exchange energy. Therefore quite accurate results can be achieved with hybrid functionals including about 50% of Hartree-Fock exchange. In contrast, in the case of [Fe(CO)$_6$]$^{2+}$, the failure is likely to be caused by the multiconfigurational character of the ground state wave-function and no suitable exchange and correlation functional has been identified.

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

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

Quantum Monte Carlo calculations of the structural properties and the B1-B2 phase transition of MgO

We report diffusion Monte Carlo (DMC) calculations on MgO in the rock-salt and CsCl structures. The calculations are based on Hartree-Fock pseudopotentials, with the single-particle orbitals entering the correlated wave function being represented by a systematically convergeable cubic-spline basis. Systematic tests are presented on system-size errors using periodically repeating cells of up to over 600 atoms. The equilibrium lattice parameter of the rock-salt structure obtained within DMC is almost identical to the Hartree-Fock result, which is close to the experimental value. The DMC result for the bulk modulus is also in good agreement with the experimental value. The B1-B2 transition pressure (between the rock-salt and CsCl structures) is predicted to be just below 600 GPa, which is beyond the experimentally accessible range, in accord with other predictions based on Hartree-Fock and density functional theories.

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
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