Jose M. Soler

The SIESTA method for ab initio order-N materials simulation

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.An efficient method is described to handle mesh indexes in multidimensional problems like numerical integration of partial differential equations, lattice model simulations, and determination of atomic neighbor lists. By creating an extended mesh, beyond the periodic unit cell, the stride in memory between equivalent pairs of mesh points is independent of their position within the cell. This allows to contract the mesh indexes of all dimensions into a single index, avoiding modulo and other implicit index operations.

Self-consistent order- N density-functional calculations for very large systems

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 {ital O}({ital 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 {ital O}({ital 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 {beta}-C{sub 3}N{sub 4}. We apply the method to solve the existing controversy about the faceting of large icosahedral fullerenes by performing dynamical simulations on C{sub 60}, C{sub 240}, and C{sub 540}. {copyright} {ital 1996 The American Physical Society.}

Ab-initio structural, elastic, and vibrational properties of carbon nanotubes

A study based on ab initio calculations is presented on the structural, elastic, and vibrational properties of single-wall carbon nanotubes with different radii and chiralities. These properties are obtained using an implementation of pseudopotential-density-functional theory which allows calculations on systems with a large number of atoms per cell. Different quantities are monitored versus tube radius. The validity of expectations based on graphite is explored down to small radii, where some deviations appear related to the curvature-induced rehibridization of the carbon orbitals. Young moduli are found to be very similar to graphite and do not exhibit a systematic variation with either the radius or the chirality. The Poisson ratio also retains graphitic values except for a possible slight reduction for small radii. It shows, however, chirality dependence. The behavior of characteristic phonon branches as the breathing mode, twistons, and high-frequency optic modes, is also studied, the latter displaying a small chirality dependence at the top of the band. The results are compared with the predictions of the simple zone-folding approximation. Except for the known defficiencies of the zone-folding procedure in the low-frequency vibrational regions, it offers quite accurate results, even for relatively small radii.

Efficient implementation of a van der Waals density functional: application to double-wall carbon nanotubes.

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(NlogN) 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.

LINEAR-SCALING AB-INITIO CALCULATIONS FOR LARGE AND COMPLEX SYSTEMS

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.

The SIESTA method; developments and applicability

Recent developments in and around the SIESTA method of first-principles simulation of condensed matter are described and reviewed, with emphasis on (i) the applicability of the method for large and varied systems, (ii) efficient basis sets for the standards of accuracy of density-functional methods, (iii) new implementations, and (iv) extensions beyond ground-state calculations.

Projection of plane-wave calculations into atomic orbitals

The projection of the eigenfunctions obtained in standard plane-wave first-principle electronic-structure calculations into atomic-orbital basis sets is proposed as a formal and practical link between the methods based on plane waves and the ones based on atomic orbitals. Given a candidate atomic basis, (i) its quality is evaluated by its projection into the plane-wave eigenfunctions, (ii) it is optimized by maximizing that projection, (iii) the associated tight-binding hamiltonian and energy bands are obtained, and (iv) population analysis is performed in a natural way. The proposed method replaces the traditional trial-and-error procedures of finding appropriate atomic bases and the fitting of bands to obtain tight-binding hamiltonians. Test calculations of some zincblende semiconductors are presented.

Density, structure, and dynamics of water: The effect of van der Waals interactions

It is known that ab initio molecular dynamics (AIMD) simulations of liquid water at ambient conditions, based on the generalized gradient approximation (GGA) to density functional theory (DFT), with commonly used functionals fail to produce structural and diffusive properties in reasonable agreement with experiment. This is true for canonical, constant temperature simulations where the density of the liquid is fixed to the experimental density. The equilibrium density, at ambient conditions, of DFT water has recently been shown by Schmidt et al. [J. Phys. Chem. B, 113, 11959 (2009)] to be underestimated by different GGA functionals for exchange and correlation, and corrected by the addition of interatomic pair potentials to describe van der Waals (vdW) interactions. In this contribution we present a DFT-AIMD study of liquid water using several GGA functionals as well as the van der Waals density functional (vdW-DF) of Dion et al. [Phys. Rev. Lett. 92, 246401 (2004)]. As expected, we find that the density of water is grossly underestimated by GGA functionals. When a vdW-DF is used, the density improves drastically and the experimental diffusivity is reproduced without the need of thermal corrections. We analyze the origin of the density differences between all the functionals. We show that the vdW-DF increases the population of non-H-bonded interstitial sites, at distances between the first and second coordination shells. However, it excessively weakens the H-bond network, collapsing the second coordination shell. This structural problem is partially associated to the choice of GGA exchange in the vdW-DF. We show that a different choice for the exchange functional is enough to achieve an overall improvement both in structure and diffusivity.

Stiff Monatomic Gold Wires with a Spinning Zigzag Geometry

Using first principles density functional calculations, gold monatomic wires are found to exhibit a zigzag shape which remains under tension, becoming linear just before breaking. At room temperature they are found to spin, what explains the extremely long apparent interatomic distances shown by electron microscopy.The zigzag structure is stable if the tension is relieved, the wire holding its chainlike shape even as a free-standing cluster. This unexpected metallic-wire stiffness stems from the transverse quantization in the wire, as shown in a simple free electron model.

Systematic generation of finite-range atomic basis sets for linear-scaling calculations

Basis sets of atomic orbitals are very efficient for density functional calculations but lack a systematic variational convergence. We present a method to optimize numerical atomic orbitals variationally, using a single parameter to control their range. The efficiency of the basis generation scheme is tested and compared with other schemes for multiple