<div>Dear all,</div><div><br></div><div> </div><div><br></div><div>When I try to run MD simulation with ensemble_DFT by cp.x,it’s hard to convergence. The ensemble-DFT (eDFT) is a method to simulate the metals in the framework of "ab-initio" molecular dynamics. As far as I know,it’s a method of BOMD. Due to my system contains Fe atom, I test the metal system with 4 Fe atoms,then the code “cg_sub: missed minimum, case 1(or 3), iteration ***”occures in the output file and the number of iteration is up to 250 ,which is the electron maxstep during each step.</div><div><br></div><div>Second, I find that the value of codes”econs”&”econt”in the *.evp file are “NaN”,what make them so large?</div><div><br></div><div>Can you tell me if some parameters are not set reasonably or some important codes was ignored.</div><div><br></div><div>Belowing is my input and output file,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-n-rrkjus_psl.0.2.4.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> </div><div><br></div><div>OUTPUT:</div><div><br></div><div>PERFORMING CONJUGATE GRADIENT MINIMIZATION OF EL. STATES</div><div><br></div><div> cg_sub: missed minimum, case 1, iteration           1</div><div> cg_sub: missed minimum, case 1, iteration          66</div><div> cg_sub: missed minimum, case 1, iteration          98</div><div> cg_sub: missed minimum, case 3, iteration         135</div><div> cg_sub: missed minimum, case 1, iteration         146</div><div><br></div><div>NOTE: eigenvalues are not computed without ortho</div><div><br></div><div>Occupations  :</div><div> 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000</div><div> 1.000000 0.999948 0.999076 0.969357 0.948499 0.153049 1.000000 1.000000 1.000000 1.000000</div><div> 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 0.999973</div><div> 0.974829 0.955270</div><div><br></div><div>  nfi     ekinc              temph  tempp    etot                 enthal               econs                econt              vnhh    xnhh0  vnhp    xnhp0</div><div><br></div><div><br></div><div>       nfi  tempp            E -T.S-mu.nbsp         +K_p #Iter</div><div><br></div><div>Step     1     11   -122.244412           NaN           NaN  250</div><div>......</div><div> cg_sub: missed minimum, case 1, iteration           6</div><div> cg_sub: missed minimum, case 1, iteration          24</div><div> cg_sub: missed minimum, case 1, iteration          46</div><div> cg_sub: missed minimum, case 3, iteration          63</div><div> cg_sub: missed minimum, case 1, iteration         102</div><div> cg_sub: missed minimum, case 1, iteration         105</div><div> cg_sub: missed minimum, case 3, iteration         112</div><div> cg_sub: missed minimum, case 1, iteration         116</div><div> cg_sub: missed minimum, case 3, iteration         120</div><div> cg_sub: missed minimum, case 1, iteration         126</div><div> cg_sub: missed minimum, case 1, iteration         140</div><div> cg_sub: missed minimum, case 1, iteration         145</div><div> cg_sub: missed minimum, case 1, iteration         147</div><div> cg_sub: missed minimum, case 3, iteration         149</div><div> cg_sub: missed minimum, case 3, iteration         154</div><div> cg_sub: missed minimum, case 3, iteration         159</div><div> cg_sub: missed minimum, case 1, iteration         163</div><div> cg_sub: missed minimum, case 3, iteration         167</div><div> cg_sub: missed minimum, case 1, iteration         179</div><div> cg_sub: missed minimum, case 3, iteration         184</div><div> cg_sub: missed minimum, case 1, iteration         188</div><div> cg_sub: missed minimum, case 3, iteration         192</div><div> cg_sub: missed minimum, case 1, iteration         198</div><div> cg_sub: missed minimum, case 3, iteration         202</div><div> cg_sub: missed minimum, case 3, iteration         205</div><div> cg_sub: missed minimum, case 1, iteration         209</div><div> cg_sub: missed minimum, case 1, iteration         212</div><div> cg_sub: missed minimum, case 3, iteration         219</div><div> cg_sub: missed minimum, case 1, iteration         222</div><div> cg_sub: missed minimum, case 1, iteration         226</div><div> cg_sub: missed minimum, case 3, iteration         243</div><div> cg_sub: missed minimum, case 3, iteration         245</div><div> cg_sub: missed minimum, case 3, iteration         248</div><div><br></div><div>NOTE: eigenvalues are not computed without ortho</div><div><br></div><div>Occupations  :</div><div><br></div><div> 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000</div><div> 1.000000 0.999924 0.993571 0.973833 0.966758 0.098386 1.000000 1.000000 1.000000 1.000000</div><div> 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 0.999703</div><div> 0.997142 0.970684 </div><div><br></div><div>       nfi  tempp            E -T.S-mu.nbsp         +K_p #Iter</div><div><br></div><div>Step     2     98   -122.228699           NaN           NaN  250</div><div><br></div><div>……</div><div><br></div><div> cg_sub: missed minimum, case 1, iteration          28</div><div> cg_sub: missed minimum, case 1, iteration          46</div><div> cg_sub: missed minimum, case 1, iteration          66</div><div> cg_sub: missed minimum, case 3, iteration         108</div><div> cg_sub: missed minimum, case 1, iteration         133</div><div> cg_sub: missed minimum, case 1, iteration         186</div><div> cg_sub: missed minimum, case 1, iteration         228</div><div> cg_sub: missed minimum, case 1, iteration         243</div><div><br></div><div>NOTE: eigenvalues are not computed without ortho</div><div><br></div><div>Occupations  :</div><div><br></div><div> 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000</div><div> 1.000000 0.999639 0.988952 0.982049 0.922040 0.208842 1.000000 1.000000 1.000000 1.000000</div><div> 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000</div><div> 0.999838 0.898639</div><div><br></div><div> </div><div><br></div><div>       nfi  tempp            E -T.S-mu.nbsp         +K_p #Iter</div><div><br></div><div>Step     3    268   -122.273098           NaN           NaN  250</div><div><br></div><div>……</div><div><br></div><div> cg_sub: missed minimum, case 1, iteration          22</div><div> cg_sub: missed minimum, case 3, iteration          57</div><div> cg_sub: missed minimum, case 1, iteration          88</div><div> cg_sub: missed minimum, case 2, iteration         131</div><div> cg_sub: missed minimum, case 3, iteration         144</div><div> cg_sub: missed minimum, case 1, iteration         148</div><div> cg_sub: missed minimum, case 3, iteration         157</div><div> cg_sub: missed minimum, case 1, iteration         160</div><div> cg_sub: missed minimum, case 1, iteration         187</div><div> cg_sub: missed minimum, case 1, iteration         198</div><div> cg_sub: missed minimum, case 1, iteration         202</div><div> cg_sub: missed minimum, case 3, iteration         220</div><div> cg_sub: missed minimum, case 1, iteration         228</div><div> cg_sub: missed minimum, case 1, iteration         242</div><div><br></div><div>NOTE: eigenvalues are not computed without ortho</div><div><br></div><div>Occupations  :</div><div><br></div><div> 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000</div><div> 1.000000 0.999856 0.991684 0.987230 0.857249 0.226278 1.000000 1.000000 1.000000 1.000000</div><div> 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000</div><div> 0.999981 0.937722 </div><div><br></div><div>       nfi  tempp            E -T.S-mu.nbsp         +K_p #Iter</div><div><br></div><div>Step     4    377   -122.252885           NaN           NaN  250</div><div><br></div><div>……</div><div><br></div><div> cg_sub: missed minimum, case 3, iteration         135</div><div> cg_sub: missed minimum, case 3, iteration         171</div><div> cg_sub: missed minimum, case 1, iteration         189</div><div> cg_sub: missed minimum, case 3, iteration         206</div><div> cg_sub: missed minimum, case 1, iteration         235</div><div> cg_sub: missed minimum, case 1, iteration         246</div><div><br></div><div>NOTE: eigenvalues are not computed without ortho</div><div><br></div><div>Occupations  :</div><div><br></div><div> 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000</div><div> 1.000000 0.999150 0.975162 0.974626 0.880859 0.240104 1.000000 1.000000 1.000000 1.000000</div><div> 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000</div><div> 0.999968 0.930132</div><div> </div><div><br></div><div>       nfi  tempp            E -T.S-mu.nbsp         +K_p #Iter</div><div><br></div><div>Step     5    458   -122.284141           NaN           NaN   250</div><div><br></div><div>……</div><div><br></div><div> cg_sub: missed minimum, case 1, iteration           4</div><div> cg_sub: missed minimum, case 1, iteration          10</div><div> cg_sub: missed minimum, case 1, iteration          18</div><div> cg_sub: missed minimum, case 1, iteration          46</div><div> cg_sub: missed minimum, case 3, iteration          72</div><div> cg_sub: missed minimum, case 1, iteration         164</div><div> cg_sub: missed minimum, case 1, iteration         166</div><div> cg_sub: missed minimum, case 1, iteration         173</div><div> cg_sub: missed minimum, case 3, iteration         179</div><div> cg_sub: missed minimum, case 1, iteration         189</div><div> cg_sub: missed minimum, case 3, iteration         211</div><div> cg_sub: missed minimum, case 3, iteration         228</div><div> cg_sub: missed minimum, case 3, iteration         238</div><div> cg_sub: missed minimum, case 3, iteration         243</div><div><br></div><div>NOTE: eigenvalues are not computed without ortho</div><div><br></div><div>Occupations  :</div><div><br></div><div> 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000</div><div> 1.000000 0.999354 0.977205 0.976795 0.849851 0.245372 1.000000 1.000000 1.000000 1.000000</div><div> 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000</div><div> 0.999993 0.951430 </div><div><br></div><div>       nfi  tempp            E -T.S-mu.nbsp         +K_p #Iter</div><div><br></div><div>Step     6    670   -122.280486           NaN           NaN  250</div><div><br></div><div>……</div><div><br></div><div> cg_sub: missed minimum, case 1, iteration           5</div><div> cg_sub: missed minimum, case 1, iteration          42</div><div> cg_sub: missed minimum, case 1, iteration          74</div><div> cg_sub: missed minimum, case 1, iteration         107</div><div> cg_sub: missed minimum, case 3, iteration         139</div><div> cg_sub: missed minimum, case 3, iteration         159</div><div> cg_sub: missed minimum, case 3, iteration         165</div><div> cg_sub: missed minimum, case 3, iteration         177</div><div> cg_sub: missed minimum, case 3, iteration         196</div><div> cg_sub: missed minimum, case 1, iteration         209</div><div><br></div><div>NOTE: eigenvalues are not computed without ortho</div><div><br></div><div>Occupations  :</div><div><br></div><div> 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000</div><div> 1.000000 0.999258 0.980387 0.979132 0.821162 0.276359 1.000000 1.000000 1.000000 1.000000</div><div> 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000</div><div> 0.999994 0.943709</div><div> </div><div><br></div><div>       nfi  tempp            E -T.S-mu.nbsp         +K_p #Iter</div><div><br></div><div>Step     7    930   -122.290608           NaN           NaN  250</div><div><br></div><div>……</div><div><br></div><div> cg_sub: missed minimum, case 1, iteration           9</div><div> cg_sub: missed minimum, case 3, iteration          30</div><div> cg_sub: missed minimum, case 1, iteration          33</div><div> cg_sub: missed minimum, case 3, iteration          52</div><div> cg_sub: missed minimum, case 1, iteration          55</div><div> cg_sub: missed minimum, case 3, iteration          92</div><div> cg_sub: missed minimum, case 1, iteration         131</div><div> cg_sub: missed minimum, case 1, iteration         145</div><div> cg_sub: missed minimum, case 1, iteration         152</div><div> cg_sub: missed minimum, case 3, iteration         169</div><div> cg_sub: missed minimum, case 3, iteration         189</div><div><br></div><div>NOTE: eigenvalues are not computed without ortho</div><div><br></div><div>Occupations  :</div><div><br></div><div> 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000</div><div> 1.000000 0.999276 0.983691 0.977505 0.811211 0.268521 1.000000 1.000000 1.000000 1.000000</div><div> 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000</div><div> 0.999999 0.959797 </div><div><br></div><div>       nfi  tempp            E -T.S-mu.nbsp         +K_p #Iter</div><div><br></div><div>Step     8   1334   -122.273479           NaN           NaN  250</div><div><br></div><div>……</div><div><br></div><div> cg_sub: missed minimum, case 3, iteration          27</div><div> cg_sub: missed minimum, case 1, iteration          73</div><div> cg_sub: missed minimum, case 1, iteration         110</div><div> cg_sub: missed minimum, case 1, iteration         195</div><div> cg_sub: missed minimum, case 1, iteration         207</div><div> cg_sub: missed minimum, case 1, iteration         209</div><div> cg_sub: missed minimum, case 3, iteration         225</div><div> cg_sub: missed minimum, case 1, iteration         227</div><div><br></div><div>NOTE: eigenvalues are not computed without ortho</div><div><br></div><div>Occupations  :</div><div><br></div><div> 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000</div><div> 1.000000 0.997041 0.978604 0.970747 0.775047 0.544445 1.000000 1.000000 1.000000 1.000000</div><div> 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000</div><div> 0.999917 0.734200</div><div> </div><div>       nfi  tempp            E -T.S-mu.nbsp         +K_p #Iter</div><div><br></div><div>Step     9   1308   -122.361095           NaN           NaN  250</div><div><br></div><div>…… </div><div><br></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>