<div dir="ltr"><div><div>Dear QE users and devalopers<br>I do a phonon calculation of pure titanium in orthorhombic phase with 4 atoms in the unit cell:<br>ibrav=  8<br>Ti      -0.000000000   0.110000000   0.250000000    0   0   0<br>Ti      -0.000000011  -0.103719384   0.749999888<br>Ti       0.500000007   0.610015023   0.250000018<br>Ti       0.499999990   0.396287113   0.749999901<br><br></div>The calculation succeed when I run it with: ecut=50 and ecutrho=500 (I use ultrasoft PP). When I increase the cutoff to ecut=50 and ecutrho=500 and don't change anything else I got this error in the from matdyn:<br>     Error in routine frc_blk (1):<br>     wrong total_weight<br><br>Relaxation to 1000 kbar work good and the stress is:<br>      total   stress  (Ry/bohr**3)                               (kbar)     P= 1000.04<br>   0.00679839   0.00000000   0.00000000       1000.08      0.00      0.00<br>   0.00000000   0.00679814   0.00000000          0.00   1000.04      0.00<br>   0.00000000   0.00000000   0.00679787          0.00      0.00   1000.00<br><br></div><div>What is the problem?<br><br></div><div>Thank you very much<br></div><div>Uri Argaman<br></div><div>Ben-Gurion university<br></div><div>Israel<br></div></div>