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
-
- 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: 4198
- 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 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/
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: 640
- 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/