BSE calculation- dense k-grid

Deals with issues related to computation of optical spectra in reciprocal space: RPA, TDDFT, local field effects.

Moderators: Davide Sangalli, andrea.ferretti, myrta gruning, andrea marini, Daniele Varsano, Conor Hogan

Post Reply
pyadav
Posts: 86
Joined: Thu Nov 26, 2020 2:56 pm
Contact:

BSE calculation- dense k-grid

Post by pyadav » Thu Jan 14, 2021 11:37 am

Dear Developers,

Can someone please help me with doing BSE calculation using a very dense k-grid, like 51x51x1 while keeping the DFT k-grid low (12x12x1)? There is a paper that has mentioned YAMBO for this purpose (PHYSICAL REVIEW B 88, 045412 (2013).

Wishes,
Pushpendra
Pushpendra Yadav
Ph.D. Research Scholar
Quantum Transport and Theory Group
Department of Physics
Indian Instittute of Technology Kanpur, India.

https://sites.google.com/site/amitkag1/

User avatar
Daniele Varsano
Posts: 3812
Joined: Tue Mar 17, 2009 2:23 pm
Contact:

Re: BSE calculation- dense k-grid

Post by Daniele Varsano » Thu Jan 14, 2021 11:47 am

Dear Pushpendra,
I may be wrong, but in the paper you mention a large grid (51x51x1) is used.
Can you explain better what you mean by keeping the DFT k-grid low (12x12x1)?

Best,
Daniele
Dr. Daniele Varsano
S3-CNR Institute of Nanoscience and MaX Center, Italy
MaX - Materials design at the Exascale
http://www.nano.cnr.it
http://www.max-centre.eu/

pyadav
Posts: 86
Joined: Thu Nov 26, 2020 2:56 pm
Contact:

Re: BSE calculation- dense k-grid

Post by pyadav » Thu Jan 14, 2021 12:38 pm

Dear Dr. Daniele,

Maybe I misunderstood the paper mentioned above but I wanted to understand that suppose we do a ground-state calculation (scf/nscf) for k-grid 12x12x1 and use these ground state wave-functions to do the G0W0 calculation. Now if I want to calculate BSE spectra for a higher k-grid (say 51x51x1). For this hogher k-grid BSE calculation do I need to go back to DFT calculation with higher k-grid or is it possible to implement higher k-grid in YAMBO?

Thank you!
Last edited by pyadav on Thu Jan 14, 2021 1:21 pm, edited 1 time in total.
Pushpendra Yadav
Ph.D. Research Scholar
Quantum Transport and Theory Group
Department of Physics
Indian Instittute of Technology Kanpur, India.

https://sites.google.com/site/amitkag1/

User avatar
Daniele Varsano
Posts: 3812
Joined: Tue Mar 17, 2009 2:23 pm
Contact:

Re: BSE calculation- dense k-grid

Post by Daniele Varsano » Thu Jan 14, 2021 12:39 pm

Dear Pushpendra,
yes, you need to recalculate the electronic structure for a denser grid, as you need it to build your BSE matrix.

Best,
Daniele
Dr. Daniele Varsano
S3-CNR Institute of Nanoscience and MaX Center, Italy
MaX - Materials design at the Exascale
http://www.nano.cnr.it
http://www.max-centre.eu/

pyadav
Posts: 86
Joined: Thu Nov 26, 2020 2:56 pm
Contact:

Re: BSE calculation- dense k-grid

Post by pyadav » Thu Jan 14, 2021 1:24 pm

Thank you, Dr. Daniele, for the clarification.
Pushpendra Yadav
Ph.D. Research Scholar
Quantum Transport and Theory Group
Department of Physics
Indian Instittute of Technology Kanpur, India.

https://sites.google.com/site/amitkag1/

Zafer Kandemir
Posts: 40
Joined: Wed May 06, 2020 9:43 am
Contact:

Re: BSE calculation- dense k-grid

Post by Zafer Kandemir » Mon Sep 05, 2022 2:21 pm

pyadav wrote: Thu Jan 14, 2021 12:38 pm Dear Dr. Daniele,

Maybe I misunderstood the paper mentioned above but I wanted to understand that suppose we do a ground-state calculation (scf/nscf) for k-grid 12x12x1 and use these ground state wave-functions to do the G0W0 calculation. Now if I want to calculate BSE spectra for a higher k-grid (say 51x51x1). For this hogher k-grid BSE calculation do I need to go back to DFT calculation with higher k-grid or is it possible to implement higher k-grid in YAMBO?

Thank you!
Dear Developers,

Is this calculation possible for the new version of the Yambo program?

Best regards,
Zafer
Dr. Zafer Kandemir
Postdoctoral Researcher
Department of Mechatranics Engineering, Sabanci University
Istanbul, Turkey

User avatar
Daniele Varsano
Posts: 3812
Joined: Tue Mar 17, 2009 2:23 pm
Contact:

Re: BSE calculation- dense k-grid

Post by Daniele Varsano » Fri Sep 09, 2022 9:05 am

Dear Zafer,

can you be more explicit on what kind of calculation you are referring?
In the previous post, the user was asking about BSE calculations using a dense grid without calculating the ground state on such grid. This is not possible.
What is possible is to speed up the calculation using a double gri method as explained here:
http://www.attaccalite.com/speed-up-die ... ith-yambo/

or via Haydock method as described here:

https://www.frontiersin.org/articles/10 ... 63946/full


Best,
Daniele
Dr. Daniele Varsano
S3-CNR Institute of Nanoscience and MaX Center, Italy
MaX - Materials design at the Exascale
http://www.nano.cnr.it
http://www.max-centre.eu/

Zafer Kandemir
Posts: 40
Joined: Wed May 06, 2020 9:43 am
Contact:

Re: BSE calculation- dense k-grid

Post by Zafer Kandemir » Mon Sep 12, 2022 8:28 am

Dear Daniele,

Thank you for replying, my question is like this: For k-grid 18x18x8, a ground state calculation of a bulk structure was made, and then GW calculation was made and ndb.QP was obtained. Now I want to calculate BSE spectra for k-grid 26x26x10. This calculation is not possible, ok.
If I want to use the double grid method, I should choose k-grid:18x18x8 for scf calc. and k-grid:26x26x10 for ncsf calc. Is it possible to calculate BSE spectra by using the ndb.QP file obtained from the previous GW calculation?

Best regards,

Zafer
Dr. Zafer Kandemir
Postdoctoral Researcher
Department of Mechatranics Engineering, Sabanci University
Istanbul, Turkey

User avatar
Daniele Varsano
Posts: 3812
Joined: Tue Mar 17, 2009 2:23 pm
Contact:

Re: BSE calculation- dense k-grid

Post by Daniele Varsano » Mon Sep 12, 2022 10:35 am

Dear Zafer,

It is not possible to use the ndb.QP directly. You can use a scissor/stretching operator with parameters extracted from the previous calculation and apply them to the calculation with denser grid.

Best,
Daniele
Dr. Daniele Varsano
S3-CNR Institute of Nanoscience and MaX Center, Italy
MaX - Materials design at the Exascale
http://www.nano.cnr.it
http://www.max-centre.eu/

Post Reply