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 "edir", 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 "eamp" 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 "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 (http://www.democritos.it/pipermail/pw_forum/2004-July/001169.html). 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>