Quantum Espresso를 이용해 NiO 같은 물질의 DFT+U+V (Dudarev) 계산을 하기위해서는 pw.x 계산 input 으로 다음과 같은 parameters 를 설정하여야 한다.
자세한 계산과정은 youtube의 DFT+U+V tutorial을 참조하도록하자.
.
.
&SYSTEM
.
.
nspin = 2
lda_plus_u = .true.
lda_plus_u_kind = 2
hubbard_v(1,1,1) = 1.0000E-10
hubbard_v(2,2,1) = 1.0000E-10
hubbard_v(3,3,1) = 1.0000E-10
starting_magnetization(1) = 1.0000E-10
starting_magnetization(2) = 1.0000E-10
u_projection_type = 'ortho-atomic'
.
.
ATOMIC_POSITIONS angstrom
Ni1 0.0000000000 0.0000000000 0.0000000000
Ni2 2.0972999043 2.0972999043 2.0972999043
O 1.0486499522 1.0486499522 1.0486499522
O 3.1459498565 3.1459498565 3.1459498565
.
.hubbard_v(na,nb,k) 에서 na,nb 는 ATOMIC_POSITIONS에 나와 있는 원자의 순서를 나타낸다.
k=1 은 na, nb 원자의 standard orbitals 간의 interaction
k=2 는 na 원자의 standard orbital과 nb 원자의 background orbital 간의 interaction
k=3 는 na, nb 원자의 background orbital 간의 interaction
k=3 는 na 원자의 background orbital, 원자의 standard orbital 간의 interaction 을 의미한다.
aiida-quantumespresso를 사용하는 경우, pw.x input parameter dict를 다음과 같이 설정하여 준다.
.
.
'SYSTEM': {
.
.
'nspin' : 2,
'starting_magnetization': {
'Ni1': 1.0,
'Ni2': -1.0,
},
'lda_plus_u': True ,
'lda_plus_u_kind' : 2,
'U_projection_type' : 'ortho-atomic',
'hubbard_v': {
(1, 1, 1, 1.0E-10), #(na,nb,k, Energy)
(2, 2, 1, 1.0E-10),
(3, 3, 1, 1.0E-10),
},
.
.형태로 parameters 를 입력해 준다.