<div dir="ltr"><div dir="ltr"><br>Dear QE users,</div><div dir="ltr">Hello</div><div dir="ltr">   <br>I hope everybody is doing well. I am using hp.x to calculate hubbard parameters for my structure. The first scf input is as follow:<br><br>&control<br> calculation='scf',<br> verbosity = 'low',<br> restart_mode = 'from_scratch',<br> prefix = 'my_file',<br> tstress = .true.,<br> tprnfor = .true.,<br> forc_conv_thr = 1.0D-4 ,<br>    pseudo_dir ='./',<br>    outdir='./tmp',<br> nstep=3000,<br> wf_collect=.true.<br> verbosity='high'<br>/ <br><br>&SYSTEM<br> ibrav=0,<br> nat=10,<br> ntyp=3,<br> ecutwfc=42,<br> ecutrho=420,<br> occupations='smearing',<br> degauss=0.001,<br> smearing='mv',<br> nspin = 2,<br> starting_magnetization(1) =  0.5,<br> starting_magnetization(2) = -0.5,<br> lda_plus_u = .true.,<br> lda_plus_u_kind = 0,<br> U_projection_type = 'atomic',<br> Hubbard_U(1) = 1.d-8<br> Hubbard_U(2) = 1.d-8<br><br>/<br>&electrons<br> mixing_beta = 0.7,<br> conv_thr =  1.0d-10,<br> electron_maxstep=600,<br>/ <br>&IONS<br> ion_dynamics  = 'bfgs',<br>/<br>&cell <br> cell_dynamics = 'bfgs',<br>/<br><br>ATOMIC_SPECIES <br> Ti1  47.867    Ti.pbe-sp-van_ak.UPF<br> Ti2  47.867    Ti.pbe-sp-van_ak.UPF<br> O   15.9994   O.pbe-van_ak.UPF<br><br>K_POINTS {automatic}<br>4 4 4 0 0 0<br><br>CELL_PARAMETERS (bohr)<br>   4.800141587   2.771363038   8.781343859<br>  -4.800141587   2.771363038   8.781343859<br>   0.000000000  -5.542726075   8.781343859<br><br>ATOMIC_POSITIONS (crystal)<br>Ti1      0.654936732   0.654936732   0.654936732<br>Ti2      0.345063271   0.345063271   0.345063271<br>Ti2      0.845063271   0.845063271   0.845063271<br>Ti1      0.154936732   0.154936732   0.154936732<br>O        0.063790525   0.436209470   0.750000001<br>O        0.936209470   0.563790525   0.250000001<br>O        0.436209470   0.750000001   0.063790525<br>O        0.563790525   0.250000001   0.936209470<br>O        0.750000001   0.063790525   0.436209470<br>O        0.250000001   0.936209470   0.563790525<br><br>and the second scf input for making the occupation fixed is as follow:<br><br><br> calculation='scf',<br> verbosity = 'low',<br> restart_mode = 'from_scratch',<br> prefix = 'my_file',<br> tstress = .true.,<br> tprnfor = .true.,<br> forc_conv_thr = 1.0D-4 ,<br>    pseudo_dir ='./',<br>    outdir='./tmp',<br> nstep=3000,<br> wf_collect=.true.<br> verbosity='high'<br>/ <br><br>&SYSTEM<br> ibrav=0,<br> nat=10,<br> ntyp=3,<br> ecutwfc=42,<br> ecutrho=420,<br> occupations='fixed',<br> nspin = 2,<br> tot_magnetization = 0.00 ,<br> lda_plus_u = .true.,<br> lda_plus_u_kind = 0,<br> U_projection_type = 'atomic',<br> Hubbard_U(1) = 1.d-8,<br> Hubbard_U(2) = 1.d-8,<br> nbnd = 50,<br>/<br>&electrons<br> mixing_beta = 0.7,<br> conv_thr =  1.0d-15,<br> electron_maxstep=600,<br> startingpot = 'file',<br> startingwfc = 'file',<br>/ <br>&IONS<br> ion_dynamics  = 'bfgs',<br>/<br>&cell <br> cell_dynamics = 'bfgs',<br>/<br><br>ATOMIC_SPECIES <br> Ti1  47.867    Ti.pbe-sp-van_ak.UPF<br> Ti2  47.867    Ti.pbe-sp-van_ak.UPF<br> O   15.9994   O.pbe-van_ak.UPF<br><br>K_POINTS {automatic}<br>4 4 4 0 0 0<br><br>CELL_PARAMETERS (bohr)<br>   4.800141587   2.771363038   8.781343859<br>  -4.800141587   2.771363038   8.781343859<br>   0.000000000  -5.542726075   8.781343859<br><br>ATOMIC_POSITIONS (crystal)<br>Ti1      0.654936732   0.654936732   0.654936732<br>Ti2      0.345063271   0.345063271   0.345063271<br>Ti2      0.845063271   0.845063271   0.845063271<br>Ti1      0.154936732   0.154936732   0.154936732<br>O        0.063790525   0.436209470   0.750000001<br>O        0.936209470   0.563790525   0.250000001<br>O        0.436209470   0.750000001   0.063790525<br>O        0.563790525   0.250000001   0.936209470<br>O        0.750000001   0.063790525   0.436209470<br>O        0.250000001   0.936209470   0.563790525<br><br><br><br>my input hp.x is as follow<br><br> &inputhp<br>  prefix = 'my_file',<br>  outdir='./tmp',<br>  nq1=2, nq2=2, nq3=2<br>  conv_thr_chi = 1.0d-8,<br>  iverbosity = 2<br> /<br><br>####################################################################################<br>After some houres running the program , I am faced with this error:<br><br>     factorization<br> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br><br>     stopping ...<br>     Error in routine broyden (3):<br>     stopping ...<br> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br><br> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br><br>     stopping ...<br><br><br>  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%      task #        69      from broyden : error #         3      factorization  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br><br><br>Could someone leave their comments on this error.<br>Any help would be appreciated.<br><br>Best<br>Mohaddeseh<br><br><br><br><br><br clear="all"><br>-- <br><div dir="ltr">---------------------------------------------------------<br>Mohaddeseh Abbasnejad, <br>Assistant Professor of Physics,<br>Faculty of Physics, <br>Shahid Bahonar University of Kerman,<br>Kerman, Iran<br>P.O. Box 76169-133<br>Tel: +98 34 31322199<br>Fax: +98 34 33257434<br>Cellphone: +98 917 731 7514<br>E-Mail:     <a href="mailto:mohaddeseh.abbasnejad@gmail.com" target="_blank">mohaddeseh.abbasnejad@gmail.com</a><br>Website:  <a href="http://academicstaff.uk.ac.ir/moabbasnejad" target="_blank">academicstaff.uk.ac.ir/moabbasnejad</a> <br>---------------------------------------------------------<br><br></div></div></div>