Francois Gygi

Towards an assessment of the accuracy of density functional theory for first principles simulations of water

A series of Car-Parrinello (CP) molecular dynamics simulations of water are presented, aimed at assessing the accuracy of density functional theory in describing the structural and dynamical properties of water at ambient conditions. We found negligible differences in structural properties obtained using the Perdew-Burke-Ernzerhof or the Becke-Lee-Yang-Parr exchange and correlation energy functionals; we also found that size effects, although not fully negligible when using 32 molecule cells, are rather small. In addition, we identified a wide range of values of the fictitious electronic mass (micro) entering the CP Lagrangian for which the electronic ground state is accurately described, yielding trajectories and average properties that are independent of the value chosen. However, care must be exercised not to carry out simulations outside this range, where structural properties may artificially depend on micro. In the case of an accurate description of the electronic ground state, and in the absence of proton quantum effects, we obtained an oxygen-oxygen correlation function that is overstructured compared to experiment, and a diffusion coefficient which is approximately ten times smaller.

Selective dispersion of high purity semiconducting single-walled carbon nanotubes with regioregular poly(3-alkylthiophene)s

Conjugated polymers, such as polyfluorene and poly(phenylene vinylene), have been used to selectively disperse semiconducting single-walled carbon nanotubes (sc-SWNTs), but these polymers have limited applications in transistors and solar cells. Regioregular poly(3-alkylthiophene)s (rr-P3ATs) are the most widely used materials for organic electronics and have been observed to wrap around SWNTs. However, no sorting of sc-SWNTs has been achieved before. Here we report the application of rr-P3ATs to sort sc-SWNTs. Through rational selection of polymers, solvent and temperature, we achieved highly selective dispersion of sc-SWNTs. Our approach enables direct film preparation after a simple centrifugation step. Using the sorted sc-SWNTs, we fabricate high-performance SWNT network transistors with observed charge-carrier mobility as high as 12 cm(2) V(-1) s(-1) and on/off ratio of >10(6). Our method offers a facile and a scalable route for separating sc-SWNTs and fabrication of 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.

The solvation of Na+ in water: First-principles simulations

First-principles molecular dynamics simulations have been performed on the solvation of Na+ in water. Consistent with the available experimental data, we find that the first solvation shell of Na+ contains on average 5.2 water molecules. A significant number of water exchanges between the first and second solvation shells are observed. However, the simulations are not long enough to reliably measure the rate of water exchange. Contrary to several previous studies, we do not find any effect of Na+ on the orientation of water molecules outside of the first solvation shell. Furthermore, the complete set of structural properties determined by first-principles molecular dynamics is not predicted by any of the known classical simulations.First-principles molecular dynamics simulations have been performed on the solvation of Na+ in water. Consistent with the available experimental data, we find that the first solvation shell of Na+ contains on average 5.2 water molecules. A significant number of water exchanges between the first and second solvation shells are observed. However, the simulations are not long enough to reliably measure the rate of water exchange. Contrary to several previous studies, we do not find any effect of Na+ on the orientation of water molecules outside of the first solvation shell. Furthermore, the complete set of structural properties determined by first-principles molecular dynamics is not predicted by any of the known classical simulations.

A first principles molecular dynamics simulation of the hydrated magnesium ion

Abstract First principles molecular dynamics has been used to investigate the solvation of Mg2+ in water. In agreement with experiment, we find that the first solvation shell around Mg2+ contains six water molecules in an octahedral arrangement. The electronic structure of first solvation shell water molecules has been examined with a localized orbital analysis. We find that water molecules tend to asymmetrically coordinate Mg2+ through one of the oxygen lone pair orbitals and that the first solvation shell dipole moments increase by 0.2 Debye relative to pure liquid water.

Structural and Vibrational Properties of Liquid Water from van der Waals Density Functionals

