<div dir="ltr"><div><div><div><div><div><div>Dear all,<br><br></div>I want to calculate phonon with soc, but when I start to calculate ph.x, the program gives me an error:<br><br>task #         1<br>     from checkallsym : error #         1<br>     some of the original symmetry operations not satisfied<br><br><br></div>this is my input file:<br> &control<br>        calculation = 'scf'<br>        prefix='$PREFIX'<br>        pseudo_dir = '/home-fs4/hp0273/si/pseudo/'<br>        outdir='$TMP_DIR/'<br>        verbosity = 'high'<br>        wf_collect=.true.<br> /<br> &system<br>        ibrav=  4<br>        celldm(1) = 6.7460000000<br>        celldm(3) = 5.4150607768<br>        ecutwfc = 36<br>        ecutrho =550<br>        nat = 4<br>        ntyp = 1<br>        occupations = 'smearing', smearing = 'mp', degauss = 0.02D0<br>        lspinorb = .true.<br>        noncolin = .true.<br> /<br> &electrons<br>        electron_maxstep = 200<br>        mixing_beta = 0.2<br>        mixing_mode = "local-TF"<br>        conv_thr    = 1.0d-7<br> /<br>ATOMIC_SPECIES<br>Pb  207.1900   Pb.rel-pbe-dn-rrkjus_psl.0.2.2.UPF<br>ATOMIC_POSITIONS {crystal}<br> Pb    0.00000000    0.00000000    0.27374760<br> Pb   -0.33333000    0.33333000    0.42458253<br> Pb    0.33333000   -0.33333000    0.57541747<br> Pb    0.00000000    0.00000000    0.72625240<br>K_POINTS automatic<br> 16 16 1   0 0 0<br></div><br><br> for phonon:<br>Phononen<br>&inputph<br>        prefix = '$PREFIX'<br>        outdir = '$TMP_DIR/'<br>        tr2_ph = 1.0d-14<br>        alpha_mix(1) = 0.3<br>        alpha_mix(10) = 0.2<br>        niter_ph = 100<br>        recover = .false.<br>        ldisp = .true.<br>        nq1 = 6, nq2 = 6, nq3 = 1<br>        amass(1) = 207.19, <br>        fildyn = '$PREFIX.dyn'<br>        reduce_io = .true.<br>/<br></div>would you please tell me what is the problem?<br><br></div>Best,<br></div>Maedeh Zahedifar<br><div><div><div><br></div></div></div></div>