第一原理分子動力学プログラム STATE Senri Wiki
開始行:
* Crystalline Si [#me90bdef]
In this example, silicon in the diamond structure (space ...
** SCF calculation [#a76d305d]
- Input file (nfinp_1)
0 0 0 0 0 0 : I_CTRL(1:6) (DUMMY)
4.00 8.00 1 2 2 : GMAX, GMAXP, NTYP...
227 2 : NUM_SPACE_GROUP T...
10.30 10.30 10.30 90.00 90.00 90.00 : A B C ALPHA BETA ...
08 08 08 1 1 1 : N1 N2 N3 M1 M2 M3
0 0 : NCORD NINV : IWEI...
0.00d0 0.00d0 0.00d0 1 1 1
0.25d0 0.25d0 0.25d0 1 1 1
14 0.50 28.09 6 1 0.2 : TYPE 1IATOMN,ALFA...
0 0 0 0 0 : ICOND INIPOS INIV...
0 1 : IPRE IPRI
20 20 0 84200.00 0 : NMD1 NMD2 LAST_IT...
6 1 : WAY_MIX MIX_WHAT
0 20 0.60 : ITER_START KBXMIX...
0.20 0.30 0.20 0.20 0.20 : DTIM1 DTIM2 DTIM3...
300.00 4 1 0.50D-09 : DTIO IMDALG IEXPL...
0.0002 0.50D+03 0 : WIDTH FORCCR ISTR...
ggapbe 1 : XCTYPE KSPIN
2.00 : DESTM
102 : NBZTYP
4 4 4 : DUMMY
4 4 4 : DUMMY
8 : KEG
1 : NEXTST
0 : DUMMY
2 : IMSD
0 : EVALUATE_EKO_DIFF
2 : NPDOSAO
1
2
-15.00 5.00 0.20 501
2.400 0.20
0.20 14
0 0.00 : SM_N DOPPING
- Execution (interactive mode)
$ mpirun -np 6 ./STATE < nfinp_1 > nfout_1
- Execution (batch job)
$ qsub qsub.sh
-- Job script (qsub.sh)
#$ -S /bin/sh
#$ -cwd
#$ -pe fillup 6
#$ -N Si
#disable OPENMP parallelism
OMP_NUM_THREADS=1
# execuable of the STATE code
ln -fs ../../src/STATE
# pseudopotential data
ln -fs ${HOME}/STATE/gncpp/pot.Si_pbe1 fort.37
# launch STATE
mpirun -np $NSLOTS ./STATE < nfinp_1 > nfout_1
- Output file (nfout_1)~
Convergence of the total energy can be monitored by using...
$ grep ETOT: nfout_1
The result is:
ETOT: 1 -6.05513096 0.6055E+01 0.3203E-02
ETOT: 2 -7.84016187 0.1785E+01 0.5187E-02
ETOT: 3 -7.87270490 0.3254E-01 0.2825E-02
ETOT: 4 -7.87351715 0.8123E-03 0.6089E-03
ETOT: 5 -7.87355245 0.3530E-04 0.1887E-03
ETOT: 6 -7.87355822 0.5769E-05 0.2560E-04
ETOT: 7 -7.87355833 0.1069E-06 0.1066E-04
ETOT: 8 -7.87355833 0.4548E-08 0.1730E-05
ETOT: 9 -7.87355833 0.1916E-09 0.4194E-06
ETOT: 10 -7.87355833 0.2984E-10 0.1669E-07
ETOT: 11 -7.87355833 0.1951E-11 0.9524E-08
Converged total energy and its components
TOTAL ENERGY AND ITS COMPONENTS
TOTAL ENERGY = -7.8735583...
FREE ENERGY = -7.8735583...
KINETIC ENERGY = 3.0192245...
HARTREE ENERGY = 0.5501422...
XC ENERGY = -2.4009866...
LOCAL ENERGY = -0.8429497...
NONLOCAL ENERGY = 0.1688528...
EWALD ENERGY = -8.3678416...
PC ENERGY = 0.0000000...
ENTROPIC ENERGY = 0.0000000...
** Band structure calculation [#q35a4f7f]
To calculate the band structure, an SCF calculation is pe...
- Input file (nfinp_scf)
0 0 0 0 0 0 : INPUT_CTRL(1:6)...
6.00 12.00 1 2 2 : GMAX GMAXP NTYP...
227 2 : num_space_group...
10.2347 10.2347 10.2347 90.0 90.0 90.0 : A B C ALPHA BET...
4 4 4 2 2 2 : N1 N2 N3 M1 M2 M3
0 0 : NCORD, NINV
0.00 0.00 0.00 1 1 1 : CPS(1,1:3) IWEI...
0.25 0.25 0.25 1 1 1 : CPS(2,1:3) IWEI...
14 0.5000 28.09 6 1 0.0 : ATOMN ALFA AMIO...
0 0 0 0 0 : ICOND INIPOS IN...
0 1 : IPRE IPRI
200 200 0 84200.00 0 : NMD1 NMD2 LAST_...
6 1 : WAY_MIX MIX_WHAT
0 20 0.60 : ITER_START KBXM...
0.20 0.30 0.20 0.20 0.20 : DTIM1 DTIM2 DTI...
300.00 4 1 0.50D-09 : DTIO IMDALG IEX...
-0.0020 0.50D+03 0 : WIDTH FORCCR IS...
ggapbe 1 : xctype,kspin
2.00 : DETSTM
101 : NBZTYP
0 0 0 : NKX NKY NKZ (DU...
0 0 0 : NKX2 NKY2 NKZ2 ...
8 : KEG
1 : NEXTST
0 : (DUMMY)
2 : IMSD
0 : EVALUATE_EKO_DIFF
0 : NPDOSAO
0 0.00 : SM_N DOPING
- Input file (nfinp_band)
0 0 0 0 0 0 : INPUT_CTRL(1:6) (...
6.00 12.00 1 2 2 : GMAX GMAXP NTYP...
227 2 : num_space_group...
10.2347 10.2347 10.2347 90.0 90.0 90.0 : A B C ALPHA BET...
4 4 4 2 2 2 : N1 N2 N3 M1 M2 M3
0 0 : NCORD, NINV
0.00 0.00 0.00 1 1 1 : CPS(1,1:3) IWEI...
0.25 0.25 0.25 1 1 1 : CPS(2,1:3) IWEI...
14 0.5000 28.09 6 1 0.0 : ATOMN ALFA AMIO...
22 0 0 0 0 : ICOND INIPOS IN...
0 1 : IPRE IPRI
200 200 0 84200.00 0 : NMD1 NMD2 LAST_...
6 1 : WAY_MIX MIX_WHAT
0 20 0.60 : ITER_START KBXM...
0.20 0.30 0.20 0.20 0.20 : DTIM1 DTIM2 DTI...
300.00 4 1 0.50D-09 : DTIO IMDALG IEX...
-0.0100 0.50D+03 0 : WIDTH FORCCR IS...
ggapbe 1 : xctype,kspin
2.00 : DETSTM
101 : NBZTYP
0 0 0 : NKX NKY NKZ (DU...
0 0 0 : NKX2 NKY2 NKZ2 ...
16 : KEG
1 : NEXTST
0 : (DUMMY)
2 : IMSD
0 : EVALUATE_EKO_DIFF
0 : NPDOSAO
0 0.00 : SM_N DOPING
&KPOINTS_BAND
NKSEG 6
KMESH 40 20 20 20 40 20
KPOINTS
0.000 0.000 0.000
0.000 0.500 0.500
0.250 0.500 0.750
0.500 0.500 0.500
0.000 0.000 0.000
0.375 0.375 0.750
0.500 0.500 1.000
&END
When the non-SCF calculation is performed, k-points and K...
Here a sample bash script, which generate the band struct...
#!/bin/sh
#
# Fermi level obtained in the SCF calculation
ef=0.0790726449
# Number of k-points
nk=161
# Number of bands
nb=16
# Divisions for each segment (in this case we have 6 seg...
ik1=40
ik2=60
ik3=80
ik4=100
ik5=140
ik6=160
# Maximum energy
emin=-12.5
# Minimum energy
emax=7.5
#
nk1=`expr $ik1 + 2`
nk2=`expr $ik2 + 2`
nk3=`expr $ik3 + 2`
nk4=`expr $ik4 + 2`
nk5=`expr $ik5 + 2`
nk6=`expr $ik6 + 2`
#
e2b=~/STATE/util/bandutil/src/energy2band
#
${e2b} << EOF
${nb}
${nb}
${nk}
${ef}
EOF
#
k1=`head -${nk1} band.data | tail -1 | awk '{print $1}'`
k2=`head -${nk2} band.data | tail -1 | awk '{print $1}'`
k3=`head -${nk3} band.data | tail -1 | awk '{print $1}'`
k4=`head -${nk4} band.data | tail -1 | awk '{print $1}'`
k5=`head -${nk5} band.data | tail -1 | awk '{print $1}'`
k6=`head -${nk6} band.data | tail -1 | awk '{print $1}'`
#
xmin=0.0
xmax=${k6}
offset=-0.5
#
pos=`echo "scale=2; $emin + $offset" | bc`
#
gnuplot<<EOF
set xrange [$xmin:$xmax]
set yrange [$emin:$emax]
set arrow from ${k1},${emin} to ${k1},${emax} nohead lt 0
set arrow from ${k2},${emin} to ${k2},${emax} nohead lt 0
set arrow from ${k3},${emin} to ${k3},${emax} nohead lt 0
set arrow from ${k4},${emin} to ${k4},${emax} nohead lt 0
set arrow from ${k5},${emin} to ${k5},${emax} nohead lt 0
unset xtics
set mytics
set ylabel 'Energy (eV)'
set label '{/Symbol G}' at 0.0, ${pos}
set label 'X' at ${k1}, ${pos}
set label 'W' at ${k2}, ${pos}
set label 'L' at ${k3}, ${pos}
set label '{/Symbol G}' at ${k4}, ${pos}
set label 'K' at ${k5}, ${pos}
set label 'X' at ${k6}, ${pos}
plot 'band.data' with lines
EOF
Note that the energy origin is set to the valence band ma...
終了行:
* Crystalline Si [#me90bdef]
In this example, silicon in the diamond structure (space ...
** SCF calculation [#a76d305d]
- Input file (nfinp_1)
0 0 0 0 0 0 : I_CTRL(1:6) (DUMMY)
4.00 8.00 1 2 2 : GMAX, GMAXP, NTYP...
227 2 : NUM_SPACE_GROUP T...
10.30 10.30 10.30 90.00 90.00 90.00 : A B C ALPHA BETA ...
08 08 08 1 1 1 : N1 N2 N3 M1 M2 M3
0 0 : NCORD NINV : IWEI...
0.00d0 0.00d0 0.00d0 1 1 1
0.25d0 0.25d0 0.25d0 1 1 1
14 0.50 28.09 6 1 0.2 : TYPE 1IATOMN,ALFA...
0 0 0 0 0 : ICOND INIPOS INIV...
0 1 : IPRE IPRI
20 20 0 84200.00 0 : NMD1 NMD2 LAST_IT...
6 1 : WAY_MIX MIX_WHAT
0 20 0.60 : ITER_START KBXMIX...
0.20 0.30 0.20 0.20 0.20 : DTIM1 DTIM2 DTIM3...
300.00 4 1 0.50D-09 : DTIO IMDALG IEXPL...
0.0002 0.50D+03 0 : WIDTH FORCCR ISTR...
ggapbe 1 : XCTYPE KSPIN
2.00 : DESTM
102 : NBZTYP
4 4 4 : DUMMY
4 4 4 : DUMMY
8 : KEG
1 : NEXTST
0 : DUMMY
2 : IMSD
0 : EVALUATE_EKO_DIFF
2 : NPDOSAO
1
2
-15.00 5.00 0.20 501
2.400 0.20
0.20 14
0 0.00 : SM_N DOPPING
- Execution (interactive mode)
$ mpirun -np 6 ./STATE < nfinp_1 > nfout_1
- Execution (batch job)
$ qsub qsub.sh
-- Job script (qsub.sh)
#$ -S /bin/sh
#$ -cwd
#$ -pe fillup 6
#$ -N Si
#disable OPENMP parallelism
OMP_NUM_THREADS=1
# execuable of the STATE code
ln -fs ../../src/STATE
# pseudopotential data
ln -fs ${HOME}/STATE/gncpp/pot.Si_pbe1 fort.37
# launch STATE
mpirun -np $NSLOTS ./STATE < nfinp_1 > nfout_1
- Output file (nfout_1)~
Convergence of the total energy can be monitored by using...
$ grep ETOT: nfout_1
The result is:
ETOT: 1 -6.05513096 0.6055E+01 0.3203E-02
ETOT: 2 -7.84016187 0.1785E+01 0.5187E-02
ETOT: 3 -7.87270490 0.3254E-01 0.2825E-02
ETOT: 4 -7.87351715 0.8123E-03 0.6089E-03
ETOT: 5 -7.87355245 0.3530E-04 0.1887E-03
ETOT: 6 -7.87355822 0.5769E-05 0.2560E-04
ETOT: 7 -7.87355833 0.1069E-06 0.1066E-04
ETOT: 8 -7.87355833 0.4548E-08 0.1730E-05
ETOT: 9 -7.87355833 0.1916E-09 0.4194E-06
ETOT: 10 -7.87355833 0.2984E-10 0.1669E-07
ETOT: 11 -7.87355833 0.1951E-11 0.9524E-08
Converged total energy and its components
TOTAL ENERGY AND ITS COMPONENTS
TOTAL ENERGY = -7.8735583...
FREE ENERGY = -7.8735583...
KINETIC ENERGY = 3.0192245...
HARTREE ENERGY = 0.5501422...
XC ENERGY = -2.4009866...
LOCAL ENERGY = -0.8429497...
NONLOCAL ENERGY = 0.1688528...
EWALD ENERGY = -8.3678416...
PC ENERGY = 0.0000000...
ENTROPIC ENERGY = 0.0000000...
** Band structure calculation [#q35a4f7f]
To calculate the band structure, an SCF calculation is pe...
- Input file (nfinp_scf)
0 0 0 0 0 0 : INPUT_CTRL(1:6)...
6.00 12.00 1 2 2 : GMAX GMAXP NTYP...
227 2 : num_space_group...
10.2347 10.2347 10.2347 90.0 90.0 90.0 : A B C ALPHA BET...
4 4 4 2 2 2 : N1 N2 N3 M1 M2 M3
0 0 : NCORD, NINV
0.00 0.00 0.00 1 1 1 : CPS(1,1:3) IWEI...
0.25 0.25 0.25 1 1 1 : CPS(2,1:3) IWEI...
14 0.5000 28.09 6 1 0.0 : ATOMN ALFA AMIO...
0 0 0 0 0 : ICOND INIPOS IN...
0 1 : IPRE IPRI
200 200 0 84200.00 0 : NMD1 NMD2 LAST_...
6 1 : WAY_MIX MIX_WHAT
0 20 0.60 : ITER_START KBXM...
0.20 0.30 0.20 0.20 0.20 : DTIM1 DTIM2 DTI...
300.00 4 1 0.50D-09 : DTIO IMDALG IEX...
-0.0020 0.50D+03 0 : WIDTH FORCCR IS...
ggapbe 1 : xctype,kspin
2.00 : DETSTM
101 : NBZTYP
0 0 0 : NKX NKY NKZ (DU...
0 0 0 : NKX2 NKY2 NKZ2 ...
8 : KEG
1 : NEXTST
0 : (DUMMY)
2 : IMSD
0 : EVALUATE_EKO_DIFF
0 : NPDOSAO
0 0.00 : SM_N DOPING
- Input file (nfinp_band)
0 0 0 0 0 0 : INPUT_CTRL(1:6) (...
6.00 12.00 1 2 2 : GMAX GMAXP NTYP...
227 2 : num_space_group...
10.2347 10.2347 10.2347 90.0 90.0 90.0 : A B C ALPHA BET...
4 4 4 2 2 2 : N1 N2 N3 M1 M2 M3
0 0 : NCORD, NINV
0.00 0.00 0.00 1 1 1 : CPS(1,1:3) IWEI...
0.25 0.25 0.25 1 1 1 : CPS(2,1:3) IWEI...
14 0.5000 28.09 6 1 0.0 : ATOMN ALFA AMIO...
22 0 0 0 0 : ICOND INIPOS IN...
0 1 : IPRE IPRI
200 200 0 84200.00 0 : NMD1 NMD2 LAST_...
6 1 : WAY_MIX MIX_WHAT
0 20 0.60 : ITER_START KBXM...
0.20 0.30 0.20 0.20 0.20 : DTIM1 DTIM2 DTI...
300.00 4 1 0.50D-09 : DTIO IMDALG IEX...
-0.0100 0.50D+03 0 : WIDTH FORCCR IS...
ggapbe 1 : xctype,kspin
2.00 : DETSTM
101 : NBZTYP
0 0 0 : NKX NKY NKZ (DU...
0 0 0 : NKX2 NKY2 NKZ2 ...
16 : KEG
1 : NEXTST
0 : (DUMMY)
2 : IMSD
0 : EVALUATE_EKO_DIFF
0 : NPDOSAO
0 0.00 : SM_N DOPING
&KPOINTS_BAND
NKSEG 6
KMESH 40 20 20 20 40 20
KPOINTS
0.000 0.000 0.000
0.000 0.500 0.500
0.250 0.500 0.750
0.500 0.500 0.500
0.000 0.000 0.000
0.375 0.375 0.750
0.500 0.500 1.000
&END
When the non-SCF calculation is performed, k-points and K...
Here a sample bash script, which generate the band struct...
#!/bin/sh
#
# Fermi level obtained in the SCF calculation
ef=0.0790726449
# Number of k-points
nk=161
# Number of bands
nb=16
# Divisions for each segment (in this case we have 6 seg...
ik1=40
ik2=60
ik3=80
ik4=100
ik5=140
ik6=160
# Maximum energy
emin=-12.5
# Minimum energy
emax=7.5
#
nk1=`expr $ik1 + 2`
nk2=`expr $ik2 + 2`
nk3=`expr $ik3 + 2`
nk4=`expr $ik4 + 2`
nk5=`expr $ik5 + 2`
nk6=`expr $ik6 + 2`
#
e2b=~/STATE/util/bandutil/src/energy2band
#
${e2b} << EOF
${nb}
${nb}
${nk}
${ef}
EOF
#
k1=`head -${nk1} band.data | tail -1 | awk '{print $1}'`
k2=`head -${nk2} band.data | tail -1 | awk '{print $1}'`
k3=`head -${nk3} band.data | tail -1 | awk '{print $1}'`
k4=`head -${nk4} band.data | tail -1 | awk '{print $1}'`
k5=`head -${nk5} band.data | tail -1 | awk '{print $1}'`
k6=`head -${nk6} band.data | tail -1 | awk '{print $1}'`
#
xmin=0.0
xmax=${k6}
offset=-0.5
#
pos=`echo "scale=2; $emin + $offset" | bc`
#
gnuplot<<EOF
set xrange [$xmin:$xmax]
set yrange [$emin:$emax]
set arrow from ${k1},${emin} to ${k1},${emax} nohead lt 0
set arrow from ${k2},${emin} to ${k2},${emax} nohead lt 0
set arrow from ${k3},${emin} to ${k3},${emax} nohead lt 0
set arrow from ${k4},${emin} to ${k4},${emax} nohead lt 0
set arrow from ${k5},${emin} to ${k5},${emax} nohead lt 0
unset xtics
set mytics
set ylabel 'Energy (eV)'
set label '{/Symbol G}' at 0.0, ${pos}
set label 'X' at ${k1}, ${pos}
set label 'W' at ${k2}, ${pos}
set label 'L' at ${k3}, ${pos}
set label '{/Symbol G}' at ${k4}, ${pos}
set label 'K' at ${k5}, ${pos}
set label 'X' at ${k6}, ${pos}
plot 'band.data' with lines
EOF
Note that the energy origin is set to the valence band ma...
ページ名: