<div dir="ltr"><div><div><div><div><div><div><div><div><div><div><div><div>Dear all, <br></div>          
 I am doing relax calculations on monolayer of graphene-like material 
(2-D) under the presence of static homogeneous electric field in 
z-direction (i.e. perpendicular to the layer). I am using a electric 
field range of 0.1-0.5 V/angstrom. The problem is that the bfgs steps 
are not converging (i.e. CASE: energy  _new > energy  _old) after 4-5
 steps, if I apply a field value > 0.2 V/angstrom (0.0055 Ry.a.u). 
Following is the relevant part of the code:<br>
...<br>/<br> &system<br>    ibrav = 4,<br>    celldm(1) = 7.5284, celldm(3) = 5.020259255 ,<br>    nat = 10, ntyp = 3,<br>    ecutwfc = 40, ecutrho = 400,<br>    london = .true.,<br>/<br> &electrons<br>    diagonalization = 'david', mixing_mode = 'local-TF',<br>


    mixing_beta = 0.4, conv_thr = 1.0d-8, electron_maxstep = 125,<br>    efield_cart(1) = 0.0,<br>    efield_cart(2) = 0.0,<br>    efield_cart(3) = .008<br>/<br> &ions<br>...<br>K_POINTS automatic<br> 12 12 1 0 0 0<br>


<br></div>            Is it because the applied field is too high? I 
have checked the paper Souza et al., PRL, 89, 117602 (2002) where they 
have put a condition :<br><br></div>            e|E.ai| < e|Ec.ai| ~ Egap/Ni<br>
</div><div>(E = magnitude of the electric field, ai = lattice vector, Ni = number of k-points in i direction, Egap = Band gap)<br></div>Is this condition true for aperiodic systems like mine where there's artificial periodicity in z-direction ? If it's true then:<br>


</div>for, Egap = 1.15 eV, N3 = 1 and a3 = 20 angstrom gives an Emax = 0.05 V/angstrom<br></div>whereas
 if I use a3 as the layer thickness which is around 5.55 angstrom then 
Emax =  0.0207 V/angstrom which is the value after which the relax 
calculations are not converging. Is this interpretation correct?<br>
<br></div>           Any kind of help is highly appreciated.<br><br></div>Thanks,<br></div>Rajdeep Banerjee<br></div><div>(Ph. D. Student)<br></div>Theoretical Sciences Unit<br></div>Jawaharlal Nehru Centre for Advanced Scientific Research<br>


</div>Bangalore, India</div>