<div dir="ltr">Hello all,<br><br>I am trying to calculate the momentum matrix elements for my material with HSE and HSE+U etc, which is not supported in QE. I already calculated the matrix elements with PBE, but I want to see if the optical selection rules change when the hybrid functional is turned on. <br><br>I set out to do this with the help of BerkeleyGW- They outline a relatively simple way that takes into account only the wavefunctions in the inner product (i.e. the dipole matrix element). This 'finite difference' method is mentioned in the BerkeleyGW manual ( <a href="https://arxiv.org/pdf/1111.4429.pdf">https://arxiv.org/pdf/1111.4429.pdf</a>, equation 46).<br><br>I tried porting my wavefunctions from QE to GW using pw2gw.x but I am encountering some errors. Please bear with me while I outline the procedure I used, and the doubts I have.<br><br>1) Use pp.x to save the wavefunction in a readable format. Here is my pp.x input. To start with, I just want to calculate the matrix element between the band edges at gamma point. The input below generates the wavefunction at VBM GAMMA. <br><br><div>&INPUTPP<br>prefix='CrPs4SCF'<br>outdir='./'<br>filplot='wfn'<br>plot_num=0<br>kpoint=1<br>kband=92<br>/<br><br><br>&PLOT<br>output_format=5<br>iflag=3<br>filepp='./plots'<br>fileout='./exportgcf'<br>/<br><br>This seems to generate the wavefunction just fine. Did not encounter an issue here.<br><br>2) To convert this wavefunction to a format that BGW can read, I used pw2gw.x. Here's the input: <br><br>&input_pw2bgw<br>   prefix = 'CrPs4SCF'<br>   outdir = './'<br>   real_or_complex = 2<br>   wfng_flag = .true.<br>   wfng_file = 'realwfnGW'<br>   wfng_kgrid = .false.<br>   wfng_nk1 = 4<br>   wfng_nk2 = 6<br>   wfng_nk3 = 1<br>/<br><br>This seems to lead to an error :   Error in routine pw2gw (1):     reading inputpp namelist.<br><br>Would someone be able to tell me where I might have gone wrong? I am also confused as to what input the pw2gw file reads - Is it the wavefunction directly from my first step i.e. the output 'plots' from my pp.x calculation? or should I copy paste the whole SCF directory that includes the .xml file and charge density data?<br><br>Best,<br>Shivesh Sivakumar<br>University of Washington-Seattle</div></div>