<div>+ ferrariruffino.fz@</div><div> </div><div>Fabrizio,</div><div> </div><div>It seems I made a mistake but cannot figure out where. Because I cannot reproduce it.</div><div> </div><div>Anyway, you are correct - both 'input_dft = "scan"' and 'input_dft = "MGGA_X_SCAN + MGGA_C_SCAN" result in the same indexes:</div><div><div>from PW output:   Exchange-correlation= MGGA_X_SCAN + MGGA_C_SCAN</div><div>                           (   0   0   0   0   0 263 267)</div><div> </div><div>However, when one takes SCAN pseudopotentials from the library on QE website, and don't use 'input_dft = 'scan' ' option, then PW should take XC from pseudopotentials. But index is different:</div><div> </div><div><div>Exchange-correlation=  SLA  PW   SCAN SCAN</div><div>                           (   1   4   0   0   0   5   0)</div><div> </div><div>Header of pseudopotential file is:</div><div>==============</div><div><div><div><PP_INFO></div><div>Generated using Vanderbilt code, version   7  3  2                              </div><div>Author: Yi Yao     Generation date:    0    0    0                              </div><div>Automatically converted from original format                                    </div><div>    0        The Pseudo was generated with a Non-Relativistic Calculation</div><div>nl pn  l   occ               Rcut            Rcut US             E pseu</div><div>0S  0  0  2.00      0.00000000000      0.00000000000      0.00000000000</div><div>0P  0  1  2.00      0.00000000000      0.00000000000      0.00000000000</div><div></PP_INFO></div><div> </div><div><PP_HEADER></div><div>   0                   Version Number</div><div>   C                   Element</div><div>   NC                  Norm - Conserving pseudopotential</div><div>    F                  Nonlinear Core Correction</div><div>SLA  PW   SCAN SCAN    SCAN Exchange-Correlation functional</div><div>    4.00000000000      Z valence</div><div>    0.00000000000      Total energy</div><div>  0.0000000  0.0000000 Suggested cutoff for wfc and rho</div><div>    0                  Max angular momentum component</div><div> 1389                  Number of points in mesh</div><div>    2    1             Number of Wavefunctions, Number of Projectors</div><div> Wavefunctions         nl  l   occ</div><div>                       0S  0  2.00</div><div>                       0P  1  2.00</div></div></div></div><div>=========</div></div><div>and this run fails with error:<div> Error in routine tau_xc (5):</div><div>     This case is not implemented</div><div> </div><div>I guess this is 'different' SCAN.</div><div> </div><div>thanks for your help,</div><div> Sergey </div><div> </div></div>