Author
José M. Soler
Other affiliations: International School for Advanced Studies, Massachusetts Institute of Technology, Harvard University
Bio: José M. Soler is an academic researcher from Autonomous University of Madrid. The author has contributed to research in topics: Density functional theory & van der Waals force. The author has an hindex of 53, co-authored 169 publications receiving 23325 citations. Previous affiliations of José M. Soler include International School for Advanced Studies & Massachusetts Institute of Technology.
Papers published on a yearly basis
Papers
More filters
••
TL;DR: In this paper, a selfconsistent density functional method using standard norm-conserving pseudopotentials and a flexible, numerical linear combination of atomic orbitals basis set, which includes multiple-zeta and polarization orbitals, was developed and implemented.
Abstract: We have developed and implemented a selfconsistent density functional method using standard norm-conserving pseudopotentials and a flexible, numerical linear combination of atomic orbitals basis set, which includes multiple-zeta and polarization orbitals. Exchange and correlation are treated with the local spin density or generalized gradient approximations. The basis functions and the electron density are projected on a real-space grid, in order to calculate the Hartree and exchange-correlation potentials and matrix elements, with a number of operations that scales linearly with the size of the system. We use a modified energy functional, whose minimization produces orthogonal wavefunctions and the same energy and density as the Kohn-Sham energy functional, without the need for an explicit orthogonalization. Additionally, using localized Wannier-like electron wavefunctions allows the computation time and memory required to minimize the energy to also scale linearly with the size of the system. Forces and stresses are also calculated efficiently and accurately, thus allowing structural relaxation and molecular dynamics simulations.
8,723 citations
••
TL;DR: A method to perform fully self-consistent density-functional calculations that scales linearly with the system size and which is well suited for very large systems is presented, using strictly localized pseudoatomic orbitals as basis functions.
Abstract: We present a method to perform fully self-consistent density-functional calculations that scales linearly with the system size and which is well suited for very large systems. It uses strictly localized pseudoatomic orbitals as basis functions. The sparse Hamiltonian and overlap matrices are calculated with an $O(N)$ effort. The long-range self-consistent potential and its matrix elements are computed in a real-space grid. The other matrix elements are directly calculated and tabulated as a function of the interatomic distances. The computation of the total energy and atomic forces is also done in $O(N)$ operations using truncated, Wannier-like localized functions to describe the occupied states, and a band-energy functional which is iteratively minimized with no orthogonality constraints. We illustrate the method with several examples, including carbon and silicon supercells with up to 1000 Si atoms and supercells of $\ensuremath{\beta}$-${\mathrm{C}}_{3}$${\mathrm{N}}_{4}$. We apply the method to solve the existing controversy about the faceting of large icosahedral fullerenes by performing dynamical simulations on ${\mathrm{C}}_{60}$, ${\mathrm{C}}_{240}$, and ${\mathrm{C}}_{540}$.
2,236 citations
••
TL;DR: In this article, a linear scaling, fully self-consistent density-functional method for performing first-principles calculations on systems with a large number of atoms, using standard norm-conserving pseudopotentials and flexible linear combinations of atomic orbitals (LCAO) basis sets, was implemented.
Abstract: We have implemented a linear scaling, fully self-consistent density-functional method for performing first-principles calculations on systems with a large number of atoms, using standard norm-conserving pseudopotentials and flexible linear combinations of atomic orbitals (LCAO) basis sets. Exchange and correlation are treated within the local-spin-density or gradient-corrected approximations. The basis functions and the electron density are projected on a real-space grid in order to calculate the Hartree and exchange–correlation potentials and matrix elements. We substitute the customary diagonalization procedure by the minimization of a modified energy functional, which gives orthogonal wave functions and the same energy and density as the Kohn–Sham energy functional, without the need of an explicit orthogonalization. The additional restriction to a finite range for the electron wave functions allows the computational effort (time and memory) to increase only linearly with the size of the system. Forces and stresses are also calculated efficiently and accurately, allowing structural relaxation and molecular dynamics simulations. We present test calculations beginning with small molecules and ending with a piece of DNA. Using double-z, polarized bases, geometries within 1% of experiments are obtained. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 65: 453–461, 1997
1,383 citations
••
TL;DR: This work presents an efficient implementation of the van der Waals density functional, which expresses the nonlocal correlation energy as a double spatial integral, and applies the method to calculate the binding energies and the barriers for relative translation and rotation in double-wall carbon nanotubes.
Abstract: We present an efficient implementation of the van der Waals density functional of Dion et al. [Phys. Rev. Lett. 92, 246401 (2004)], which expresses the nonlocal correlation energy as a double spatial integral. We factorize the integration kernel and use fast Fourier transforms to evaluate the self-consistent potential, total energy, and atomic forces, in $O(N\mathrm{log} N)$ operations. The resulting overhead, for medium and large systems, is a small fraction of the total computational cost, representing a dramatic speedup over the $O({N}^{2})$ evaluation of the double integral. This opens the realm of first-principles simulations to the large systems of interest in soft matter and biomolecular problems. We apply the method to calculate the binding energies and the barriers for relative translation and rotation in double-wall carbon nanotubes.
1,353 citations
••
TL;DR: A brief review of the Siesta project is presented in the context of linear-scaling density-functional methods for electronic-structure calculations and molecular-dynamics simulations of systems with a large number of atoms as discussed by the authors.
Abstract: A brief review of the Siesta project is presented in the context of linear-scaling density-functional methods for electronic-structure calculations and molecular-dynamics simulations of systems with a large number of atoms. Applications of the method to different systems are reviewed, including carbon nanotubes, gold nanostructures, adsorbates on silicon surfaces, and nucleic acids. Also, progress in atomic-orbital bases adapted to linear-scaling methodology is presented.
894 citations
Cited by
More filters
••
TL;DR: In this paper, a selfconsistent density functional method using standard norm-conserving pseudopotentials and a flexible, numerical linear combination of atomic orbitals basis set, which includes multiple-zeta and polarization orbitals, was developed and implemented.
Abstract: We have developed and implemented a selfconsistent density functional method using standard norm-conserving pseudopotentials and a flexible, numerical linear combination of atomic orbitals basis set, which includes multiple-zeta and polarization orbitals. Exchange and correlation are treated with the local spin density or generalized gradient approximations. The basis functions and the electron density are projected on a real-space grid, in order to calculate the Hartree and exchange-correlation potentials and matrix elements, with a number of operations that scales linearly with the size of the system. We use a modified energy functional, whose minimization produces orthogonal wavefunctions and the same energy and density as the Kohn-Sham energy functional, without the need for an explicit orthogonalization. Additionally, using localized Wannier-like electron wavefunctions allows the computation time and memory required to minimize the energy to also scale linearly with the size of the system. Forces and stresses are also calculated efficiently and accurately, thus allowing structural relaxation and molecular dynamics simulations.
8,723 citations
••
TL;DR: The most relevant features of WSXM, a freeware scanning probe microscopy software based on MS-Windows, are described and some relevant procedures of the software are carried out.
Abstract: In this work we briefly describe the most relevant features of WSXM, a freeware scanning probe microscopy software based on MS-Windows. The article is structured in three different sections: The introduction is a perspective on the importance of software on scanning probe microscopy. The second section is devoted to describe the general structure of the application; in this section the capabilities of WSXM to read third party files are stressed. Finally, a detailed discussion of some relevant procedures of the software is carried out.
6,996 citations
••
6,213 citations
••
TL;DR: In this paper, the 2D counterpart of layered black phosphorus, which is called phosphorene, is introduced as an unexplored p-type semiconducting material and the authors find that the band gap is direct, depends on the number of layers and the in-layer strain, and significantly larger than the bulk value of 0.31-0.36 eV.
Abstract: We introduce the 2D counterpart of layered black phosphorus, which we call phosphorene, as an unexplored p-type semiconducting material. Same as graphene and MoS2, single-layer phosphorene is flexible and can be mechanically exfoliated. We find phosphorene to be stable and, unlike graphene, to have an inherent, direct, and appreciable band gap. Our ab initio calculations indicate that the band gap is direct, depends on the number of layers and the in-layer strain, and is significantly larger than the bulk value of 0.31–0.36 eV. The observed photoluminescence peak of single-layer phosphorene in the visible optical range confirms that the band gap is larger than that of the bulk system. Our transport studies indicate a hole mobility that reflects the structural anisotropy of phosphorene and complements n-type MoS2. At room temperature, our few-layer phosphorene field-effect transistors with 1.0 μm channel length display a high on-current of 194 mA/mm, a high hole field-effect mobility of 286 cm2/V·s, and an...
5,233 citations
••
TL;DR: It is shown how derivatives of the GPW energy functional, namely ionic forces and the Kohn–Sham matrix, can be computed in a consistent way and the computational cost is scaling linearly with the system size, even for condensed phase systems of just a few tens of atoms.
4,047 citations