<div dir="ltr"><p class="MsoNormal"><span lang="EN-US">Dear Marcella,</span></p>

<p class="MsoNormal"><span lang="EN-US">Thanks for your reply!</span></p>

<p class="MsoNormal"><span lang="EN-US">Yes, I need the thermal electron
temperature is equal to the atom temperature, because I simulate the argon
system in the high temperature and high density range. The argon atom will
gradually be ionized and the system become plasma. I hope the system reach
thermodynamic equilibrium. </span></p><p class="MsoNormal"> </p>

<p class="MsoNormal"><span lang="EN-US">I have done the same thing by VASP code,
but the calculation using VASP need more time. So I want to use CP2K do this
simulation.</span></p>

<p class="MsoNormal"><span lang="EN-US"> </span></p>

<p class="MsoNormal"><span lang="EN-US">Here is vasp INCAR</span></p>

<p class="MsoNormal"><span lang="EN-US">##############################<wbr>############</span></p>

<p class="MsoNormal"><span lang="EN-US"> 
SYSTEM = 64Ar-1.93g-17200k</span></p>

<p class="MsoNormal" style="text-indent:10.5pt"><span lang="EN-US">PREC=Accurate</span></p>

<p class="MsoNormal"><span lang="EN-US"> 
ENCUT = 600.0</span></p>

<p class="MsoNormal"><span lang="EN-US"> 
IBRION = 0</span></p>

<p class="MsoNormal"><span lang="EN-US"> 
ISTART = 0</span></p>

<p class="MsoNormal"><span lang="EN-US"> 
IALGO = 48</span></p>

<p class="MsoNormal"><span lang="EN-US">  NSW
= 10000</span></p>

<p class="MsoNormal"><span lang="EN-US"> 
NBLOCK = 1; KBLOCK = 50</span></p>

<p class="MsoNormal"><span lang="EN-US"> 
ISYM = 0</span></p>

<p class="MsoNormal"><span lang="EN-US"> 
SMASS = 1</span></p>

<p class="MsoNormal"><span lang="EN-US"> 
POTIM = 1.0</span></p>

<p class="MsoNormal"><span lang="EN-US"> 
TEBEG = 17200  #temperature</span></p>

<p class="MsoNormal"><span lang="EN-US"> 
TEEND = 17200</span></p>

<p class="MsoNormal"><span lang="EN-US"> 
ISMEAR = -1; SIGMA = 1.48218 #unit eV 1.48218 eV=17200 Kelvin</span></p>

<p class="MsoNormal"><span lang="EN-US"> 
LREAL = A</span></p>

<p class="MsoNormal"><span lang="EN-US"> 
ISIF = 2</span></p>

<p class="MsoNormal"><span lang="EN-US"> 
NELMIN = 6</span></p>

<p class="MsoNormal"><span lang="EN-US"> 
NBANDS = 500</span></p>

<p class="MsoNormal"><span lang="EN-US">##############################<wbr>############</span></p>

<p class="MsoNormal"><span lang="EN-US">ISMEAR determines how the partial
occupancies are set for each orbital. For the finite temperature LDA SIGMA determines
the width of the smearing in eV.</span></p>

<p class="MsoNormal"><span lang="EN-US">ISMEAR=-1 is Fermi smearing. </span></p>

<p class="MsoNormal"><span lang="EN-US"> </span></p>

<p class="MsoNormal"><span lang="EN-US">I think the keyword of CP2K (SMEAR and
ELECTRONIC_TEMPERTURE) have the same meanings.</span></p>

<p class="MsoNormal"><span lang="EN-US"> </span></p>

<p class="MsoNormal"><span lang="EN-US">Is that right about my understanding?</span></p>

<p class="MsoNormal"><span lang="EN-US"> </span></p>

<p class="MsoNormal"><span lang="EN-US">Best wishes,</span></p>

<p class="MsoNormal"><span lang="EN-US">jianzhifu</span></p></div>