Расчет ширины запр. зоны Eg в Gaussian с PBC.

[IndexOutOfBounds]

Я сделал попытку посчитать запрещенную зону у бесконечного графенового слоя. Для этого задается простая графеновая элементарная ячейка и делается расчет.

%chk=Graphene STO.chk
%mem=60MW
%nproc=1
#p pbepbe/sto-3g*/auto

Graphene

0 1
C -0.00000000 0.00000000 0.00000000
C 2.45951210 1.41999997 -0.00000000
Tv 2.45951210 0.00000000 0.00000000
Tv 1.22975605 2.12999996 0.00000000

В файле out я получаю такие данные:

HOCO at k-point 260 kxyz= 7 32 1 -3.219791 eV
LUCO at k-point 260 kxyz= 7 32 1 -2.940898 eV
Indirect gap 0.278893 eV
Max direct gap at k-point 511 kxyz= 22 25 1 14.808817 eV
Min direct gap at k-point 260 kxyz= 7 32 1 0.278893 eV


Вопрос, почему здесь есть gap ведь известно, что Eg в графене равна 0? Если считать в PBEPBE (что по умолчанию), то Eg порядка 0.3-0.4.

[Darth Vasya]

Очевидно, потому, что вы не попали в точку K зоны Бриллюэна:

> Min direct gap at k-point 260 kxyz= 7 32 1 0.278893 eV

[Darth Vasya]

Кстати, это тоже странно:

> Max direct gap at k-point 511 kxyz= 22 25 1 14.808817 eV

---
Вообще, перестаньте использовать гауссиан в задачах, для которых он не приспособлен.

[IndexOutOfBounds]

А вы мне не подскажите программу для таких рассчетов? Еще интересно что за координаты kxyz он пишет. Как оси в обратном пространстве направлены и масштаб применительно к зоне Бриллюэна. А я-то как раз не против попасть в точку K, но Гауссиан неправильно представляет себе или дисперсию энергии или зону Бриллюэна.

[Himera]

По поводу программ -- посмотрите вот в этой теме

[Darth Vasya]

Для таких расчётов подойдёт много программ. Вы лучше сразу опишите свою "большую" задачу целиком, чтобы можно было относительно этого уже что-то советовать.

[fixatorx]

Я использую SIESTA для таких вещей. А вы можете либо ее, либо Abinit.

[dhvol]

Помогите мне! Где можно найти псевдопотенциал атома B (Boron) для расчета в SIESTA . В файле test of Siesta бывает папка born , внутри этой папке есть file born.pseudos --> Si ( Si.psf) . Я взял это файл для расчета, но результат расчета неправильно. Где ошибка ? thanks you very much.