IBZ point weight in r_setup of 5.0

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
Quxiao
Posts: 96
Joined: Fri Mar 26, 2021 11:27 am

IBZ point weight in r_setup of 5.0

Post by Quxiao » Fri Sep 15, 2023 2:44 am

Dear Daniele:
I wrote some code to process some data in yambo 4.4, but I found that in yambo 5.0 and above, I can't find the information about the weight of IBZ points in FBZ(full brillouin zone) in r_setup, how do I get this weight information? (Of course, I can get all the points through the fatlog, although their format has changed, but I still want to get the weight column).As shown in the figure, r_setup of 4.4 and 5.0, what should I do to obtain the weight column in the 5.0 output?
Thank you very much for your continuous help!
Best wishes!
Quxiao
You do not have the required permissions to view the files attached to this post.
Quxiao in BIT,calculate the exciton

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

Re: IBZ point weight in r_setup of 5.0

Post by Daniele Varsano » Mon Sep 18, 2023 8:50 am

Dear Quxiao,
to analyse the k point grid, you can use the ypp post-processing.

Code: Select all

>ypp -k k

to generate the input file, uncommenting ListPts and ExpandPts
then run ypp.

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: 614
Joined: Tue May 29, 2012 4:49 pm
Location: Via Salaria Km 29.3, CP 10, 00016, Monterotondo Stazione, Italy
Contact:

Re: IBZ point weight in r_setup of 5.0

Post by Davide Sangalli » Mon Sep 18, 2023 2:18 pm

Dear Quxiao,
we will restore this.

For now you can try using the following patch>

Code: Select all

diff --git a/src/communicate/REPORT_Energies.F b/src/communicate/REPORT_Energies.F
index 64c8af1f8f..cec09386e6 100644
--- a/src/communicate/REPORT_Energies.F
+++ b/src/communicate/REPORT_Energies.F
@@ -81,6 +81,7 @@ subroutine REPORT_Energies(E,k,Xk,n_bands_,k_range,what,l_kpoints,l_what)
      call msg('r',trim(dumb_ch1),k_al,'[rlu]')
      if (fat_log) call msg('r','      :',k%pt(ik,:),'[iku]')
      if (fat_log) call msg('r','      :',k_cc,'[cc]')
+     call msg('r','      : weight ',k%weights(ik))
      if (trim(what)/="Q") then
        call msg('r','      : wf components',wf_nc_k(ik) )
      endif
Best,
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/

Post Reply