VASP+phonopy (Fe2VAl-dfpt)

------------------------------------------------------------------------------
■ Youtube動画:
基本的な使い方:https://www.youtube.com/watch?v=bKwjCIrvfSk
エラーが表示された場合:https://www.youtube.com/watch?v=5Y6tLWvlLJE
------------------------------------------------------------------------------
■ phonon calculation for Fe2VAl (Primitive cell) by DFPT method
Need: POTCAR and POSCAR-unitcell
(Not need: KPOINT (for vasp verasion 5.3)

1. mkdir Fe2VAl-dfpt
2. cd Fe2*
3. gedit POSCAR-unitcell
-----PSOCAR-unitcell----
Generated by cif2cell 1.2.10.  :  Failed to get author information, No journal information.
 1.0
  0.000000000000000   2.882550000000000   2.882550000000000
  2.882550000000000   0.000000000000000   2.882550000000000
  2.882550000000000   2.882550000000000   0.000000000000000
  Al  Fe   V
   1   2   1
Direct
  0.500000000000000   0.500000000000000   0.500000000000000
  0.250000000000000   0.250000000000000   0.250000000000000
  0.750000000000000   0.750000000000000   0.750000000000000
  0.000000000000000   0.000000000000000   0.000000000000000
-------------
4. gedit INCAR
-----INCAR----
   PREC = Accurate
#  ENCUT = 500
 IBRION = 8
  EDIFF = 1.0e-08
  IALGO = 38
 ISMEAR = 0; SIGMA = 0.1
  LREAL = .FALSE.
ADDGRID = .TRUE.
  LWAVE = .FALSE.
 LCHARG = .FALSE.
-------------
5. phonopy -d --dim="2 2 2" -c POSCAR-unitcell
6. mv SPOSCAR POSCAR
7. mpirun -np 8 /home/wien2k3/vasp/vasp.5.3/vasp
8. phonopy --fc vasprun.xml

# band dispersion
1. gedit band.conf
-----band.conf----G-X-K-G-L
DIM = 2 2 2
PRIMITIVE_AXIS = 1.0 0.0 0.0  0.0 1.0 0.0  0.0 0.0 1.0
BAND = 0.0 0.0 0.0  0.0 0.5 0.5  0.375 0.75 0.375  0.0 0.0 0.0  0.5 0.5 0.5
FORCE_CONSTANTS = READ
-------------
* BAND is half as for 2x2x2 as for 1x1x1
2. phonopy --dim="2 2 2" -c POSCAR-unitcell band.conf
3. cp $HOME/phonopy*/scripts/bandplot ./
4. bandplot


# other bandplot option
◆ bandplot -h
◆ bandplot --xlabel="" --ylabel="Frequency (THz)" --title="Fe2VAl" --fmin=0 --fmax=12
◆ bandplot --gnuplot
◆ bandplot --xlabel="" --ylabel="Frequency (THz)" --title="Fe2VAl" --fmin=0 --fmax=12 --dos=total_dos.dat
◆ bandplot --xlabel="" --ylabel="Frequency (THz)" --title="Fe2VAl" --fmin=0 --fmax=12 --dos=partial_dos.dat

Fig. band dispersion for G-X-K-G-L line with "bandplot --xlabel="" --ylabel="Frequency (THz)" --title="Fe2VAl" --fmin=0 --fmax=12 --dos=partial_dos.dat".

# TDOS and Cv (Specific heat)
1. gedit mesh.conf
-----mesh.conf-----
DIM = 2 2 2
MP = 16 16 16
FORCE_CONSTANTS = READ
-------------
2. phonopy --dim="2 2 2" -c POSCAR-unitcell -p mesh.conf

3. phonopy --dim="2 2 2" -c POSCAR-unitcell -t -p mesh.conf


# PDOS
1. gedit pdos.conf
-----pdos.conf----plot pdos (Al, Fe, V, Total)-----
DIM = 2 2 2
MP = 16 16 16
PDOS = 1, 2 3, 4, 1 2 3 4
FORCE_CONSTANTS = READ
-------------
Maybe, relate POSCAR-unitcell
-----POSCAR-unitcell----
  Al  Fe   V
   1   2   1
-------------
1=Al, 2=Fe, 3=Fe, V=4
2. phonopy --dim="2 2 2" -c POSCAR-unitcell -p pdos.conf

------------------------------------------------------------------------------
■ benchmark
-----------
Fe2VAl(primitive cell) DFPT calculation
Openmpi + intel compiler + MKL
mpirun -machinefile $HOME/machinefile -genv I_MPI_DEVICE sock -n 4 /home/student/vasp/vasp.5.3/vasp -r ssh
Intel(R) Core(TM) i7-4770 (8M Cache, 3.40 GHz, up to 3.90 GHz)
(cpu value $303.00 - $312.00)
Elapsed time (sec):      1171.669 = 19.5 min
Maximum memory used (kb):      248236. = 0.25 GB
-----------
Hyper threading OFF (on BIOS)
※ Hyper threading OFFの有無を調べたが、計算速度に違いはない。加えて、PC安定に動作する(ハングアップしたりシャットダウンしたりする頻度が低下するわ)けではなかった(PCがハングアップしたりシャットダウンしたりしなくなるようにするにはそれよりも並列コア数を最大の半分くらいに指定したりHDDやメモリをクリーンにする方が効果がある)。
Fe2VAl(primitive cell) DFPT calculation
Openmpi + intel compiler + MKL
mpirun -machinefile $HOME/machinefile -genv I_MPI_DEVICE sock -n 4 /home/student/vasp/vasp.5.3/vasp -r ssh
Intel(R) Core(TM) i7-4770 (8M Cache, 3.40 GHz, up to 3.90 GHz)
(cpu value $303.00 - $312.00)
Elapsed time (sec):      1138.405 = 18.9 min
Maximum memory used (kb):      246372. = 0.25 GB
-----------
Fe2VAl(primitive cell) DFPT calculation
intelmpi + intel compiler + MKL
mpirun -machinefile $HOME/machinefile -genv I_MPI_DEVICE sock -n 4 /home/student/vasp/vasp.5.3/vasp -r ssh
Intel(R) Core(TM) i7-4770S (8M Cache, 3.10 GHz, up to 3.90 GHz)
(cpu value $303.00 - $305.00)
Elapsed time (sec):      931.538 = 15.6 min
Maximum memory used (kb):      288624. = 0.29 GB
-----------
Fe2VAl(primitive cell) DFPT calculation
intelmpi + intel compiler + MKL
mpirun -machinefile $HOME/machinefile -genv I_MPI_DEVICE sock -n 4 /home/student/vasp/vasp.5.3/vasp -r ssh
Intel(R) Xeon(R) E5-2640 v3 (20M Cache, 2.60 GHz, up to 3.40 GHz)
(cpu value $939.00 - $944.00)
Elapsed time (sec):      489.297 = 8.15 min
Maximum memory used (kb):      206484. = 0.21 GB
----------Results----------
Sorry, I didnt memorize second time (hh:mm:ss)
3.1 GB
start main + 3 part
Sat Nov 19 9:48:XX JST 2016, DAV, main, 17+3 cycle
Sat Nov 19 9:49:XX JST 2016, RMM, 1/3, 18 cycle
Sat Nov 19 9:53:XX JST 2016, RMM, 2/3, 18 cycle
Sat Nov 19 9:55:XX JST 2016, RMM, 3/3, 17 cycle
Sat Nov 19 9:57:XX JST 2016, END
primitive cell / about 9 min
------------------------------
-----------
------------------------------------------------------------------------------
QRコード
携帯用QRコード
アクセス数
ページビュー数
[無料でホームページを作成] [通報・削除依頼]