<div dir="ltr">Hi,<div><br></div><div><br></div><div>i completely agree with Thomas. As I understand this should not be a metal and you should not encounter troubles integrating around the Fermi level. So I would not specify smearing and degauss (which should be around 0.01 anyway if you really have to). Moreover, I am not entirely sure about placing the starting magnetization on C atom type...i would rather place it on iron and (most important) would start with a 0.6 value in order to avoid sticking in magnetic troubling states. Finally, how did you generate the structure? Where does it come from? And your box specification? Cut off is low, i agree, but should not be the problem...however i would set it at least around 40 because you are far from convergency anyway!Very last thing...you should start with electronic threshold set to the default.</div>
<div><br></div><div>Best regards</div><div><br></div><div>Vincenzo</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Sep 1, 2014 at 1:57 PM, 琨陶 <span dir="ltr"><<a href="mailto:taokun76@gmail.com" target="_blank">taokun76@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div>Dear all, <br><br></div>I just start to learn Quantum Espresso, and want to optimize the geometry of an FePc molecule. But, it doesn't converge even after 500 steps. Could you give me some suggestions? Any suggestions are welcome!<br>
<br></div><div>Regards,<br></div><div>Kun Tao<br><br><br></div>Attached below is the input file:<br><br>&CONTROL<br> calculation = 'scf'<br> restart_mode='from_scratch',<br> prefix='FePc',<br>
pseudo_dir = './',<br> outdir='./temp'<br> etot_conv_thr=1.D-4<br> forc_conv_thr=1.D-2<br> /<br> &SYSTEM<br> ibrav=6,<br> celldm(1)=41.58,<br> celldm(3)=0.454545,<br> nat= 57,<br>
ntyp= 4,<br> ecutwfc = 27.0,<br> ecutrho = 250.0,<br> & occupations='fixed',<br> smearing='marzari-vanderbilt',<br> degauss=0.1,<br> nspin=2,<br> assume_isolated = 'martyna-tuckerman'<br>
starting_magnetization(1)= 1,<br> & lda_plus_u = .true.<br> & Hubbard_U(1)= 1.d-20<br> & Hubbard_U(2)= 1.d-20<br> & nbnd = 380<br> /<br> &ELECTRONS<br> diagonalization='cg',<br> conv_thr = 1.D-4,<br>
mixing_beta = 0.05D0,<br> electron_maxstep = 500,<br> /<br>&IONS<br> & ion_dynamics = 'bfgs'<br> bfgs_ndim = 3<br>/<br>ATOMIC_SPECIES<br> C 12.0 C.pw91-van_ak.UPF<br> H 1.0 H.pw91-van_ak.UPF<br>
N 14.0 N.pw91-van_ak.UPF<br> Fe 55.8 Fe.pw91-sp-van_ak.UPF<br>ATOMIC_POSITIONS {Angstrom}<br> C 12.15539372 13.73909330 5.04996061<br> C 13.73649156 12.16178989 5.05578041<br> C 13.43098617 14.43384445 5.05163491<br>
C 14.42783630 13.43936718 5.06136000<br> C 13.76218534 15.79065418 5.04836023<br> C 15.78316009 13.77488291 5.07699788<br> C 13.73909330 9.84460628 5.04996061<br> C 12.16178989 8.26350933 5.05578041<br>
C 14.43384445 8.56901383 5.05163491<br> C 13.43936718 7.57216460 5.06136000<br> C 15.79065418 8.23781467 5.04836023<br> C 13.77488291 6.21683902 5.07699788<br> C 9.84460628 8.26090580 5.04996061<br>
C 8.26350933 9.83821011 5.05578041<br> C 8.56901383 7.56615555 5.05163491<br> C 7.57216460 8.56063193 5.06136000<br> C 8.23781467 6.20934671 5.04836023<br> C 6.21683902 8.22511798 5.07699788<br>
C 8.26090580 12.15539372 5.04996061<br> C 9.83821011 13.73649156 5.05578041<br> C 7.56615555 13.43098617 5.05163491<br> C 8.56063193 14.42783630 5.06136000<br> C 6.20934671 13.76218534 5.04836023<br>
C 8.22511798 15.78316009 5.07699788<br> C 6.87043524 16.11305392 5.07808685<br> C 5.87526149 15.11545813 5.06072640<br> C 5.88694608 6.87043524 5.07808685<br> C 6.88454098 5.87526149 5.06072640<br>
C 15.11545813 16.12473762 5.06072640<br> C 16.11305392 15.12956476 5.07808685<br> C 16.12473762 6.88454098 5.06072640<br> C 15.12956476 5.88694608 5.07808685<br> H 12.98551416 16.55329919 5.03687501<br>
H 16.54878414 13.00105667 5.08993745<br> H 16.55329919 9.01448494 5.03687501<br> H 13.00105667 5.45121586 5.08993745<br> H 9.01448494 5.44669992 5.03687501<br> H 5.45121586 8.99894244 5.08993745<br>
H 5.44669992 12.98551416 5.03687501<br> H 8.99894244 16.54878414 5.08993745<br> H 6.57077259 17.16044581 5.09259462<br> H 4.82672697 15.41213369 5.05869567<br> H 4.83955419 6.57077259 5.09259462<br>
H 6.58786541 4.82672697 5.05869567<br> H 15.41213369 17.17327213 5.05869567<br> H 17.16044581 15.42922831 5.09259462<br> H 17.17327213 6.58786541 5.05869567<br> H 15.42922831 4.83955419 5.09259462<br>
N 12.36901915 9.63413960 5.05092382<br> N 9.63413960 9.63097996 5.05092382<br> N 9.63097996 12.36585951 5.05092382<br> N 14.37677419 11.00406444 5.05336285<br> N 11.00406444 7.62322491 5.05336285<br>
N 7.62322491 10.99593645 5.05336285<br> N 10.99593645 14.37677419 5.05336285<br> N 12.36585951 12.36901915 5.05092382<br> Fe 11.00000000 11.00000000 5.04961908<br>K_POINTS Gamma<br></div>
<br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div><br></div>