<div dir="ltr"><div dir="ltr"><div><div><font color="#000000">Dear Giuseppe,</font></div><div><font color="#000000"><br></font></div><div><font color="#000000">I can not thank you enough for your detailed suggestions and observations.</font></div><div><font color="#000000"><br></font></div><div><font color="#000000">For reference, the surface I am working on is based on the material presented in Energy Environ. Sci., 2013, 6, 148 (DOI: <a href="https://doi.org/10.1039/C2EE23355J" target="_blank">https://doi.org/10.1039/C2EE23355J</a>). I took the structural details from the Supporting information (<a href="https://doi.org/10.1039/C2EE23355J" target="_blank">https://doi.org/10.1039/C2EE23355J</a>) (Table 3, Page 4). </font></div><span class="gmail-im" style=""><font color="#000000"><div>I have constructed surface slabs from this structure, and comparable surface studies on solid-state electrolytes have been reported in publications such as Nanomaterials 2022, 12, 2795 (DOI: <a href="https://doi.org/10.3390/nano12162795" target="_blank">https://doi.org/10.3390/nano12162795</a>), Mater. Adv., 2022, 3, 3143−3150 (DOI: 10.1039/D1MA01246K), ACS Appl. Energy Mater. 2020, 3, 11061−11072 (DOI: <a href="https://dx.doi.org/10.1021/acsaem.0c02033" target="_blank">https://dx.doi.org/10.1021/acsaem.0c02033</a>), ACS Appl. Mater. Interfaces 2020, 12, 16350−16358 (DOI: <a href="https://dx.doi.org/10.1021/acsami.9b23019" target="_blank">https://dx.doi.org/10.1021/acsami.9b23019</a>), and Chem. Mater. 2018, 30, 3019−3027 (DOI: <a href="https://doi.org/10.1021/acs.chemmater.8b00649" target="_blank">https://doi.org/10.1021/acs.chemmater.8b00649</a>).</div><div><br></div><div>Despite implementing your recommendations, </div></font></span><div><font color="#000000">(a) including simplifying the setup by removing dipole corrections and focusing on the gamma point, I continue to face convergence issues. </font></div><div><font color="#000000">(b) I am using pseudopotentials from PseudoDojo, generated with the ONCVPSP code (version 3.3.0, 08/16/2017, D. R. Hamann), and I have raised the ecutwfc to 90 Ry, as per your suggestion. However, I still have not achieved convergence. </font></div><div><font color="#000000">(c) Additionally, I have tried lowering electron_maxstep and using scf_must_converge=.false., but without success, </font><span style="color:rgb(0,0,0)">and etot_conv_thr is not reached, </span><font color="#000000" style="">I am attaching a typical total energy convergence trend for ionic minimization using the BFGS algorithm. </font></div><div><font color="#000000" style="">(d) Additionally, I have tried using PAW and pseudopotentials from the SSSP library without any luck! convergence is not achieved in 200 scf steps.</font></div><span class="gmail-im" style=""><font color="#000000"><div style="">I greatly appreciate your continued guidance and would be grateful for any further insights you may have to help resolve this issue.</div><div style=""><br></div><div style="">Thank you again for your time and valuable assistance.</div><div style=""><br></div><div style="">Best regards,</div><div style="">Vinay</div></font></span></div><div><br></div><div><br style="color:rgb(0,0,0);font-family:arial,helvetica,sans-serif;font-size:16px"></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Sep 5, 2024 at 6:11 PM Giuseppe Mattioli <<a href="mailto:giuseppe.mattioli@ism.cnr.it">giuseppe.mattioli@ism.cnr.it</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"><br>
Dear Vinay Maithani<br>
Ok, let us say "unconventional amorphous structure" as you are sure  <br>
that it is meaningful, for the sake of discussion, even if I don't see  <br>
any pattern or repeated unit in the structure and I see strange P-S  <br>
"molecules" and direct Al-Li and Li-Li bonds which are unconventional.  <br>
Only for my curiosity, can you post a link to some publication where  <br>
this or a similar structure is investigated?<br>
<br>
However, the structure is so unconventional that it is possible that  <br>
you aren't able to obtain convergence, no matter what you try. To give  <br>
it a try I suggest that you:<br>
<br>
1) start it simpler (no dipole correction, gamma point as the  <br>
structure doesn't seem periodic).<br>
<br>
2) If you are using NC pseudopotentials (which ones?), then raise  <br>
ecutwfc as there are no Li pp in my experience that is well converged  <br>
at 50 Ry. You don't need to indicate ecutrho, instead, when you use NC  <br>
PPs.<br>
<br>
3) Don't specify nbnd when using occupation='smearing', the code  <br>
provised all the orbitals you need. Try Gaussian smearing, first, as  <br>
it is more stable.<br>
<br>
4) Try to reduce electron_maxstep to a much smaller value (if it  <br>
doesn't converge after 50-100 step it will probably not converge after  <br>
2000), and use this smaller value together with  <br>
scf_must_converge=.false. This will force the bfgs step on an  <br>
unconverged wavefunction, and can be very helpful to move your system  <br>
out of an unlucky initial structure with a few BFGS steps. However,  <br>
your structure is so unconventional that even such trick may fail.<br>
<br>
HTH<br>
Giuseppe<br>
<br>
<br>
Quoting Vinay Maithani <<a href="mailto:vinaymaithani21@kgpian.iitkgp.ac.in" target="_blank">vinaymaithani21@kgpian.iitkgp.ac.in</a>>:<br>
<br>
> Dear Giuseppe,<br>
><br>
> Thank you for your feedback. I would like to clarify that this is an<br>
> *asymmetric<br>
> polar slab* structure of a *solid-state electrolyte*. I have thoroughly<br>
> checked the structure, and it appears correct based on my verification.<br>
><br>
> Could you please explain why you believe the slab is wrong? Your insights<br>
> would be helpful as I continue to address the convergence issue.<br>
><br>
> On Wed, Sep 4, 2024 at 3:26 PM Vinay Maithani <<br>
> <a href="mailto:vinaymaithani21@kgpian.iitkgp.ac.in" target="_blank">vinaymaithani21@kgpian.iitkgp.ac.in</a>> wrote:<br>
><br>
>> Dear Quantum ESPRESSO community,<br>
>><br>
>> I hope this email finds you well.<br>
>><br>
>> I am currently working on large band gap semiconductors and have<br>
>> encountered persistent difficulties in achieving convergence while<br>
>> calculating the surface properties of these materials. Below are the<br>
>> strategies I have attempted so far:<br>
>><br>
>>    - *Mixing Beta:* Adjusted up to 0.1<br>
>>    - *Mixdim:* Increased to 24<br>
>>    - *Smearing:* Gaussian with degauss = 0.2<br>
>>    - *Nbnd:* Added a minimum of 20% more than the number of Kohn-Sham<br>
>>    states<br>
>>    - *Ecutwfc:* Increased from 50 Ry to 70 Ry<br>
>>    - *Ecutrho:* Set to 8 times Ecutwfc<br>
>>    - *Electron Max Step:* Extended to 2000 steps<br>
>><br>
>> Despite these efforts, my calculations still fail to converge, returning<br>
>> the error:<br>
>> *“Convergence NOT achieved after *** iterations: stopping.”*<br>
>><br>
>> I am using norm-conserving pseudopotentials and have applied dipole<br>
>> corrections, but the issue persists. I have attached the input and output<br>
>> files, along with the plot of potential data obtained via pp.x<br>
>> (plot_num=11) and averaged in the xy-plane where a discontinuity occurs in<br>
>> the vacuum region, available here: Google Drive link<br>
>> <<a href="https://drive.google.com/drive/folders/1HoSFk11TMN4qAo8jZoSvyEymK92xGu2r?usp=drive_link" rel="noreferrer" target="_blank">https://drive.google.com/drive/folders/1HoSFk11TMN4qAo8jZoSvyEymK92xGu2r?usp=drive_link</a>><br>
>> .<br>
>><br>
>> I would greatly appreciate any insights or suggestions you may have to<br>
>> help resolve this problem.<br>
>><br>
>> Thank you for your time and guidance.<br>
>><br>
>> Best regards,<br>
>> Vinay Maithani<br>
>> PhD Candidate, IIT Kharagpur<br>
>> <a href="mailto:vinaymaithani21@kgpian.iitkgp.ac.in" target="_blank">vinaymaithani21@kgpian.iitkgp.ac.in</a><br>
>><br>
<br>
<br>
<br>
GIUSEPPE MATTIOLI<br>
CNR - ISTITUTO DI STRUTTURA DELLA MATERIA<br>
Via Salaria Km 29,300 - C.P. 10<br>
I-00015 - Monterotondo Scalo (RM)<br>
Mob (*preferred*) +39 373 7305625<br>
Tel + 39 06 90672342 - Fax +39 06 90672316<br>
E-mail: <<a href="mailto:giuseppe.mattioli@ism.cnr.it" target="_blank">giuseppe.mattioli@ism.cnr.it</a>><br>
<br>
_______________________________________________<br>
The Quantum ESPRESSO community stands by the Ukrainian<br>
people and expresses its concerns about the devastating<br>
effects that the Russian military offensive has on their<br>
country and on the free and peaceful scientific, cultural,<br>
and economic cooperation amongst peoples<br>
_______________________________________________<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></div>