How does the scissor correction work

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
Hsiao-Yi-Chen
Posts: 4
Joined: Sun Mar 29, 2020 3:30 am

How does the scissor correction work

Post by Hsiao-Yi-Chen » Tue Mar 15, 2022 10:13 am

Dear all,

I am now trying to compute linear response on semi-metal.
And I want to use the XfnQP_E to do some corrections on the band edge.

Code: Select all

% XfnQP_E
 aaaaa | bbbbb | ccccc|
%
I know for semiconductor/ insulator, aaaaa will open the gap while bbbbb and ccccc will stretch the conduction/valance band respectively with the band edge maintained.

However, what if I work on the semi-metal?
There is no well-defined conduction/valance band.
Further stretching and moving the metallic bands may change the occupation for each k-point.

So how does yambo do scissor correction for semi-metal?
And is there a way for me to find the number of the correction?
or a way to plot the scissor corrected band out?

Thank you.

Best
Hsiao-Yi Chen

RIKEN Center for Emergent Matter Science (CEMS),
Wako, Saitama, Japan

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

Re: How does the scissor correction work

Post by Daniele Varsano » Tue Mar 15, 2022 11:35 am

Dear Hsiao-Yi Chen,

that variable is meant for semiconductors, you can give it a try using it for semimetal but I'm afraid that it can produce wrong results.

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/

User avatar
Davide Sangalli
Posts: 610
Joined: Tue May 29, 2012 4:49 pm
Location: Via Salaria Km 29.3, CP 10, 00016, Monterotondo Stazione, Italy
Contact:

Re: How does the scissor correction work

Post by Davide Sangalli » Tue Mar 15, 2022 11:51 am

I think you can anyway try to use it. Especially the stretching parameters.
They are defined with respect to the conduction band minimum and valence band maximum which are reported in the report file

For example here

Code: Select all

  [X] Conduction Band Min                           :  1.776788 [eV]
  [X] Valence Band Max                              :  0.000000 [eV]
In a semi-metal I would expect both to be the same number.
The occupations will be influenced as you say.

You can plot the QP corrected bands using ypp and adding the QP correction variables in input

Code: Select all

ypp -electron b -V qp -F bands.in
and changing

Code: Select all

% GfnQP_E
 0.000000 | 1.000000 | 1.000000 |        # [EXTQP G] E parameters  (c/v) eV|adim|adim
%
in input

It is not tested, but it may work just fine.

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/

Hsiao-Yi-Chen
Posts: 4
Joined: Sun Mar 29, 2020 3:30 am

Re: How does the scissor correction work

Post by Hsiao-Yi-Chen » Wed Mar 16, 2022 2:54 am

Ciao Daniele, Ciao Davide,

I give a test and plot them.
The operation looks based on the Fermi surface.
So it will cut the metallic bands along Fermi level for rigidly opening a gap.
For the stretching, it also treats the points above fermi as conduction while below fermi as a valance.
So if you use different stretching values to v/c, the result will be very weird.
But if the parameters are the same, it would work pretty well.

BTW,
I remember you can generate QP manually by ypp -q g.
But it was reported with bugs a long time ago.
Is that fixed in some version?


Best
Hsiao-Yi Chen

RIKEN Center for Emergent Matter Science (CEMS),
Wako, Saitama, Japan

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

Re: How does the scissor correction work

Post by Daniele Varsano » Wed Mar 16, 2022 3:54 pm

Dear Hsiao,

the QP database generator (ypp -q g) it is not much tested and used, anyway it should work.
You can give a try and report if something does not work.
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