<div>Dear all,</div><div><br></div><div><br></div><div>When I try to run MD simulation with ensemble_DFT by cp.x,the convergence procedure is different due to the difference between USPP file. I test </div><div>the metal system with 4 Fe atoms.</div><div><br></div><div>If the USPP file is Fe.pbe-spn-rrkjus_psl.0.2.1.UPF,then the code “cg_sub: missed minimum, case 1(or 3), iteration ***”occures in the output </div><div>file and the number of iteration is up to 250 ,which is the maximum number of conjugate gradient iterations for conjugate gradient minimizations of electronic states, I find that the value of codes'econs' & 'econt' in the *.evp file are “NaN”. </div><div><br></div><div>Then I test the system with the USPP file Fe.pbe-sp-van_mit.UPF,it's ok to convergence,while I'm not sure the USPP file is effective or not.</div><div><br></div><div>Can you tell me what's the difference between the two USPP file.</div><div>Belowing is my input file and USPP files,any suggestions would be appreciated.</div><div><br></div><div><br></div><div>INPUT:</div><div><br></div><div>&control</div><div><br></div><div>   calculation  = 'cp',</div><div>   title = ' MD Simulation' ,</div><div>   restart_mode = 'from_scratch',</div><div>   ndr =90,ndw =91,</div><div>   nstep  = 10,iprint = 1,isave  = 1,</div><div>   dt = 125.D0,</div><div>   tstress = .TRUE.,</div><div>   tprnfor = .TRUE.,</div><div>   prefix = 'Fe',</div><div>   pseudo_dir = '/public/home/---/pseudo/ ',</div><div>   outdir = '/public/home/---/tmp/'  </div><div><br></div><div>/</div><div><br></div><div>&system</div><div><br></div><div> ibrav = 0,</div><div> nat = 4, ntyp = 1,nspin=2,</div><div> nr1b = 24,nr2b = 24,nr3b = 24,</div><div> ecutwfc = 37.0, ecutrho = 235.0,</div><div> occupations= 'ensemble',</div><div> smearing='fd',</div><div> degauss=0.018,</div><div> vdw_corr='DFT-D',    </div><div> london_s6 = 0.75</div><div><br></div><div>/</div><div><br></div><div>&electrons</div><div><br></div><div> orthogonalization ='Gram-Schmidt',</div><div> startingwfc='random',</div><div> ampre = 0.02,</div><div> tcg = .true.,</div><div> passop= 0.3,</div><div> maxiter = 250,</div><div> emass = 300, emass_cutoff = 2.5,</div><div> conv_thr=1.d-6</div><div> n_inner = 2,</div><div> lambda_cold = 0.03,</div><div> niter_cold_restart = 2,</div><div><br></div><div>/</div><div><br></div><div>&ions</div><div><br></div><div>ion_dynamics = 'damp',</div><div>ion_damping = 0.05,</div><div><br></div><div>/</div><div><br></div><div>CELL_PARAMETERS  angstrom</div><div><br></div><div>       2.866400000000000       0.000000000000000       0.000000000000000</div><div>       0.000000000000000       4.053701755186239       0.000000000000000</div><div>       0.000000000000000       0.000000000000000      12.026850877593123</div><div><br></div><div>ATOMIC_SPECIES</div><div>  Fe  55.847  Fe.pbe-spn-rrkjus_psl.0.2.1.UPF/Fe.pbe-sp-van_mit.UPF</div><div>ATOMIC_POSITIONS  crystal</div><div><div>Fe   -0.0000000000000002   0.5000000000000000    0.1662945704038070 </div><div>Fe   0.4999999999999990    0.0000000000000002    0.1662945704038070 </div><div>Fe   0.0000000000000000    0.0000000000000000    0.3348217183847710 </div><div>Fe   0.5000000000000000    0.5000000000000000    0.3348217183847710 </div></div><div><br></div><div><b>Fe.pbe-spn-rrkjus_psl.0.2.1.UPF</b></div><div><div><UPF version="2.0.1"></div><div>  <PP_INFO></div><div>    Generated using "atomic" code by A. Dal Corso  v.5.0.99 svn rev. 10869</div><div>    Author: ADC</div><div>    Generation date: 10Apr2014</div><div>    Pseudopotential type: USPP</div><div>    Element: Fe</div><div>    Functional: PBE</div><div><br></div><div>    Suggested minimum cutoff for wavefunctions:  64. Ry</div><div>    Suggested minimum cutoff for charge density: 782. Ry</div><div>    The Pseudo was generated with a Scalar-Relativistic Calculation</div><div>    Local Potential according to Troullier-Martins recipe, cutoff radius:   1.8000</div><div><br></div><div>    Valence configuration: </div><div>    nl pn  l   occ       Rcut    Rcut US       E pseu</div><div>    3S  1  0  2.00      1.100      1.300    -6.910118</div><div>    4S  2  0  2.00      0.800      1.400    -0.388932</div><div>    3P  2  1  6.00      1.000      1.300    -4.413014</div><div>    4P  3  1  0.00      1.000      1.600    -0.097406</div><div>    3D  3  2  6.00      1.400      2.000    -0.551556</div><div>    Generation configuration:</div><div>    3S  1  0  2.00      1.100      1.300    -6.910160</div><div>    4S  2  0  2.00      0.800      1.400    -0.388938</div><div>    3P  2  1  6.00      1.000      1.300    -4.413054</div><div>    4P  3  1  0.00      1.000      1.600    -0.097411</div><div>    3D  3  2  6.00      1.400      2.000    -0.551588</div><div>    3D  3  2  0.00      1.400      2.000    -0.400000</div><div><br></div><div>    Pseudization used: troullier-martins</div><div><PP_INPUTFILE></div><div> &input</div><div>   title='Fe',</div><div>   zed=26.,</div><div>   rel=1,</div><div>   config='[Ar] 4s2.0 3d6.0 4p0',</div><div>   iswitch=3,</div><div>   dft='PBE'</div><div> /</div><div> &inputp</div><div>   pseudotype=3,</div><div>   file_pseudopw='Fe.pbe-spn-rrkjus_psl.0.2.1.UPF',</div><div>   author='ADC',</div><div>   lloc=-2,</div><div>   rcloc=1.8,</div><div>   which_augfun='PSQ',</div><div>   rmatch_augfun=0.68,</div><div>   nlcc=.true.,</div><div>   new_core_ps=.true.,</div><div>   rcore=0.6,</div><div>   tm=.true.</div><div> /</div><div>6</div><div>3S  1  0  2.00  0.00  1.10  1.30  0.0</div><div>4S  2  0  2.00  0.00  0.80  1.40  0.0</div><div>3P  2  1  6.00  0.00  1.00  1.30  0.0</div><div>4P  3  1  0.00  0.00  1.00  1.60  0.0</div><div>3D  3  2  6.00  0.00  1.40  2.00  0.0</div><div>3D  3  2  0.00 -0.40  1.40  2.00  0.0</div><div></PP_INPUTFILE></div><div>  </PP_INFO></div><div>  <!--                               --></div><div>  <!-- END OF HUMAN READABLE SECTION --></div><div>  <!--                               --></div><div>  <PP_HEADER generated='Generated using "atomic" code by A. Dal Corso  v.5.0.99 svn rev. 10869'</div><div>             author="ADC"</div><div>             date="10Apr2014"</div><div>             comment=""</div><div>             element="Fe"</div><div>             pseudo_type="USPP"</div><div>             relativistic="scalar"</div><div>             is_ultrasoft="T"</div><div>             is_paw="F"</div><div>             is_coulomb="F"</div><div>             has_so="F"</div><div>             has_wfc="F"</div><div>             has_gipaw="F"</div><div>             paw_as_gipaw="F"</div><div>             core_correction="T"</div><div>             functional="PBE"</div><div>             z_valence="1.600000000000000E+001"</div><div>             total_psenergy="-2.537079270394615E+002"</div><div>             wfc_cutoff="6.392614047675073E+001"</div><div>             rho_cutoff="7.817369406582422E+002"</div><div>             l_max="2"</div><div>             l_max_rho="4"</div><div>             l_local="-2"</div><div>             mesh_size="1191"</div><div>             number_of_wfc="5"</div><div>             number_of_proj="6"/></div><div>  <PP_MESH dx="1.250000000000000E-002" mesh="1191" xmin="-7.000000000000000E+000" rmax="1.000000000000000E+002"</div><div>zmesh="2.600000000000000E+001"></div><div>    <PP_R type="real" size="1191" columns="4"></div><div> 3.507238329055831E-005  3.551353956419316E-005  3.596024490063778E-005  3.641256909850977E-005 </div><div> 3.687058283438538E-005  3.733435767384265E-005  3.780396608264379E-005  3.827948143805812E-005 </div></div><div>......</div><div><br></div><div><b>Fe.pbe-sp-van_mit.UPF</b></div><div><div><UPF version="2.0.1"></div><div><br></div><div><PP_INFO></div><div>Generated by new atomic code, or converted to UPF format</div><div>Author:</div><div>Generation date:</div><div>Pseudopotential type: US</div><div>Element: Fe</div><div>Functional:  SLA  PW   PBE  PBE</div><div>Suggested minimum cutoff for wavefunctions:   0. Ry</div><div>Suggested minimum cutoff for charge density:   0. Ry</div><div>The Pseudo was generated with a Scalar-Relativistic Calculation</div><div>L component and cutoff radius for Local Potential:  0   0.0000</div><div>Valence configuration:</div><div>nl pn  l   occ       Rcut    Rcut US       E pseu</div><div>3S  0  0  2.00     10.000      2.000     0.000000</div><div>3P  0  1  6.00     10.000      2.000     0.000000</div><div>3D  0  2  5.50     10.000      2.000     0.000000</div><div>4S  0  0  0.40     10.000      2.000     0.000000</div><div>4P  0  1  0.10     10.000      2.000     0.000000</div><div>Generation configuration: not available.</div><div></PP_INFO></div><div><!--</div><div>END OF HUMAN READABLE SECTION</div><div>--></div><div><PP_HEADER generated="Generated by new atomic code, or converted to UPF format" author="" date="" comment="" element="Fe" pseudo_type="US" relativistic="scalar" is_ultrasoft="true" is_paw="false" is_coulomb="false" has_so="false" has_wfc="false" has_gipaw="false" paw_as_gipaw="false" core_correction="false" functional=" SLA  PW   PBE  PBE" z_valence="1.600000000000e1" total_psenergy="-2.505735262170e2" wfc_cutoff="0.000000000000e0" rho_cutoff="0.000000000000e0" l_max="2" l_max_rho="0" l_local="0" mesh_size="861" number_of_wfc="5" number_of_proj="6"/></div><div><PP_MESH mesh="861" rmax="2.040110545010e2" zmesh="2.600000000000e1"></div><div><PP_R></div><div>0.000000000000000e0 1.629646503080000e-6 3.287149538759999e-6 4.972985275609999e-6</div><div>6.687638021610001e-6 8.431600363309999e-6 1.020537330730000e-5 1.200946642430000e-5</div></div><div><br></div><div><br></div><div><div>Thank you!</div><div><br></div><div>Best regards,</div><div><br></div><div><br></div><div>------</div><div> </div><div><br></div><div>LU Hailin, Graduate Student</div><div><br></div><div>State Key Laboratory of Mechanical Transmission, Chongqing University</div><div> </div><div>Chongqing 400044, China</div><div><br></div><div>Tel: (+86) 18362360286</div></div>