## 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
NWChem packages.
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.

** Fiesta has now evolved into the beDeft (beyond-DFT) package, ** a major C++/MPI restructuration including
a newly developed analytic continuation approach [JCTC 16, 1742], cubic-scaling "space-time" RPA and GW [JCTC 17, 2383],
and a subsystem-based GW formalism.
On-going efforts to implement excited-states BSE gradients within the Z-vector formalism [JCP 159, 024116].
The Fiesta and beDeft codes have been so far shared with a few selected academic partners.
** A github open-access version is under active preparation. **

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