ABINIT Зонная структура

[JackOo]

Всем привет, хочу попросить о помощи. Готов взять консультацию и заплатить за неё.

Нужно рассчитать в абините зонную структуру Бора, гранецентрированная кубическая решётка, а = 2.3284 ангстрем(параметр решётки).

Цель: сделать адекватный и правильный входной файл .in для моей задачи, чтобы получить выходной файл .out, как из .out получить график я знаю.

Брал из туториала файл для силикона и менял его, но всё равно не получил зонной структуры.

Вот файл, #red отмечено где я внёс изменения.

rystalline silicon

#

# Computation of the band structure.

# First, a SCF density computation, then a non-SCF band structure calculation.

 

ndtset 2

 

#Dataset 1 : usual self-consistent calculation

kptopt1 1 # Option for the automatic generation of k points,

# taking into account the symmetry

nshiftk1 4

shiftk1 0.5 0.5 0.5 # These shifts will be the same for all grids

0.5 0.0 0.0

0.0 0.5 0.0

0.0 0.0 0.5

ngkpt1 4 4 4

prtden1 1 # Print the density, for use by dataset 2

toldfe1 1.0d-6 # This value is way too large for most realistic studies of materials

 

#Dataset 2 : the band structure

iscf2 -2

getden2 -1

kptopt2 -3

nband2 8

ndivk2 10 12 17 # 10, 12 and 17 divisions of the 3 segments, delimited

# by 4 points.

kptbounds2 0.5 0.0 0.0 # L point

0.0 0.0 0.0 # Gamma point

0.0 0.5 0.5 # X point

1.0 1.0 1.0 # Gamma point in another cell.

tolwfr2 1.0d-12

enunit2 1 # Will output the eigenenergies in eV

 

 

#Definition of the unit cell

acell 3*2.3284 #red #My params # This is equivalent to 10.217 10.217 10.217

rprim 0.0 0.5 0.5 # FCC primitive vectors (to be scaled by acell)

0.5 0.0 0.5

0.5 0.5 0.0

 

#Definition of the atom types

ntypat 1 # There is only one type of atom

znucl 14 #red # The keyword "znucl" refers to the atomic number of the

# possible type(s) of atom. The pseudopotential(s)

# mentioned in the "files" file must correspond

# to the type(s) of atom. Here, the only type is Silicon.

 

 

#Definition of the atoms

natom 1 #red #One atom # There are two atoms

typat 1 #red #Bor # They both are of type 1, that is, Silicon.

xred #red # This keyword indicate that the location of the atoms

# will follow, one triplet of number for each atom

0.0 0.0 0.0 # Triplet giving the REDUCED coordinate of atom 1.

# Triplet giving the REDUCED coordinate of atom 2.

 

#Definition of the planewave basis set

ecut 8.0 # Maximal kinetic energy cut-off, in Hartree

 

#Definition of the SCF procedure

nstep 10 # Maximal number of SCF cycles

diemac 12.0 # Although this is not mandatory, it is worth to

# precondition the SCF cycle. The model dielectric

# function used as the standard preconditioner

# is described in the "dielng" input variable section.

# Here, we follow the prescription for bulk silicon.

# add to conserve old < 6.7.2 behavior for calculating forces at each SCF step

optforces 1

 

 

#%%

#%% [setup]

#%% executable = abinit

#%% [files]

#%% files_to_test =

#%% tbase3_5.out, tolnlines= 0, tolabs= 0.000e+00, tolrel= 0.000e+00

#%% psp_files = 14si.pspnc

#%% [paral_info]

#%% max_nprocs = 1

#%% [extra_info]

#%% author =

#%% keywords =

#%% description =

#%%

 

 

Был бы очень признателен, если бы вы смогли дать контакты человека, который может решить мою задачу.

Изменено 6 Марта, 2013 в 10:23 пользователем JackOo

[Скорпион]

А в на каком програмном языке это написано?

[aversun]

Вероятно на внутреннем скриптовом программы.

[JackOo]

Сам ABINIT на фортране, но к делу это отношения не имеет.

aversun прав, это внутренний язык ABINIT.

[JackOo]

Кто-нибудь может посоветовать специалиста по моему вопросу?

[mypucm]

Вы меня заинтересовали. Я не специалист, но очень хочется попробовать потыркаться и посчитать структуру некого гипотетического кристаллического вещества, которое еще никем не было получено. Конкретно - K[Fe(Fe(CN)₆)]*H₂O (не верьте никому, кто скажет, что это вещество получить - раз плюнуть, старые учебники нагло врут!). Начинаю качать ABINIT. Может, сам тут в этой теме потом буду вопить о помощи, может, и сам разберусь, и Вам помогу!

Изменено 7 Марта, 2013 в 08:57 пользователем mypucm

[JackOo]

Всем привет, турист, есть успехи?

[JackOo]

И ещё вопрос, каким образом сегодня считают зонную структуру? У меня есть подозрения что ABINIT устарел и плохо поддерживается.

[mypucm]

JackOo, я в этих вещах ньюб, и сейчас потихоньку изучаю книжку, рекомендованную авторами программы. А она что-то объемистая оказалась.

[JackOo]

Преподаватель сказал мне, что в моём случае в элементарной ячейке будет 4 атома, на http://rruff.geo.arizona.edu/AMS/result.php , я нашёл только 5 и 2 атома.

И откуда вообще берут координаты атомов в таких случаях?

 

 

#Definition of the unit cell

acell 3*4.400 #red #My params # This is equivalent to 10.217 10.217 10.217

\\

rprim

0.0 0.5 0.5 # FCC primitive vectors (to be scaled by acell)

0.5 0.0 0.5

0.5 0.5 0.0

 

Здесь описана ГЦК решётка с параметром 4.40 а.е.

\\\

 

Мне нужно изменить значения znucl c 14 на 5

natom на 4(или 5, ещё вопрос где взять координаты этих 4-х атомов)

в typat внести (1 1 1 1) или (1 1 1 1 1)

и в xred записать 4-5 векторов.