Базисный набор UGBS

[Droog_Andrey]

Оказалось, это просто набор экспонент вида exp(0.12+0.672k), где k принимает целые значения. Для каждого элемента и типа симметрии базисных функций (s, p, d, f, g, ...) зашит определённый диапазон значений k. Например, для нобелия в стандартном UGBS используются следующие диапазоны:

s: от -5 до 26
p: от -4 до 20
d: от -2 до 16
f: от -3 до 11

Если указать поляризационные функции (UGBS1P), то диапазон p-функций расширяется до диапазона s-функций стандартного базиса и т.д.

Повозившись немного, я пришёл к выводу, что если нужен полностью универсальный (т.е. одинаковый вообще для всех атомов) базисный набор, то можно взять следующие диапазоны значений k:

s: от -7 до 28
p: от -7 до 24
d: от -7 до 20
f: от -7 до 16
g: от -7 до 12

Таким образом, получается набор 36s 32p 28d 24f 20g ёмкостью 36 + 32*3 + 28*5 + 24*7 + 20*9 = 620 базисных функций. Расчёт энергии функционалом PBE1PBE для одного атома в этом базисе занимает время от нескольких минут до часа в зависимости от числа электронов.

Вот пример расчёта для атома бория (Z=107) в основном состоянии: http://www.primefan.ru/stuff/chem/bh.rar

[VTur]

Это сжатый базис или полный?

[Droog_Andrey]

В наборе UGBS не применяли сжатие, чтобы не терять универсальность. Просто выбрали оптимальные диапазоны k для каждого элемента.

[Droog_Andrey]

А ещё проще взять тупо степени двойки в качестве экспонент:

viewtopic.php?f=71&t=83489