<div dir="ltr"><div>Thank You for the response. I have already optimised the structure and removed top layer, and running for only the bottom part . So all the remaining atoms needs to be frozen or only ones near the interface will do the job ?</div><div><br></div><div>Its a 3d structure having vacuum along z-axix.</div><div><br></div><div>Thank You.</div><div><br></div><div>Satya</div><div>IIT Madras<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, May 31, 2023 at 12:26 PM Giovanni Cantele <<a href="mailto:giovanni.cantele@spin.cnr.it">giovanni.cantele@spin.cnr.it</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">If you want to visualize the charge density variations induced by the formation of the interface, once you have the optimized structure of the heterostructure and<div>its charge density rho_tot, you have to make, at frozen atomic positions two more scf calculations, obtained from the former by removing one layer at a time, thus obtaining</div><div>let's say rho_1 and rho_2. What you're looking for is then easily computed from rho_tot - rho_1 - rho_2.</div><div><br></div><div>In this way, you just visualize charge density reconstruction induced by the interface, getting rid of contributions related to geometry changes induced by the presence of</div><div>the interface itself. This is because if you change the structure, the result would be harder visualized.</div><div><br></div><div>Giovanni</div><div><br></div><div><div><div><div dir="ltr" class="gmail_signature"><div dir="ltr"><span style="color:rgb(0,0,0)">-- <br><br>Giovanni Cantele, PhD<br>CNR-SPIN<br>c/o Dipartimento di Fisica<br>Universita' di Napoli "Federico II"<br>Complesso Universitario M. S. Angelo - Ed. 6<br>Via Cintia, I-80126, Napoli, Italy<br><a href="mailto:giovanni.cantele@spin.cnr.it" style="color:rgb(17,85,204)" target="_blank">e-mail: giovanni.cantele@spin.cnr.it</a><br>Phone: +39 081 676910<br>Skype contact: giocan74<br><br>ResearcherID: <a href="http://www.researcherid.com/rid/A-1951-2009" style="color:rgb(17,85,204)" target="_blank">http://www.researcherid.com/rid/A-1951-2009</a><br>Web page: </span><a href="https://sites.google.com/view/giovanni-cantele/home" style="color:rgb(17,85,204)" target="_blank">https://sites.google.com/view/giovanni-cantele/home</a><br></div></div></div><br></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Il giorno mar 30 mag 2023 alle ore 22:19 Satyasiban Dash ph19d005 <<a href="mailto:ph19d005@smail.iitm.ac.in" target="_blank">ph19d005@smail.iitm.ac.in</a>> ha scritto:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="auto">Dear users,<div dir="auto"><br></div><div dir="auto">I am planning to calculate charge density difference at interface of a heterostructure. So I underwent relaxation of cell ,then scf of hetero structure, so here are my small doubts</div><div dir="auto"><br></div><div dir="auto">1. Do I need to delete one layer and calculate scf with same lattice cell and atomic position and do the same for other layer later ?</div><div dir="auto"><br></div><div dir="auto">2. Instead of above approach do I need to make a complete different structure of upper or lower layer and take scf then rho ?</div><div dir="auto"><br></div><div dir="auto"><br></div><div dir="auto">Thank You.</div><div dir="auto"><br></div><div dir="auto"><br></div><div dir="auto">Satya</div><div dir="auto">PhD Scholar</div><div dir="auto">IIT Madras</div></div>
_______________________________________________<br>
The Quantum ESPRESSO community stands by the Ukrainian<br>
people and expresses its concerns about the devastating<br>
effects that the Russian military offensive has on their<br>
country and on the free and peaceful scientific, cultural,<br>
and economic cooperation amongst peoples<br>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer" target="_blank">www.max-centre.eu</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div>
_______________________________________________<br>
The Quantum ESPRESSO community stands by the Ukrainian<br>
people and expresses its concerns about the devastating<br>
effects that the Russian military offensive has on their<br>
country and on the free and peaceful scientific, cultural,<br>
and economic cooperation amongst peoples<br>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer" target="_blank">www.max-centre.eu</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div>