<div dir="ltr">Dear Prof. Blaha,<div><br></div><div>Thanks so much for your suggestion.</div><div><br></div><div>I tried bulk Al supercell with 39ML without vacuum with the same k-mesh as used in 39 ML thin film supercell. In fact I get the same results for plasma frequency and dielectric function as those from Al unit cell with only one Al atom. I think the k-mesh of 61x61x1 I used in my calculation is dense enough to give a precise result. </div>
<div><br></div><div>The main difference for the dielectric function between thin film geometry and bulk geometry is at the low energy range (<1.2eV). I researched some paper for studying the anisotropic surface reflectance in semiconductor surface, say, GaAs(110). Even 15 atomic layers are used in the LDA calculation but still some difference around the band gap regime for the dielectric function is found between surface calculation and bulk calculation. I think the difference I encountered for teh dielectric function between slab Al(111) and bulk Al might be similar to the case in semiconductor system. However, from the physical point of view, it's hard to understand why there is still appreciable difference out there even though very thick film is used. Physically the dielectric function of the very thick slab should converge to that in the bulk counterpart.</div>
<div><br></div><div>Thank you so much for sharing any understanding about this,</div><div><br></div><div>Wenmei</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">2013/11/24 Peter Blaha <span dir="ltr"><<a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">pblaha@theochem.tuwien.ac.at</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">As you probably know, the dielectric function of Al converges VERY slowly<br>
with respect to the k-mesh.<br>
<br>
When you do slab calculations, you include the surface effect, but you also replace<br>
the periodicity in k-z (and thus the k-mesh in k-z) to a backfoldung according to<br>
your slab. Even a 39 ML slab corresponds probably not to a very large k-z mesh and<br>
in addition the integration over k-z is limited to a "root"-sampling instead of the<br>
tetrahedron method. I could even imagine large numerical problems in this 2-D integration<br>
using a 3-D algorithm in joint due to large degeneracy of the tetrahedra.<br>
<br>
At least you could differentiate between "integration problems" and surface effects<br>
by using a 39-layer bulk structure (i.e. remove the vacuum in your supercell, so that<br>
you get 3D-Al again, but restrict yourself to 1-k point in k-z) and compare the<br>
resulting eps to bulk Al (with 1 atom/cell and good k-meshes.<br>
<br>
Am 23.11.2013 16:54, schrieb phlhj phlhj:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5">
Dear all,<br>
<br>
I was trying to calculate the optical properties of Al(111) slab. For the bulk FCC Al, I can reproduce the dielectric functions and plasma frequency very precisely reported<br>
in literature before. However, I did find some difference between the slab dielectric functions and the corresponding bulk values.<br>
<br>
Especially even though I used very thick slab, say 39MLs, in the low photo energy range (<1eV), the imaginary part is much larger than the bulk. I doubt this may be related<br>
to the band-folding and symmetry reduction in the direction normal to the surface.<br>
<br>
Also, I found the plasma frequency of the slab is smaller than the bulk plasma frequency.<br>
<br>
Mathematically, this behavior of the imaginary parts of the interband and intraband transitions contributions seems to be able to be understood from the f-sum rule.<br>
<br>
1) However, physically it's hard to believe, because when the slab thickness is very thick for example the 39MLs used in my test calculation, the slab wavefunctions should<br>
be very very close to the bulk wavefunction except in the very thin slab/vacuum interface region. This should give us the dielectric functions for the slab which are very<br>
very close to the bulk values. This argument should be also true for the slab plasma frequency.<br>
<br>
2) If the different values are because of the surface slab structure we used in the calculation, which indeed breaks the translational symmetry in the normal direction.<br>
Then the question is that in real experiment because the sample always is finite with the boundary surface, how can we get the dielectric information really for the ideal<br>
bulk rather than the slab similar as that mentioned above. Or in calculating dielectric function, when should we use bulk geometry? when should we use slab geometry?<br>
<br>
Thanks a lot for any idea.<br>
<br>
Wenmei<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br></div></div>
______________________________<u></u>_________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.<u></u>at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">http://zeus.theochem.tuwien.<u></u>ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at: <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" target="_blank">http://www.mail-archive.com/<u></u>wien@zeus.theochem.tuwien.ac.<u></u>at/index.html</a><br>
<br>
</blockquote>
<br>
-- <br>
------------------------------<u></u>-----------<br>
Peter Blaha<br>
Inst. Materials Chemistry, TU Vienna<br>
Getreidemarkt 9, A-1060 Vienna, Austria<br>
Tel: <a href="tel:%2B43-1-5880115671" value="+4315880115671" target="_blank">+43-1-5880115671</a><br>
Fax: <a href="tel:%2B43-1-5880115698" value="+4315880115698" target="_blank">+43-1-5880115698</a><br>
email: <a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">pblaha@theochem.tuwien.ac.at</a><br>
------------------------------<u></u>-----------<br>
______________________________<u></u>_________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.<u></u>at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">http://zeus.theochem.tuwien.<u></u>ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at: <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" target="_blank">http://www.mail-archive.com/<u></u>wien@zeus.theochem.tuwien.ac.<u></u>at/index.html</a><br>
</blockquote></div><br></div>