Physics Laws

Можно спорить о формулировках законов Эйнштейна или термодинамики, однако, я думаю, в подобной форме они лучше запоминаются и значительно понятнее. Мне особенно интересно было узнать о третьем законе Эйнштейна и четвёртом законе термодинамики.

eko

Я продолжаю делать расчёты. Пока не удаётся смоделировать адсорбцию и в связи с этим я осваиваю теорию. А также делаю новогодние игрушки и мастерю сайт eko: исправил скрипт, подправил css и добавил link. Еще не хватает фоновых картинок.

Гаусс 03 (the next day)

Обнаружил, что оптимизация геометрии очень сильно влияет на время расчета. Так для вчерашних попыток получил новые данные:

число атомов	время		энергия, a.u.
	вчера	сегодня	вчера	сегодня
1	6 секунд	–	−48.0796518	–
7	12 минут 14 секунд	5 минут 24 секунды	−48.06444983	−48.06769137
13	99 минут 5 секунд	45 минут 12 секунд	−48.06769137	−48.06769243
20	–	353 минуты 15 секунд	–	−48.06943069

Кристалл кадмия задаётся значения: α = β = 90°, γ = 120°, a = b = 2.9793 Å, c = 5.6196 Å. Данные из CDC Handbook of Chemistry and Physics, 85th Edition и Smithell's Metals Reference Book, 7th Edition. Довольно долго искал эту информацию. Впредь буду ею пользоваться.

Гаусс 03

Обнаружил ошибку во вчерашних расчетах. Я не правильно моделировал кристаллическую решетку кадмия.

Сегодня тестирую "Гаусса 03". Для расчета кластера (cluster, est – klaster) из 7 атомов при помощи DFT (B3LYP) на основе базиса (LanL2DZ) потребовалось:

1. 1 атом – 6 секунд (−48.0796518 a.u.);
   
2. 7 атомов – 12 минут 14 секунд (−48.06444983 a.u.);
3. 13 атомов – 99 минут 5 секунд (−48.06769137 a.u.).

ADF computations

Вчерашний вечер и сегодняшнее утро прошло у компьютера - я изучаю работу ADF! Tutorial Band GUI был сразу пройден и, вот, уже сегодня я создаю собственные модели.

• На расчет пластины (slab, est. – plaat) из атомов Cd с использованием TZ2P (triple Z, 2 polarization function, without core) потребовалось:
  • 1 слой – 3 минуты 54 секунды (−28.8891 kcal/mol, −3.7554 eV);
  • 2 слоя – 43 минуты 43 секунды (−29.3295 kcal/mol, −4.1003 eV);
  • 3 слоя – 83 минуты 57 секунд (−29.5215 kcal/mol, −4.2531 eV);
  • 4 слоя – 249 минуты 13 секунд (−29.6199 kcal/mol, −4.4081 eV);
    
• На расчет кристалла (bulk, est. – kogu, ehk tahkis) из атомов Cd с использованием той же базы понадобилось 46 минут 53 секунды (−29.9198 kcal/mol, −4.8985 eV).

В скобках указана энергия образования (эксп. энтальния образования (CDC Handbook) −26,13 kcal/mol) и через запятиую – энергия Ферми.

К сожалению, мне так и не удалось задефинировать I⁻-ион. Мне также не удалось пока запустить .run файл для расчета адсорбции атомов йода на поверхности кадмиевой пластины.

Все это подталкнуло меня опробовать еще раз Gaussian. Особенно интересно сравнить базисы. Интересно, что в Gaussian 03 я могу использовать: STO-3G (P, –), 3-21G (P, D), SDDAll, LanL2DZ (–, –) и DGDZVP (–, –). Видимо, базисы от ADF всё же выигрывают.