Hi everyone,
I did two BSE calculations using a different number of RL's in the screening calculation.
To my surprise the FFT grid changed when I increased the number of RL for the screening.
I thought the FFT grid size is solely related the the number of plane-waves of the wave function.
Can you guide me to the subroutine where yambo decides the FFT grid size and could you explain why it depends on the screening G's?
Thanks a lot in advance.
Best,
Franz
FFT grid in BSE calculation
Moderators: Davide Sangalli, andrea.ferretti, myrta gruning, andrea marini, Daniele Varsano, Conor Hogan, Nicola Spallanzani
-
- Posts: 43
- Joined: Wed Jul 20, 2022 9:36 am
FFT grid in BSE calculation
Franz Fischer
PhD student / IMPRS-UFAST fellow
Institute of Physical Chemistry
University of Hamburg
PhD student / IMPRS-UFAST fellow
Institute of Physical Chemistry
University of Hamburg
- Davide Sangalli
- Posts: 621
- Joined: Tue May 29, 2012 4:49 pm
- Location: Via Salaria Km 29.3, CP 10, 00016, Monterotondo Stazione, Italy
- Contact:
Re: FFT grid in BSE calculation
Ciao,
the FFT grid size in real space is related to the cut-off in G space.
Setting a cutoff of the screening selects a G-max on which W_{G,G'} is defined.
Accordingly the sum over G done in this step runs up to G-max
with (below FT is the Fourier Transform)
(e.g. he code prefers to do a Fourier transform rather than a convolution in G-space)
This means that it first has to load the WFs in real space.
The FFT grid used in real space is related to the needed G-max.
Best,
D.
the FFT grid size in real space is related to the cut-off in G space.
Setting a cutoff of the screening selects a G-max on which W_{G,G'} is defined.
Accordingly the sum over G done in this step runs up to G-max
Code: Select all
< cvk | W | c'v'k'> = sum_{G,G'} (psi_{ck}*psi_{vk})(G) W_{G,G'} (psi_{ck}*psi_{vk})(G')
Code: Select all
(psi_{ck}*psi_{vk})(G)= FT [ psi_{ck}(r)*psi_{vk}(r) ]
This means that it first has to load the WFs in real space.
The FFT grid used in real space is related to the needed G-max.
Best,
D.
Davide Sangalli, PhD
CNR-ISM, Division of Ultrafast Processes in Materials (FLASHit) and MaX Centre
https://sites.google.com/view/davidesangalli
http://www.max-centre.eu/
CNR-ISM, Division of Ultrafast Processes in Materials (FLASHit) and MaX Centre
https://sites.google.com/view/davidesangalli
http://www.max-centre.eu/