Tamm-Dancoff approximation

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
zchen
Posts: 41
Joined: Thu Apr 01, 2010 2:23 pm

Tamm-Dancoff approximation

Post by zchen » Wed Apr 07, 2010 2:30 pm

Hi, everyone:
In GPL version of Yambo Bethe-Sapeter equation can only be solved by Tammm-Dancoff approximation. My calculation of absorption for ce@C82 disagrees with the experiment. How can I solve BSE without TDA?
If I am only interested in the excitation from the bands near HOMO and final states are near LOMO what should I do?
Third question is how can I calculate the photoionization (electron goes to the continuum state) using Yambo?

Zhifan Chen
Dept of Physics
Clark Atlanta University
Atlanta GA 30314

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

Re: Tamm-Dancoff approximation

Post by Daniele Varsano » Wed Apr 07, 2010 2:48 pm

Dear Zchen,
the solution of the Bethe Salpeter equation beyond the Tamm-Dancoff approximation, unfortunately even if implemented
in Yambo, at the moment it is still not present in the GPL version. We hope to include it soon as GPL
in the next release of the code.

Regards,

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/

User avatar
andrea marini
Posts: 325
Joined: Mon Mar 16, 2009 4:27 pm
Contact:

Re: Tamm-Dancoff approximation

Post by andrea marini » Wed Apr 07, 2010 2:55 pm

zchen wrote: My calculation of absorption for ce@C82 disagrees with the experiment. How can I solve BSE without TDA?
Dear Zhifan Chen,

why do you think that the problem is related to the TDA ? What do you mean as disagreement ?

zchen wrote:
Third question is how can I calculate the photoionization (electron goes to the continuum state) using Yambo?
You cannot calculate directly photoemission (PE) as you do not know the shape of the outgoing electron wavefunction. With Yambo, however, you can calculate the peaks of the PE spectra and you can try to code yourself an expression for the intensities assuming some expression for the outgoing states (plane waves, for example). The Peaks of the PE are simply the quasiparticle energies.
Andrea MARINI
Istituto di Struttura della Materia, CNR, (Italy)

zchen
Posts: 41
Joined: Thu Apr 01, 2010 2:23 pm

Re: Tamm-Dancoff approximation

Post by zchen » Thu Apr 08, 2010 2:45 pm

Thanks. If I am only interested in the excitation from near HOMO states what should I do? This will save lots of computing time.

User avatar
andrea marini
Posts: 325
Joined: Mon Mar 16, 2009 4:27 pm
Contact:

Re: Tamm-Dancoff approximation

Post by andrea marini » Thu Apr 08, 2010 2:49 pm

zchen wrote:Thanks. If I am only interested in the excitation from near HOMO states what should I do? This will save lots of computing time.
What observable are you referring to ? Photemission, absorption ? And what are you exactly asking help for ? Please provide as more informations as possible so that we can try to help you.
Andrea MARINI
Istituto di Struttura della Materia, CNR, (Italy)

zchen
Posts: 41
Joined: Thu Apr 01, 2010 2:23 pm

Re: Tamm-Dancoff approximation

Post by zchen » Fri Apr 09, 2010 4:31 pm

I try to use Bethe-Sapeter equation to calculate the absorption of the Ce@C82. The calculation needs to include 177 occupied bands and some empty bands. Acturely I am only interested in the excitation between bands 175-177 to 178-180 etc.

Thanks

User avatar
andrea marini
Posts: 325
Joined: Mon Mar 16, 2009 4:27 pm
Contact:

Re: Tamm-Dancoff approximation

Post by andrea marini » Fri Apr 09, 2010 4:57 pm

zchen wrote:I try to use Bethe-Sapeter equation to calculate the absorption of the Ce@C82. The calculation needs to include 177 occupied bands and some empty bands. Acturely I am only interested in the excitation between bands 175-177 to 178-180 etc.
Dear zchen,

please edit your signature to include your complete affiliation and name so that it will appear in all your posts.

Regarding your question, can you post the input file of your calculation and any report file you already have so that we can understand what exactly you're running ? Without some more informations it is really difficult for us to give you any help.

From the report file we will se home electron, G-vectors, k-points you are using. The BSE hamiltonian size depends on many factors.
Andrea MARINI
Istituto di Struttura della Materia, CNR, (Italy)

Post Reply