第一原理分子動力学プログラム STATE Senri Wiki
開始行:
* Hydrogen terminated Au(111) surface [#a734259b]
This tutorial explains how to analyze the electronic stru...
#ref(http://www-cp.prec.eng.osaka-u.ac.jp/puki_state/grap...
** SCF calculation [#tb8cd26e]
To obtain the self-consistent charge density, electronic ...
- Input file (nfinp_scf)
0 0 0 0 0 0 : I_CTROL (D...
6.00 20.00 2 6 6 : GMAX GMAXP...
1 0 : NUM_SPAC_G...
Cartesian
5.498778336527 0.000000000000 0.00000000...
-2.749389168264 4.762081729212 0.00000000...
0.000000000000 0.000000000000 85.30494051...
12 12 1 1 1 1 : NK1 NK2 NK...
1 0 : NCORD NINV...
0.000000000000 0.000000000000 8.97946742...
0.000000000000 3.174721152808 4.48973371...
0.000000000000 -3.174721152808 0.00000000...
0.000000000000 0.000000000000 -4.48973371...
0.000000000000 3.174721152808 -8.97946742...
0.000000026956 -3.174721183174 -10.61015812...
79.00 0.50 12.0107 1 1 0.0 : IATOMN ALF...
1.00 0.50 12.0107 1 1 0.0 : IATOMN ALF...
0 0 0 0 0 : ICOND INIP...
0 1 : IPRE IPRI
200 200 0 80000.00 0 : NMD1 NMD2 ...
6 1 : WAY_MIX MI...
0 20 0.60 : ITER_START...
0.20 0.30 0.20 0.20 0.20 : DTIM1 DTIM...
200.00 2 1 1.0D-12 : DTIO IMDAL...
-0.0010 1.00D-03 0 : WIDTH FORC...
rev-vdW-DF2 1 : XCTYPE KSPIN
2.00 : DESTM
101 : NBZTYP
4 4 4 : DUMMY
4 4 4 : DUMMY
48 : KEG
1 : NEXTST
0 : DUMMY
2 : IMSD
0 : EVALUATE_E...
0 : NPDOSAO
0 0.000 : SM_N DOPING
- Execution~
In the command line or in the job script, execute
> mpirun -np 24 < nfinp_scf > nfout_scf
** PDOS (AOLDOS) calculation [#jb76cbfe]
Then densities of state projected onto the atomic orbital...
- Input file (nfinp_aoldos)
0 0 0 0 0 0 : I_CTROL (D...
6.00 20.00 2 6 6 : GMAX GMAXP...
1 0 : NUM_SPAC_G...
Cartesian
5.498778336527 0.000000000000 0.00000000...
-2.749389168264 4.762081729212 0.00000000...
0.000000000000 0.000000000000 85.30494051...
12 12 1 1 1 1 : NK1 NK2 NK...
1 0 : NCORD NINV...
0.000000000000 0.000000000000 8.97946742...
0.000000000000 3.174721152808 4.48973371...
0.000000000000 -3.174721152808 0.00000000...
0.000000000000 0.000000000000 -4.48973371...
0.000000000000 3.174721152808 -8.97946742...
0.000000026956 -3.174721183174 -10.61015812...
79.00 0.50 12.0107 1 1 0.0 : IATOMN ALF...
1.00 0.50 12.0107 1 1 0.0 : IATOMN ALF...
14 0 0 0 0 : ICOND INIP...
0 1 : IPRE IPRI
200 200 0 80000.00 0 : NMD1 NMD2 ...
6 1 : WAY_MIX MI...
0 20 0.60 : ITER_START...
0.20 0.30 0.20 0.20 0.20 : DTIM1 DTIM...
200.00 2 1 1.0D-12 : DTIO IMDAL...
-0.0010 1.00D-03 0 : WIDTH FORC...
rev-vdW-DF2 1 : XCTYPE KSPIN
2.00 : DESTM
101 : NBZTYP
4 4 4 : DUMMY
4 4 4 : DUMMY
48 : KEG
1 : NEXTST
0 : DUMMY
2 : IMSD
0 : EVALUATE_E...
6 : NPDOSAO
1
2
3
4
5
6
-15.00 5.00 0.20 501 : EPDOS(1) E...
2.2 0.3 : RAD WIDTH ...
0.7 0.3 : RAD WIDTH ...
0.2 12 : DR NR
0 0.000 : SM_N DOPING
Here to calculate PDOS/AOLDOS by a postprocessing, ICOND ...
14 0 0 0 0 : ICOND INIPO...
NPDOSAO (Number of Atomic orbitals for PDOS) > 0
6 : NPDOSAO
followed by the atomic indices for which PDOS are calcula...
1
2
3
4
5
6
Then minimum energy (EPDOS(1)), maximum energy (EPDOS(2))...
-15.00 5.00 0.20 501 : EPDOS(1) E...
Cutoff radii (RAD) for the atomic orbitals and smearing w...
2.2 0.3 : RAD WIDTH ...
0.7 0.3 : RAD WIDTH ...
Width (DR) and number of grids (NR) for the radial integr...
0.2 12 : DR NR
Note that EPSDOS(3) is used only when
&OTHERS
GAUSSDOS
&END
is written in the input file for the calculation of densi...
- Execution
> mpirun -np 24 < nfinp_aoldos > nfout_aoldos
- Analysis
PDOS/AOLDOS is printed by the keyword "AO_LDOS" in the ou...
AO_LDOS: atm spn energy s px py pz dzz...
AO_LDOS: 1 1 -15.00 0.000 0.000 0.000 0.000 0.000...
AO_LDOS: 1 1 -14.96 0.000 0.000 0.000 0.000 0.000...
AO_LDOS: 1 1 -14.92 0.000 0.000 0.000 0.000 0.000...
...
PDOS/AOLDOS can be extracted by using the state2pdos.pl s...
state2pdos.pl nfout_aoldos
then the following messages are printed
make dos data file :: 0001
make dos data file :: 0002
make dos data file :: 0003
make dos data file :: 0004
make dos data file :: 0005
make dos data file :: 0006
meaning that the dos data are printed to files 0001, 0002...
NOTE: with the latest version of the code, they will be p...
The data are given in the above files in the following or...
Energy s px py pz dzz dxx-yy dxy ...
Finally plot the data using for e.g., gnuplot or xmgrace.
Alternatively, you can perform SCF and AOLDOS calculation...
0 0 0 0 0 0 : I_CTROL (D...
6.00 20.00 2 6 6 : GMAX GMAXP...
1 0 : NUM_SPAC_G...
Cartesian
5.498778336527 0.000000000000 0.00000000...
-2.749389168264 4.762081729212 0.00000000...
0.000000000000 0.000000000000 85.30494051...
12 12 1 1 1 1 : NK1 NK2 NK...
1 0 : NCORD NINV...
0.000000000000 0.000000000000 8.97946742...
0.000000000000 3.174721152808 4.48973371...
0.000000000000 -3.174721152808 0.00000000...
0.000000000000 0.000000000000 -4.48973371...
0.000000000000 3.174721152808 -8.97946742...
0.000000026956 -3.174721183174 -10.61015812...
79.00 0.50 12.0107 1 1 0.0 : IATOMN ALF...
1.00 0.50 12.0107 1 1 0.0 : IATOMN ALF...
0 0 0 0 0 : ICOND INIPO...
0 1 : IPRE IPRI
200 200 0 80000.00 0 : NMD1 NMD2 ...
6 1 : WAY_MIX MI...
0 20 0.60 : ITER_START...
0.20 0.30 0.20 0.20 0.20 : DTIM1 DTIM...
200.00 2 1 1.0D-12 : DTIO IMDAL...
-0.0010 1.00D-03 0 : WIDTH FORC...
rev-vdW-DF2 1 : XCTYPE KSPIN
2.00 : DESTM
101 : NBZTYP
4 4 4 : DUMMY
4 4 4 : DUMMY
48 : KEG
1 : NEXTST
0 : DUMMY
2 : IMSD
0 : EVALUATE_E...
6 : NPDOSAO
1
2
3
4
5
6
-15.00 5.00 0.20 501 : EPDOS(1) E...
2.2 0.3 : RAD WIDTH ...
0.7 0.3 : RAD WIDTH ...
0.2 12 : DR NR
0 0.000 : SM_N DOPING
-Visualization
Let us visualize PDOSs for the bottom Au and H atoms usin...
Here is an example gnuplot script to generate an eps file...
set terminal postscript eps color 'Helvetica' 20
set output 'pdos.eps'
xmin=-11.0
xmax=4.0
ymin=0.0
ymax=2.0
set xrange [xmin:xmax]
set yrange [ymin:ymax]
set yzeroaxis
set xlabel '{/Helvetica-Oblique E}-{/Helvetica-Oblique E...
set ylabel 'PDOS (Arb. unit)'
plot 'pdos_0005.dat' using ($1):($11) title 'Au#5 s' wit...
'pdos_0005.dat' using ($1):($12) title 'Au#5 p' wit...
'pdos_0005.dat' using ($1):($13) title 'Au#5 d' wit...
'pdos_0006.dat' using ($1):($2) title 'H s' with li...
Type
> gnuplot pdos.gp
to generate the eps file called "pdos.eps" as shown below
#ref(http://www-cp.prec.eng.osaka-u.ac.jp/puki_state/grap...
** Band structure calculation [#l6de7d56]
The band structure calculation can be performed after the...
- Input file (nfinp_band)
0 0 0 0 0 0 : I_CTROL (D...
6.00 20.00 2 6 6 : GMAX GMAXP...
1 0 : NUM_SPAC_G...
Cartesian
5.498778336527 0.000000000000 0.00000000...
-2.749389168264 4.762081729212 0.00000000...
0.000000000000 0.000000000000 85.30494051...
12 12 1 1 1 1 : NK1 NK2 NK...
1 0 : NCORD NINV...
0.000000000000 0.000000000000 8.97946742...
0.000000000000 3.174721152808 4.48973371...
0.000000000000 -3.174721152808 0.00000000...
0.000000000000 0.000000000000 -4.48973371...
0.000000000000 3.174721152808 -8.97946742...
0.000000026956 -3.174721183174 -10.61015812...
79.00 0.50 12.0107 1 1 0.0 : IATOMN ALF...
1.00 0.50 12.0107 1 1 0.0 : IATOMN ALF...
22 0 0 0 0 : ICOND INIP...
0 1 : IPRE IPRI
200 200 0 80000.00 0 : NMD1 NMD2 ...
6 1 : WAY_MIX MI...
0 20 0.60 : ITER_START...
0.20 0.30 0.20 0.20 0.20 : DTIM1 DTIM...
200.00 2 1 1.0D-12 : DTIO IMDAL...
-0.0010 1.00D-03 0 : WIDTH FORC...
rev-vdW-DF2 1 : XCTYPE KSPIN
2.00 : DESTM
101 : NBZTYP
4 4 4 : DUMMY
4 4 4 : DUMMY
48 : KEG
1 : NEXTST
0 : DUMMY
2 : IMSD
0 : EVALUATE_E...
0 : NPDOSAO
0 0.000 : SM_N DOPIN...
!
&KPOINTS_BAND
NKSEG 3
KMESH 40 20 40
KPOINTS
0.0 0.0 0.0
0.5 0.0 0.0
0.33333333 0.33333333 0.0
0.0 0.0 0.0
&END
For the band structure calculation, ICOND is set to 22, a...
&KPOINTS_BAND
NKSEG 3
KMESH 40 20 40
KPOINTS
0.0 0.0 0.0
0.5 0.0 0.0
0.33333333 0.33333333 0.0
0.0 0.0 0.0
&END
where the variables in the &KPOINTS_BAND ... &END block are
- the number of k-point segment (NKSEG)
- the k-point mesh for each segment (KMESH)
- the k-points which define the symmetry k-points (# of k...
K-points are given in the unit of the reciprocal lattice ...
After the band structure run, the file "energy.data" is g...
To convert the "energy.data" into plottable one, "energy2...
Type
> energy2band
and input number of bands, number of bands to be plotted,...
Here is an example script to run energy2band:
#!/bin/sh
E2B=energy2band
scfout_file='nfout_scf'
bandout_file='nfout_band'
ef=`grep 'FERMI' ${scfout_file} | awk '{print $4}'`
nb=`grep 'NUMBER OF BANDS CONSIDERED' ${bandout_file} | ...
nk=`grep 'NUMBER OF K-POINTS' ${bandout_file} | grep -v...
${E2B} << EOF
${nb}
${nb}
${nk}
${ef}
EOF
By using the following gnuplot script, the band structure...
#ref(http://www-cp.prec.eng.osaka-u.ac.jp/puki_state/grap...
The following gnuplot script can be used to generate the ...
emin=-12
emax=4
xmin=0.0000
xmax=1.80236106
G=0.00000000
M=0.65970994
K=1.04059365
G2=1.80236106
offset=-0.5
set xrange [xmin:xmax]
set yrange [emin:emax]
set arrow from M,emin to M,emax nohead lt 1
set arrow from K,emin to K,emax nohead lt 1
set xzeroaxis
unset xtics
unset key
set ylabel '{/Helvetica-Oblique E-E}_F (eV)'
set label '{/Symbol G}' at G,emin+offset center
set label 'M' at M,emin+offset center
set label 'K' at K,emin+offset center
set label '{/Symbol G}' at G2,emin+offset center
#set terminal postscript eps 'Helvetica' 20
#set output 'band.eps'
plot 'band.data' with lines lw 3
pause -1
終了行:
* Hydrogen terminated Au(111) surface [#a734259b]
This tutorial explains how to analyze the electronic stru...
#ref(http://www-cp.prec.eng.osaka-u.ac.jp/puki_state/grap...
** SCF calculation [#tb8cd26e]
To obtain the self-consistent charge density, electronic ...
- Input file (nfinp_scf)
0 0 0 0 0 0 : I_CTROL (D...
6.00 20.00 2 6 6 : GMAX GMAXP...
1 0 : NUM_SPAC_G...
Cartesian
5.498778336527 0.000000000000 0.00000000...
-2.749389168264 4.762081729212 0.00000000...
0.000000000000 0.000000000000 85.30494051...
12 12 1 1 1 1 : NK1 NK2 NK...
1 0 : NCORD NINV...
0.000000000000 0.000000000000 8.97946742...
0.000000000000 3.174721152808 4.48973371...
0.000000000000 -3.174721152808 0.00000000...
0.000000000000 0.000000000000 -4.48973371...
0.000000000000 3.174721152808 -8.97946742...
0.000000026956 -3.174721183174 -10.61015812...
79.00 0.50 12.0107 1 1 0.0 : IATOMN ALF...
1.00 0.50 12.0107 1 1 0.0 : IATOMN ALF...
0 0 0 0 0 : ICOND INIP...
0 1 : IPRE IPRI
200 200 0 80000.00 0 : NMD1 NMD2 ...
6 1 : WAY_MIX MI...
0 20 0.60 : ITER_START...
0.20 0.30 0.20 0.20 0.20 : DTIM1 DTIM...
200.00 2 1 1.0D-12 : DTIO IMDAL...
-0.0010 1.00D-03 0 : WIDTH FORC...
rev-vdW-DF2 1 : XCTYPE KSPIN
2.00 : DESTM
101 : NBZTYP
4 4 4 : DUMMY
4 4 4 : DUMMY
48 : KEG
1 : NEXTST
0 : DUMMY
2 : IMSD
0 : EVALUATE_E...
0 : NPDOSAO
0 0.000 : SM_N DOPING
- Execution~
In the command line or in the job script, execute
> mpirun -np 24 < nfinp_scf > nfout_scf
** PDOS (AOLDOS) calculation [#jb76cbfe]
Then densities of state projected onto the atomic orbital...
- Input file (nfinp_aoldos)
0 0 0 0 0 0 : I_CTROL (D...
6.00 20.00 2 6 6 : GMAX GMAXP...
1 0 : NUM_SPAC_G...
Cartesian
5.498778336527 0.000000000000 0.00000000...
-2.749389168264 4.762081729212 0.00000000...
0.000000000000 0.000000000000 85.30494051...
12 12 1 1 1 1 : NK1 NK2 NK...
1 0 : NCORD NINV...
0.000000000000 0.000000000000 8.97946742...
0.000000000000 3.174721152808 4.48973371...
0.000000000000 -3.174721152808 0.00000000...
0.000000000000 0.000000000000 -4.48973371...
0.000000000000 3.174721152808 -8.97946742...
0.000000026956 -3.174721183174 -10.61015812...
79.00 0.50 12.0107 1 1 0.0 : IATOMN ALF...
1.00 0.50 12.0107 1 1 0.0 : IATOMN ALF...
14 0 0 0 0 : ICOND INIP...
0 1 : IPRE IPRI
200 200 0 80000.00 0 : NMD1 NMD2 ...
6 1 : WAY_MIX MI...
0 20 0.60 : ITER_START...
0.20 0.30 0.20 0.20 0.20 : DTIM1 DTIM...
200.00 2 1 1.0D-12 : DTIO IMDAL...
-0.0010 1.00D-03 0 : WIDTH FORC...
rev-vdW-DF2 1 : XCTYPE KSPIN
2.00 : DESTM
101 : NBZTYP
4 4 4 : DUMMY
4 4 4 : DUMMY
48 : KEG
1 : NEXTST
0 : DUMMY
2 : IMSD
0 : EVALUATE_E...
6 : NPDOSAO
1
2
3
4
5
6
-15.00 5.00 0.20 501 : EPDOS(1) E...
2.2 0.3 : RAD WIDTH ...
0.7 0.3 : RAD WIDTH ...
0.2 12 : DR NR
0 0.000 : SM_N DOPING
Here to calculate PDOS/AOLDOS by a postprocessing, ICOND ...
14 0 0 0 0 : ICOND INIPO...
NPDOSAO (Number of Atomic orbitals for PDOS) > 0
6 : NPDOSAO
followed by the atomic indices for which PDOS are calcula...
1
2
3
4
5
6
Then minimum energy (EPDOS(1)), maximum energy (EPDOS(2))...
-15.00 5.00 0.20 501 : EPDOS(1) E...
Cutoff radii (RAD) for the atomic orbitals and smearing w...
2.2 0.3 : RAD WIDTH ...
0.7 0.3 : RAD WIDTH ...
Width (DR) and number of grids (NR) for the radial integr...
0.2 12 : DR NR
Note that EPSDOS(3) is used only when
&OTHERS
GAUSSDOS
&END
is written in the input file for the calculation of densi...
- Execution
> mpirun -np 24 < nfinp_aoldos > nfout_aoldos
- Analysis
PDOS/AOLDOS is printed by the keyword "AO_LDOS" in the ou...
AO_LDOS: atm spn energy s px py pz dzz...
AO_LDOS: 1 1 -15.00 0.000 0.000 0.000 0.000 0.000...
AO_LDOS: 1 1 -14.96 0.000 0.000 0.000 0.000 0.000...
AO_LDOS: 1 1 -14.92 0.000 0.000 0.000 0.000 0.000...
...
PDOS/AOLDOS can be extracted by using the state2pdos.pl s...
state2pdos.pl nfout_aoldos
then the following messages are printed
make dos data file :: 0001
make dos data file :: 0002
make dos data file :: 0003
make dos data file :: 0004
make dos data file :: 0005
make dos data file :: 0006
meaning that the dos data are printed to files 0001, 0002...
NOTE: with the latest version of the code, they will be p...
The data are given in the above files in the following or...
Energy s px py pz dzz dxx-yy dxy ...
Finally plot the data using for e.g., gnuplot or xmgrace.
Alternatively, you can perform SCF and AOLDOS calculation...
0 0 0 0 0 0 : I_CTROL (D...
6.00 20.00 2 6 6 : GMAX GMAXP...
1 0 : NUM_SPAC_G...
Cartesian
5.498778336527 0.000000000000 0.00000000...
-2.749389168264 4.762081729212 0.00000000...
0.000000000000 0.000000000000 85.30494051...
12 12 1 1 1 1 : NK1 NK2 NK...
1 0 : NCORD NINV...
0.000000000000 0.000000000000 8.97946742...
0.000000000000 3.174721152808 4.48973371...
0.000000000000 -3.174721152808 0.00000000...
0.000000000000 0.000000000000 -4.48973371...
0.000000000000 3.174721152808 -8.97946742...
0.000000026956 -3.174721183174 -10.61015812...
79.00 0.50 12.0107 1 1 0.0 : IATOMN ALF...
1.00 0.50 12.0107 1 1 0.0 : IATOMN ALF...
0 0 0 0 0 : ICOND INIPO...
0 1 : IPRE IPRI
200 200 0 80000.00 0 : NMD1 NMD2 ...
6 1 : WAY_MIX MI...
0 20 0.60 : ITER_START...
0.20 0.30 0.20 0.20 0.20 : DTIM1 DTIM...
200.00 2 1 1.0D-12 : DTIO IMDAL...
-0.0010 1.00D-03 0 : WIDTH FORC...
rev-vdW-DF2 1 : XCTYPE KSPIN
2.00 : DESTM
101 : NBZTYP
4 4 4 : DUMMY
4 4 4 : DUMMY
48 : KEG
1 : NEXTST
0 : DUMMY
2 : IMSD
0 : EVALUATE_E...
6 : NPDOSAO
1
2
3
4
5
6
-15.00 5.00 0.20 501 : EPDOS(1) E...
2.2 0.3 : RAD WIDTH ...
0.7 0.3 : RAD WIDTH ...
0.2 12 : DR NR
0 0.000 : SM_N DOPING
-Visualization
Let us visualize PDOSs for the bottom Au and H atoms usin...
Here is an example gnuplot script to generate an eps file...
set terminal postscript eps color 'Helvetica' 20
set output 'pdos.eps'
xmin=-11.0
xmax=4.0
ymin=0.0
ymax=2.0
set xrange [xmin:xmax]
set yrange [ymin:ymax]
set yzeroaxis
set xlabel '{/Helvetica-Oblique E}-{/Helvetica-Oblique E...
set ylabel 'PDOS (Arb. unit)'
plot 'pdos_0005.dat' using ($1):($11) title 'Au#5 s' wit...
'pdos_0005.dat' using ($1):($12) title 'Au#5 p' wit...
'pdos_0005.dat' using ($1):($13) title 'Au#5 d' wit...
'pdos_0006.dat' using ($1):($2) title 'H s' with li...
Type
> gnuplot pdos.gp
to generate the eps file called "pdos.eps" as shown below
#ref(http://www-cp.prec.eng.osaka-u.ac.jp/puki_state/grap...
** Band structure calculation [#l6de7d56]
The band structure calculation can be performed after the...
- Input file (nfinp_band)
0 0 0 0 0 0 : I_CTROL (D...
6.00 20.00 2 6 6 : GMAX GMAXP...
1 0 : NUM_SPAC_G...
Cartesian
5.498778336527 0.000000000000 0.00000000...
-2.749389168264 4.762081729212 0.00000000...
0.000000000000 0.000000000000 85.30494051...
12 12 1 1 1 1 : NK1 NK2 NK...
1 0 : NCORD NINV...
0.000000000000 0.000000000000 8.97946742...
0.000000000000 3.174721152808 4.48973371...
0.000000000000 -3.174721152808 0.00000000...
0.000000000000 0.000000000000 -4.48973371...
0.000000000000 3.174721152808 -8.97946742...
0.000000026956 -3.174721183174 -10.61015812...
79.00 0.50 12.0107 1 1 0.0 : IATOMN ALF...
1.00 0.50 12.0107 1 1 0.0 : IATOMN ALF...
22 0 0 0 0 : ICOND INIP...
0 1 : IPRE IPRI
200 200 0 80000.00 0 : NMD1 NMD2 ...
6 1 : WAY_MIX MI...
0 20 0.60 : ITER_START...
0.20 0.30 0.20 0.20 0.20 : DTIM1 DTIM...
200.00 2 1 1.0D-12 : DTIO IMDAL...
-0.0010 1.00D-03 0 : WIDTH FORC...
rev-vdW-DF2 1 : XCTYPE KSPIN
2.00 : DESTM
101 : NBZTYP
4 4 4 : DUMMY
4 4 4 : DUMMY
48 : KEG
1 : NEXTST
0 : DUMMY
2 : IMSD
0 : EVALUATE_E...
0 : NPDOSAO
0 0.000 : SM_N DOPIN...
!
&KPOINTS_BAND
NKSEG 3
KMESH 40 20 40
KPOINTS
0.0 0.0 0.0
0.5 0.0 0.0
0.33333333 0.33333333 0.0
0.0 0.0 0.0
&END
For the band structure calculation, ICOND is set to 22, a...
&KPOINTS_BAND
NKSEG 3
KMESH 40 20 40
KPOINTS
0.0 0.0 0.0
0.5 0.0 0.0
0.33333333 0.33333333 0.0
0.0 0.0 0.0
&END
where the variables in the &KPOINTS_BAND ... &END block are
- the number of k-point segment (NKSEG)
- the k-point mesh for each segment (KMESH)
- the k-points which define the symmetry k-points (# of k...
K-points are given in the unit of the reciprocal lattice ...
After the band structure run, the file "energy.data" is g...
To convert the "energy.data" into plottable one, "energy2...
Type
> energy2band
and input number of bands, number of bands to be plotted,...
Here is an example script to run energy2band:
#!/bin/sh
E2B=energy2band
scfout_file='nfout_scf'
bandout_file='nfout_band'
ef=`grep 'FERMI' ${scfout_file} | awk '{print $4}'`
nb=`grep 'NUMBER OF BANDS CONSIDERED' ${bandout_file} | ...
nk=`grep 'NUMBER OF K-POINTS' ${bandout_file} | grep -v...
${E2B} << EOF
${nb}
${nb}
${nk}
${ef}
EOF
By using the following gnuplot script, the band structure...
#ref(http://www-cp.prec.eng.osaka-u.ac.jp/puki_state/grap...
The following gnuplot script can be used to generate the ...
emin=-12
emax=4
xmin=0.0000
xmax=1.80236106
G=0.00000000
M=0.65970994
K=1.04059365
G2=1.80236106
offset=-0.5
set xrange [xmin:xmax]
set yrange [emin:emax]
set arrow from M,emin to M,emax nohead lt 1
set arrow from K,emin to K,emax nohead lt 1
set xzeroaxis
unset xtics
unset key
set ylabel '{/Helvetica-Oblique E-E}_F (eV)'
set label '{/Symbol G}' at G,emin+offset center
set label 'M' at M,emin+offset center
set label 'K' at K,emin+offset center
set label '{/Symbol G}' at G2,emin+offset center
#set terminal postscript eps 'Helvetica' 20
#set output 'band.eps'
plot 'band.data' with lines lw 3
pause -1
ページ名: