I found that I used the unit of system as [unit_system = 'a.u.'] instead of [unit_system = 'A_eV_fs']
I re-calculated ground state of CrI3 with [unit_system = 'A_eV_fs'], and the result was changed. However the band gap is still negative.
*****************************************
Inho Lee, Ph.D
Theoretical Division, T-4
Los Alamos National Laboratory
Los Alamos, New Mexico 87545
Dear Salmon developers and users,
I am Inho Lee, postdoc in Los Alamos National Lab.
I performed the spin-polarized LDA calculation of CrI3 monolayer, however it shows unphysical results with negative band gap.
(Moreover, DOS is not correct.)
Therefore, I tried calculation with larger r-grid, however the band gap is still negative.
I used r-grid with ~0.19 Angstrom (a,b-axis) and ~0.205 Angstrom(c-axis).
Also I used 8*8*1 kgrid.
Is r-grid with 0.2 Ang not enough? or Is there inappropriate option in my input file?
(I checked the atomic coordinate, and crystal structure is correct)
I attached input/pseudopotential/output(w/o intermediate iteration data) files.
I welcome any comments and suggestions,
Thank you,
Sincerely,
Inho Lee
*****************************************
Inho Lee, Ph.D
Theoretical Division, T-4
Los Alamos National Laboratory
Los Alamos, New Mexico 87545