Matt Probert

First-principles simulation: ideas, illustrations and the CASTEP code

First-principles simulation, meaning density-functional theory calculations with plane waves and pseudopotentials, has become a prized technique in condensed-matter theory. Here I look at the basics of the suject, give a brief review of the theory, examining the strengths and weaknesses of its implementation, and illustrating some of the ways simulators approach problems through a small case study. I also discuss why and how modern software design methods have been used in writing a completely new modular version of the CASTEP code.

First principles methods using CASTEP

Abstract The CASTEP code for first principles electronic structure calculations will be described. A brief, non-technical overview will be given and some of the features and capabilities highlighted. Some features which are unique to CASTEP will be described and near-future development plans outlined.

Exploring atomic defects in molybdenum disulphide monolayers

Defects usually play an important role in tailoring various properties of two-dimensional materials. Defects in two-dimensional monolayer molybdenum disulphide may be responsible for large variation of electric and optical properties. Here we present a comprehensive joint experiment–theory investigation of point defects in monolayer molybdenum disulphide prepared by mechanical exfoliation, physical and chemical vapour deposition. Defect species are systematically identified and their concentrations determined by aberration-corrected scanning transmission electron microscopy, and also studied by ab-initio calculation. Defect density up to 3.5 × 1013 cm−2 is found and the dominant category of defects changes from sulphur vacancy in mechanical exfoliation and chemical vapour deposition samples to molybdenum antisite in physical vapour deposition samples. Influence of defects on electronic structure and charge-carrier mobility are predicted by calculation and observed by electric transport measurement. In light of these results, the growth of ultra-high-quality monolayer molybdenum disulphide appears a primary task for the community pursuing high-performance electronic devices.

Reproducibility in density functional theory calculations of solids

A comparison of DFT methods Density functional theory (DFT) is now routinely used for simulating material properties. Many software packages are available, which makes it challenging to know which are the best to use for a specific calculation. Lejaeghere et al. compared the calculated values for the equation of states for 71 elemental crystals from 15 different widely used DFT codes employing 40 different potentials (see the Perspective by Skylaris). Although there were variations in the calculated values, most recent codes and methods converged toward a single value, with errors comparable to those of experiment. Science, this issue p. 10.1126/science.aad3000; see also p. 1394 A survey of recent density functional theory methods shows a convergence to more accurate property calculations. [Also see Perspective by Skylaris] INTRODUCTION The reproducibility of results is one of the underlying principles of science. An observation can only be accepted by the scientific community when it can be confirmed by independent studies. However, reproducibility does not come easily. Recent works have painfully exposed cases where previous conclusions were not upheld. The scrutiny of the scientific community has also turned to research involving computer programs, finding that reproducibility depends more strongly on implementation than commonly thought. These problems are especially relevant for property predictions of crystals and molecules, which hinge on precise computer implementations of the governing equation of quantum physics. RATIONALE This work focuses on density functional theory (DFT), a particularly popular quantum method for both academic and industrial applications. More than 15,000 DFT papers are published each year, and DFT is now increasingly used in an automated fashion to build large databases or apply multiscale techniques with limited human supervision. Therefore, the reproducibility of DFT results underlies the scientific credibility of a substantial fraction of current work in the natural and engineering sciences. A plethora of DFT computer codes are available, many of them differing considerably in their details of implementation, and each yielding a certain “precision” relative to other codes. How is one to decide for more than a few simple cases which code predicts the correct result, and which does not? We devised a procedure to assess the precision of DFT methods and used this to demonstrate reproducibility among many of the most widely used DFT codes. The essential part of this assessment is a pairwise comparison of a wide range of methods with respect to their predictions of the equations of state of the elemental crystals. This effort required the combined expertise of a large group of code developers and expert users. RESULTS We calculated equation-of-state data for four classes of DFT implementations, totaling 40 methods. Most codes agree very well, with pairwise differences that are comparable to those between different high-precision experiments. Even in the case of pseudization approaches, which largely depend on the atomic potentials used, a similar precision can be obtained as when using the full potential. The remaining deviations are due to subtle effects, such as specific numerical implementations or the treatment of relativistic terms. CONCLUSION Our work demonstrates that the precision of DFT implementations can be determined, even in the absence of one absolute reference code. Although this was not the case 5 to 10 years ago, most of the commonly used codes and methods are now found to predict essentially identical results. The established precision of DFT codes not only ensures the reproducibility of DFT predictions but also puts several past and future developments on a firmer footing. Any newly developed methodology can now be tested against the benchmark to verify whether it reaches the same level of precision. New DFT applications can be shown to have used a sufficiently precise method. Moreover, high-precision DFT calculations are essential for developing improvements to DFT methodology, such as new density functionals, which may further increase the predictive power of the simulations. Recent DFT methods yield reproducible results. Whereas older DFT implementations predict different values (red darts), codes have now evolved to mutual agreement (green darts). The scoreboard illustrates the good pairwise agreement of four classes of DFT implementations (horizontal direction) with all-electron results (vertical direction). Each number reflects the average difference between the equations of state for a given pair of methods, with the green-to-red color scheme showing the range from the best to the poorest agreement. The widespread popularity of density functional theory has given rise to an extensive range of dedicated codes for predicting molecular and crystalline properties. However, each code implements the formalism in a different way, raising questions about the reproducibility of such predictions. We report the results of a community-wide effort that compared 15 solid-state codes, using 40 different potentials or basis set types, to assess the quality of the Perdew-Burke-Ernzerhof equations of state for 71 elemental crystals. We conclude that predictions from recent codes and pseudopotentials agree very well, with pairwise differences that are comparable to those between different high-precision experiments. Older methods, however, have less precise agreement. Our benchmark provides a framework for users and developers to document the precision of new applications and methodological improvements.

