<div dir="ltr"><div>Dear QE users,<br><br></div><div>I just noticed that pw.x was giving me different energy values. I don't know if it's just due to machine accuracy or something strange happens.<br><br>1) I did the scf calculation with k points generated automatically and got energy_1 (-362.64960323 Ry). The input file is attached at the end of this email.<br></div><div>2) I copy the k points and weights spit out by step 1) to the input file and redo the calculation, I got exactly the same value as energy_1.<br></div><div>3) I add nosym=.true. in the input and redo the calculation. I got energy_3 (-362.64951123 Ry) which is different but not far from energy_1.<br>The difference is 0.00009200 Ry<br>Is this normal?<br><br></div><div>Then I tried to do a calculation with hybrid functional (HSE). q-mesh 1x1x1.<br>the same step 1) I got energy energy_hybrid_1 (-362.47024611 Ry)_<br></div><div>the same step 2) I got the same energy value as energy_hybrid_1<br></div><div>the same step 3) I got energy_hybrid_3 (-362.46997237 Ry)<br></div><div>Now the difference is 0.00027374 Ry about 3 times as the previous one.<br></div><div>Is this normal?<br><br></div><div>For my understanding, nosym affects the k point generation. <br>I checked the output of 1) 2) 3) and they have the same kpoints and weights.<br>In principle, I thought step 2) and 3) should yield the same energy. Only machine accuracy plays a little bit. If nosym is doing more than K points, please bear my limited knowledge and explain me where else symmetry is also involved and which might be relevant to my calculation.<br><br></div><div>Let me know your thoughts. Thanks!<br></div><div><br></div><div>Ye<br><br> &control<br>    calculation = 'scf'<br>    restart_mode='from_scratch',<br>    prefix='TiO2-rutile',<br>    tstress = .true.<br>    tprnfor = .true.<br>    pseudo_dir = '/home/yeluo/pw_pseudo/',<br>    outdir='./out',<br> /<br> &system<br>    ibrav     = 0,<br>    nat       = 6,<br>    ntyp      = 2,<br>    ecutwfc   = 300,<br>    input_dft = 'pbe'<br>    lda_plus_u = .true.,<br>    Hubbard_U(1) = 4.86041<br> /<br> &electrons<br>    mixing_beta = 0.7<br>    conv_thr =  1.0d-10<br> /<br>CELL_PARAMETERS bohr<br>8.681023246 0 0<br>0 8.681023246 0<br>0 0 5.59094331<br>ATOMIC_SPECIES<br> Ti  1.0  Ti.opt.upf<br> O   1.0  O.ncpp<br>ATOMIC_POSITIONS crystal<br> Ti  0.0000000   0.0000000   0.0000000<br> Ti  0.5000000   0.5000000   0.5000000<br> O   0.3050990   0.3050990   0.0000000<br> O   0.6949010   0.6949010   0.0000000<br> O   0.1949010   0.8050990   0.5000000<br> O   0.8050990   0.1949010   0.5000000<br>K_POINTS automatic<br>8 8 8 1 1 1<br><br></div><div><br></div><div><div><div><div><div dir="ltr">===================<br>
Ye Luo, Ph.D.<br>
Leadership Computing Facility<br>
Argonne National Laboratory</div></div></div>
</div></div></div>