<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none"><!--P{margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr" style="font-size:12pt;color:#000000;background-color:#FFFFFF;font-family:Calibri,Arial,Helvetica,sans-serif;">
<p>Dear Prem Sen,</p>
<p><br>
</p>
<p>> How to check the convergence for both the codes ??</p>
<p><br>
</p>
<p>First of all you need to make convergence tests for the ground-state calculation (with respect to the supercell size, kinetic energy cutoff for the wavefunctions (and density/potentials in the case of USPP)). In the turboLanczos calculation, check the convergence
of the spectrum with respect to the number of Lanczos iterations. In the turboDavidson calculation, check that you compute enough eigenvalues of the Liouvillian (num_eign) to obtain the spectrum in the energy range of interest. I suggest you to read papers
about these codes (please see TDDFPT/README), if you have not done this so.</p>
<p><br>
</p>
<p>> 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. </p>
<p><br>
</p>
<p>Did you compare the examples "CH4" and "CH4-DAV"? These examples do not have converged parameters. First of all, in both calculations (Lanczos and Davidson) you have to start from the same ground state calculation: note that in "CH4" and "CH4-DAV" the input
files are not the same. </p>
<p>Did you perform convergence tests in these examples in the linear-response part? Note that with Lanczos you will obtain an absorption spectrum in a wide energy range, whereas with Davidson you will obtain just a few peaks (depends on how many eigenvalues
do you compute). Hence those peaks which you compute with Davidson should coincide with the peaks in the much more extended spectrum computed with Lanczos. </p>
<p><br>
</p>
<p>> 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.</p>
<p><br>
</p>
<p>Before studying your system, try to solve the issue mentioned above (i.e. reach the agreement between spectra computed with two codes). Try first more simple XC functionals (LDA or GGA), and see if the spectra agree. Then try hybrid functionals. Try also
another pseudopotential (without the word "gipaw"). </p>
<p><br>
</p>
<p>You use different broadening parameters in Lorentzian smearing: In Davidson you have 0.005 Ry and in Lanczos you have 0.01 Ry. If you want the spectra to agree, you must use the same broadening parameter.</p>
<p><br>
</p>
<p>When you use the turbo_spectrum.x code after the turbo_lanczos.x code, set</p>
<p> itermax0 = 2000<br>
itermax = 10000</p>
<p> extrapolation="osc"</p>
<p>Please see the paper turboTDDFT (see TDDFPT/README) in order to understand what does this mean.</p>
<p><br>
</p>
<p>HTH</p>
<p><br>
</p>
<p>Regards,</p>
<p>Iurii<br>
</p>
<p> <br>
</p>
<div id="Signature">
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<font size="3" face="'Times New Roman', Times, serif" color="808080">--<br>
Dr. Iurii Timrov<br>
Postdoctoral Researcher<br>
Swiss Federal Institute of Technology Lausanne (EPFL<font color="808080"><font face="'Times New Roman', Times, serif">)</font></font>
</font></div>
<font color="808080"></font>
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<font size="3" face="'Times New Roman', Times, serif" color="808080">Laboratory of Theory and Simulation of Materials
<font color="808080"><font face="'Times New Roman', Times, serif">(THEOS)</font></font></font><font size="3" face="'Times New Roman', Times, serif" color="808080"><br>
</font></div>
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<font size="3" face="'Times New Roman', Times, serif" color="808080">CH-1015 Lausanne, Switzerland<br>
</font><br>
</div>
</div>
</div>
<div style="color: rgb(33, 33, 33);">
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" face="Calibri, sans-serif" color="#000000"><b>From:</b> pw_forum-bounces@pwscf.org <pw_forum-bounces@pwscf.org> on behalf of prem sen <premsen1502@gmail.com><br>
<b>Sent:</b> Tuesday, January 24, 2017 11:53 AM<br>
<b>To:</b> pw_forum@pwscf.org<br>
<b>Subject:</b> Re: [Pw_forum] Learning TDDFT calculation in Atomic label (CH4 molecule) using turbo_lanczos.x and turbo_davidson.x</font>
<div> </div>
</div>
<div>
<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>
</div>
</div>
</body>
</html>