第一原理分子動力学プログラム STATE Senri Wiki
開始行:
*Crystalline Ni [#wa9506a3]
In this example, SCF calculation of Ni in the fcc structu...
-Input file (nfinp_1)
0 0 0 0 0 0
5.00 15.00 1 1 1 : GMAX GMAXP NTYP N...
221 2 : num_space_group, ...
6.70 6.70 6.70 90.00 90.00 90.00 : a,b,c,alpha,beta,...
12 12 12 1 1 1 : K_mesh
0 0 : NCORD NINV
0.00 0.00 0.00 1 0 1 : CPS(1:3) IWEI IMD...
28 0.50 58.69 6 1 0.1 : IATOMN,ALFA,AMION...
0 0 0 0 0 : ICOND INIPOS INIV...
0 1 : ipre, ipri
30 30 0 84200.00 0 : nmd1,nmd2,last_it...
6 1 : WAY_MIX MIX_WHAT
0 20 0.30 : iter_start, KBXMI...
0.20 0.30 0.20 0.20 0.20 : dtim1,dtim2,dtim3...
300.00 4 1 0.50D-09 : dtio ,imdalg, iex...
-0.002 0.50D+03 0 : width,forccr,istr...
ggapbe 2 : xctype,kspin
2.00 : destm
101 : nbztyp
4 4 4 : dummy
4 4 4 : dummy
10 : keg
1 : nextst
0 : dummy
2 : imsd
0 : evaluate_eko_diff
1 : npdosao
1
-15.00 5.00 0.20 501
2.40000 0.200000
0.20000 14
0 0.00 : sm_n, dopping
- Job script (qsub.sh)
#$ -S /bin/sh
#$ -cwd
#$ -pe fillup 6
#$ -N Ni
#disable OPENMP parallelism
OMP_NUM_THREADS=1
# execuable of the STATE code
ln -fs ../../src/STATE
# pseudopotential data
ln -fs ${HOME}/STATE/gncpp/pot.Ni_pbe4 fort.37
# launch STATE
mpirun -np $NSLOTS ./STATE < nfinp_1 > nfout_1
- Output file (nfout_1)~
Monitor the convergence of the SCF calculation: Note that...
$ grep FTOT: nfout_1
The result is:
FTOT: 1 -48.03422182 0.4803E+02 0.1740E-01
FTOT: 2 -48.09019134 0.5597E-01 0.2913E-01
FTOT: 3 -48.35446197 0.2643E+00 0.4213E-02
FTOT: 4 -48.35583377 0.1372E-02 0.3191E-02
FTOT: 5 -48.35596669 0.1329E-03 0.2878E-02
FTOT: 6 -48.35602991 0.6322E-04 0.3500E-02
FTOT: 7 -48.35607610 0.4619E-04 0.2013E-02
FTOT: 8 -48.35607759 0.1490E-05 0.1434E-02
FTOT: 9 -48.35608861 0.1102E-04 0.3847E-03
FTOT: 10 -48.35609075 0.2137E-05 0.1860E-03
FTOT: 11 -48.35609115 0.4076E-06 0.6846E-04
FTOT: 12 -48.35609122 0.6508E-07 0.4412E-05
FTOT: 13 -48.35609122 0.4191E-08 0.4316E-05
FTOT: 14 -48.35609122 0.4136E-09 0.7494E-06
FTOT: 15 -48.35609122 0.2350E-10 0.2593E-06
FTOT: 16 -48.35609122 0.3081E-10 0.5970E-07
FTOT: 17 -48.35609122 0.1107E-10 0.3456E-07
FTOT: 18 -48.35609122 0.7724E-11 0.2828E-07
FTOT: 19 -48.35609122 0.8050E-11 0.9248E-08
Converged total energy and its component
TOTAL ENERGY AND ITS COMPONENTS
TOTAL ENERGY = -48.3560927...
FREE ENERGY = -48.3560912...
KINETIC ENERGY = 6.2464447...
HARTREE ENERGY = 10.9194438...
XC ENERGY = -18.0779981...
LOCAL ENERGY = -8.9288357...
NONLOCAL ENERGY = -4.2997589...
EWALD ENERGY = -34.2153886...
PC ENERGY = 0.0000000...
ENTROPIC ENERGY = 0.0000015...
We can see that (electronic) entropic contribution ("ENTR...
終了行:
*Crystalline Ni [#wa9506a3]
In this example, SCF calculation of Ni in the fcc structu...
-Input file (nfinp_1)
0 0 0 0 0 0
5.00 15.00 1 1 1 : GMAX GMAXP NTYP N...
221 2 : num_space_group, ...
6.70 6.70 6.70 90.00 90.00 90.00 : a,b,c,alpha,beta,...
12 12 12 1 1 1 : K_mesh
0 0 : NCORD NINV
0.00 0.00 0.00 1 0 1 : CPS(1:3) IWEI IMD...
28 0.50 58.69 6 1 0.1 : IATOMN,ALFA,AMION...
0 0 0 0 0 : ICOND INIPOS INIV...
0 1 : ipre, ipri
30 30 0 84200.00 0 : nmd1,nmd2,last_it...
6 1 : WAY_MIX MIX_WHAT
0 20 0.30 : iter_start, KBXMI...
0.20 0.30 0.20 0.20 0.20 : dtim1,dtim2,dtim3...
300.00 4 1 0.50D-09 : dtio ,imdalg, iex...
-0.002 0.50D+03 0 : width,forccr,istr...
ggapbe 2 : xctype,kspin
2.00 : destm
101 : nbztyp
4 4 4 : dummy
4 4 4 : dummy
10 : keg
1 : nextst
0 : dummy
2 : imsd
0 : evaluate_eko_diff
1 : npdosao
1
-15.00 5.00 0.20 501
2.40000 0.200000
0.20000 14
0 0.00 : sm_n, dopping
- Job script (qsub.sh)
#$ -S /bin/sh
#$ -cwd
#$ -pe fillup 6
#$ -N Ni
#disable OPENMP parallelism
OMP_NUM_THREADS=1
# execuable of the STATE code
ln -fs ../../src/STATE
# pseudopotential data
ln -fs ${HOME}/STATE/gncpp/pot.Ni_pbe4 fort.37
# launch STATE
mpirun -np $NSLOTS ./STATE < nfinp_1 > nfout_1
- Output file (nfout_1)~
Monitor the convergence of the SCF calculation: Note that...
$ grep FTOT: nfout_1
The result is:
FTOT: 1 -48.03422182 0.4803E+02 0.1740E-01
FTOT: 2 -48.09019134 0.5597E-01 0.2913E-01
FTOT: 3 -48.35446197 0.2643E+00 0.4213E-02
FTOT: 4 -48.35583377 0.1372E-02 0.3191E-02
FTOT: 5 -48.35596669 0.1329E-03 0.2878E-02
FTOT: 6 -48.35602991 0.6322E-04 0.3500E-02
FTOT: 7 -48.35607610 0.4619E-04 0.2013E-02
FTOT: 8 -48.35607759 0.1490E-05 0.1434E-02
FTOT: 9 -48.35608861 0.1102E-04 0.3847E-03
FTOT: 10 -48.35609075 0.2137E-05 0.1860E-03
FTOT: 11 -48.35609115 0.4076E-06 0.6846E-04
FTOT: 12 -48.35609122 0.6508E-07 0.4412E-05
FTOT: 13 -48.35609122 0.4191E-08 0.4316E-05
FTOT: 14 -48.35609122 0.4136E-09 0.7494E-06
FTOT: 15 -48.35609122 0.2350E-10 0.2593E-06
FTOT: 16 -48.35609122 0.3081E-10 0.5970E-07
FTOT: 17 -48.35609122 0.1107E-10 0.3456E-07
FTOT: 18 -48.35609122 0.7724E-11 0.2828E-07
FTOT: 19 -48.35609122 0.8050E-11 0.9248E-08
Converged total energy and its component
TOTAL ENERGY AND ITS COMPONENTS
TOTAL ENERGY = -48.3560927...
FREE ENERGY = -48.3560912...
KINETIC ENERGY = 6.2464447...
HARTREE ENERGY = 10.9194438...
XC ENERGY = -18.0779981...
LOCAL ENERGY = -8.9288357...
NONLOCAL ENERGY = -4.2997589...
EWALD ENERGY = -34.2153886...
PC ENERGY = 0.0000000...
ENTROPIC ENERGY = 0.0000015...
We can see that (electronic) entropic contribution ("ENTR...
ページ名: