<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Arial,Helvetica,sans-serif;" dir="ltr">
<p>Hi there,</p>
<p><br>
</p>
<p> I am trying to perform a berry phase calculation with a homogeneous electric field (lelfield=.true.) however I receive the "wrong k-strings?" in c_phase error while running in parallel with OpenMPI on 4 processors:</p>
<p><br>
</p>
<p>From the documentation, if K_POINTS is automatic, then field strength should be specified in cartesian coordinates (efield_cart(i)) and 'efield' is not used, however reading the berry phase documentation says that the direction of the efield is taken along
gdir when a homogeneous electric field is applied.</p>
<p><br>
</p>
<p>Must I still specify efield_cart values when gdir is specified (it would seem so)? The code doesn't complain until after bandstructure computation. I have tried efield calculation with the following specification (see below for input file used):</p>
<p><br>
</p>
<p>efield_cart(1) = 0.0001,</p>
<p>efield_cart(2) = 0,</p>
<p>efield_cart(3) = 0,</p>
<p><br>
</p>
<p>which results in the electric field existing in two directions in crystal reference system, but the error is the same whether I do or do not include these fields and leave efield=0.0001<br>
</p>
<p><br>
</p>
<p>My process is as follows:</p>
<p> - scf calculation of trigonal alpha-Al2O3 cell (input crystal is correct)</p>
<p> - band structure calculation</p>
<p> - Berry phase calculation *without* electric field (this works and output is correct) with input:</p>
<p><br>
</p>
<p></p>
<div> &CONTROL<br>
title = 'alpha-Alumina' ,<br>
calculation = 'nscf' ,<br>
outdir = '${PWD}' ,<br>
pseudo_dir = '${PWD}/pseudo/' ,<br>
verbosity = 'high' ,<br>
prefix = 'run_${RUN_NUMBER}',<br>
lberry = .true.,<br>
gdir = 1,<br>
nppstr = 10,<br>
lelfield = .false.,<br>
nberrycyc = 1<br>
/<br>
&SYSTEM<br>
ibrav = 5,<br>
A = 5.136,<br>
cosAB = 0.56956647115,<br>
nat = 10,<br>
ntyp = 2,<br>
occupations = 'fixed',<br>
ecutwfc = 50,<br>
ecutrho = 350,<br>
nbnd = 48<br>
/<br>
&ELECTRONS<br>
conv_thr = 1d-8,<br>
diagonalization = 'cg',<br>
mixing_beta = 0.7,<br>
electron_maxstep = 200,<br>
efield = 0,<br>
efield_phase = 'write'<br>
/<br>
ATOMIC_SPECIES<br>
O 15.99900 O.pbe-van_ak.UPF <br>
Al 26.98200 Al.pbe-n-van.UPF <br>
ATOMIC_POSITIONS crystal <br>
O 0.556000000 0.944000000 0.250000000 <br>
O 0.444000000 0.056000000 0.750000000 <br>
O 0.250000000 0.556000000 0.944000000 <br>
O 0.750000000 0.444000000 0.056000000 <br>
O 0.944000000 0.250000000 0.556000000 <br>
O 0.056000000 0.750000000 0.444000000 <br>
Al 0.352000000 0.352000000 0.352000000 <br>
Al 0.648000000 0.648000000 0.648000000 <br>
Al 0.148000000 0.148000000 0.148000000 <br>
Al 0.852000000 0.852000000 0.852000000 <br>
K_POINTS automatic<br>
10 4 4 0 0 0 <br>
<br>
</div>
- Then Berry phase calculation *with* electric field:
<p></p>
<p><br>
</p>
<p></p>
<div> &CONTROL<br>
title = 'alpha-Alumina' ,<br>
calculation = 'nscf' ,<br>
outdir = '${PWD}' ,<br>
pseudo_dir = '${PWD}/pseudo/' ,<br>
verbosity = 'high' ,<br>
prefix = 'run_${RUN_NUMBER}',<br>
lberry = .true.,<br>
gdir = 1,<br>
nppstr = 10,<br>
lelfield = .true.,<br>
nberrycyc = 1<br>
/<br>
&SYSTEM<br>
ibrav = 5,<br>
A = 5.136,<br>
cosAB = 0.56956647115,<br>
nat = 10,<br>
ntyp = 2,<br>
occupations = 'fixed',<br>
ecutwfc = 50,<br>
ecutrho = 350,<br>
nbnd = 48,<br>
nosym = .false.<br>
/<br>
&ELECTRONS<br>
conv_thr = 1d-8,<br>
diagonalization = 'cg',<br>
mixing_beta = 0.7,<br>
electron_maxstep = 200,<br>
efield = 0.0001,<br>
efield_phase = 'read'<br>
/<br>
ATOMIC_SPECIES<br>
O 15.99900 O.pbe-van_ak.UPF <br>
Al 26.98200 Al.pbe-n-van.UPF <br>
ATOMIC_POSITIONS crystal <br>
O 0.556000000 0.944000000 0.250000000 <br>
O 0.444000000 0.056000000 0.750000000 <br>
O 0.250000000 0.556000000 0.944000000 <br>
O 0.750000000 0.444000000 0.056000000 <br>
O 0.944000000 0.250000000 0.556000000 <br>
O 0.056000000 0.750000000 0.444000000 <br>
Al 0.352000000 0.352000000 0.352000000 <br>
Al 0.648000000 0.648000000 0.648000000 <br>
Al 0.148000000 0.148000000 0.148000000 <br>
Al 0.852000000 0.852000000 0.852000000 <br>
K_POINTS automatic<br>
10 4 4 0 0 0 </div>
<br>
My last thought is that when the electric field is specified in cartesian coordinates as above, the fact that the field in the crystal reference system is along two directions causes the problem, I have struggled to get the field to exist along one direction
in the crystal reference system (does anyone know a reference I could use to do this properly?). Maybe the sensible thing to do is not use automatic K_POINTS?<br>
<br>
Kindest regards,<br>
Louis<br>
<p></p>
</div>
</body>
</html>