Dear users,<br><br>I was reading the documentation for the saw tooth like potential in pwscf.<br>However, I still did not understood it in full. Can anybody gives me more details about how the potential is implemented.<br>
I need to apply an electric field to a supercell of graphene and BN layers with both ends in vacuum. eamp = 0.002. <pre><font size="4"><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". <span style="background-color: rgb(255, 255, 0);">Important: the change of slope of this</span><br style="background-color: rgb(255, 255, 0);">
<span style="background-color: rgb(255, 255, 0);">potential must be located in the empty region, or else</span><br style="background-color: rgb(255, 255, 0);"><span style="background-color: rgb(255, 255, 0);">unphysical forces will result. Used only if tefield is .TRUE.</span></font><br>
</pre>The following keywords that need to be used are still not clear. In particular the position of the slope change.Does it need to be in vacuum?<br>emaxpos=0.5 (Position of the maximum of the sawlike potential along crystal <font size="2">axis "edir"</font>); eopreg=0.5 (Zone in the unit cell where the sawlike potential decreases.)?<br>
How the sawlike potential is modeled.<br>Thank you in advance for your help and comment.<br><br>Here is the system I want to apply electric field to it:<br><br>&control<br> title = 'GphBN'<br> calculation = 'relax'<br>
outdir = '/pwscf/pwscftemp/relaxGphBN'<br> prefix = 'GphBN'<br> pseudo_dir = /pwscf/input/Gph-BN/pseudopot-C-B-N'<br> verbosity = 'high'<br> tprnfor = .t.<br>
restart_mode = 'from_scratch'<br><br> / <br> &system<br> ibrav = 0,<br> celldm(1) = 1.8897261<br> nat = 72,<br> ntyp = 3,<br> ecutwfc = 25.0<br> ecutrho = 250.0<br>
occupations = 'smearing'<br> smearing = 'gaussian'<br> degauss = 0.003675<br><br> /<br> &electrons<br> mixing_beta = 0.05<br> diagonalization = 'david'<br><br> /<br>
&ions<br> trust_radius_ini = 0.10<br><br> /<br> &cell<br> <br>ATOMIC_SPECIES<br> B 10.81100 B.pz-vbc.UPF<br> C 12.01070 C.pz-vbc.UPF<br> N 14.00674 N.pz-vbc.UPF<br>ATOMIC_POSITIONS angstrom<br> B 0.710000 8.607000 11.306000<br>
B 0.710000 11.066000 11.306000<br> B 2.840000 7.377000 11.306000<br> B 0.710000 6.148000 11.306000<br> B 6.390000 6.148000 8.000000<br> B 6.390000 8.607000 8.000000<br>
B 6.390000 11.066000 8.000000<br> B 2.840000 9.836000 11.306000<br> B 7.100000 7.377000 11.306000<br> B 7.100000 9.836000 11.306000<br> B 7.100000 12.295000 11.306000<br>
B 4.970000 11.066000 11.306000<br> B 2.840000 12.295000 11.306000<br> B 4.970000 6.148000 11.306000<br> B 4.970000 8.607000 11.306000<br> B 4.260000 7.377000 8.000000<br>
B 2.130000 11.066000 8.000000<br> B 2.130000 8.607000 8.000000<br> B 2.130000 6.148000 8.000000<br> B 4.260000 9.836000 8.000000<br> B 0.000000 7.377000 8.000000<br>
B 0.000000 9.836000 8.000000<br> B 4.260000 12.295000 8.000000<br> B 0.000000 12.295000 8.000000<br> C 2.130000 6.148000 14.612000<br> C 2.130000 8.607000 14.612000<br>
C 2.130000 11.066000 14.612000<br> C 0.000000 12.295000 14.612000<br> C 0.000000 9.836000 14.612000<br> C 0.000000 7.377000 14.612000<br> C 0.710000 11.066000 14.612000<br>
C 0.710000 8.607000 14.612000<br> C 0.710000 6.148000 14.612000<br> C 6.390000 8.607000 14.612000<br> C 6.390000 6.148000 14.612000<br> C 4.970000 11.066000 14.612000<br>
C 6.390000 11.066000 14.612000<br> C 7.100000 12.295000 14.612000<br> C 7.100000 9.836000 14.612000<br> C 7.100000 7.377000 14.612000<br> C 4.970000 8.607000 14.612000<br>
C 2.840000 12.295000 14.612000<br> C 2.840000 9.836000 14.612000<br> C 2.840000 7.377000 14.612000<br> C 4.260000 7.377000 14.612000<br> C 4.970000 6.148000 14.612000<br>
C 4.260000 12.295000 14.612000<br> C 4.260000 9.836000 14.612000<br> N 0.000000 12.295000 11.306000<br> N 0.000000 9.836000 11.306000<br> N 0.000000 7.377000 11.306000<br>
N 2.840000 7.377000 8.000000<br> N 2.840000 9.836000 8.000000<br> N 2.840000 12.295000 8.000000<br> N 4.970000 11.066000 8.000000<br> N 4.970000 8.607000 8.000000<br>
N 4.970000 6.148000 8.000000<br> N 7.100000 12.295000 8.000000<br> N 7.100000 9.836000 8.000000<br> N 7.100000 7.377000 8.000000<br> N 0.710000 6.148000 8.000000<br>
N 0.710000 8.607000 8.000000<br> N 0.710000 11.066000 8.000000<br> N 6.390000 11.066000 11.306000<br> N 6.390000 8.607000 11.306000<br> N 6.390000 6.148000 11.306000<br>
N 2.130000 11.066000 11.306000<br> N 2.130000 8.607000 11.306000<br> N 2.130000 6.148000 11.306000<br> N 4.260000 12.295000 11.306000<br> N 4.260000 9.836000 11.306000<br>
N 4.260000 7.377000 11.306000<br>K_POINTS automatic<br>10 10 1 1 1 0 <br>CELL_PARAMETERS<br>8.51980 0.00000 0.00000<br>0.00000 7.37600 0.00000<br>0.00000 0.00000 22.6120<br><br><br>