A periodic genetic algorithm with real-space representation for crystal structure and polymorph prediction

A genetic algorithm is described that is suitable for determining the global minimum energy configurations of crystal structures and which can also be used as a polymorph search technique. This algorithm requires no prior assumptions about unit cell size, shape, or symmetry, nor about the ionic configuration within the unit cell. This therefore enables true ab initio crystal structure and polymorph prediction. Our algorithm uses a real-space representation of the population members, and makes use of a periodic cut for the crossover operation. Results on large Lennard-Jones systems with fcc- and hcp-commensurate cells show robust convergence to the bulk structure from a random initial assignment and an ability to successfully discriminate between competing low enthalpy configurations. Results from an ab initio carbon polymorph search show the spontaneous emergence of both Lonsdaleite and graphite-like structures.

Stiffness and thermal expansion of ZrB2 : an ab initio study

The stiffness and thermal expansion coefficient of ZrB2 are calculated within the density functional theory formalism. The stiffness tensor obtained here using the static finite strain technique is in good agreement with the results of resonant ultrasonic measurements and points to a possible misinterpretation of the experimentally obtained compression data. The methodology of evaluating thermal expansion coefficients from molecular dynamics simulations for small unit cells is validated for a number of systems: metals, semiconductors and insulators.

Improving the convergence of defect calculations in supercells: An ab initio study of the neutral silicon vacancy

We present a systematic methodology for the accurate calculation of defect structures in supercells, which we illustrate with a study of the neutral vacancy in silicon. This is a prototypical defect which has been studied extensively using ab initio methods, yet remarkably there is still no consensus about the energy or structure of this defect, or even whether the nearest-neighbor atoms relax inwards or outwards. In this paper, we show that the differences between previous calculations can be attributed to supercell convergence errors, and we demonstrate how to systematically reduce each such source of error. The various sources of scatter in previous theoretical studies are discussed and a different effect, that of supercell symmetry, is identified. It is shown that a consistent treatment of this effect is crucial in understanding the systematic effects of increasing the supercell size. This work therefore also presents the best converged ab initio study of the neutral silicon vacancy to date.

Density functional theory in the solid state.

Density functional theory (DFT) has been used in many fields of the physical sciences, but none so successfully as in the solid state. From its origins in condensed matter physics, it has expanded into materials science, high-pressure physics and mineralogy, solid-state chemistry and more, powering entire computational subdisciplines. Modern DFT simulation codes can calculate a vast range of structural, chemical, optical, spectroscopic, elastic, vibrational and thermodynamic phenomena. The ability to predict structure–property relationships has revolutionized experimental fields, such as vibrational and solid-state NMR spectroscopy, where it is the primary method to analyse and interpret experimental spectra. In semiconductor physics, great progress has been made in the electronic structure of bulk and defect states despite the severe challenges presented by the description of excited states. Studies are no longer restricted to known crystallographic structures. DFT is increasingly used as an exploratory tool for materials discovery and computational experiments, culminating in ex nihilo crystal structure prediction, which addresses the long-standing difficult problem of how to predict crystal structure polymorphs from nothing but a specified chemical composition. We present an overview of the capabilities of solid-state DFT simulations in all of these topics, illustrated with recent examples using the CASTEP computer program.

Langevin Dynamics in Constant Pressure Extended Systems

The advantages of performing Langevin dynamics in extended systems are discussed. A simple Langevin dynamics scheme for producing the canonical ensemble is reviewed, and is then extended to the Hoover ensemble. We show that the resulting equations of motion generate the isobaric-isothermal ensemble. The Parrinello-Rahman ensemble is then discussed and we show that despite the presence of intrinsic probability gradients in this system, a Langevin dynamics approach samples the extended phase space in the correct fashion. The implementation of these methods in the ab initio plane wave density functional theory code CASTEP [M. D. Segall, P. L. D. Lindan, M. J. Probert, C. J. Pickard, P. J. Hasnip, S. J. Clarke, and M. C. Payne, J. Phys.: Condens. Matter 14, 2717 (2003)] is demonstrated.

Quantum simulation of low-temperature metallic liquid hydrogen

The melting temperature of solid hydrogen drops with pressure above ~65 GPa, suggesting that a liquid state might exist at low temperatures. It has also been suggested that this low-temperature liquid state might be non-molecular and metallic, although evidence for such behaviour is lacking. Here we report results for hydrogen at high pressures using ab initio methods, which include a description of the quantum motion of the protons. We determine the melting temperature as a function of pressure and find an atomic solid phase from 500 to 800 GPa, which melts at <200 K. Beyond this and up to 1,200 GPa, a metallic atomic liquid is stable at temperatures as low as 50 K. The quantum motion of the protons is critical to the low melting temperature reported, as simulations with classical nuclei lead to considerably higher melting temperatures of ~300 K across the entire pressure range considered.