Dear all,
I want to read the many-body wave functions and charge density. Specifically, I want to get the quasiparticle corrected wave function and substitute the wave function files in the 'prefix.save' directory. Which file should I read and how to transfer it into Quantum-ESPRESSO format?
Best regards,
Jason
How to get the GW wavefunction
Moderators: Davide Sangalli, andrea.ferretti, myrta gruning, andrea marini, Daniele Varsano, Conor Hogan, Nicola Spallanzani
- jasonhan0710
- Posts: 58
- Joined: Wed Dec 23, 2020 6:48 am
- Location: China
How to get the GW wavefunction
Jason Han
Department of Physics
National University of Defense Technology
Hunan, China
Department of Physics
National University of Defense Technology
Hunan, China
- Daniele Varsano
- Posts: 3609
- Joined: Tue Mar 17, 2009 2:23 pm
- Contact:
Re: How to get the GW wavefunction
Dear Jason,
in order to have GW wavefunction a sort of self-consistency is needed, at the moment yambo calculate G0W0 or eigenvalue self-consistent GW which is a diagonal approximation so the wfs are not changed. There is the possibility to perform a coshex self-consistency, but it is not properly tested so it is not distributed with the GPL release.
Best,
Daniele
in order to have GW wavefunction a sort of self-consistency is needed, at the moment yambo calculate G0W0 or eigenvalue self-consistent GW which is a diagonal approximation so the wfs are not changed. There is the possibility to perform a coshex self-consistency, but it is not properly tested so it is not distributed with the GPL release.
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/
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/
- jasonhan0710
- Posts: 58
- Joined: Wed Dec 23, 2020 6:48 am
- Location: China
Re: How to get the GW wavefunction
Dear Daniele,
Thanks for your reply.
Also, I have another question. What is the output in "ndb.QP" file? How to read it and get the quasiparticle energy?
Best regards,
Jason
Thanks for your reply.
What kind of self-consistency do we need? Can we just get the self-energy with Yambo and put it into the SCF calculation, just like a DFT+U calculation? Do you think it's a reasonable approximation?in order to have GW wavefunction a sort of self-consistency is needed
If possible, would you please send me a copy of this version?There is the possibility to perform a coshex self-consistency
Also, I have another question. What is the output in "ndb.QP" file? How to read it and get the quasiparticle energy?
Best regards,
Jason
Jason Han
Department of Physics
National University of Defense Technology
Hunan, China
Department of Physics
National University of Defense Technology
Hunan, China
- Daniele Varsano
- Posts: 3609
- Joined: Tue Mar 17, 2009 2:23 pm
- Contact:
Re: How to get the GW wavefunction
Daer Jason,
this is implemented in the quest all code: https://www.questaal.org/get/
Anyway, are you sure you need to perform quasi-particle self consistent calculations?
Best,
Daniele
Self energy are non-local object and QE is not meant to do that. If you want to perform a quasi-particle self consistent calculation you may have a look to the works of of Mark van Schilfgaarde: M van Schilfgaarde, T Kotani, S Faleev Physical review letters 96 (22), 226402What kind of self-consistency do we need? Can we just get the self-energy with Yambo and put it into the SCF calculation, just like a DFT+U calculation? Do you think it's a reasonable approximation?
this is implemented in the quest all code: https://www.questaal.org/get/
Anyway, are you sure you need to perform quasi-particle self consistent calculations?
This is a develop version and it s not yet ready to be distributed.If possible, would you please send me a copy of this version?
This is the yambo database in netcdf format containing the quasi-particle corrections. QP energy can be read both in the report and in the output o.qp of the calculation.Also, I have another question. What is the output in "ndb.QP" file? How to read it and get the quasiparticle energy?
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/
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/