<div dir="ltr"><p>Dear everyone,</p><p>I have a convergency problem facing the pw.x after version 5.0.2</p><p>The example input is given in the last section.</p><p>The structure is a bcc-like supercell containing Fe and Cu atoms (relaxed) (with GBRV uspp).<br>The structure can be run properly with smearing, the problem comes from "tetrahedra" occupation.</p><p>If I use a non-offset k point mesh, as given in the input, the convergency would get stuck around 1E-6~1E-8 Ry, never get better like most other supercells of Fe and Cu.</p><p>The problem can be solved by offset the k mesh. However, my local mentors do not suggest me to do so. Because we have lots of different supercells, the offseted k mesh would give very different centering k point in different supercells, it would appear to be trouble if I compare their energies. (btw, what is your opinions on this?)</p><p>The problem can also be solved by increasing the size of k mesh, but still, not cheap and absolute necessary for the calculation of lots supercells.</p><p>However, in version 5.0.2, there is no such problem, and convergency can be achieved quickly (if I rerun [in post 5.0.2 versions] starting from the converged charge density [by v 5.0.2], convergency still cannot be achieved ). That's why I'm asking for help here, before I turn to other solutions.</p><p>So far, I tried to change number of bands, mixing beta, and parallel running parameter...all those appear to be irrelevant. Only k mesh would change the result.</p><p>Your suggestions and help is appreciated.</p><p>Yi Wang<br>Nanjing university of Science and Technology</p><p>&control<br>   calculation = 'scf' ,<br>   prefix = 'pwscf' ,<br>   outdir = './temp-p_t/' ,<br>   pseudo_dir = '/home/sim/pot/pwscf/' ,<br>   restart_mode = from_scratch ,<br>   disk_io= 'none'<br>/<br>&system<br>   ibrav = 0 ,<br>   ecutwfc = 40 ,<br>   ecutrho = 400 ,<br>   occupations = 'tetrahedra' ,<br>   nspin=2,<br>   nbnd= 120,<br>   nat = 12 ,<br>   ntyp = 2 , <br>   starting_magnetization( 1 )= 0.01,<br>   starting_magnetization( 2 )= 0.26,<br>   use_all_frac=.true.<br>/<br>&electrons<br>   conv_thr = 1d-9 ,<br>   diagonalization = 'david' ,<br>   mixing_mode = 'local-TF' ,<br>   startingpot = 'atomic' ,<br>   startingwfc = 'atomic+random'<br>   mixing_beta = 0.12 ,<br>   mixing_ndim =  16,<br>   diago_david_ndim=4<br>/<br>CELL_PARAMETERS bohr<br>5.689675841  0.000000000  0.000000000<br>0.000000000  -10.558955237  0.000000000<br>0.000000000  0.000000000  -15.996363568<br>ATOMIC_SPECIES<br>  Cu 63.546 cu_pbe_v1.2.uspp.F.UPF<br>  Fe 55.845 fe_pbe_v1.5.uspp.F.UPF<br>ATOMIC_POSITIONS crystal<br>Fe       1.000000000   0.498189104   0.666054422<br>Fe       0.500000000   0.251812190   0.500611021<br>Fe       1.000000000   0.003971761   0.664566351<br>Fe       0.500000000   0.746028153   0.502102015<br>Fe       1.000000000   0.500799770   0.333333330<br>Fe       0.500000000   0.249203627   0.833333320<br>Cu       1.000000000   0.998064252   0.333333150<br>Cu       0.500000000   0.751929936   0.833333540<br>Fe       0.500000000   0.251812380   0.166055619<br>Fe       0.500000000   0.746028003   0.164564485<br>Fe       1.000000000   0.498188874   0.000612278<br>Fe       1.000000000   0.003971951   0.002100469<br>K_POINTS automatic<br>   12 7 4 0 0 0<span><span><br></span></span></p></div>