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
			
			
						IBZ point weight in r_setup of 5.0
Moderators: Davide Sangalli, andrea marini, Daniele Varsano
- 
				
				Quxiao
 - Posts: 102
 - Joined: Fri Mar 26, 2021 11:27 am
 
IBZ point weight in r_setup of 5.0
You do not have the required permissions to view the files attached to this post.
			
									Quxiao in BIT,calculate the exciton
						- Daniele Varsano
 - Posts: 4283
 - Joined: Tue Mar 17, 2009 2:23 pm
 - Contact:
 
Re: IBZ point weight in r_setup of 5.0
Dear Quxiao, 
to analyse the k point grid, you can use the ypp post-processing. 
to generate the input file, uncommenting ListPts and ExpandPts
then run ypp.
Best,
Daniele
			
			
									
									to analyse the k point grid, you can use the ypp post-processing.
Code: Select all
>ypp -k kto 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/
						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/
- Davide Sangalli
 - Posts: 649
 - 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
Dear Quxiao,
we will restore this.
For now you can try using the following patch>
Best,
D.
			
			
									
									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
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/
						CNR-ISM, Division of Ultrafast Processes in Materials (FLASHit) and MaX Centre
https://sites.google.com/view/davidesangalli
http://www.max-centre.eu/