<div dir="ltr"><div>Tetrahedra should work properly for any symmetry, but may be subject to the usual problem of "quasi-symmetry": if the system is close to a given symmetry but not really symmetric, one may encounter all kinds of funny failures.</div><div><br></div><div>Paolo<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Feb 5, 2021 at 4:28 PM Lorenzo Paulatto <<a href="mailto:paulatz@gmail.com">paulatz@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div>I'm not 100% sure this is the cause, but tetraheadra may not work properly for non-cubic systems, I would try occupations="smearing" unless you really want tetraheadra for some reason.</div><br><div><u></u>--<br><span style="font-size:0.9em;color:gray">Lorenzo Paulatto - Paris<span></span></span><u></u></div><div>On Feb 5 2021, at 3:54 pm, Mayuri Bora <<a href="mailto:mayu219@tezu.ernet.in" target="_blank">mayu219@tezu.ernet.in</a>> wrote:</div><blockquote><div><div>Dear all,</div><br><div>I have been trying the calculation metal oxide system where scf have been</div><div>converged but while calculating the nscf calculation i am getting error as</div><div>follows-</div><div>Subspace diagonalization in iterative solution of the eigenvalue problem:</div><div>a serial algorithm will be used</div><br><div>%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</div><div>Error in routine tetrahedra (366):</div><div>cannot remap grid on k-point list</div><div>%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</div><br><div>However, the input is as follows-</div><div>&CONTROL</div><div>calculation='nscf',</div><div>outdir='/global/home/pritam/mayuri/MNO/new',</div><div>prefix='MnO',</div><div>pseudo_dir='/global/home/pritam/mayuri/MNO/new',</div><div>verbosity='high',</div><div>tstress=.true.,</div><div>tprnfor=.true.,</div><div>/</div><br><div>&SYSTEM</div><div>ibrav=14,</div><div>celldm(1)=11.9917635700d0, celldm(2)=1.0000000000d0,</div><div>celldm(3)=1.0000000000d0,</div><div>celldm(4)=-0.5719136861d0, celldm(5)=-0.5719136861d0,</div><div>celldm(6)=0.1438273409d0,</div><div>nat=7,</div><div>ntyp=2,</div><div>nspin=1,</div><div>ecutwfc=30,</div><div>ecutrho=120,</div><div>input_dft='PBE',</div><div>occupations='tetrahedra',</div><div>smearing='gaussian',</div><div>starting_magnetization(1)=0.16,</div><div>degauss=0.005d0,</div><div>/</div><br><div>&ELECTRONS</div><div>diagonalization= 'david',</div><div>electron_maxstep= 200,</div><div>conv_thr=1d-06,</div><div>mixing_beta=0.7d0,</div><div>/</div><br><div>ATOMIC_SPECIES</div><div>Mn 54.938000d0 Mn.pbe-spn-kjpaw_psl.0.3.1.UPF</div><div>O 15.999400d0 O.pbe-n-kjpaw_psl.1.0.0.UPF</div><br><div>&CONTROL</div><div>calculation='nscf',</div><div>outdir='/global/home/pritam/mayuri/MNO/new',</div><div>prefix='MnO',</div><div>pseudo_dir='/global/home/pritam/mayuri/MNO/new',</div><div>verbosity='high',</div><div>tstress=.true.,</div><div>tprnfor=.true.,</div><div>/</div><br><div>&SYSTEM</div><div>ibrav=14,</div><div>celldm(1)=11.9917635700d0, celldm(2)=1.0000000000d0,</div><div>celldm(3)=1.0000000000d0,</div><div>celldm(4)=-0.5719136861d0, celldm(5)=-0.5719136861d0,</div><div>celldm(6)=0.1438273409d0,</div><div>nat=7,</div><div>ntyp=2,</div><div>nspin=1,</div><div>ecutwfc=30,</div><div>ecutrho=120,</div><div>input_dft='PBE',</div><div>occupations='tetrahedra',</div><div>smearing='gaussian',</div><div>starting_magnetization(1)=0.16,</div><div>degauss=0.005d0,</div><div>/</div><br><div>&ELECTRONS</div><div>diagonalization= 'david',</div><div>electron_maxstep= 200,</div><div>conv_thr=1d-06,</div><div>mixing_beta=0.7d0,</div><div>/</div><br><div>ATOMIC_SPECIES</div><div>Mn 54.938000d0 Mn.pbe-spn-kjpaw_psl.0.3.1.UPF</div><div>O 15.999400d0 O.pbe-n-kjpaw_psl.1.0.0.UPF</div><br><div>ATOMIC_POSITIONS {crystal}</div><div>Mn 0.2615056310d0 0.7591129683d0 0.5458583547d0</div><div>Mn 0.5729298910d0 0.4364554765d0 0.0011162978d0</div><div>Mn 0.8659373932d0 0.6702645038d0 0.7259815947d0</div><div>O 0.6299195554d0 0.8901722878d0 0.7607555811d0</div><div>O 0.0627574912d0 0.4004316641d0 0.1974996665d0</div><div>O 0.6173271969d0 0.3686682914d0 0.2785206979d0</div><div>O 0.1558038413d0 0.8323428081d0 0.7989998074d0</div><br><div>K_POINTS {automatic}</div><div>27 27 1 0 0 0</div><br><div>regards</div><div>Mayuri</div><br><br><br><div>> Send users mailing list submissions to</div><div>> <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a></div><div>></div><div>> To subscribe or unsubscribe via the World Wide Web, visit</div><div>> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></div><div>> or, via email, send a message with subject or body 'help' to</div><div>> <a href="mailto:users-request@lists.quantum-espresso.org" target="_blank">users-request@lists.quantum-espresso.org</a></div><div>></div><div>> You can reach the person managing the list at</div><div>> <a href="mailto:users-owner@lists.quantum-espresso.org" target="_blank">users-owner@lists.quantum-espresso.org</a></div><div>></div><div>> When replying, please edit your Subject line so it is more specific</div><div>> than "Re: Contents of users digest..."</div><div>></div><div>></div><div>> Today's Topics:</div><div>></div><div>> 1. Iteration, nstep and multiple convergence</div><div>> (Hemant Verma (M19PH015))</div><div>> 2. Re: Relax calculation of strontium hexaferrite using QE</div><div>> (ABDELHAMID AIT M\'HID)</div><div>> 3. Re: Iteration, nstep and multiple convergence (Lorenzo Paulatto)</div><div>> 4. Re: Relax calculation of strontium hexaferrite using QE</div><div>> (Paolo Giannozzi)</div><div>> 5. Re: Relax calculation of strontium hexaferrite using QE</div><div>> (Lorenzo Paulatto)</div><div>> 6. Re: Relax calculation of strontium hexaferrite using QE</div><div>> (Dr. SUNIL KUMAR)</div><div>> 7. Re: Iteration, nstep and multiple convergence</div><div>> (Hemant Verma (M19PH015))</div><div>> 8. Re: Iteration, nstep and multiple convergence (Lorenzo Paulatto)</div><div>> 9. Converging SCF for adsorbed molecules with EXX (Guido Fratesi)</div><div>></div><div>></div><div>> ----------------------------------------------------------------------</div><div>></div><div>> Message: 1</div><div>> Date: Thu, 4 Feb 2021 19:49:35 +0530</div><div>> From: "Hemant Verma (M19PH015)" <<a href="mailto:verma.12@iitj.ac.in" target="_blank">verma.12@iitj.ac.in</a>></div><div>> To: Quantum ESPRESSO users Forum <<a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a>></div><div>> Subject: [QE-users] Iteration, nstep and multiple convergence</div><div>> Message-ID:</div><div>> <CAJoHbNf38q=uJ-RV=<a href="mailto:E_fJuLdEfAGnwz2Y6RHZjyUM0pzZCvbCA@mail.gmail.com" target="_blank">E_fJuLdEfAGnwz2Y6RHZjyUM0pzZCvbCA@mail.gmail.com</a>></div><div>> Content-Type: text/plain; charset="utf-8"</div><div>></div><div>> Dear Developers and users</div><div>></div><div>> I have a query regarding the iteration to achieve convergence and the</div><div>> option "nstep" in any calculation using QE. Are they same or different</div><div>> in their role?</div><div>> If they are different, please elaborate.</div><div>> In 'vc-relax' calculation we observe, from the output file, that the code</div><div>> runs several iterations and achieves convergence. However, it again starts</div><div>> iterations and achieves the convergence again. I am unable to understand</div><div>> why it keeps on achieving the convergence multiple times.</div><div>> Thank you in advance.</div><div>></div><div>> Regards</div><div>> Hemant Verma</div><div>> IIT Jodhpur</div><div>> -------------- next part --------------</div><div>> An HTML attachment was scrubbed...</div><div>> URL:</div><div>> <<a href="http://lists.quantum-espresso.org/pipermail/users/attachments/20210204/a8840a7d/attachment-0001.html" target="_blank">http://lists.quantum-espresso.org/pipermail/users/attachments/20210204/a8840a7d/attachment-0001.html</a>></div><div>></div><div>> ------------------------------</div><div>></div><div>> Message: 2</div><div>> Date: Thu, 4 Feb 2021 15:55:42 +0100</div><div>> From: "ABDELHAMID AIT M\\'HID" <<a href="mailto:a.aitmhid@edu.umi.ac.ma" target="_blank">a.aitmhid@edu.umi.ac.ma</a>></div><div>> To: Quantum ESPRESSO users Forum <<a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a>></div><div>> Subject: Re: [QE-users] Relax calculation of strontium hexaferrite</div><div>> using QE</div><div>> Message-ID:</div><div>> <CAC1VTbRgqT9beX1JSrsSmeXyciSO=<a href="mailto:7HOLhdMK5cTLxyXaZG35g@mail.gmail.com" target="_blank">7HOLhdMK5cTLxyXaZG35g@mail.gmail.com</a>></div><div>> Content-Type: text/plain; charset="utf-8"</div><div>></div><div>> Dear SUNIL,</div><div>> convergence issues are often related to the BAD CHOICE of ecutcfc and/Or</div><div>> K_POINTS. I suggest you increase ecutwfc it seems too small and for</div><div>> ecutrho</div><div>> let it be around ecutwfc* 8~12 since you use ultrasoft pseudopotential.</div><div>> Also, increase K-POINTS try 4*4*3. another remarque is that since ibrav=4,</div><div>> a and c seems to be bad description of cell parameters try to use</div><div>> celldm(1)</div><div>> (in bohr) and celldm(3) instead (read pw.x input description)</div><div>> I hope this solves the problem</div><div>> best wishes.</div><div>> Abdelhamid</div><div>></div><div>></div><div>></div><div>> Le mar. 2 f?vr. 2021 ? 10:55, Dr. SUNIL KUMAR <<a href="mailto:suniliitd14@gmail.com" target="_blank">suniliitd14@gmail.com</a>> a</div><div>> ?crit :</div><div>></div><div>>> Dear Developers and users of QE.</div><div>>> I am struggling to optimize the structure of strontium hexaferrite durin</div><div>>> relax calculation of DFT using QE as given in the following script and</div><div>>> attached snapshot. It is taking too much time and send some unexpected</div><div>>> error as</div><div>>></div><div>>> convergence NOT achieved after 100 iterations: stoppingv</div><div>>> or</div><div>>> convergence NOT achieved after 200 iterations: stopping</div><div>>></div><div>>> Can anyone suggest me any modification in the following script.</div><div>>> I will be too much grateful.</div><div>>> Thanks</div><div>>> With regards</div><div>>> SUNIL</div><div>>></div><div>>> &CONTROL</div><div>>> calculation = "relax"</div><div>>> prefix = "espresso"</div><div>>> tprnfor = .TRUE.</div><div>>> tstress = .TRUE.</div><div>>> wf_collect = .TRUE.</div><div>>> verbosity = 'high'</div><div>>> disk_io = 'high'</div><div>>> forc_conv_thr = 1.00000e-03</div><div>>> nstep = 200</div><div>>> pseudo_dir = "/home/pseudopot"</div><div>>> outdir = "./outdir"</div><div>>> wfcdir = "./wfcdir"</div><div>>> /</div><div>>></div><div>>> &SYSTEM</div><div>>> a = 5.99774e+00</div><div>>> c = 2.33482e+01</div><div>>> degauss = 1.00000e-02</div><div>>> ecutrho = 250</div><div>>> ecutwfc = 25</div><div>>> ibrav = 4</div><div>>> nat = 64</div><div>>> nspin = 2</div><div>>> ntyp = 3</div><div>>> occupations = "smearing"</div><div>>> smearing = "gaussian"</div><div>>> starting_magnetization(1) = 0.00000e+00</div><div>>> starting_magnetization(2) = 0.2</div><div>>> starting_magnetization(3) = 0.00000e+00</div><div>>> lda_plus_u = .true.</div><div>>> lda_plus_u_kind = 0</div><div>>> Hubbard_U(2) = 4</div><div>>> Hubbard_J(1,1) = 0</div><div>>> Hubbard_J(2,2) = 1</div><div>>></div><div>>> /</div><div>>></div><div>>> &ELECTRONS</div><div>>> conv_thr = 1.00000e-06</div><div>>> electron_maxstep = 500</div><div>>> mixing_beta = 4.00000e-01</div><div>>> startingpot = "atomic"</div><div>>> startingwfc = "atomic+random"</div><div>>></div><div>>> /</div><div>>></div><div>>> &IONS</div><div>>> ion_dynamics = "bfgs"</div><div>>> /</div><div>>></div><div>>> &CELL</div><div>>> /</div><div>>></div><div>>></div><div>>> K_POINTS {automatic}</div><div>>> 3 3 1 0 0 0</div><div>>></div><div>>> ATOMIC_SPECIES</div><div>>> Sr 87.62000 Sr.pbe-nsp-van.UPF</div><div>>> Fe 55.84500 Fe.pbe-nd-rrkjus.UPF</div><div>>> O 15.99940 O.pbe-rrkjus.UPF</div><div>>></div><div>>> ATOMIC_POSITIONS {angstrom}</div><div>>> Sr -0.000000 3.462796 17.511161</div><div>>> Sr 2.998869 1.731398 5.837054</div><div>>> Fe 0.000000 0.000000 0.000000</div><div>>> Fe 0.000000 0.000000 11.674107</div><div>>> Fe 0.000000 0.000000 17.511161</div><div>>> Fe 0.000000 0.000000 5.837054</div><div>>> Fe 2.998869 1.731398 22.706489</div><div>>> Fe -0.000000 3.462796 0.641726</div><div>>> Fe -0.000000 3.462796 11.032382</div><div>>> Fe 2.998869 1.731398 12.315833</div><div>>> Fe 2.998869 1.731398 18.887702</div><div>>> Fe -0.000000 3.462796 4.460513</div><div>>> Fe -0.000000 3.462796 7.213594</div><div>>> Fe 2.998869 1.731398 16.134621</div><div>>> Fe 2.998869 3.447345 2.560692</div><div>>> Fe 1.512816 0.873424 2.560692</div><div>>> Fe 4.484923 0.873424 2.560692</div><div>>> Fe 0.000000 1.746849 20.787523</div><div>>> Fe 1.486054 4.320769 20.787523</div><div>>> Fe -1.486054 4.320769 20.787523</div><div>>> Fe 0.000000 1.746849 14.234800</div><div>>> Fe 1.486054 4.320769 14.234800</div><div>>> Fe -1.486054 4.320769 14.234800</div><div>>> Fe 2.998869 3.447345 9.113415</div><div>>> Fe 1.512816 0.873424 9.113415</div><div>>> Fe 4.484923 0.873424 9.113415</div><div>>> O 0.000000 0.000000 19.782055</div><div>>> O 0.000000 0.000000 3.566160</div><div>>> O 0.000000 0.000000 8.107948</div><div>>> O 0.000000 0.000000 15.240267</div><div>>> O 2.998869 1.731398 1.300869</div><div>>> O -0.000000 3.462796 22.047346</div><div>>> O -0.000000 3.462796 12.974977</div><div>>> O 2.998869 1.731398 10.373238</div><div>>> O 2.998869 3.294521 17.511161</div><div>>> O 1.645165 0.949836 17.511161</div><div>>> O 4.352574 0.949836 17.511161</div><div>>> O 0.000000 1.899672 5.837054</div><div>>> O 1.353705 4.244357 5.837054</div><div>>> O -1.353705 4.244357 5.837054</div><div>>> O 2.998869 3.580534 22.111857</div><div>>> O 1.397470 0.806830 22.111857</div><div>>> O 4.600268 0.806830 22.111857</div><div>>> O 0.000000 1.613659 1.236358</div><div>>> O 1.601399 4.387364 1.236358</div><div>>> O -1.601399 4.387364 1.236358</div><div>>> O 0.000000 1.613659 10.437749</div><div>>> O 1.601399 4.387364 10.437749</div><div>>> O -1.601399 4.387364 10.437749</div><div>>> O 2.998869 3.580534 12.910466</div><div>>> O 1.397470 0.806830 12.910466</div><div>>> O 4.600268 0.806830 12.910466</div><div>>> O 0.000000 5.134086 19.813669</div><div>>> O -1.447380 2.627150 19.813669</div><div>>> O 1.447380 2.627150 19.813669</div><div>>> O 2.998869 0.060107 3.534546</div><div>>> O 4.446249 2.567043 3.534546</div><div>>> O 1.551489 2.567043 3.534546</div><div>>> O 2.998869 0.060107 8.139561</div><div>>> O 4.446249 2.567043 8.139561</div><div>>> O 1.551489 2.567043 8.139561</div><div>>> O 0.000000 5.134086 15.208654</div><div>>> O -1.447380 2.627150 15.208654</div><div>>> O 1.447380 2.627150 15.208654</div><div>>></div><div>>></div><div>>> Dr. Sunil Kumar</div><div>>> Ph.D (Chemical Engg. IIT Delhi)</div><div>>> M.Tech (Chemical Engg. IIT Delhi)</div><div>>> B.Tech (Chemical Engg. IET-CSJMU Kanpur)</div><div>>> Scientist-C and Assistant Professor</div><div>>> CSIR-National Metallurgical Laboratory Jamshedpur-831007</div><div>>> <a href="http://www.nmlindia.org/" target="_blank">http://www.nmlindia.org/</a></div><div>>> <a href="https://scholar.google.co.in/citations?user=OchYqugAAAAJ&hl=en&oi=sra" target="_blank">https://scholar.google.co.in/citations?user=OchYqugAAAAJ&hl=en&oi=sra</a></div><div>>></div><div>>> _______________________________________________</div><div>>> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" target="_blank">www.max-centre.eu</a>)</div><div>>> users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a></div><div>>> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></div><div>> -------------- next part --------------</div><div>> An HTML attachment was scrubbed...</div><div>> URL:</div><div>> <<a href="http://lists.quantum-espresso.org/pipermail/users/attachments/20210204/f53e7f73/attachment-0001.html" target="_blank">http://lists.quantum-espresso.org/pipermail/users/attachments/20210204/f53e7f73/attachment-0001.html</a>></div><div>></div><div>> ------------------------------</div><div>></div><div>> Message: 3</div><div>> Date: Thu, 4 Feb 2021 16:01:18 +0100</div><div>> From: Lorenzo Paulatto <<a href="mailto:paulatz@gmail.com" target="_blank">paulatz@gmail.com</a>></div><div>> To: Quantum ESPRESSO users Forum <<a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a>></div><div>> Subject: Re: [QE-users] Iteration, nstep and multiple convergence</div><div>> Message-ID: <<a href="mailto:BB1C6931-95FB-4051-8E75-6A63400DD4D3@getmailspring.com" target="_blank">BB1C6931-95FB-4051-8E75-6A63400DD4D3@getmailspring.com</a>></div><div>> Content-Type: text/plain; charset="utf-8"</div><div>></div><div>> Are you asking , why, after printing</div><div>> A final scf calculation at the relaxed structure.</div><div>> The G-vectors are recalculated for the final unit cell</div><div>> Results may differ from those at the preceding step.</div><div>> the code proceeds to do a final scf of the relaxed structure,</div><div>> recalculating the G-vectors for the final unit cell? It is because the</div><div>> results may differ from those at the preceding step.</div><div>></div><div>> regards</div><div>> --</div><div>> Lorenzo Paulatto - Paris</div><div>> On Feb 4 2021, at 3:19 pm, Hemant Verma (M19PH015) <<a href="mailto:verma.12@iitj.ac.in" target="_blank">verma.12@iitj.ac.in</a>></div><div>> wrote:</div><div>>> Dear Developers and users</div><div>>></div><div>>> I have a query regarding the iteration to achieve convergence and the</div><div>>> option "nstep" in any calculation using QE. Are they same or different</div><div>>> in their role?</div><div>>> If they are different, please elaborate.</div><div>>> In 'vc-relax' calculation we observe, from the output file, that the</div><div>>> code runs several iterations and achieves convergence. However, it again</div><div>>> starts iterations and achieves the convergence again. I am unable to</div><div>>> understand why it keeps on achieving the convergence multiple times.</div><div>>> Thank you in advance.</div><div>>></div><div>>> Regards</div><div>>> Hemant Verma</div><div>>> IIT Jodhpur</div><div>>></div><div>>> _______________________________________________</div><div>>> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" target="_blank">www.max-centre.eu</a>)</div><div>>> users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a></div><div>>> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></div><div>></div><div>> -------------- next part --------------</div><div>> An HTML attachment was scrubbed...</div><div>> URL:</div><div>> <<a href="http://lists.quantum-espresso.org/pipermail/users/attachments/20210204/9ece23f9/attachment-0001.html" target="_blank">http://lists.quantum-espresso.org/pipermail/users/attachments/20210204/9ece23f9/attachment-0001.html</a>></div><div>></div><div>> ------------------------------</div><div>></div><div>> Message: 4</div><div>> Date: Thu, 4 Feb 2021 16:01:23 +0100</div><div>> From: Paolo Giannozzi <<a href="mailto:p.giannozzi@gmail.com" target="_blank">p.giannozzi@gmail.com</a>></div><div>> To: Quantum ESPRESSO users Forum <<a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a>></div><div>> Subject: Re: [QE-users] Relax calculation of strontium hexaferrite</div><div>> using QE</div><div>> Message-ID:</div><div>> <<a href="mailto:CAPMgbCtrwNwauPcdu3An_f38VoqRUY08kyAKCHr3_a_XgSSanQ@mail.gmail.com" target="_blank">CAPMgbCtrwNwauPcdu3An_f38VoqRUY08kyAKCHr3_a_XgSSanQ@mail.gmail.com</a>></div><div>> Content-Type: text/plain; charset="utf-8"</div><div>></div><div>> On Thu, Feb 4, 2021 at 3:56 PM ABDELHAMID AIT M\'HID <</div><div>> <a href="mailto:a.aitmhid@edu.umi.ac.ma" target="_blank">a.aitmhid@edu.umi.ac.ma</a>> wrote:</div><div>></div><div>></div><div>>> since ibrav=4, a and c seems to be bad description of cell parameters</div><div>>> try</div><div>>> to use celldm(1) (in bohr) and celldm(3) instead (read pw.x input</div><div>>> description)</div><div>>></div><div>></div><div>> using a and c is 100% equivalent to using celldm(1) and celldm(3), as long</div><div>> as their values are correctly set in the correct units</div><div>></div><div>> Paolo</div><div>> --</div><div>> Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,</div><div>> Univ. Udine, via delle Scienze 206, 33100 Udine, Italy</div><div>> Phone +39-0432-558216, fax +39-0432-558222</div><div>> -------------- next part --------------</div><div>> An HTML attachment was scrubbed...</div><div>> URL:</div><div>> <<a href="http://lists.quantum-espresso.org/pipermail/users/attachments/20210204/dd59fa02/attachment-0001.html" target="_blank">http://lists.quantum-espresso.org/pipermail/users/attachments/20210204/dd59fa02/attachment-0001.html</a>></div><div>></div><div>> ------------------------------</div><div>></div><div>> Message: 5</div><div>> Date: Thu, 4 Feb 2021 16:34:16 +0100</div><div>> From: Lorenzo Paulatto <<a href="mailto:paulatz@gmail.com" target="_blank">paulatz@gmail.com</a>></div><div>> To: Quantum ESPRESSO users Forum <<a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a>></div><div>> Subject: Re: [QE-users] Relax calculation of strontium hexaferrite</div><div>> using QE</div><div>> Message-ID: <<a href="mailto:D6C4DD37-4274-4EE5-8BCF-32CAD7339162@getmailspring.com" target="_blank">D6C4DD37-4274-4EE5-8BCF-32CAD7339162@getmailspring.com</a>></div><div>> Content-Type: text/plain; charset="utf-8"</div><div>></div><div>> Converging the self-consistent field with DFT+U can be difficult. To</div><div>> start, reduce mixing_beta (see manual), and also remove disk_io='high'</div><div>> unless you have a reason to use it. Also, your values of the U parameters</div><div>> seem to be arbitrary, if they are absurd it is very likely that SCF will</div><div>> never converge, try to relax without the +U to begin, then find a first</div><div>> approximation of U (by first principles or fitting some physical</div><div>> parameter) and relax again</div><div>></div><div>> --</div><div>> Lorenzo Paulatto - Paris</div><div>> On Feb 2 2021, at 10:54 am, Dr. SUNIL KUMAR <<a href="mailto:suniliitd14@gmail.com" target="_blank">suniliitd14@gmail.com</a>> wrote:</div><div>>> Dear Developers and users of QE.</div><div>>> I am struggling to optimize the structure of strontium hexaferrite durin</div><div>>> relax calculation of DFT using QE as given in the following script and</div><div>>> attached snapshot. It is taking too much time and send some unexpected</div><div>>> error as</div><div>>></div><div>>> convergence NOT achieved after 100 iterations: stoppingv</div><div>>> or</div><div>>> convergence NOT achieved after 200 iterations: stopping</div><div>>></div><div>>></div><div>>> Can anyone suggest me any modification in the following script.</div><div>>> I will be too much grateful.</div><div>>> Thanks</div><div>>> With regards</div><div>>> SUNIL</div><div>>></div><div>>> &CONTROL</div><div>>> calculation = "relax"</div><div>>> prefix = "espresso"</div><div>>> tprnfor = .TRUE.</div><div>>> tstress = .TRUE.</div><div>>> wf_collect = .TRUE.</div><div>>> verbosity = 'high'</div><div>>> disk_io = 'high'</div><div>>> forc_conv_thr = 1.00000e-03</div><div>>> nstep = 200</div><div>>> pseudo_dir = "/home/pseudopot"</div><div>>> outdir = "./outdir"</div><div>>> wfcdir = "./wfcdir"</div><div>>> /</div><div>>></div><div>>> &SYSTEM</div><div>>> a = 5.99774e+00</div><div>>> c = 2.33482e+01</div><div>>> degauss = 1.00000e-02</div><div>>> ecutrho = 250</div><div>>> ecutwfc = 25</div><div>>> ibrav = 4</div><div>>> nat = 64</div><div>>> nspin = 2</div><div>>> ntyp = 3</div><div>>> occupations = "smearing"</div><div>>> smearing = "gaussian"</div><div>>> starting_magnetization(1) = 0.00000e+00</div><div>>> starting_magnetization(2) = 0.2</div><div>>> starting_magnetization(3) = 0.00000e+00</div><div>>> lda_plus_u = .true.</div><div>>> lda_plus_u_kind = 0</div><div>>> Hubbard_U(2) = 4</div><div>>> Hubbard_J(1,1) = 0</div><div>>> Hubbard_J(2,2) = 1</div><div>>></div><div>>> /</div><div>>> &ELECTRONS</div><div>>> conv_thr = 1.00000e-06</div><div>>> electron_maxstep = 500</div><div>>> mixing_beta = 4.00000e-01</div><div>>> startingpot = "atomic"</div><div>>> startingwfc = "atomic+random"</div><div>>></div><div>>> /</div><div>>> &IONS</div><div>>> ion_dynamics = "bfgs"</div><div>>> /</div><div>>></div><div>>> &CELL</div><div>>> /</div><div>>></div><div>>></div><div>>> K_POINTS {automatic}</div><div>>> 3 3 1 0 0 0</div><div>>></div><div>>> ATOMIC_SPECIES</div><div>>> Sr 87.62000 Sr.pbe-nsp-van.UPF</div><div>>> Fe 55.84500 Fe.pbe-nd-rrkjus.UPF</div><div>>> O 15.99940 O.pbe-rrkjus.UPF</div><div>>></div><div>>> ATOMIC_POSITIONS {angstrom}</div><div>>> Sr -0.000000 3.462796 17.511161</div><div>>> Sr 2.998869 1.731398 5.837054</div><div>>> Fe 0.000000 0.000000 0.000000</div><div>>> Fe 0.000000 0.000000 11.674107</div><div>>> Fe 0.000000 0.000000 17.511161</div><div>>> Fe 0.000000 0.000000 5.837054</div><div>>> Fe 2.998869 1.731398 22.706489</div><div>>> Fe -0.000000 3.462796 0.641726</div><div>>> Fe -0.000000 3.462796 11.032382</div><div>>> Fe 2.998869 1.731398 12.315833</div><div>>> Fe 2.998869 1.731398 18.887702</div><div>>> Fe -0.000000 3.462796 4.460513</div><div>>> Fe -0.000000 3.462796 7.213594</div><div>>> Fe 2.998869 1.731398 16.134621</div><div>>> Fe 2.998869 3.447345 2.560692</div><div>>> Fe 1.512816 0.873424 2.560692</div><div>>> Fe 4.484923 0.873424 2.560692</div><div>>> Fe 0.000000 1.746849 20.787523</div><div>>> Fe 1.486054 4.320769 20.787523</div><div>>> Fe -1.486054 4.320769 20.787523</div><div>>> Fe 0.000000 1.746849 14.234800</div><div>>> Fe 1.486054 4.320769 14.234800</div><div>>> Fe -1.486054 4.320769 14.234800</div><div>>> Fe 2.998869 3.447345 9.113415</div><div>>> Fe 1.512816 0.873424 9.113415</div><div>>> Fe 4.484923 0.873424 9.113415</div><div>>> O 0.000000 0.000000 19.782055</div><div>>> O 0.000000 0.000000 3.566160</div><div>>> O 0.000000 0.000000 8.107948</div><div>>> O 0.000000 0.000000 15.240267</div><div>>> O 2.998869 1.731398 1.300869</div><div>>> O -0.000000 3.462796 22.047346</div><div>>> O -0.000000 3.462796 12.974977</div><div>>> O 2.998869 1.731398 10.373238</div><div>>> O 2.998869 3.294521 17.511161</div><div>>> O 1.645165 0.949836 17.511161</div><div>>> O 4.352574 0.949836 17.511161</div><div>>> O 0.000000 1.899672 5.837054</div><div>>> O 1.353705 4.244357 5.837054</div><div>>> O -1.353705 4.244357 5.837054</div><div>>> O 2.998869 3.580534 22.111857</div><div>>> O 1.397470 0.806830 22.111857</div><div>>> O 4.600268 0.806830 22.111857</div><div>>> O 0.000000 1.613659 1.236358</div><div>>> O 1.601399 4.387364 1.236358</div><div>>> O -1.601399 4.387364 1.236358</div><div>>> O 0.000000 1.613659 10.437749</div><div>>> O 1.601399 4.387364 10.437749</div><div>>> O -1.601399 4.387364 10.437749</div><div>>> O 2.998869 3.580534 12.910466</div><div>>> O 1.397470 0.806830 12.910466</div><div>>> O 4.600268 0.806830 12.910466</div><div>>> O 0.000000 5.134086 19.813669</div><div>>> O -1.447380 2.627150 19.813669</div><div>>> O 1.447380 2.627150 19.813669</div><div>>> O 2.998869 0.060107 3.534546</div><div>>> O 4.446249 2.567043 3.534546</div><div>>> O 1.551489 2.567043 3.534546</div><div>>> O 2.998869 0.060107 8.139561</div><div>>> O 4.446249 2.567043 8.139561</div><div>>> O 1.551489 2.567043 8.139561</div><div>>> O 0.000000 5.134086 15.208654</div><div>>> O -1.447380 2.627150 15.208654</div><div>>> O 1.447380 2.627150 15.208654</div><div>>></div><div>>></div><div>>> Dr. Sunil Kumar</div><div>>> Ph.D (Chemical Engg. IIT Delhi)</div><div>>> M.Tech (Chemical Engg. IIT Delhi)</div><div>>></div><div>>> B.Tech (Chemical Engg. IET-CSJMU Kanpur)</div><div>>> Scientist-C and Assistant Professor</div><div>>> CSIR-National Metallurgical Laboratory Jamshedpur-831007</div><div>>> <a href="http://www.nmlindia.org/" target="_blank">http://www.nmlindia.org/</a></div><div>>></div><div>>> <a href="https://scholar.google.co.in/citations?user=OchYqugAAAAJ&hl=en&oi=sra" target="_blank">https://scholar.google.co.in/citations?user=OchYqugAAAAJ&hl=en&oi=sra</a></div><div>>></div><div>>></div><div>>> _______________________________________________</div><div>>> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" target="_blank">www.max-centre.eu</a>)</div><div>>> users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a></div><div>>> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></div><div>></div><div>> -------------- next part --------------</div><div>> An HTML attachment was scrubbed...</div><div>> URL:</div><div>> <<a href="http://lists.quantum-espresso.org/pipermail/users/attachments/20210204/5f745cf5/attachment-0001.html" target="_blank">http://lists.quantum-espresso.org/pipermail/users/attachments/20210204/5f745cf5/attachment-0001.html</a>></div><div>></div><div>> ------------------------------</div><div>></div><div>> Message: 6</div><div>> Date: Thu, 4 Feb 2021 21:08:39 +0530</div><div>> From: "Dr. SUNIL KUMAR" <<a href="mailto:suniliitd14@gmail.com" target="_blank">suniliitd14@gmail.com</a>></div><div>> To: Quantum ESPRESSO users Forum <<a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a>></div><div>> Subject: Re: [QE-users] Relax calculation of strontium hexaferrite</div><div>> using QE</div><div>> Message-ID:</div><div>> <CADe9O7+4XugcwhX=<a href="mailto:QQwEqQ2EZGckJSHZfQ0BS_WktL_W-QOsiA@mail.gmail.com" target="_blank">QQwEqQ2EZGckJSHZfQ0BS_WktL_W-QOsiA@mail.gmail.com</a>></div><div>> Content-Type: text/plain; charset="utf-8"</div><div>></div><div>> Thankyou for your suggestions..</div><div>> i am too much grateful to you all</div><div>></div><div>> Dr. Sunil Kumar</div><div>> Ph.D (Chemical Engg. IIT Delhi)</div><div>> M.Tech (Chemical Engg. IIT Delhi)</div><div>> B.Tech (Chemical Engg. IET-CSJMU Kanpur)</div><div>> Scientist-C and Assistant Professor</div><div>> CSIR-National Metallurgical Laboratory Jamshedpur-831007</div><div>> <a href="http://www.nmlindia.org/" target="_blank">http://www.nmlindia.org/</a></div><div>> <a href="https://scholar.google.co.in/citations?user=OchYqugAAAAJ&hl=en&oi=sra" target="_blank">https://scholar.google.co.in/citations?user=OchYqugAAAAJ&hl=en&oi=sra</a></div><div>></div><div>></div><div>></div><div>> On Thu, Feb 4, 2021 at 9:04 PM Lorenzo Paulatto <<a href="mailto:paulatz@gmail.com" target="_blank">paulatz@gmail.com</a>> wrote:</div><div>></div><div>>> Converging the self-consistent field with DFT+U can be difficult. To</div><div>>> start, reduce mixing_beta (see manual), and also remove disk_io='high'</div><div>>> unless you have a reason to use it. Also, your values of the U</div><div>>> parameters</div><div>>> seem to be arbitrary, if they are absurd it is very likely that SCF will</div><div>>> never converge, try to relax without the +U to begin, then find a first</div><div>>> approximation of U (by first principles or fitting some physical</div><div>>> parameter)</div><div>>> and relax again</div><div>>></div><div>>> --</div><div>>> Lorenzo Paulatto - Paris</div><div>>> On Feb 2 2021, at 10:54 am, Dr. SUNIL KUMAR <<a href="mailto:suniliitd14@gmail.com" target="_blank">suniliitd14@gmail.com</a>></div><div>>> wrote:</div><div>>></div><div>>> Dear Developers and users of QE.</div><div>>> I am struggling to optimize the structure of strontium hexaferrite durin</div><div>>> relax calculation of DFT using QE as given in the following script and</div><div>>> attached snapshot. It is taking too much time and send some unexpected</div><div>>> error as</div><div>>></div><div>>> convergence NOT achieved after 100 iterations: stoppingv</div><div>>> or</div><div>>> convergence NOT achieved after 200 iterations: stopping</div><div>>></div><div>>> Can anyone suggest me any modification in the following script.</div><div>>> I will be too much grateful.</div><div>>> Thanks</div><div>>> With regards</div><div>>> SUNIL</div><div>>></div><div>>> &CONTROL</div><div>>> calculation = "relax"</div><div>>> prefix = "espresso"</div><div>>> tprnfor = .TRUE.</div><div>>> tstress = .TRUE.</div><div>>> wf_collect = .TRUE.</div><div>>> verbosity = 'high'</div><div>>> disk_io = 'high'</div><div>>> forc_conv_thr = 1.00000e-03</div><div>>> nstep = 200</div><div>>> pseudo_dir = "/home/pseudopot"</div><div>>> outdir = "./outdir"</div><div>>> wfcdir = "./wfcdir"</div><div>>> /</div><div>>></div><div>>> &SYSTEM</div><div>>> a = 5.99774e+00</div><div>>> c = 2.33482e+01</div><div>>> degauss = 1.00000e-02</div><div>>> ecutrho = 250</div><div>>> ecutwfc = 25</div><div>>> ibrav = 4</div><div>>> nat = 64</div><div>>> nspin = 2</div><div>>> ntyp = 3</div><div>>> occupations = "smearing"</div><div>>> smearing = "gaussian"</div><div>>> starting_magnetization(1) = 0.00000e+00</div><div>>> starting_magnetization(2) = 0.2</div><div>>> starting_magnetization(3) = 0.00000e+00</div><div>>> lda_plus_u = .true.</div><div>>> lda_plus_u_kind = 0</div><div>>> Hubbard_U(2) = 4</div><div>>> Hubbard_J(1,1) = 0</div><div>>> Hubbard_J(2,2) = 1</div><div>>></div><div>>> /</div><div>>></div><div>>> &ELECTRONS</div><div>>> conv_thr = 1.00000e-06</div><div>>> electron_maxstep = 500</div><div>>> mixing_beta = 4.00000e-01</div><div>>> startingpot = "atomic"</div><div>>> startingwfc = "atomic+random"</div><div>>></div><div>>> /</div><div>>></div><div>>> &IONS</div><div>>> ion_dynamics = "bfgs"</div><div>>> /</div><div>>></div><div>>> &CELL</div><div>>> /</div><div>>></div><div>>></div><div>>> K_POINTS {automatic}</div><div>>> 3 3 1 0 0 0</div><div>>></div><div>>> ATOMIC_SPECIES</div><div>>> Sr 87.62000 Sr.pbe-nsp-van.UPF</div><div>>> Fe 55.84500 Fe.pbe-nd-rrkjus.UPF</div><div>>> O 15.99940 O.pbe-rrkjus.UPF</div><div>>></div><div>>> ATOMIC_POSITIONS {angstrom}</div><div>>> Sr -0.000000 3.462796 17.511161</div><div>>> Sr 2.998869 1.731398 5.837054</div><div>>> Fe 0.000000 0.000000 0.000000</div><div>>> Fe 0.000000 0.000000 11.674107</div><div>>> Fe 0.000000 0.000000 17.511161</div><div>>> Fe 0.000000 0.000000 5.837054</div><div>>> Fe 2.998869 1.731398 22.706489</div><div>>> Fe -0.000000 3.462796 0.641726</div><div>>> Fe -0.000000 3.462796 11.032382</div><div>>> Fe 2.998869 1.731398 12.315833</div><div>>> Fe 2.998869 1.731398 18.887702</div><div>>> Fe -0.000000 3.462796 4.460513</div><div>>> Fe -0.000000 3.462796 7.213594</div><div>>> Fe 2.998869 1.731398 16.134621</div><div>>> Fe 2.998869 3.447345 2.560692</div><div>>> Fe 1.512816 0.873424 2.560692</div><div>>> Fe 4.484923 0.873424 2.560692</div><div>>> Fe 0.000000 1.746849 20.787523</div><div>>> Fe 1.486054 4.320769 20.787523</div><div>>> Fe -1.486054 4.320769 20.787523</div><div>>> Fe 0.000000 1.746849 14.234800</div><div>>> Fe 1.486054 4.320769 14.234800</div><div>>> Fe -1.486054 4.320769 14.234800</div><div>>> Fe 2.998869 3.447345 9.113415</div><div>>> Fe 1.512816 0.873424 9.113415</div><div>>> Fe 4.484923 0.873424 9.113415</div><div>>> O 0.000000 0.000000 19.782055</div><div>>> O 0.000000 0.000000 3.566160</div><div>>> O 0.000000 0.000000 8.107948</div><div>>> O 0.000000 0.000000 15.240267</div><div>>> O 2.998869 1.731398 1.300869</div><div>>> O -0.000000 3.462796 22.047346</div><div>>> O -0.000000 3.462796 12.974977</div><div>>> O 2.998869 1.731398 10.373238</div><div>>> O 2.998869 3.294521 17.511161</div><div>>> O 1.645165 0.949836 17.511161</div><div>>> O 4.352574 0.949836 17.511161</div><div>>> O 0.000000 1.899672 5.837054</div><div>>> O 1.353705 4.244357 5.837054</div><div>>> O -1.353705 4.244357 5.837054</div><div>>> O 2.998869 3.580534 22.111857</div><div>>> O 1.397470 0.806830 22.111857</div><div>>> O 4.600268 0.806830 22.111857</div><div>>> O 0.000000 1.613659 1.236358</div><div>>> O 1.601399 4.387364 1.236358</div><div>>> O -1.601399 4.387364 1.236358</div><div>>> O 0.000000 1.613659 10.437749</div><div>>> O 1.601399 4.387364 10.437749</div><div>>> O -1.601399 4.387364 10.437749</div><div>>> O 2.998869 3.580534 12.910466</div><div>>> O 1.397470 0.806830 12.910466</div><div>>> O 4.600268 0.806830 12.910466</div><div>>> O 0.000000 5.134086 19.813669</div><div>>> O -1.447380 2.627150 19.813669</div><div>>> O 1.447380 2.627150 19.813669</div><div>>> O 2.998869 0.060107 3.534546</div><div>>> O 4.446249 2.567043 3.534546</div><div>>> O 1.551489 2.567043 3.534546</div><div>>> O 2.998869 0.060107 8.139561</div><div>>> O 4.446249 2.567043 8.139561</div><div>>> O 1.551489 2.567043 8.139561</div><div>>> O 0.000000 5.134086 15.208654</div><div>>> O -1.447380 2.627150 15.208654</div><div>>> O 1.447380 2.627150 15.208654</div><div>>></div><div>>></div><div>>> Dr. Sunil Kumar</div><div>>> Ph.D (Chemical Engg. IIT Delhi)</div><div>>> M.Tech (Chemical Engg. IIT Delhi)</div><div>>> B.Tech (Chemical Engg. IET-CSJMU Kanpur)</div><div>>> Scientist-C and Assistant Professor</div><div>>> CSIR-National Metallurgical Laboratory Jamshedpur-831007</div><div>>> <a href="http://www.nmlindia.org/" target="_blank">http://www.nmlindia.org/</a></div><div>>> <a href="https://scholar.google.co.in/citations?user=OchYqugAAAAJ&hl=en&oi=sra" target="_blank">https://scholar.google.co.in/citations?user=OchYqugAAAAJ&hl=en&oi=sra</a></div><div>>></div><div>>> _______________________________________________</div><div>>> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" target="_blank">www.max-centre.eu</a>)</div><div>>> users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a></div><div>>> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></div><div>>></div><div>>> _______________________________________________</div><div>>> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" target="_blank">www.max-centre.eu</a>)</div><div>>> users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a></div><div>>> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></div><div>> -------------- next part --------------</div><div>> An HTML attachment was scrubbed...</div><div>> URL:</div><div>> <<a href="http://lists.quantum-espresso.org/pipermail/users/attachments/20210204/7fa597a8/attachment-0001.html" target="_blank">http://lists.quantum-espresso.org/pipermail/users/attachments/20210204/7fa597a8/attachment-0001.html</a>></div><div>></div><div>> ------------------------------</div><div>></div><div>> Message: 7</div><div>> Date: Thu, 4 Feb 2021 22:06:44 +0530</div><div>> From: "Hemant Verma (M19PH015)" <<a href="mailto:verma.12@iitj.ac.in" target="_blank">verma.12@iitj.ac.in</a>></div><div>> To: Quantum ESPRESSO users Forum <<a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a>></div><div>> Subject: Re: [QE-users] Iteration, nstep and multiple convergence</div><div>> Message-ID:</div><div>> <<a href="mailto:CAJoHbNeaSQmhnhpa0Sefb9qNwEE41fR0gCWJsX%2BOe89P%2BAUmmA@mail.gmail.com" target="_blank">CAJoHbNeaSQmhnhpa0Sefb9qNwEE41fR0gCWJsX+Oe89P+AUmmA@mail.gmail.com</a>></div><div>> Content-Type: text/plain; charset="utf-8"</div><div>></div><div>> Dear Dr. Paulatto,</div><div>></div><div>> Thank you for your prompt response. Can you suggest any reference material</div><div>> to read in detail regarding this?</div><div>></div><div>> Regards</div><div>> Hemant</div><div>></div><div>> On Thu, Feb 4, 2021 at 8:31 PM Lorenzo Paulatto <<a href="mailto:paulatz@gmail.com" target="_blank">paulatz@gmail.com</a>> wrote:</div><div>></div><div>>> Are you asking , why, after printing</div><div>>> A final scf calculation at the relaxed structure.</div><div>>> The G-vectors are recalculated for the final unit cell</div><div>>> Results may differ from those at the preceding step.</div><div>>> the code proceeds to do a final scf of the relaxed structure,</div><div>>> recalculating the G-vectors for the final unit cell? It is because the</div><div>>> results may differ from those at the preceding step.</div><div>>></div><div>>> regards</div><div>>> --</div><div>>> Lorenzo Paulatto - Paris</div><div>>> On Feb 4 2021, at 3:19 pm, Hemant Verma (M19PH015) <<a href="mailto:verma.12@iitj.ac.in" target="_blank">verma.12@iitj.ac.in</a>></div><div>>> wrote:</div><div>>></div><div>>> Dear Developers and users</div><div>>></div><div>>> I have a query regarding the iteration to achieve convergence and the</div><div>>> option "nstep" in any calculation using QE. Are they same or different</div><div>>> in their role?</div><div>>> If they are different, please elaborate.</div><div>>> In 'vc-relax' calculation we observe, from the output file, that the</div><div>>> code</div><div>>> runs several iterations and achieves convergence. However, it again</div><div>>> starts</div><div>>> iterations and achieves the convergence again. I am unable to understand</div><div>>> why it keeps on achieving the convergence multiple times.</div><div>>> Thank you in advance.</div><div>>></div><div>>> Regards</div><div>>> Hemant Verma</div><div>>> IIT Jodhpur</div><div>>> _______________________________________________</div><div>>> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" target="_blank">www.max-centre.eu</a>)</div><div>>> users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a></div><div>>> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></div><div>>></div><div>>> _______________________________________________</div><div>>> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" target="_blank">www.max-centre.eu</a>)</div><div>>> users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a></div><div>>> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></div><div>> -------------- next part --------------</div><div>> An HTML attachment was scrubbed...</div><div>> URL:</div><div>> <<a href="http://lists.quantum-espresso.org/pipermail/users/attachments/20210204/c3df2aec/attachment-0001.html" target="_blank">http://lists.quantum-espresso.org/pipermail/users/attachments/20210204/c3df2aec/attachment-0001.html</a>></div><div>></div><div>> ------------------------------</div><div>></div><div>> Message: 8</div><div>> Date: Thu, 4 Feb 2021 19:58:52 +0100</div><div>> From: Lorenzo Paulatto <<a href="mailto:paulatz@gmail.com" target="_blank">paulatz@gmail.com</a>></div><div>> To: Quantum ESPRESSO users Forum <<a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a>></div><div>> Subject: Re: [QE-users] Iteration, nstep and multiple convergence</div><div>> Message-ID: <<a href="mailto:2D86ED9B-17E5-4437-BDD5-2F546859606B@getmailspring.com" target="_blank">2D86ED9B-17E5-4437-BDD5-2F546859606B@getmailspring.com</a>></div><div>> Content-Type: text/plain; charset="utf-8"</div><div>></div><div>> It is more or less explained here:</div><div>> <a href="https://www.quantum-espresso.org/resources/faq/self-consistency#6.11" target="_blank">https://www.quantum-espresso.org/resources/faq/self-consistency#6.11</a></div><div>></div><div>> --</div><div>> Lorenzo Paulatto - Paris</div><div>> On Feb 4 2021, at 5:36 pm, Hemant Verma (M19PH015) <<a href="mailto:verma.12@iitj.ac.in" target="_blank">verma.12@iitj.ac.in</a>></div><div>> wrote:</div><div>>> Dear Dr. Paulatto,</div><div>>></div><div>>> Thank you for your prompt response. Can you suggest any reference</div><div>>> material to read in detail regarding this?</div><div>>></div><div>>> Regards</div><div>>> Hemant</div><div>>></div><div>>></div><div>>></div><div>>> On Thu, Feb 4, 2021 at 8:31 PM Lorenzo Paulatto <<a href="mailto:paulatz@gmail.com" target="_blank">paulatz@gmail.com</a></div><div>>> (mailto:<a href="mailto:paulatz@gmail.com" target="_blank">paulatz@gmail.com</a>)> wrote:</div><div>>> > Are you asking , why, after printing</div><div>>> > A final scf calculation at the relaxed structure.</div><div>>> > The G-vectors are recalculated for the final unit cell</div><div>>> > Results may differ from those at the preceding step.</div><div>>> > the code proceeds to do a final scf of the relaxed structure,</div><div>>> recalculating the G-vectors for the final unit cell? It is because the</div><div>>> results may differ from those at the preceding step.</div><div>>> ></div><div>>> > regards</div><div>>> > --</div><div>>> > Lorenzo Paulatto - Paris</div><div>>> ></div><div>>> > On Feb 4 2021, at 3:19 pm, Hemant Verma (M19PH015)</div><div>>> <<a href="mailto:verma.12@iitj.ac.in" target="_blank">verma.12@iitj.ac.in</a> (mailto:<a href="mailto:verma.12@iitj.ac.in" target="_blank">verma.12@iitj.ac.in</a>)> wrote:</div><div>>> > > Dear Developers and users</div><div>>> > ></div><div>>> > > I have a query regarding the iteration to achieve convergence and</div><div>>> the option "nstep" in any calculation using QE. Are they same or</div><div>>> different in their role?</div><div>>> > > If they are different, please elaborate.</div><div>>> > > In 'vc-relax' calculation we observe, from the output file, that the</div><div>>> code runs several iterations and achieves convergence. However, it</div><div>>> again starts iterations and achieves the convergence again. I am</div><div>>> unable to understand why it keeps on achieving the convergence</div><div>>> multiple times.</div><div>>> > > Thank you in advance.</div><div>>> > ></div><div>>> > > Regards</div><div>>> > > Hemant Verma</div><div>>> > > IIT Jodhpur</div><div>>> > ></div><div>>> > > _______________________________________________</div><div>>> > > Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" target="_blank">www.max-centre.eu</a></div><div>>> (<a href="http://www.max-centre.eu" target="_blank">http://www.max-centre.eu</a>))</div><div>>> > > users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a></div><div>>> (mailto:<a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a>)</div><div>>> > > <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></div><div>>> ></div><div>>> > _______________________________________________</div><div>>> > Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" target="_blank">www.max-centre.eu</a></div><div>>> (<a href="http://www.max-centre.eu" target="_blank">http://www.max-centre.eu</a>))</div><div>>> > users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a></div><div>>> (mailto:<a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a>)</div><div>>> > <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></div><div>>></div><div>>></div><div>>> _______________________________________________</div><div>>> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" target="_blank">www.max-centre.eu</a>)</div><div>>> users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a></div><div>>> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></div><div>></div><div>> -------------- next part --------------</div><div>> An HTML attachment was scrubbed...</div><div>> URL:</div><div>> <<a href="http://lists.quantum-espresso.org/pipermail/users/attachments/20210204/d0cb392a/attachment-0001.html" target="_blank">http://lists.quantum-espresso.org/pipermail/users/attachments/20210204/d0cb392a/attachment-0001.html</a>></div><div>></div><div>> ------------------------------</div><div>></div><div>> Message: 9</div><div>> Date: Fri, 05 Feb 2021 10:33:15 +0100</div><div>> From: Guido Fratesi <<a href="mailto:guido.fratesi@unimi.it" target="_blank">guido.fratesi@unimi.it</a>></div><div>> To: <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a></div><div>> Subject: [QE-users] Converging SCF for adsorbed molecules with EXX</div><div>> Message-ID: <<a href="mailto:5830de7f-07a7-279b-3cb0-d964664c8565@unimi.it" target="_blank">5830de7f-07a7-279b-3cb0-d964664c8565@unimi.it</a>></div><div>> Content-Type: text/plain; charset=utf-8; format=flowed</div><div>></div><div>> Dear all,</div><div>></div><div>> I'm running a calculation for phthalocyanines on Al surface. On top of</div><div>> the standard PBE calculation (that exposed no difficulty), I would like</div><div>> to run a SCF with PBE0.</div><div>></div><div>> For the molecule in the gas phase, the PBE0-SCF also went smoothly and I</div><div>> recover literature results for the energy levels.</div><div>></div><div>> However, for the molecule adsorbed on 5 layers of Al(001), the SCF loop</div><div>> (phi-psi) is not converging. At least, "dexx" fluctuates around 0.002Ry</div><div>> for more than 10 iterations. In another run with 3 layers, no progress</div><div>> is made after 30 iterations, and sometimes (say 1/10) also the phi-scf</div><div>> loop is not converging after several iterations.</div><div>></div><div>> I switched to NCPP for the EXX calculation and kept the "m-p" smearing</div><div>> as used for the surface, with nqx=1 for EXX. Part of the input is below,</div><div>> full I/O can be retrieved at:</div><div>> <a href="https://drive.google.com/drive/folders/1yTYWGrfWFUTcSkl_pXYzbOVYT9dsrXim?usp=sharing" target="_blank">https://drive.google.com/drive/folders/1yTYWGrfWFUTcSkl_pXYzbOVYT9dsrXim?usp=sharing</a></div><div>></div><div>> To be noticed that the molecule (here, AlPh) may be spin-polarized (one</div><div>> electron in a pi orbital, in the gas phase) or not depending on the</div><div>> interaction with the substrate... Here I'm looking first to a</div><div>> spin-compensated case.</div><div>></div><div>> Thank you for your help,</div><div>></div><div>> Guido</div><div>></div><div>></div><div>> &CONTROL</div><div>> ? calculation = 'scf' ...</div><div>> /</div><div>> &SYSTEM</div><div>> ? ibrav = 6</div><div>> ? a=14.2751509299088549285</div><div>> ? c=22.0</div><div>> ? nat = 182</div><div>> ? ntyp = 4</div><div>> ? ecutwfc = 70</div><div>> ? ecutrho = 280</div><div>> ? occupations = "smearing",</div><div>> ? smearing??? = "m-p",</div><div>> ? degauss???? = 0.02D0,</div><div>> ? vdw_corr='grimme-d3'</div><div>> ? input_dft='PBE0'</div><div>> ? nqx1=1, ? nqx2=1, ? nqx3=1</div><div>> /</div><div>> &ELECTRONS</div><div>> ? conv_thr??? = 1.0d-7</div><div>> ??? mixing_mode = "local-TF"</div><div>> ??? mixing_beta = 0.7</div><div>> /</div><div>> ATOMIC_SPECIES</div><div>> C 12.0 C.0.0.UPF</div><div>> H? 1.0 H.nc2.UPF</div><div>> N 14.0 N.nc2.0.0.UPF</div><div>> Al 26? Al.vbc.UPF</div><div>> ATOMIC_POSITIONS {angstrom}</div><div>></div><div>></div><div>> --</div><div>> Guido Fratesi</div><div>></div><div>> Dipartimento di Fisica</div><div>> Universita` degli Studi di Milano</div><div>> Via Celoria 16, 20133 Milano, Italy</div><div>></div><div>></div><div>></div><div>> ------------------------------</div><div>></div><div>> Subject: Digest Footer</div><div>></div><div>> _______________________________________________</div><div>> users mailing list</div><div>> <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a></div><div>> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></div><div>></div><div>> ------------------------------</div><div>></div><div>> End of users Digest, Vol 163, Issue 5</div><div>> *************************************</div><div>></div><br><br><div>Mayuri Bora</div><div>INSPIRE Fellow</div><div>Advanced Functional Material Laboratory</div><div>Tezpur University</div><div>Napaam</div><div><a href="http://www.tezu.ernet.in/afml/" target="_blank">http://www.tezu.ernet.in/afml/</a></div><br><br><div>* * * D I S C L A I M E R * * *</div><div>This e-mail may contain privileged information and is intended solely for the individual named. If you are not the named addressee you should not disseminate, distribute or copy this e-mail. Please notify the sender immediately by e-mail if you have received this e-mail in error and destroy it from your system. Though considerable effort has been made to deliver error free e-mail messages but it can not be guaranteed to be secure or error-free as information could be intercepted, corrupted, lost, destroyed, delayed, or may contain viruses. The recipient must verify the integrity of this e-mail message.</div><br><div>_______________________________________________</div><div>Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" target="_blank">www.max-centre.eu</a>)</div><div>users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a></div><div><a href="https://lists.quantum-espresso.org/mailman/listinfo/users" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></div></div></blockquote>_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer" target="_blank">www.max-centre.eu</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div><br clear="all"><br>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 206, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div>