Fiesta: a Gaussian-basis
GW and Bethe-Salpeter code
This is a very minimal Fiesta website. Detailed features, benchmark calculations, examples, can be found in the related papers (see Publications list).
Fiesta ad minima.
The Fiesta code implements the GW and Bethe-Salpeter formalisms using Gaussian bases and
resolution-of-the-identity techniques (RI-SVS density and RI-V Coulomb metric). Dynamical
screening contribution to the self-energy is explicitely accounted for through a contour
deformation approach. Self-consistency on the wavefunctions is implemented at the static
COHSEX level. Tamm-Dancoff approximation (TDA) or full Bethe-Salpeter
calculations can be performed. The code presently reads input Kohn-Sham eigenstates
from the Orca and
Any DFT code dumping all Kohn-Sham eigenstates (occupied/unoccupied) expressed on a Gaussian basis set can be branched rather straighforwardly onto the Fiesta code.
The Fiesta code implements continuous polarizable models (PCM) and is merged with the Mescal Discrete Polarizable Model (DPM) [see: D'Avino et al., JCTC 2014; Li, D'Avino et al. JPC Lett. 2016] so as to provide embedded GW and Bethe-Salpeter QM/MM formalisms.
The Fiesta code is evolving towards a cubic scaling "space-time" approach following the methodology implemented at the RPA level (see Duchemin, Blase, J. Chem. Phys. 2019, 150, 174120).
Not started yet: periodic boundary conditions. As it stands, Fiesta allows only finite size systems (clusters, molecules) calculations.
Fiesta is not yet an open-source package freely distributed, even though it is shared by a few academic friendly partners.
The Fiesta code is written in F90/MPI. It was initially developed having in mind to run GW/BSE
calculations for moderate size systems (~100 atoms) on a standard desk PC (e.g. 8 cores/32 Go)
with a reasonable basis (e.g. 6-311G**). If you happen to have access to larger computers, the code
has been shown to sustain excellent scaling up to several thousand cores, allowing GW/BSE
calculations on several hundred atoms with e.g. aug-cc-pvtz quality basis.
Left image: Performance of one GW iteration for one or two fullerenes (RI-SVS metric,
DZP basis). The processor grid size varies from 16 to 61440 cores. 138 TeraFlops were reached
for the largest grid (European PRACE project, Ivan Duchemin, L_sim/INAC/CEA/Grenoble).