We present results for the structural and vibrational properties of the water molecule, water dimer, and liquid water at the experimental equilibrium density, as obtained with several van der Waals density functionals. The functional form originally proposed by Dion et al. [ Phys. Rev. Lett. 2004 , 92 , 246401 ], with an appropriately chosen local exchange functional, yields a description of the liquid superior to that of the semilocal functional PBE. In particular, a specific choice of the local exchange functional (optB88) fitted to quantum chemistry calculations yields the best agreement with experimental results for pair correlation functions although it is slightly inferior to other van der Waals functionals in describing infrared spectra. When using optB88, liquid water displays a hydrogen-bonded network less tightly bound than when using the PBE approximation. The performance of optB88 is definitely inferior to that of the PBE0 hybrid functional for the isolated molecule but only moderately so for the liquid. However, the computational cost of optB88 is much less than that of hybrid functionals; therefore the use of optB88 appears to be a sensible alternative to calculations implying the evaluation of the Fock operator, in cases when simulations of large systems are required.

Density functional theory for efficient ab initio molecular dynamics simulations in solution

We present a density functional for first‐principles molecular dynamics simulations that includes the electrostatic effects of a continuous dielectric medium. It allows for numerical simulations of molecules in solution in a model polar solvent. We propose a smooth dielectric model function to model solvation into water and demonstrate its good numerical properties for total energy calculations and constant energy molecular dynamics.

Optimization algorithm for the generation of ONCV pseudopotentials

Abstract We present an optimization algorithm to construct pseudopotentials and use it to generate a set of Optimized Norm-Conserving Vanderbilt (ONCV) pseudopotentials for elements up to Z = 83 (Bi) (excluding Lanthanides). We introduce a quality function that assesses the agreement of a pseudopotential calculation with all-electron FLAPW results, and the necessary plane-wave energy cutoff. This quality function allows us to use a Nelder–Mead optimization algorithm on a training set of materials to optimize the input parameters of the pseudopotential construction for most of the periodic table. We control the accuracy of the resulting pseudopotentials on a test set of materials independent of the training set. We find that the automatically constructed pseudopotentials ( http://www.quantum-simulation.org ) provide a good agreement with the all-electron results obtained using the FLEUR code with a plane-wave energy cutoff of approximately 60 Ry.

REAL-SPACE ADAPTIVE-COORDINATE ELECTRONIC-STRUCTURE CALCULATIONS

We present a real-space adaptive-coordinate method, which combines the advantages of the finite-difference approach with the accuracy and flexibility of the adaptive coordinate method. The discretized Kohn-Sham equations are written in generalized curvilinear coordinates and solved self-consistently by means of an iterative approach. The Poisson equation is solved in real space using the Multigrid algorithm. We implemented the method on a massively parallel computer, and applied it to the calculation of the equilibrium geometry and harmonic vibrational frequencies of the CO_2, CO, N_2 and F_2 molecules, yielding excellent agreement with the results of accurate quantum chemistry and Local Density Functional calculations.

Early chemistry in hot and dense nitromethane: Molecular dynamics simulations

We report density functional molecular dynamic simulations to determine the early chemical events of hot (T=3000 K) and dense (rho=1.97 g/cm(3), V/V(0)=0.68) nitromethane (CH(3)NO(2)). The first step in the decomposition process is an intermolecular proton abstraction mechanism that leads to the formation of CH(3)NO(2)H(+) and the aci ion H(2)CNO(2) (-). This event is also confirmed to occur in a fast annealing simulation to a final temperature of 4000 K at rho=2.20 g/cm(3). An intramolecular hydrogen transfer that transforms nitromethane into the aci acid form, CH(2)NO(2)H, accompanies this event. To our knowledge, this is the first confirmation of chemical reactivity with bond selectivity for an energetic material near the Chapman-Jouget state of the fully reacted material. We also report the decomposition mechanism followed up to the formation of H(2)O as the first stable product. We note that similarities in the global features of reactants, intermediates, and products of the reacting fluid seem to indicate a threshold for similar chemistry in the range of high densities and temperatures reported herein.