Official Resources
- Homepage: https://siesta-project.org/siesta/
- Documentation: https://docs.siesta-project.org/
- Source Repository: https://gitlab.com/siesta-project/siesta
- License: GNU General Public License v3.0
Overview
SIESTA (Spanish Initiative for Electronic Simulations with Thousands of Atoms) is an efficient DFT code using strictly localized numerical atomic orbital basis sets. It excels at large-scale calculations with linear-scaling capabilities and is particularly strong for low-dimensional systems, molecules, and quantum transport calculations via TranSIESTA.
Scientific domain: Large-scale materials, nanostructures, surfaces, molecules, quantum transport
Target user community: Researchers needing efficient DFT for large systems (1000+ atoms) and electronic transport
Theoretical Methods
- Density Functional Theory (DFT)
- Numerical atomic orbital (NAO) basis sets
- Strictly localized basis functions
- Norm-conserving pseudopotentials
- LDA, GGA functionals
- van der Waals corrections (DFT-D, VDW-DF)
- DFT+U for correlated systems
- Spin-orbit coupling (full implementation)
- Non-collinear magnetism
- Time-Dependent DFT (Real-Time TDDFT)
- Linear-scaling O(N) DFT
Capabilities (CRITICAL)
- Ground-state electronic structure
- Geometry optimization and MD (NVE, NVT, NPT)
- Large-scale systems (1000+ atoms)
- Linear-scaling DFT for very large systems
- Band structure and DOS
- Forces and stress tensors
- Phonon calculations via finite differences
- Molecular dynamics
- Quantum transport (TranSIESTA)
- Non-equilibrium Green's function (NEGF) for transport
- Spinor quantum transport (with SOC)
- STM image simulation
- Optical properties
- Electric polarization
- Wannier functions
- Constrained DFT
- Thermostats and barostats
- Variable cell dynamics
Sources: Official SIESTA documentation, cited in 7/7 source lists
Key Strengths
Linear-Scaling DFT:
- O(N) algorithms for large systems
- Strict locality of basis functions
- Efficient sparse matrix operations
- Systems up to millions of atoms
TranSIESTA Transport:
- Non-Equilibrium Green's Function (NEGF)
- Ballistic electron transport
- Zero-bias and finite-bias I-V curves
- Multi-electrode calculations
- Semi-infinite electrode treatment
Spinor Quantum Transport:
- Full spinor wave functions in transport
- Spin-orbit coupling effects
- Topological material transport
- Non-collinear spin transport
- Ultra-low-energy electronics applications
Computational Efficiency:
- Strictly localized basis sets
- Sparse matrix storage
- Efficient pseudopotential handling
- MPI, OpenMP, GPU parallelization
Open-Source Ecosystem:
- GPL v3 license
- Active GitLab development
- Large user community
- Extensive third-party tools
Inputs & Outputs
-
Input formats:
- fdf input files (Flexible Data Format)
- XV, XYZ coordinate files
- Pseudopotential files (.psf, .vps, .psml)
-
Output data types:
- Standard output with energies, forces
- XV files (structures)
- Density matrices
- DOS and band structure files
- LDOS, PDOS files
- Molecular dynamics trajectories
- Transmission coefficients
Interfaces & Ecosystem
-
Framework integrations:
- ASE - calculator interface
- Phonopy - phonon calculations
- pymatgen - structure I/O
- Deneb - GUI and workflow
- AiiDA-SIESTA - automated workflows
-
Transport calculations:
- TranSIESTA - quantum transport (integrated)
- TBtrans - post-processing transport data
- Smeagol - transport properties
- inelastica - inelastic transport
-
Utilities:
- Util/ directory with analysis tools
- sisl - Python interface to SIESTA
- ATOM - pseudopotential generation
-
Post-processing:
- Denchar - charge density plotting
- grid2cube - grid file conversion
- fatbands - fat band analysis
Advanced Features
SIESTA 5.0 (May 2024):
- CMake-based building framework
- Real-Time TDDFT implementation
- D3 dispersion corrections
- PSML pseudopotential format support
- Improved parallelization
Multi-Electrode Transport:
- N-electrode TranSIESTA calculations
- Complex device geometries
- Improved inversion algorithms
- TBtrans rewrite for flexibility
Spin-Orbit Coupling:
- Full SOC implementation
- Spin texture calculations
- Topological material studies
- Interface with sisl for analysis
Workflow Automation:
- MPI parallelization
- OpenMP threading
- GPU offloading
- High-throughput ready
Performance Characteristics
- Speed: Very efficient for O(N) calculations
- Accuracy: Good for localized basis systems
- System size: Up to millions of atoms (O(N))
- Memory: Efficient sparse storage
- Parallelization: Excellent MPI/OpenMP/GPU
Computational Cost
- Standard DFT: Efficient with NAO basis
- Large systems: Linear scaling available
- Transport: Moderate cost for NEGF
- MD: Production-ready speeds
- Typical: Competitive performance
Limitations & Known Constraints
- Basis sets: NAO basis sets require careful convergence testing
- Basis completeness: Strictly localized basis less complete than plane-waves
- Pseudopotentials: Limited to norm-conserving; quality varies
- Accuracy: Generally less accurate than plane-wave codes for same computational cost
- Overlap matrix: Can become ill-conditioned for small basis cutoffs
- Documentation: Comprehensive but can be overwhelming
Comparison with Other Codes
- vs VASP/QE: SIESTA faster for large systems, less accurate per atom
- vs FHI-aims: SIESTA pseudopotential, FHI-aims all-electron
- vs CONQUEST: Both O(N), different basis approaches
- vs OpenMX: Similar capabilities, different communities
- Unique strength: TranSIESTA for quantum transport, large-scale O(N), open-source
Application Areas
Nanoscale Transport:
- Molecular electronics
- Nanojunctions
- 2D material devices
- Spintronic devices
Large-Scale Materials:
- Nanostructures
- Amorphous materials
- Complex interfaces
- Biological systems
Surface Science:
- Adsorption studies
- Surface reconstruction
- Catalysis
- STM simulations
2D Materials:
- Graphene and derivatives
- Transition metal dichalcogenides
- Van der Waals heterostructures
- Topological systems
Best Practices
Basis Set Selection:
- SZ (single-zeta) for quick tests
- DZP (double-zeta polarized) for production
- TZP for high accuracy
- Optimize PAO.EnergyShift
Pseudopotentials:
- Use quality-tested PSML sets
- Check transferability
- Test against reference calculations
- Include semicore if needed
Transport Calculations:
- Converge electrode calculations first
- Check contact self-energies
- Test k-point sampling
- Verify buffer regions
Convergence:
- MeshCutoff for grid spacing
- k-point sampling for periodic
- PAO cutoff radius
- SCF tolerance
Community and Support
- Open-source GPL v3
- Active GitLab development
- Mailing lists for support
- Annual SIESTA schools
- Large international community
Verification & Sources
Primary sources:
- Official website: https://siesta-project.org/siesta/
- Documentation: https://docs.siesta-project.org/
- GitLab repository: https://gitlab.com/siesta-project/siesta
- J. M. Soler et al., J. Phys. Condens. Matter 14, 2745 (2002) - SIESTA method
- E. Artacho et al., Phys. Status Solidi B 215, 809 (1999) - Linear-scaling
- A. García et al., J. Chem. Phys. 152, 204108 (2020) - Recent developments
Secondary sources:
- SIESTA manual and tutorials
- Published DFT studies using SIESTA (>10,000 citations)
- Workshop materials
- Confirmed in 7/7 source lists (claude, g, gr, k, m, q, z)
Confidence: CONFIRMED - Appears in ALL 7 independent source lists
Verification status: ✅ VERIFIED
- Official homepage: ACCESSIBLE
- Documentation: COMPREHENSIVE and ACCESSIBLE
- Source code: OPEN (GitLab, GPL v3)
- Community support: Active mailing list, workshops
- Academic citations: >10,000
- Active development: Regular releases, large community
- Specialized strength: O(N) methods, TranSIESTA quantum transport, large systems, open-source