Dear Forum,<br><br>I am using PWscf for interface structural optimization.<br>The geometry relaxation on interface (LiCoO2/sulfides) seems very difficult to converge so far.<br>Can anyone give me help?<br><br>Here is my input file:
<br><br>&amp;CONTROL<br>&nbsp; calculation = &quot;relax&quot;,<br>&nbsp; prefix&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = &quot;intf66a_rlx&quot;<br>&nbsp; pseudo_dir&nbsp; = &quot;./&quot;,<br>&nbsp; outdir&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = &quot;./&quot;,<br>/<br>&amp;SYSTEM<br>&nbsp; ibrav&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 14,<br>&nbsp; celldm(1)&nbsp;&nbsp; = 
61.18796,<br>&nbsp; celldm(2)&nbsp;&nbsp; = 0.18633,<br>&nbsp; celldm(3)&nbsp;&nbsp; = 0.38711,<br>&nbsp; celldm(4)&nbsp;&nbsp; = -0.22231431,<br>&nbsp; celldm(5)&nbsp;&nbsp; = -0.39398,<br>&nbsp; celldm(6)&nbsp;&nbsp; = 0.267793,<br>&nbsp; nat&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 66,<br>&nbsp; ntyp&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 5,<br>&nbsp; ecutwfc&nbsp;&nbsp;&nbsp;&nbsp; = 42.D0
,<br>&nbsp; ecutrho&nbsp;&nbsp;&nbsp;&nbsp; = 320.D0,<br>&nbsp; occupations = &quot;smearing&quot;,<br>&nbsp; smearing&nbsp;&nbsp;&nbsp; = &quot;methfessel-paxton&quot;,<br>&nbsp; degauss&nbsp;&nbsp;&nbsp;&nbsp; = 0.05D0<br>/<br>&amp;ELECTRONS<br>&nbsp; conv_thr&nbsp;&nbsp;&nbsp; = 1.D-6,<br>&nbsp; mixing_beta = 0.1D0,<br>
&nbsp; mixing_mode = &quot;local-TF&quot;,<br>/<br>&amp;IONS<br>&nbsp; ion_dynamics&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = &quot;bfgs&quot;,<br>&nbsp; opt_scheme&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = &quot;broyden&quot;,<br>&nbsp; bfgs_ndim&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 3,<br>&nbsp; trust_radius_max&nbsp; = 0.8D0,<br>&nbsp; pot_extrapolation = &quot;second_order&quot;,
<br>&nbsp; wfc_extrapolation = &quot;second_order&quot;,<br>/<br>ATOMIC_SPECIES<br>P&nbsp;&nbsp; 30.9737&nbsp; P.pbe-n-van.UPF<br>Li&nbsp; 6.941&nbsp;&nbsp;&nbsp; Li.pbe-n-van.UPF<br>S&nbsp;&nbsp; 32.066&nbsp;&nbsp; S.pbe-van_bm.UPF<br>Co&nbsp; 58.9332&nbsp; Co.pbe-nd-rrkjus.UPF<br>O&nbsp;&nbsp; 15.9994
&nbsp; O.pbe-rrkjus.UPF<br>ATOMIC_POSITIONS {angstrom}<br>........<br>K_POINTS crystal <br>&nbsp;1&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.2914500000&nbsp;&nbsp;&nbsp; 0.5294200000&nbsp;&nbsp; -0.2781600000&nbsp;&nbsp;&nbsp; 1.0<br><br>===================================================<br><br>
The output is as following:<br>&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 1.677255&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.001300<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.607898&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.002820<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.448659&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 
0.002098<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.352568&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.001579<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.291404&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.001432<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.254533&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 
0.001601<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.233841&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.001331<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.223115&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.001884<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.222679&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 
0.001979<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.230224&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.001892<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.241234&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.001275<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.252683&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 
0.001494<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.266314&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.001429<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.284588&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.001093<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.309129&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 
0.001386<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.312575&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.001656<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.306147&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.001511<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.280490&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 
0.001399<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.261216&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000735<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.244781&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.001154<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.230320&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 
0.001218<br><br>Thank you in advance<br><br>Sei Ou<br>Simulatio Co.<br>Yokohama-shi, Japan<br><br><br>