<html><head><base href="x-msg://10/"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div><br></div>You need to consider a slab with equally spaced atomic planes (at least in the inner part; you may allow to relax the outern planes, I think) and to choose an appropriate width for the window of the macroscopic average. If you do this&nbsp;correctly, you should not see these oscillations in the slab interior.<div><br></div><div><br></div><div>HTH</div><div><br></div><div>GS<br><div><br></div><div><br><div><div>Il giorno 29/apr/2011, alle ore 23.33, Izaak Williamson ha scritto:</div><br class="Apple-interchange-newline"><blockquote type="cite"><span class="Apple-style-span" style="border-collapse: separate; font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; font-size: medium; "><div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font-family: 'times new roman', 'new york', times, serif; font-size: 12pt; "><font size="2" style="font-family: arial, helvetica, sans-serif; ">Dear all,<br></font><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font-family: 'times new roman', 'new york', times, serif; font-size: 12pt; "><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font-family: 'times new roman', 'new york', times, serif; font-size: 12pt; "><div class="gmail_quote" style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font size="2"><font face="georgia,serif"><font size="2" style="font-family: arial, helvetica, sans-serif; "><br style="font-family: arial, helvetica, sans-serif; ">I am trying to calculate the work function for&nbsp;Cu(100) surface with and without bulk reference. The equations used are as follows:<br><br>work function w/o bulkref&nbsp; --&gt;&nbsp; wf1 = V<font size="1">vac</font><span class="Apple-converted-space">&nbsp;</span>- E<font size="1">F,slab</font><span class="Apple-converted-space">&nbsp;</span><br>work function with bulkref&nbsp; --&gt;&nbsp; wf2 = V<font size="1">vac</font><span class="Apple-converted-space">&nbsp;</span>- V<font size="1">slab,int</font><span class="Apple-converted-space">&nbsp;</span>+ V<font size="1">bulk,macro</font><span class="Apple-converted-space">&nbsp;</span>- E<font size="1">F,bulk</font><br><br>The problem is that the macroscopic average of VH+Vbare for&nbsp;the slab is oscillating as a function of z (see 05.png). I do not understand why the macroscopic average is oscillating. We want to&nbsp;see the convergence of the work function with the number of atomic layers in the slab. Due to the oscillation in macroscopic average,&nbsp;&nbsp;the work function dramatically changes when we go from 5 to 6 layers (4.86641 eV to 3.62728 eV, respectively). This makes it hard to see the convergence of wf with number of layers.<br><br>Would it be&nbsp;correct to use the average&nbsp;of the macroscopic average&nbsp;of the peaks and troughs in calculating wf2?&nbsp;Even if we calculate the average of upper oscillations (05.png), I don't know what to do with&nbsp;the trough that comes up just before the vacuum region?</font><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font-family: arial, helvetica, sans-serif; ">&nbsp;</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font size="4"><font size="2"><span style="font-family: arial, helvetica, sans-serif; ">Thank you&nbsp;for any help.</span></font></font><br></div></font></font></div><br>--<span class="Apple-converted-space">&nbsp;</span><br><font size="3">Izaak Williamson</font><br>Research Assistant<br>Physics Department<br>Boise State University<br><br></div></div></div><span>&lt;05.png&gt;</span><span>&lt;Al100.wf.png&gt;</span>_______________________________________________<br>Pw_forum mailing list<br><a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br><a href="http://www.democritos.it/mailman/listinfo/pw_forum">http://www.democritos.it/mailman/listinfo/pw_forum</a><br></div></span></blockquote></div><br><div>
<span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div><span class="Apple-style-span" style="color: rgb(126, 126, 126); font-size: 16px; font-style: italic; "><br class="Apple-interchange-newline">§ Gabriele Sclauzero,&nbsp;EPFL SB ITP CSEA</span></div><div><font class="Apple-style-span" color="#7E7E7E"><i>&nbsp;&nbsp; PH H2 462, Station 3,&nbsp;CH-1015 Lausanne</i></font></div></span>
</div>
<br></div></div></body></html>