<div dir="ltr"><div><div><div><div><div><div>Dear Dr. Iurii Timrov,<br></div>                             Thanks for the reply. <br>"The two codes (Lanczos and Davidson) must give you the same absorption spectrum (when both are converged)." <br><br></div><div>How to check the convergence for both the codes ??<br></div><br>I was checking the optical absorption spectra for CH4 molecule (using the given example for the TDDFPT package) using both the code. They are <span style="color:rgb(255,0,0)">not</span> matching with each other. <br><br><br></div>Also I was trying to do for Carbon-dimer.  In this case also, the spectrum using both the codes are not matching with each other.  The input files are given below. <br><br>I have checked the ecut & ecutfock convergence and also relaxed the structure.<br><br></div><div><font size="4"><b>For turbo_davidson.x ::</b></font><br></div><div><br></div><b><font size="2"><span style="color:rgb(255,0,0)"><span style="background-color:rgb(255,255,255)">PW.X</span></span>:</font><br><br></b>&CONTROL<br>  calculation  = "scf",<br>  restart_mode = "from_scratch"<br>  prefix       = "C2",<br>  pseudo_dir   = "../pseudo/",<br>  outdir       = "./tmp/",<br>/<br>&SYSTEM<br>  ibrav     = 0,<br>  nat       = 2,<br>  ntyp      = 1,<br>  ecutwfc   = 75,<br>  ecutfock  = 75,<br>  nosym     = .true.,<br>  input_dft = 'PBE0',<br>  x_gamma_extrapolation = .false.,<br>  exxdiv_treatment = 'vcut_spherical',<br>/<br>&ELECTRONS<br>  conv_thr    = 1.D-10,<br>  mixing_beta = 0.6,<br>  adaptive_thr = .true.,<br>/<br>CELL_PARAMETERS angstrom<br>10.0  0.0  0.0<br> 0.0 10.0  0.0<br> 0.0  0.0 10.0<br>ATOMIC_SPECIES<br>C  12.01  C.pbe-mt_gipaw.UPF<br>ATOMIC_POSITIONS (angstrom)<br>C        1.248029552   0.000000000   0.000000000    <br>C        0.000000000   0.000000000   0.000000000   <br>K_POINTS {Gamma}<b><br><br></b></div><b><font size="2"><span style="color:rgb(255,0,0)">Turbo_davidson.x</span>:</font><br></b>&lr_input<br>    prefix="C2",<br>    outdir="./tmp/",<br>/<br>&lr_dav<br>    ecutfock  = 75,<br>    num_eign=30,<br>    num_init=60,<br>    num_basis_max = 300,<br>    residue_conv_thr = 1.0E-4,<br>    start = 0.0<br>    finish = 1.0,<br>    step = 0.0005,<br>    broadening = 0.005,<br>    reference = 0.0,<br>    d0psi_rs = .true.,<br>    lshift_d0psi = .true.,<br>    p_nbnd_virt = 20<br>/<b><br><br></b></div><b><font size="4">For turbo_lanczos.x ::</font><br></b><div><br></div><div>I used the same file for pw.x as given above.<br></div><div><br><br><b><font size="2"><b><span style="color:rgb(255,0,0)">Turbo_lanczos.x</span>:</b></font><br><br></b>&lr_input<br>    prefix="C2",<br>    outdir="./tmp/",<br>    restart_step=500,<br>/<br>&lr_control<br>    itermax=2000,<br>    ipol=4,<br>    ecutfock  = 75,<br>    pseudo_hermitian = .true.,<br>    d0psi_rs = .true.,<br>/<b><br><br></b><br><b><b><font size="2"><b><span style="color:rgb(255,0,0)">Turbo_spectrum.x</span>:</b></font></b><br><br></b>&lr_input<br>  prefix='C2',<br>  outdir="./tmp/",<br>  itermax0 = 410<br>  itermax  = 10000<br>  epsil=0.01<br>  start=0.0d0<br>  end=1.50d0<br>  increment=0.001d0<br>  ipol=4<br>/<b><br><br><br><br></b></div><div>Thanks & regards,<br></div><div>Prem Sen<br></div><div>Ph.D. Student,<br></div><div>IIT Bombay<b>,<br></b></div><div>Mumbai,India<b>.<br></b></div></div>