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

Главная Проекты Семестры


Вычисление точечных зарядов и VdW параметров для молекулярной механики



    1. С помощью babel был сделан pdb файл этана, из результатов оптимизации из предыдущего практикума. Добавим путь к скриптам в системный путь:

      export PATH=${PATH}:/home/preps/golovin/progs/bin
      С помощью скрипта Ante_RED.pl подготовлен pdb файл et_out.pdb, p2n файл переименован в Mol_red1.p2n, запущен RED
       RED-vIII.4.pl
      .
      В итоге получен файл Mol_m1-o1.mol2 с координатами атомов и зарядами.
    2. Создание файла описания молекулы в формате пакета программ GROMACS.
      Единица измерения расстояния в GROMACS нанометр.
      Содержание итогового файла et.top
      Первые две строчки:

       
      [ defaults ]
      ; nbfunc        comb-rule       gen-pairs       fudgeLJ fudgeQQ
      1               2               yes              0.5     0.8333
      

      Дальше мы задаём типы атомов и собственно параметры для функции Ленорда-Джонса. Будем считать, что в случае этана Ван-дер-Ваальсовое взаимодействие между атомами углерода разных молекул минимально, так как углероды почти полностью экранированы атомами водорода. Поэтому поставим для углерода некоторые параметры. Ван-дер-Ваальсовый радиус водорода, т.е. сигма нам известен из многих источников, см. webelements.com. Итак у нас получилось, что в этом разделе всего одна переменная это epsilon для водорода.
       
      [ atomtypes ]
      ; name      at.num  mass     charge ptype  sigma      epsilon
      H          1        1.008    0.0000  A   1.06908e-01  1.00000e-00
      C          6        12.01    0.0000  A   3.39967e-01  3.59824e-01
      

      Дальше переходим непосредственно к описанию молекулы. Здесь мы описываем имя и указываем, что соседи через три связи не учитываются при расчете Ван-дер-Ваальсовых взаимодействий. Это верно так, как мы включаем это взаимодействие в торсионные углы:
       
      [ moleculetype ]
      ; Name            nrexcl
      et            3
      

      Добавим атомы этана. В первом столбце идёт номер атома. На него мы будем ссылаться при описании связей. Остальные столбца самоочевидные.
       
      [ atoms ]
      ;   nr  type  resnr  residue  atom   cgnr     charge       mass
           1   C      1    ETH      C1      1     -0.0189        12.01
           2   C      1    ETH      C2      2     -0.0155        12.01
           3   H      1    ETH      H1      3     0.0059        1.008
           4   H      1    ETH      H2      4     0.0059        1.008
           5   H      1    ETH      H3      5     0.0059        1.008
           6   H      1    ETH      H4      6     0.0056        1.008
           7   H      1    ETH      H5      7     0.0056        1.008
           8   H      1    ETH      H6      8     0.0056        1.008
      
      Переходим к описанию связей. Константу жесткости и длина связи С-С взяты из занятия 4.
       
      [ bonds ]
      [ bonds ]
      ;  ai    aj funct  b0       kb
           1   2   1  0.1525   147975.2
           1   3   1  0.1085   300000.0
           1   4   1  0.1085   300000.0
           1   5   1  0.1085   300000.0          
           2   6   1  0.1085   300000.0     
           2   7   1  0.1085   300000.0
           2   8   1  0.1085   300000.0
      
      Описание углов.
       
      [ angles ]
      ;  ai    aj    ak funct  phi0   kphi
      ;around c1
          3     1     4     1  109.500    200.400
          3     1     5     1  109.500    200.400
          4     1     5     1  109.500    200.400
          2     1     3     1  109.500    400.400
          2     1     4     1  109.500    400.400
          2     1     5     1  109.500    400.400
      
      ;around c2
          1     2     6     1   109.500    400.400
          1     2     7     1   109.500    400.400
          1     2     8     1   109.500    400.400
          6     2     7     1   109.500    200.400
          6     2     8     1   109.500    200.400
          7     2     8     1   109.500    200.400
      
      
      Торсионные углы:
       
      [ dihedrals ]
      ;  ai    aj    ak    al funct  t0           kt      mult
          3    1     2     6      1  0.0      0.62760     3
          3    1     2     7      1  0.0      0.62760     3
          3    1     2     8      1  0.0      0.62760     3
          4    1     2     6      1  0.0      0.62760     3
          4    1     2     7      1  0.0      0.62760     3
          4    1     2     8      1  0.0      0.62760     3
          5    1     2     6      1  0.0      0.62760     3
          5    1     2     7      1  0.0      0.62760     3
          5    1     2     8      1  0.0      0.62760     3
      
      Список пар атомов которые не должны считаться при расчете VdW:
       
      [ pairs ]
      ;  ai    aj funct
         3  6
         3  7
         3  8
         4  6
         4  7
         4  8
         5  6
         5  7
         5  8
      
      Описание системы. Есть уже готовые координаты с 38 молекулами этана. Укажем это в описании:
       
      [ System ]
      ; any text here
      first one
      [ molecules ]
      ;Name count
       et    38
      
    3. Моделирование испарения этана.
      Уже подготовлены два состояния системы, первое соответствует газовой фазе, где расстояния между молекулами порядка 50 ангстрем; вторая система имеет такую же плотность как и жидкий этан.
      Наша задача провести короткое моделирование динамики каждой из этих систем о определить разницу в энергии VdW взаимодействий между системами. И сравнить эту разнице с энтальпией испарения этана. При Т=25 это значение равно 5.4 кДж/моль. Вспомним, что epsilon для водорода нам не известна. По аналогии с занятием 4 создадим 7 топологий с разными значениями epsilon.
      проведем для каждой системы молекулярную динамику с каждым файлом топологии (используем md.mdp - файл с настройкам для динамики). Траекторию trr конвертируем в pdb:

       
      trjconv_d -f v_1 -s v_1 -o v_1.pdb
      
      Во всех полученных файлах, кроме v_1, молекулы сразу же разлетаются.В v_1 куб из молекул частично сохраняет свою структуру.
      Посчитаем сами значения энергий, для этого воспользуемся утилитой g_energy.
      Все вышеперечисленные действия выполняются в скрипте et.bash.
    4. На основе полученных txt файлов установлены средния значения энергии для каждого значения epsilon водорода: жидкость, газ Для жидкости и газа разница средней энергии VdW в среднем 4 порядка, поэтому энтальпия испарения примерно равна энергиии жидкости. Энергия кулоновских взаимодействий мала по сравнению с VdW. Epsilon водорода должна лежать в диапозоне 0.03703-0.01562 , чтобы воспроизводилась энтальпия испарения этана.



© Smirnova Victoriya, 2011