<div dir="ltr"><div>Hi, <br></div><div>I wonder why is the band calculation much longer than a scf calcuation.<br></div><div>I did a calculation for a 64-atoms supercell using 20 cpu cores. <br></div><div>In summary, I did<br></div><div>1) SCF calculation with only gamma point, run in two minutes. <br></div><div>2) Non-self-consistent band calculation with 27 k-points. <br><br></div><div>In the nscf calculation each k-point takes 15-20 minutes and the total calculation lasts 7 hours. I understand that for the k-points along the line the wfc are complex, the symmetry is smaller, but even so the time scale is quite different.  <br>The machine has to keep all the wfc in memory, but it was not swapping. <br>Even the gamma point took 15 minutes, why does it takes longer than the scf calculation?<br><br></div><div>Please, let me know if  I am doing something wrong. Below are the inputs for the scf and the bnd calculations. <br><br></div><div>SCF: <br><br>&CONTROL<br><br>calculation = 'scf'<br><br>restart_mode='from_scratch',<br><br>pseudo_dir = "./",<br><br>outdir="./tmp/",<br><br>tstress = .true.,<br><br>tprnfor = .true. ,<br><br>prefix= 'cdte64',<br><br>disk_io = 'low',<br><br>etot_conv_thr = 1.0D-3,<br><br>wf_collect=.TRUE.,<br><br>verbosity='high'<br><br>/<br><br>&SYSTEM<br><br>ibrav= 1,<br><br>celldm(1) = 25.02,<br><br>nat = 64,<br><br>ntyp = 2,<br><br>ecutwfc = 28.67,<br><br>ecutrho = 180,<br><br>degauss = 0.001D0,<br><br>occupations = "smearing",<br><br>smearing = "gauss",<br><br>tot_charge = 0,<br><br>nspin = 1 , ! set to 2 for spin polarized<br><br>! tot_magnetization = 0 , ! change to 1<br><br>/<br><br>&ELECTRONS<br><br>diagonalization='david'<br><br>conv_thr = 1.D-6,<br><br>mixing_beta = 0.4D0,<br><br>diago_david_ndim= 2,<br><br>/<br><br>&IONS<br><br>pot_extrapolation='second_order',<br><br>trust_radius_max = 0.7,<br><br>trust_radius_ini = 0.7,<br><br>/<br><br>ATOMIC_SPECIES<br><br>Cd 112.411 cd_pbe_v1.uspp.F.UPF<br><br>Te 127.600 te_pbe_v1.uspp.F.UPF<br><br><br>ATOMIC_POSITIONS (alat)<br><br>Cd 0.5000000000000000 0.5000000000000000 0.5000000000000000<br><br>Cd 0.0000000000000000 0.0000000000000000 0.0000000000000000<br><br>Cd 0.2500000000000000 0.2500000000000000 0.0000000000000000<br><br>….<br><br><br>K_POINTS gamma<br><br><span style="color:rgb(0,0,255)"><br>The NSCF band calculation input is:</span><br><br>&CONTROL<br><br>calculation = 'bands'<br><br>! restart_mode='from_scratch',<br><br>pseudo_dir = "./",<br><br>outdir="./tmp/",<br><br>tstress = .true.,<br><br>tprnfor = .true. ,<br><br>prefix= 'cdte64',<br><br>disk_io = 'low',<br><br>etot_conv_thr = 1.0D-3,<br><br>verbosity='high'<br><br>/<br><br>&SYSTEM<br><br>ibrav= 1,<br><br>celldm(1) = 25.02,<br><br>nat = 64,<br><br>ntyp = 2,<br><br>ecutwfc = 28.67,<br><br>ecutrho = 180,<br><br>degauss = 0.001D0,<br><br>occupations = "smearing",<br><br>smearing = "gauss",<br><br>tot_charge = 0,<br><br>nspin = 1 , ! set to 2 for spin polarized<br><br>! tot_magnetization = 0 , ! change to 1<br><br>/<br><br>&ELECTRONS<br><br>diagonalization='david'<br><br>conv_thr = 1.D-6,<br><br>mixing_beta = 0.4D0,<br><br>diago_david_ndim= 2,<br><br>/<br><br>&IONS<br><br>pot_extrapolation='second_order',<br><br>trust_radius_max = 0.7,<br><br>trust_radius_ini = 0.7,<br><br>/<br><br>ATOMIC_SPECIES<br><br>Cd 112.411 cd_pbe_v1.uspp.F.UPF<br><br>Te 127.600 te_pbe_v1.uspp.F.UPF<br><br><br>ATOMIC_POSITIONS (alat)<br><br>Cd 0.5000000000000000 0.5000000000000000 0.5000000000000000<br><br>Cd 0.0000000000000000 0.0000000000000000 0.0000000000000000<br><br>Cd 0.2500000000000000 0.2500000000000000 0.0000000000000000<br><br>….<br><br>K_POINTS crystal_b<br><br>3<br><br>0.500 0.500 0.500 17<br><br>0.000 0.000 0.000 10<br><br>0.500 0.000 0.000 1<p style="margin-bottom:0cm"><br>
</p>Thank you. <br></div><div><br clear="all"><div><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr">
<div>Eduardo Menendez Proupin<br>Departamento de Fisica, Facultad de Ciencias, Universidad de Chile <br>URL: <a href="http://www.gnm.cl/emenendez" target="_blank">http://www.gnm.cl/emenendez</a></div><div><br></div></div></div></div></div></div></div></div></div></div>
</div></div>