Hello, <div><br></div><div>I have recently tried running an ionic relaxation with a cell using the effective screening medium (ESM) method, but have run into a hanging issue. </div><div><br></div><div>Whenever I try running the relaxation, the program hangs after the first scf minimization converges and before the forces on the atoms are calculated. The output stops immediately after the ESM charge and potential are printed. What's interesting is this error also happens when I set the esm boundary conditions to "pbc" which should remove any implementation of ESM should it not? The relaxation runs without problems when all esm related variables are removed. </div>
<div><br></div><div>This issue is reproducible on 4.3.1 and 4.3.2. Is ESM just not yet able to do 'relax' runs? </div><div><br></div><div>Thank you for any help you may provide.</div><div><br></div><div>My input file is attached below:<div>
<br></div><div>--------------------------------------------</div><div><div>&control</div><div> calculation='relax',</div><div> restart_mode='from_scratch',</div><div> prefix='qe-test',</div><div>
 pseudo_dir='/usr/global/qe/4.3.2/pseudo',</div><div> outdir='/gpfs/home/ams751/scratch/test',</div><div> tprnfor = .TRUE.</div><div>/</div><div>&system</div><div>    ibrav = 0, celldm(1) = 1.889725989, </div>
<div>    tot_charge=0,</div><div>    nat= 2, ntyp= 2,</div><div>    ecutwfc = 22.05,</div><div>    starting_magnetization(2) = 1.0</div><div>    nspin = 2,</div><div>    nosym=.TRUE.</div><div>    occupations='smearing', smearing='gauss', degauss=0.004</div>
<div>    assume_isolated='esm', esm_bc='pbc'</div><div>/</div><div>&electrons</div><div>    mixing_beta = 0.7</div><div>/</div><div>&ions</div><div>/</div><div>ATOMIC_SPECIES</div><div> C   12.000000   C.pbe-rrkjus.UPF</div>
<div> H    1.000000   H.pbe-rrkjus.UPF</div><div>CELL_PARAMETERS</div><div> 12.3000000000000000    0.0000000000000000    0.0000000000000000</div><div> -6.1500000000000000   10.6521120000000000    0.0000000000000000</div><div>
  0.0000000000000000    0.0000000000000000   20.0000000000000000</div><div>ATOMIC_POSITIONS (crystal)</div><div>C 0.46667 0.53333 0.00000</div><div>H 0.46667 0.53333 0.05000</div><div><br></div><div>K_POINTS automatic</div>
<div> 1 1 1 0 0 0</div><div><br></div><br><br>Alejandro Suarez<br>Graduate Student<br>Department of Physics<br>Penn State University<br>
</div></div>