졸논 D+2,3 / silicon 1

1. silicon 0의 결과를 기준으로 최적의 k 값 결정.

Monk k 값이 6-10 사이일때
에너지와 구조가 수렴 했고 + 시간도 급격하게 늘어나지 않았다.
let’s
k = 9

2. ENCUT 결정.

POTCAR file 안에

– ENMAX  =  245.345

INCAR file 안에

ENCUT을 250 ~ 650 으로 바꾸어 보면서
최적화를 해 본후 -> 값들이 수렴 시작  + 시간 최대한 적게 걸리는 ENCUT   선택

그래프 넣기!

그런데 이상하게도 연속적인 그래프가 나온게 아니고 450에서 최적화된 lattice constant가 이상하리만큼 큰 변화가 나옴. 그럼에도 불구하고 에너지는 연속적으로 특정 값에 수렴함

그래도 끝까지 해 보기 위해

ENCUT 500으로 선택!

3. 이번엔 원자중심은 안움직인 채로 최적화
( silicon 0이나 위의 2번 항목에서의 최적화는 다음의 과정을 반복
:POSCAR에 저장된 전자 정보와 원자 중심과의  전자 최적화 -> 전자 밀도 최적화 정보로 원자 중심 or lattice위치 조정 -> 전자 원자 관계이용하여 전자 밀도 재조정 -> 재조정 정보로 원자 중심 재조정 ->)

그런데, 이 과정을 중심이 움직이는게 마지막,

좀 더 정확한 전자 밀도 정모를 알고 싶으면 중심을 고정 해 놓고 최적화를 한 번 더 해야 한다.

CONTCAR -> POSCAR로 새로 지정하고

INCAR file
– NSW = 0

으로 해 놓으면 중심움직이지 말라는 이야기

INCAR

ENCUT = 500    #위에서 찾은 값
ENAUG = 1000

NSW = 0             # 0  은 중심 움직이지 말라는 명령
IBRION = 2

LWAVE = .T.        # 지금 까지는 이 정보들 안 받았는데 이제는 받을 거라는 뜻! TRUE
LCHARG = .T.     # CHG, CHGCAR  가 이제 파일 크기가 0이 아닐 거야

ISMEAR = 0
SIGMA = 0.05

ISIF = 3
LREAL = A

$grep E-fermi OUTCAR 로 페르미 에너지도 볼 수 있음.

sub_vasp.pl 로 최적화!

4. 밴드 구조 그릴 정보 받기

INCAR file

ENCUT = 500
ENAUG = 1000

NSW = 0
IBRION = 2

LWAVE = .F.      # 이미 구해 놨으니까 이제는 다시 받을 필요 없다
LCHARG = .F.   # 마찬가지

ISMEAR = 0
SIGMA = 0.05

ISIF = 3
LREAL = A

ICHARG = 11  # 1 (CHGCAR 정보를 써라 + 10 최적화 실행 하지 말아라, CHGCAR 정보를 바꾸지 말아라)

KPOINTS file

Line
20                           #
Line
reciprocal               # reciprocal 좌표들은 쓴다.
0.5    0.5    0.5 ! L
0.0    0.0    0.0 ! G

0.0    0.0    0.0 !G
0.0    0.5    0.5 !X

0.0    0.5    0.5 !X
0.25    0.75    0.5 !W

0.25    0.75    0.5 !W
0.475    0.75    0.475 !K

0.475    0.75    0.475 !K
0.0    0.0    0.0 !G

두 개 씩 짝지어진 좌표들을 잇는 선을 따라 eigenvalue를 구하도로 명령하는 것
L G X W K 는 나름 유명한 reciprocal 좌표,

위의
CHGCAR :최적화 된 charge density 정보
INCAR : 어떤걸 어느 범위안에서구해 볼 지
POTCAR : Si 정보 provided by vasp 사실 먼지 잘 몰라
POSCAR : 최적화된 lattice 정보
KPOINTS : 공간을 연속적으로 계산 할 수 없으니까, 계산 할 points 들을 결정해 주는 grid 정보

를 이용해서
sub_vasp.pl 을 돌리면!!

EIGENVAL file에
energy eigenvalue값 저장 되어 있다
그걸로 밴드 그래 보면 됨!

5.  silicon band 를 그려보자

Leave a Reply