<div dir="ltr">Hi QE Experts,<br><br>To find out the dielectric function (imaginary and real parts) of silicon, I have run <a href="http://scf.in">scf.in</a>, <a href="http://nscf.in">nscf.in</a>, <a href="http://epsilon.in">epsilon.in</a> (I have shared results in the attachment). I have compared the imaginary part of the dielectric function with one study and the results are in good agreement. However, the problem is in the real part of the dielectric function. It contains  several negative values and due to this absorption coefficient I found, has not matched with the study.  I am very confused why the real part shows absurd value, if the imaginary part matches with the study. Would you please share your opinions?<br>Here is the <a href="http://scf.in">scf.in</a>:<br> &control<br>  calculation = 'scf'<br>  restart_mode = 'from_scratch'<br>  prefix = 'ZG-si'<br>  pseudo_dir = './'<br>  outdir = './out'<br>/<br>&system<br>  ibrav = 0<br>  nat =   16<br>  ntyp =    1<br>  ecutwfc =  90.00<br>/<br>&electrons<br>  diagonalization = 'david'<br>  mixing_mode= 'plain'<br>  mixing_beta = 0.70<br>  conv_thr =   1D-06<br>/<br>ATOMIC_SPECIES<br>   Si   28.085 Si.pz-hgh.UPF<br>K_POINTS automatic<br>   8   8   8   0   0   0<br>CELL_PARAMETERS (angstrom)<br>      5.38154150      5.38154150      0.00000000<br>      0.00000000      5.38154150      5.38154150<br>      5.38154150      0.00000000      5.38154150<br>ATOMIC_POSITIONS (angstrom)<br>%%%%%%<br><a href="http://nscf.in">nscf.in</a><br>&control<br>  calculation = 'nscf'<br>  restart_mode = 'from_scratch'<br>  prefix = 'ZG-si'<br>  pseudo_dir = './'<br>  outdir = './out'<br>/<br>&system<br>  ibrav = 0<br>  nat =   16<br>  ntyp =    1<br>  ecutwfc =  90.00<br>  nbnd = 56<br>  nosym = .true.<br>  noinv = .true.<br>/<br>&electrons<br>  diagonalization = 'david'<br>  mixing_mode= 'plain'<br>  mixing_beta = 0.70<br>  conv_thr =   1D-06<br>/<br>ATOMIC_SPECIES<br>   Si   28.085 Si.pz-hgh.UPF<br><br>CELL_PARAMETERS (angstrom)<br>      5.38154150      5.38154150      0.00000000<br>      0.00000000      5.38154150      5.38154150<br>      5.38154150      0.00000000      5.38154150<br>ATOMIC_POSITIONS (angstrom)<br>K_POINTS crystal<br>%%%%%%<br>Here is the <a href="http://epsilon.in">epsilon.in</a>:<br>&inputpp<br>outdir      = './out',<br>prefix      = 'ZG-si',<br>calculation = 'eps'<br>/<br>&energy_grid<br>smeartype   = 'gauss',<br>intersmear  = 0.03d0,<br>wmax        = 4.5d0,<br>wmin        = 0.2d0,<br>nw          = 600,<br>shift       = 0.0d0,<br>/<br>%%%%%%<br><br>Thanks<br><br>Best Regards<br>Md Jahid Hasan<br>PhD student<br>Mechanical Engineering<br>University of Maine</div>