Ab initio molecular dynamics for liquid metals.
TL;DR: In this paper, the authors present an ab initio quantum-mechanical molecular-dynamics calculations based on the calculation of the electronic ground state and of the Hellmann-Feynman forces in the local density approximation.
Abstract: We present ab initio quantum-mechanical molecular-dynamics calculations based on the calculation of the electronic ground state and of the Hellmann-Feynman forces in the local-density approximation at each molecular-dynamics step. This is possible using conjugate-gradient techniques for energy minimization, and predicting the wave functions for new ionic positions using subspace alignment. This approach avoids the instabilities inherent in quantum-mechanical molecular-dynamics calculations for metals based on the use of a fictitious Newtonian dynamics for the electronic degrees of freedom. This method gives perfect control of the adiabaticity and allows us to perform simulations over several picoseconds.
Citations
More filters
TL;DR: A detailed description and comparison of algorithms for performing ab-initio quantum-mechanical calculations using pseudopotentials and a plane-wave basis set is presented in this article. But this is not a comparison of our algorithm with the one presented in this paper.
47,666 citations
TL;DR: The simulation allows us to study in detail the changes in the structure-property relationship through the metal-semiconductor transition, and a detailed analysis of the local structural properties and their changes induced by an annealing process is reported.
Abstract: We present ab initio quantum-mechanical molecular-dynamics simulations of the liquid-metal--amorphous-semiconductor transition in Ge. Our simulations are based on (a) finite-temperature density-functional theory of the one-electron states, (b) exact energy minimization and hence calculation of the exact Hellmann-Feynman forces after each molecular-dynamics step using preconditioned conjugate-gradient techniques, (c) accurate nonlocal pseudopotentials, and (d) Nos\'e dynamics for generating a canonical ensemble. This method gives perfect control of the adiabaticity of the electron-ion ensemble and allows us to perform simulations over more than 30 ps. The computer-generated ensemble describes the structural, dynamic, and electronic properties of liquid and amorphous Ge in very good agreement with experiment. The simulation allows us to study in detail the changes in the structure-property relationship through the metal-semiconductor transition. We report a detailed analysis of the local structural properties and their changes induced by an annealing process. The geometrical, bonding, and spectral properties of defects in the disordered tetrahedral network are investigated and compared with experiment.
16,744 citations
TL;DR: In this article, an algorithm for decomposition of electronic charge density into atomic contributions is presented. But instead of explicitly finding and representing the dividing surfaces, which is a challenging task, the algorithm assigns each point on a regular (x,y,z) grid to one of the regions by following a steepest ascent path on the grid.
7,231 citations
TL;DR: Titanium dioxide is the most investigated single-crystalline system in the surface science of metal oxides, and the literature on rutile (1.1) and anatase surfaces is reviewed in this paper.
7,056 citations
TL;DR: This paper describes how accurate off-lattice ascent paths can be represented with respect to the grid points, and maintains the efficient linear scaling of an earlier version of the algorithm, and eliminates a tendency for the Bader surfaces to be aligned along the grid directions.
Abstract: A computational method for partitioning a charge density grid into Bader volumes is presented which is efficient, robust, and scales linearly with the number of grid points. The partitioning algorithm follows the steepest ascent paths along the charge density gradient from grid point to grid point until a charge density maximum is reached. In this paper, we describe how accurate off-lattice ascent paths can be represented with respect to the grid points. This improvement maintains the efficient linear scaling of an earlier version of the algorithm, and eliminates a tendency for the Bader surfaces to be aligned along the grid directions. As the algorithm assigns grid points to charge density maxima, subsequent paths are terminated when they reach previously assigned grid points. It is this grid-based approach which gives the algorithm its efficiency, and allows for the analysis of the large grids generated from plane-wave-based density functional theory calculations.
5,417 citations
Cites methods from "Ab initio molecular dynamics for li..."
...The charge density was generated using the VASP planewave-based DFT code [22], using the PW91 generalized gradient functional [23]....
[...]