<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none"><!-- P { margin-top: 0px; margin-bottom: 0px; }--></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,<br>
</p>
<p><br>
</p>
<p>Yes, you are right, both codes (turbo_lanczos.x and turbo_davidson.x) currently cannot be used for metallic systems (using smearing). If a smearing is needed, then it must be implemented (<span style="font-size: 12pt;">playing with</span><span style="font-size: 12pt;">
the parameters </span><span style="font-size: 12pt;">mixing_mode, mixing_beta, etc. will not help without the smearing, I think).</span></p>
<p><span style="font-size: 12pt;"><br>
</span></p>
<p><font size="3">There is also another code turbo_eels.x which can be used with the smearing. But this code was created to compute electron energy loss spectra (EELS) for solids. You may try to use it for your system in a supercell with k=0 (do not use K_POINTS{gamma},
because it is not supported). The code will produce a file with the <span style="font-family: Calibri, Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255);">
dielectric function</span> and inverse dielectric function for a specified transferred momentum q. You may try to set q very small (but not exactly zero), perform the calculation, and then plot the </font>imaginary<font size="3"> part of the dielectric function,
and see if this it what you need. Be aware that I did not test the turbo_eels.x code for finite systems, so there is no warranty that the results will be correct. Please report what you get, if you decide to try it out. And use the latest version of QE (i.e.
6.0). Thanks!</font></p>
<p><font size="3"><br>
</font></p>
<p><font size="3">HTH</font></p>
<p><font size="3"><br>
</font></p>
<p><font size="3">Regards,</font></p>
<p><font size="3">Iurii</font></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 color="808080" size="3" face="'Times New Roman', Times, serif">--<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 color="808080" size="3" face="'Times New Roman', Times, serif">Laboratory of Theory and Simulation of Materials
<font color="808080"><font face="'Times New Roman', Times, serif">(THEOS)</font></font></font><font color="808080" size="3" face="'Times New Roman', Times, serif"><br>
</font></div>
<font color="808080"></font>
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<font color="808080" size="3" face="'Times New Roman', Times, serif">CH-1015 Lausanne, Switzerland<br>
</font></div>
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<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 face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><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> Sunday, January 29, 2017 11:58 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>Dear Dr. Iurii Timrov,<br>
</div>
Thanks again for your valuable reply.<br>
</div>
I knew, how to check the convergence using both the codes but to cross verify myself, I asked you the question. Now I get the same absorption spectra for CH4 molecule using both the code (turbo_lanczos.x & turbo_davidson.x).
<br>
<br>
</div>
I was also trying for Carbon-dimer. But I am facing a basic problem. Without using, occupations = "smearing", this system is not getting converged. If I use "smearing", turbo_davidson.x or turbo_lanczs.x is not possible. I have tried with different kinds of
PPs, "mixing_mode" and "mixing_beta". Is there any solution for this problem ??<br>
<br>
</div>
<div>For Example I have used ::<br>
<br>
&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>
! celldm(1) = 20,<br>
nat = 2,<br>
ntyp = 1,<br>
ecutwfc = 25,<br>
ecutrho = 250,<br>
nosym = .true.,<br>
! occupations = "smearing",<br>
! smearing = "mv",<br>
! degauss = 0.001,<br>
/<br>
&ELECTRONS<br>
diagonalization='david'<br>
mixing_mode = 'plain'<br>
mixing_beta = 0.6,<br>
conv_thr = 1.0d-8<br>
/<br>
CELL_PARAMETERS angstrom<br>
15.0 0.0 0.0<br>
0.0 15.0 0.0<br>
0.0 0.0 15.0<br>
ATOMIC_SPECIES<br>
C 12.01 C.pbe-rrkjus.UPF<br>
ATOMIC_POSITIONS (angstrom)<br>
C 1.255833278 0.000000000 0.000000000<br>
C 0.000000000 0.000000000 0.000000000<br>
K_POINTS {Gamma}<br>
<br>
</div>
<div><br>
</div>
<font size="-1">Thanks & regards,<br>
Prem Sen<br>
Ph.D. Student,<br>
IIT Bombay,<br>
Mumbai,India.</font></div>
</div>
</div>
</body>
</html>