<div dir="ltr"><p class="MsoNormal"><span lang="EN-US">Thank Samuel and Matt for your reply.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">May be the system is difficult for
simulating, which is in high temperature and high density, and the system will
gradually become a high density and temperature plasma, so I need the nuclear
temperature= electronic temperature, we test, but some warnings about molecular
orbital.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> </span></p>
<p class="MsoNormal"><span lang="EN-US">*** 19:54:49 WARNING in
qs_mo_occupation:set_mo_occupation_1 :: ***<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> ***
Fermi-Dirac smearing includes the first MO ***<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> </span></p>
<p class="MsoNormal"><span lang="EN-US"> </span></p>
<p class="MsoNormal"><span lang="EN-US">By the way, when temperature is lower and about
300 K, the pressure is negative, why?<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">The input file <o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">##########################################<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">&GLOBAL<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
PROJECT 64O2-MD<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
RUN_TYPE MD<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
PRINT_LEVEL LOW<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">&END GLOBAL<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> </span></p>
<p class="MsoNormal"><span lang="EN-US">&MOTION<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&MD<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
ENSEMBLE NVT<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
STEPS 10000<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
TIMESTEP 0.5<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
TEMPERATURE 10000<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&THERMOSTAT<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
REGION GLOBAL<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
TYPE NOSE<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&NOSE<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
LENGTH 3<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
TIMECON 1.0<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&END<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&END<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&END MD<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&PRINT<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&STRESS SILENT<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&END STRESS<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&END PRINT<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">&END MOTION<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> </span></p>
<p class="MsoNormal"><span lang="EN-US">&FORCE_EVAL<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
METHOD Quickstep<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
STRESS_TENSOR ANALYTICAL <o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&DFT<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
BASIS_SET_FILE_NAME
/home/myu1410/work/cp2k/261/cp2k-2.6.1/data/BASIS_SET<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
POTENTIAL_FILE_NAME
/home/myu1410/work/cp2k/261/cp2k-2.6.1/data/POTENTIAL<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
MULTIPLICITY 1<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&MGRID<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
CUTOFF 400<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
NGRIDS 4<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> REL_CUTOFF 30<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&END MGRID<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&QS<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
EPS_DEFAULT 1.0E-12<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
EPS_GVG 1.0E-6<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
EPS_PGF_ORB 1.0E-6<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
MAP_CONSISTENT<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&END QS<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&SCF<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
MAX_SCF 200<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
SCF_GUESS atomic<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
EPS_SCF 1.00E-04<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
CHOLESKY INVERSE<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&DIAGONALIZATION<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&END DIAGONALIZATION<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">#
&OT<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">#
&END OT<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
ADDED_MOS 800 800<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&SMEAR<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
ELECTRONIC_TEMPERATURE 10000<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
METHOD FERMI_DIRAC<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&END SMEAR<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&MIXING<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
METHOD BROYDEN_MIXING<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
ALPHA 0.6<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
BETA 1.0<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
NBROYDEN 15<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&END<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&END SCF<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&XC<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&XC_GRID<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
XC_SMOOTH_RHO NN10<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
XC_DERIV SPLINE2<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&END XC_GRID<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&XC_FUNCTIONAL Pade<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&END XC_FUNCTIONAL<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&END XC<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&END DFT<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&SUBSYS<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&CELL<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
ABC 10.1007 10.1007 10.1007<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&END CELL<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&COORD<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
SCALED<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
@INCLUDE './64o2.xyz0'<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&END COORD <o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&KIND O <o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
BASIS_SET DZVP-GTH-PADE <o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
POTENTIAL GTH-PADE-q6 <o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&END KIND <o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
&END SUBSYS <o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">&END FORCE_EVAL <o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">##########################################<o:p></o:p></span></p><br><br></div>