SIESTA (Spanish Initiative for Electronic Simulations with Thousands of Atoms) is an original method and its computer program implementation, to efficiently perform electronic structure calculations and ab initiomolecular dynamics simulations of molecules and solids. SIESTA uses strictly localized basis sets and the implementation of linear-scaling algorithms. Accuracy and speed can be set in a wide range, from quick exploratory calculations to highly accurate simulations matching the quality of other approaches, such as the plane-wave and all-electron methods.
SIESTA's backronym is the Spanish Initiative for Electronic Simulations with Thousands of Atoms.
Since 13 May 2016, with the 4.0 version announcement, SIESTA is released under the terms of the GPL open-source license. Source packages and access to the development versions can be obtained from the DevOps platform on GitLab.[2] The latest version Siesta-4.1.5 was released on 4 February 2021.
It uses norm-conserving pseudopotentials in their fully non-local (Kleinman-Bylander) form.
It uses atomic orbitals as a basis set, allowing unlimited multiple-zeta and angular momenta, polarization, and off-site orbitals. The radial shape of every orbital is numerical, and any shape can be used and provided by the user, with the only condition that it has to be of finite support, i.e., it has to be strictly zero beyond a user-provided distance from the corresponding nucleus. Finite-support basis sets are the key to calculating the Hamiltonian and overlap matrices in O(N) operations.
Projects the electron wave functions and density onto a real-space grid to calculate the Hartree and exchange-correlation potentials and their matrix elements.
Besides the standard Rayleigh-Ritz eigenstate method, it allows the use of localized linear combinations of the occupied orbitals (valence-bond or Wannier-like functions), making the computer time and memory scale linearly with the number of atoms. Simulations with several hundred atoms are feasible with modest workstations.
It is written in Fortran 95 and memory is allocated dynamically.
It may be compiled for serial or parallel execution (under MPI).
Ballistic electron transport under non-equilibrium (through TranSIESTA)
Strengths of SIESTA
SIESTA's main strengths are:
Flexible accuracy and speed.
It can tackle computationally demanding systems (systems currently out of the reach of plane-wave codes).[citation needed]
Efficient parallelization.
The use of a linear combination of numerical atomic orbitals makes SIESTA a DFT code. SIESTA can produce very fast calculations with small basis sets, allowing the computation of systems with thousands of atoms. Alternatively, the use of more complete and accurate bases achieves accuracies comparable to those of standard plane wave calculations, with competitive performance.
Implemented Solutions
SIESTA is in continuous development since it was implemented in 1996. The main solutions implemented in the current version are:
Collinear and non-collinear spin-polarised calculations
Efficient implementation of Van der Waals functional
Several post-processing tools for SIESTA have been developed. These programs process SIESTA output or provide additional features.
Applications
Since its implementation, SIESTA has been used by researchers in geosciences, biology, and engineering (extending beyond materials physics and chemistry) and has been applied to a large variety of systems including surfaces, adsorbates, nanotubes, nanoclusters, biological molecules, amorphous semiconductors, ferroelectric films, low-dimensional metals, etc.[3][4][5]
García, Alberto; Papior, Nick; Akhtar, Arsalan; Artacho, Emilio; Blum, Volker; Bosoni, Emanuele; Brandimarte, Pedro; Brandbyge, Mads; Cerdá, J.I.; Corsetti, Fabiano; Cuadrado, Ramón; Dikan, Vladimir; Ferrer, Jaime; Gale, Julian; García-Fernández, Pablo; García-Suárez, V.M.; García, Sandra; Huhs, Georg; Illera, Sergio; Korytár, Richard; Koval, Peter; Lebedeva, Irina; Lin, Lin; López-Tarifa, Pablo; G. Mayo, Sara; Mohr, Stephan; Ordejón, Pablo; Postnikov, Andrei; Pouillon, Yann; Pruneda, Miguel; Robles, Roberto; Sánchez-Portal, Daniel; Soler, Jose M.; Ullah, Rafi; Yu, Victor Wen-zhe; Junquera, Javier (2020). "Siesta: Recent developments and applications". Journal of Chemical Physics. 152 (20): 204108. doi:10.1063/5.0005077. hdl:10902/20680. PMID32486661. S2CID219179270. Postprint is available at hdl:10261/213028.
Izquierdo, J.; Vega, A.; Balbás, L.; Sánchez-Portal, Daniel; Junquera, Javier; Artacho, Emilio; Soler, Jose; Ordejón, Pablo (2000). "Systematic ab initio study of the electronic and magnetic properties of different pure and mixed iron systems". Physical Review B. 61 (20): 13639. Bibcode:2000PhRvB..6113639I. doi:10.1103/PhysRevB.61.13639.
^Mashaghi A et al. Hydration strongly affects the molecular and electronic structure of membrane phospholipids J. Chem. Phys. 136, 114709 (2012) [1]
^Mashaghi A et al. Interfacial Water Facilitates Energy Transfer by Inducing Extended Vibrations in Membrane Lipids, J. Phys. Chem. B, 2012, 116 (22), pp 6455–6460 [2]
^Mashaghi A et al. Enhanced Autoionization of Water at Phospholipid Interfaces. J. Phys. Chem. C, 2013, 117 (1), pp 510–514 [3]
