<div dir="ltr"><div dir="ltr"><div dir="ltr"><span style="font-size:16px;color:rgb(29,28,29);font-family:Slack-Lato,appleLogo,sans-serif">Dear QE Community and Developers</span><span class="gmail-c-mrkdwn__br" style="font-size:16px;box-sizing:inherit;display:block;height:8px;color:rgb(29,28,29);font-family:Slack-Lato,appleLogo,sans-serif"></span><span style="font-size:16px;color:rgb(29,28,29);font-family:Slack-Lato,appleLogo,sans-serif">My colleagues and I are currently trying to find a way to output 'free electron wavefunctions' from QE. More specifically, we wish to simulate a 1D array of silicon atoms. We were thinking to do this by deleting the cell-dependent part of the Bloch wavefunctions in the directions transverse to the length of the 1D array keeping only the kpoint-dependent e^ikx term in those directions. What would be a reasonable way to accomplish this in QE? We are totally fine with modifying the code. Our end goal is to have the modified wavefunctions in the same format as QE normally writes production wavefunctions to files.</span><span class="gmail-c-mrkdwn__br" style="font-size:16px;box-sizing:inherit;display:block;height:8px;color:rgb(29,28,29);font-family:Slack-Lato,appleLogo,sans-serif"></span><span style="font-size:16px;color:rgb(29,28,29);font-family:Slack-Lato,appleLogo,sans-serif">Best,</span><br></div><div><span style="font-size:16px;color:rgb(29,28,29);font-family:Slack-Lato,appleLogo,sans-serif">Qinxue Cao</span></div><div><font face="arial, sans-serif"><span style="color:rgb(82,88,98);background-color:rgb(250,250,250)">Department of Electrical and Computer Engineering</span><br></font></div><div><span style="color:rgb(82,88,98);background-color:rgb(250,250,250)"><font face="arial, sans-serif">Boston University</font></span></div></div></div>