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.

The performance of basis sets made of numerical atomic orbitals is explored in density-functional calculations of solids and molecules. With the aim of optimizing basis quality while maintaining strict localization of the orbitals, as needed for linear-scaling calculations, several schemes have been tried. The best performance is obtained for the basis sets generated according to a new scheme presented here, a flexibilization of previous proposals. The basis sets are tested versus converged plane-wave calculations on a significant variety of systems, including covalent, ionic and metallic. Satisfactory convergence (deviations significantly smaller than the accuracy of the underlying theory) is obtained for reasonably small basis sizes, with a clear improvement over previous schemes. The transferability of the obtained basis sets is tested in several cases and it is found to be satisfactory as well.

We have implemented a linear scaling, fully self-consistent densityfunctional 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.

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.

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.

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.

Electromagnetic field localization in nanoantennas is one of the leitmotivs that drives the development of plasmonics. The near-fields in these plasmonic nanoantennas are commonly addressed theoretically within classical frameworks that neglect atomic-scale features. This approach is often appropriate since the irregularities produced at the atomic scale are typically hidden in far-field optical spectroscopies. However, a variety of physical and chemical processes rely on the fine distribution of the local fields at this ultraconfined scale. We use time-dependent density functional theory and perform atomistic quantum mechanical calculations of the optical response of plasmonic nanoparticles, and their dimers, characterized by the presence of crystallographic planes, facets, vertices, and steps. Using sodium clusters as an example, we show that the atomistic details of the nanoparticles morphologies determine the presence of subnanometric near-field hot spots that are further enhanced by the action of the underlying nanometric plasmonic fields. This situation is analogue to a self-similar nanoantenna cascade effect, scaled down to atomic dimensions, and it provides new insights into the limits of field enhancement and confinement, with important implications in the optical resolution of field-enhanced spectroscopies and microscopies.

No abstract