<div dir="ltr">Hello everyone, <div>Firstly I am grateful to the Quantum Espresso Team for building this free software and maintaining this forum, all without any costs. You have truly taken a step forward for growing education without any possible fees.</div><div><br></div><div><u><b>Coming to my questio</b></u><b><u>n</u></b>, I am trying to calculate the vacuum potential of a hexagonal unit cell of graphene of  two atoms. I proceed in the following way:</div><div><div><br></div><div><b><u>What I did:</u></b></div><div>1) scf with pw.x</div><div>2) <a href="http://pp.in">pp.in</a> with pp.x</div><div>3) using a code to find the planar average electrostatic potential in the z (vacuum) direction. As I could not figure out how to use the average.x tool.</div><div>4) identify the vacuum potential, which in my case is not zero (should be zero as written in literature <a href="https://doi.org/10.1016/j.ijhydene.2022.09.249">https://doi.org/10.1016/j.ijhydene.2022.09.249</a>).</div><div><br></div><div>all the files are attached. </div><div><br></div><div><b><u>My question is that:</u></b></div><div>1) is my method correct ?</div><div>2) if yes then why is the vacuum potential non zero? What can I do to make it zero ?</div><div>3) I am trying to find the vacuum potential so that I can calculate work function and vbm cbm values with respect to vacuum potential for photocatalytic application.</div><div><br></div><div>Kindly help me. I have been stuck at this for a long time.</div><div><br></div><div><br></div><div><a href="https://drive.google.com/drive/folders/1knK2PNlKFH26ri84sLJGKM2LVRN789BC?usp=sharing">google drive link</a> . I have attached all the files in the google drive link. If you can't read it I'll send in whichever manner needed.<br></div><div><br></div><div><br></div><div>Ashley Cooper</div><div>Phd</div><div>Assam University Silchar</div><div>Assam, India</div><div><br></div><div><br></div></div></div>