Перейти к содержанию
Форум химиков на XuMuK.ru
β

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


JackOo

Рекомендуемые сообщения

🚑 Решение задач, контроши, рефераты, курсовые и другое! Онлайн сервис помощи учащимся. Цены в 2-3 раза ниже!

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

Нужно рассчитать в абините зонную структуру Бора, гранецентрированная кубическая решётка, а = 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 =

#%%

 

 

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

Изменено пользователем JackOo
Ссылка на комментарий

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

Изменено пользователем mypucm
Ссылка на комментарий
  • 2 недели спустя...

Преподаватель сказал мне, что в моём случае в элементарной ячейке будет 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 векторов.

Ссылка на комментарий

Для публикации сообщений создайте учётную запись или авторизуйтесь

Вы должны быть пользователем, чтобы оставить комментарий

Создать аккаунт

Зарегистрируйте новый аккаунт в нашем сообществе. Это очень просто!

Регистрация нового пользователя

Войти

Уже есть аккаунт? Войти в систему.

Войти
  • Последние посетители   0 пользователей онлайн

    • Ни одного зарегистрированного пользователя не просматривает данную страницу
×
×
  • Создать...