you can draw a saw potential by hand, which will help you understand.<br>&nbsp;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,  &lt;<a href="mailto:Mohamed.Majdoub@mail.uh.edu">Mohamed.Majdoub@mail.uh.edu</a>&gt; 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 &amp;system:<br>
------------------------------------------------------------------------<br>edir&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; INTEGER<br><br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; The direction of the electric field or dipole correction is <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; parallel to the bg(:,edir) reciprocal lattice vector, so the<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; potential is constant in planes defined by FFT grid points;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; edir = 1, 2 or 3. Used only if tefield is .TRUE.<br><br>emaxpos&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; REAL ( default = 0.5D0 )<br><br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Position of the maximum of the sawlike potential along crystal<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; axis &quot;edir&quot;, within the&nbsp; unit cell (see below), 0 &lt; emaxpos &lt; 1<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Used only if tefield is .TRUE.<br><br>eopreg&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; REAL( default = 0.1D0 )<br><br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Zone in the unit cell where the sawlike potential decreases.<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ( see below, 0 &lt; eopreg &lt; 1 ). Used only if tefield is .TRUE.<br><br>eamp&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; REAL ( default = 0.001 a.u. ) 0.000001 a.u<br><br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Amplitude of the electric field (in a.u. = 51.44 10^10 V/m )<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; The sawlike potential increases with slope &quot;eamp&quot; in the<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; region from (emaxpos+eopreg-1) to (emaxpos), then decreases<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; to 0 until (emaxpos+eopreg), in units of the crystal<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; vector &quot;edir&quot;. 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&lt; x&lt;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>