<div dir="ltr"><div><div><div><div><div><div>Dear Prof. Marzari,<br></div>                           thank you very much for all your suggestions and concerns. I'm really overwhelmed ! <br>                           I've looked at work function calculation (thanks to the links you provided). Now, I have a monolayer of crystalline non-polar semiconductor material and I can calculate the planar and macroscopic averages for unstrained and strained cases. But I'm still a little bit confused about fixing the reference. Should I take macroscopic average of Hartree potential at vacuum as the reference ?<br>
<br></div>I've thought about the following steps (in light of your previous explanation):<br><br></div>1. Do a bulk calculation and find out the macroscopic average <br></div>2. Do the strained monolayer and see where the macroscopic average sits with respect to vacuum level.<br>
</div>   But then I don't understand why do we need a bulk calculation if I can get the vacuum level (in sigle-layer) and can align it for all strain cases and refer the HOMO-LUMO directly to that !<br><br></div>                              I know I have a lot of doubts and I truly appreciate your patience and all your helpful suggestions.<br>
<div><div><div><div><div><br></div><div>Thanks,<br></div><div>Rajdeep Banerjee<br></div><div>JNCASR<br>Bangalore, India<br></div></div></div></div></div></div>