Учебный сайт Смирновой Виктории

1. Оптимизированная структура этана в виде z-matrix dj входном файле для GAMESS с шапкой для DFT: et.inp (вместо значений длин и углов связей стоят переменные).
   

   C помощью скрипта были созданы файлы с различными значениями длины связи с шагом 0.02 ангстрема. По файлам был запущен рассчет GAMESS, из полученных log-файлов были извлечены длины связей и рассчитанные для них значения энергии, эти данные были сохранены в файле bond.
   
   По зависимости энергии молекулы от длины одной связи был построен график в Gnuplot, для чего проводилась подгонка коэффициентов функции

   f(x)=a + k*x*x - 2*k*x*b + k*b*b 

   Полученные значения коэффициентов:

   Final set of parameters            Asymptotic Standard Error
   =======================            ==========================
   
   a               = -79.7652         +/- 0.0004522    (0.000567%)
   k               = 0.563608         +/- 0.02335      (4.142%)
   b               = 1.55432          +/- 0.002455     (0.1579%)
   

   График:
   
   
   

2. Аналогичные операции были проделаны для валентного угла HCH, его значения изменялись от 109.2 до 113.2. Углы и рассчитанные для них значения энергии были сохранены в файле angle. Полученные значения коэффициентов:

   Final set of parameters            Asymptotic Standard Error
   =======================            ==========================
   
   a               = -79.7647         +/- 1.21e-08     (1.517e-08%)
   k               = 3.56076e-05      +/- 6.229e-09    (0.01749%)
   b               = 111.38           +/- 9.954e-05    (8.937e-05%)
   
   

   График:
   
   
   

3. Аналогичные операции были проделаны для торсионного угла d3, его значения изменялись от -180 до 180 c шагом 12. Углы и рассчитанные для них значения энергии были сохранены в файле torsion. Подгонка не проводилась. Количество минимумов функции: 4 (значения угла -180, -60, 60, 180), но для угла значения -180 и 180 совпадают. График: