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