you can draw a saw potential by hand, which will help you understand.<br> for you question , ［0，1］ of coz define the whole supercell under your examination.<br><br><br><br><div class="gmail_quote">On Thu, Jun 5, 2008 at 5:23 AM,  <<a href="mailto:Mohamed.Majdoub@mail.uh.edu">Mohamed.Majdoub@mail.uh.edu</a>> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Dear users,<br><br>I really find difficulties in understanding and ambiguities on how the saw like potential is implemented in pwscf. I need to apply an external electric filed to a slab structure in supercell (including vacuum region) in the z direction (edir =3) of amplitude eamp.According to the manual, besides setting tefiled and dipfield (dipole correction for slab) to be true, we need to specify the following variables in &system:<br>
------------------------------------------------------------------------<br>edir           INTEGER<br><br>               The direction of the electric field or dipole correction is <br>               parallel to the bg(:,edir) reciprocal lattice vector, so the<br>
               potential is constant in planes defined by FFT grid points;<br>               edir = 1, 2 or 3. Used only if tefield is .TRUE.<br><br>emaxpos        REAL ( default = 0.5D0 )<br><br>               Position of the maximum of the sawlike potential along crystal<br>
               axis "edir", within the  unit cell (see below), 0 < emaxpos < 1<br>               Used only if tefield is .TRUE.<br><br>eopreg         REAL( default = 0.1D0 )<br><br>               Zone in the unit cell where the sawlike potential decreases.<br>
               ( see below, 0 < eopreg < 1 ). Used only if tefield is .TRUE.<br><br>eamp           REAL ( default = 0.001 a.u. ) 0.000001 a.u<br><br>               Amplitude of the electric field (in a.u. = 51.44 10^10 V/m )<br>
               The sawlike potential increases with slope "eamp" in the<br>               region from (emaxpos+eopreg-1) to (emaxpos), then decreases<br>               to 0 until (emaxpos+eopreg), in units of the crystal<br>
               vector "edir". Used only if tefield is .TRUE.<br>---------------------------------------------------------------------<br>Since it is external electric field it shoud be only applied to the vacuum region in the supercell as it was pointed out before (<a href="http://www.democritos.it/pipermail/pw_forum/2004-July/001169.html" target="_blank">http://www.democritos.it/pipermail/pw_forum/2004-July/001169.html</a>). However, it is not clear if the 0< x<1 interval is the vacuum interval or the full length of the supercell (structure + vacuum). How the decrease in the sawlike potential is defined? Is it in the slab or the vaccum?<br>
Can any one familiar with this kind of calculation explain how it is implemented?<br><br>Thank you in advance!<br><br>Mohamed Sabri Majdoub<br>PhD mechanical engineering<br>University of Houston<br>Houston, Tx, USA<br>
<br>_______________________________________________<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" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br>Hai-Ping Lan <br>Department of Electronics ,<br>Peking University , Bejing, 100871<br><a href="mailto:lanhaiping@gmail.com">lanhaiping@gmail.com</a>, <a href="mailto:hplan@pku.edu.cn">hplan@pku.edu.cn</a>