about "states"

Anything regarding the post-processing utility (e.g. excitonic wavefunction analysis) is dealt with in this forum.

Moderators: Davide Sangalli, andrea marini, Daniele Varsano

Post Reply
xixi
Posts: 30
Joined: Tue Mar 24, 2009 1:08 pm

about "states"

Post by xixi » Thu Oct 22, 2009 3:39 am

Dear developers,
what's the meaning of the two "1" in

States= "1 - 1"


if I want to plot the 11th states in o.exc_E_sorted file,
I use "11-1" or "1-11"?

Thanks !

Xi

Nanyang Technological University
Singapore

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

Re: about "states"

Post by andrea marini » Thu Oct 22, 2009 8:50 am

xixi wrote:Dear developers,
what's the meaning of the two "1" in

States= "1 - 1"


if I want to plot the 11th states in o.exc_E_sorted file,
I use "11-1" or "1-11"?
The variable States follows the syntax " 1st_state - last_state" with the index referring to the energy ordered states (as written in o.exc_E_sorted). So to plot the 11th state you need to set States= " 11 - 11 " . Be careful to leave spaces otherwise the parser will get confused.

Please also note that ypp will add to the state you asked for all other states that will be energy degenerate. This in order to respect spatial symmetries.
Andrea MARINI
Istituto di Struttura della Materia, CNR, (Italy)

xixi
Posts: 30
Joined: Tue Mar 24, 2009 1:08 pm

Re: about "states"

Post by xixi » Wed Jan 06, 2010 11:38 am

Dear andrea;

thanks for your reply,
but how can I identify the transition energy of E11, E12 in the o.exc_E_sorted file?
I didn't understand the meaning of the "1,2,3" ordered in the o.exc_E_sorted file.

Thanks!


Xi Zhu
School of Materials Science And Engineering,
Nanyang Technological University,
Blk N4.1, Nanyang Avenue,
Singapore 639798

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

Re: about "states"

Post by Daniele Varsano » Wed Jan 06, 2010 1:09 pm

Dear Xi Zhu,
could you be more explicit? What do you want to calculate/visualze?
In the o.exc_E_sorted you have the excitation energies sorted in energy with this structure:
# E [ev] Strength Index
#
6.495275 0.389558 1.000000
6.500834 0.102615 2.000000
6.505801 1.000000 3.000000
7.390 0.3106E-6 4.000

First column: Excitation energy
Second column : Oscillator strength
Third Column: Index

If you look the file sorted by intensity the index it is not monotone.

This is the index you should consider in the input file for plotting
an excitonic wave function referred to excitation you are interested in..

If you want to look at the KS transition that the more participate on each excitation consider
the ypp -e a option,
indicate the index of the excitation you want to look at (in you case 11 or 12 ) and you will have this information in the o.exc_weights_ file.

But I'm not sure on what you are meaning for E11 , E12
Hope it helps,

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/

xixi
Posts: 30
Joined: Tue Mar 24, 2009 1:08 pm

Re: about "states"

Post by xixi » Sun Jan 31, 2010 4:10 pm

Dear Varsano;
my problem is
I want to get the transition energy between valence band 1 (v1)and conduction band 1 (c1) ,
but I can't get this information from ypp,
Thanks!

Xi zhu

School of Materials Science And Engineering,
Nanyang Technological University,
Blk N4.1, Nanyang Avenue,
Singapore 63979

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

Re: about "states"

Post by Daniele Varsano » Sun Jan 31, 2010 4:35 pm

Dear Xi Zhu,

when you solve the Bethe-Salpeter equation you are not in a single-particle picture!!
The excitation energy you get comes from the diagonalization of an excitonic Hamiltonian
that is written in the valence-conduction wave function basis.
The wave function corresponding to the excitation can be written as a
\sum_ij A_ij phi_i phi_j, where i runs over occupied states and j over unoccupied states, so
you can look at the weight of each single-particle transition that "participate" in the
excitation energy (as I told you in the last post).
I want to get the transition energy between valence band 1 (v1)and conduction band 1 (c1)
If you want to look at the single-particle transition energy, you don't need ypp, and neither yambo,
you can look at the differences between the KS eigenvalues that comes from your DFT calculation.
This information is also reported in the BS spectrum file in the last columns, you will find
peaks in correspondence of the single-particle energy differences with the height that correspond
to their oscillator strength.

I hope this helps to clarify.

Cheers,

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