Помогите рассчитать энергию Гиббса
Помогите рассчитать энергию Гиббса
Здравствуйте. У меня 2 реакции, которые могут конкурировать между собой. Хочу узнать, какой процесс термодинамически более вероятен. Для этого, насколько я знаю нужно рассчитать энергию Гиббса. Что я делаю: рассчитываю гессиан оптимизированных молекул, участвующих в реакции. Вычисляю тепловые эффекты как разницу продуктов + колебательные поправки и поправки к энтальпии (H total) и реагентов. Также вычисляю энтропию реакции. Затем по формуле G= H - TS вычиляю (температура в кельвинах, энтропия и энтальпия в кдж/моль) энергию Гиббса и получаю, то процесс запрещен, хотя он легко идет идет на практике.
Скажите, что я упустил и что значат поправки к энергии Гиббса (G total), и к E total которые Gamess выдает при расчете гессиана, может их тоже надо как-то использовать?
И еще: я во многих работах видел, что что для изучения того, какой процесс термодинамически вероятен более, люди не считают энергию Гиббса, а довольствуются тепловым эффектом. Неужели, этого достаточно?
Скажите, что я упустил и что значат поправки к энергии Гиббса (G total), и к E total которые Gamess выдает при расчете гессиана, может их тоже надо как-то использовать?
И еще: я во многих работах видел, что что для изучения того, какой процесс термодинамически вероятен более, люди не считают энергию Гиббса, а довольствуются тепловым эффектом. Неужели, этого достаточно?
Re: Помогите рассчитать энергию Гиббса
H total и G total после расчета гессиана - это лишь термодинамические поправки к TOTAL ENERGY, которая появляется на этапе сведения SCF. Т.е. нужно вычислить (TOTAL ENERGY)*627.5 + (G total) для продуков (627.5 - коэффициент перевода атомных единиц энергии в ккал), вычесть то же самое для реагентов.
Действительно, часто G total пренебрегают, используют только TOTAL ENERGY. Одна из причин - не нужно вычислять гессиан (это долго). Но главное, почему это оправдано - как правило ошибка TOTAL ENERGY в "обычном", рутинном расчете существенно превосходит термодинамические поправки. Даже для реакции в газовой фазе, чтобы расчет можно было использовать для термохимии, нужно применять специально подобранные высокоуровневые методы и базисы. К тому же, для родственных соединений термодинамические поправки бывают близки и при вычитании взаимно сокращаются.
Действительно, часто G total пренебрегают, используют только TOTAL ENERGY. Одна из причин - не нужно вычислять гессиан (это долго). Но главное, почему это оправдано - как правило ошибка TOTAL ENERGY в "обычном", рутинном расчете существенно превосходит термодинамические поправки. Даже для реакции в газовой фазе, чтобы расчет можно было использовать для термохимии, нужно применять специально подобранные высокоуровневые методы и базисы. К тому же, для родственных соединений термодинамические поправки бывают близки и при вычитании взаимно сокращаются.
Re: Помогите рассчитать энергию Гиббса
А что у вас за молекулы, если не секрет?
А я вот паровоз поднимал... Но не поднял.
Re: Помогите рассчитать энергию Гиббса
А тогда у меня еще такая проблема: если сравнить Etotal в хартри одного продукта и другого, то между ними малая разница, но энергия одной молекулы все же ниже, чем другой. При этом тепловой эффект реакции мне упорно показывает то, что термодинамически более выгодно образование менее устойчивой молекулы? это вообще нормально?
Для термохимии я использую метод функционала плотности с b3lyp, базис 6-31 d,p Gamess.
тепловой эффект реакции рассчитываю вот так:H(T) = E + ZPE + H(total) вроде все верно, а получается не совсем то, что нужно...может еще, что посоветуете для решения данной проблемы?
Я рассчитываю органику, различные гетероциклы.
Для термохимии я использую метод функционала плотности с b3lyp, базис 6-31 d,p Gamess.
тепловой эффект реакции рассчитываю вот так:H(T) = E + ZPE + H(total) вроде все верно, а получается не совсем то, что нужно...может еще, что посоветуете для решения данной проблемы?
Я рассчитываю органику, различные гетероциклы.
Re: Помогите рассчитать энергию Гиббса
Etotal которая в той же табличке, что и H total, G total? Тогда это не та величина. Нужно TOTAL ENERGY, которая ДО гессиана.
Re: Помогите рассчитать энергию Гиббса
Я использую TOTAL ENERGY взятую из файла оптимизации геометрииamge писал(а):Etotal которая в той же табличке, что и H total, G total? Тогда это не та величина. Нужно TOTAL ENERGY, которая ДО гессиана.
Re: Помогите рассчитать энергию Гиббса
Не правильно. ZPE уже включена в H(total). Т.е. H(T) = E + H(total), где E = TOTAL ENERGY умноженное на 627.5Nicolass писал(а):тепловой эффект реакции рассчитываю вот так:H(T) = E + ZPE + H(total) вроде все верно
Еще мне кажется, что при использованном Вами методе расчета нельзя требовать от него термохимической точности. Ошибка при расчете TOTAL ENERGY в пару килокалорий представляется не слишком большой. Но для термохомии пара килокалорий - это огромная величина. А если к тому же экспериментальные данные для раствора, а расчет для газовой фазы...
Re: Помогите рассчитать энергию Гиббса
Спасибо. А какой метод и базис вы бы посоветовали взять для моего случая?amge писал(а): Т.е. H(T) = E + H(total), где E = TOTAL ENERGY умноженное на 627.5
Re: Помогите рассчитать энергию Гиббса
Обратите внимание вот на это:
Это важно. Для раствора можно применять модель PCM (polarizable continuum model), причем там сразу выдается Free energy, и не надо никаких гессианов. Посмотрите по соседним темам, тут уже обсуждали кое-какие вопросы и проблемы, связанные с PCM расчетами.amge писал(а):А если к тому же экспериментальные данные для раствора, а расчет для газовой фазы...
Вот и вся моя работа. Стеречь ребят над пропастью во ржи. (Дж. Д. Сэлинджер)
Re: Помогите рассчитать энергию Гиббса
Если имеется в виду термохомия, то ничего не посоветую, ибо не специалист здесь. Кроме и так очевидного - посмотреть литературу или самому попробовать подобрать метод/базис на представительном массиве родственных реакций, для которых есть экспериментальные данные.Nicolass писал(а):А какой метод и базис вы бы посоветовали взять для моего случая?
Re: Помогите рассчитать энергию Гиббса
про PCM я в курсе. Только вот (простите за глупость и не знание ) что за Free energy и как мне ее применить для того, чтобы узнать какой процесс вероятнее? сумма Free energy продуктов - сумма Free energyреагентов что ли?sanya1024 писал(а):Обратите внимание вот на это:Это важно. Для раствора можно применять модель PCM (polarizable continuum model), причем там сразу выдается Free energy, и не надо никаких гессианов. Посмотрите по соседним темам, тут уже обсуждали кое-какие вопросы и проблемы, связанные с PCM расчетами.amge писал(а):А если к тому же экспериментальные данные для раствора, а расчет для газовой фазы...
Re: Помогите рассчитать энергию Гиббса
А гетероциклы структурно жесткие, или там возможны конформации, вращения? Это существенно может повлиять на оценку энтропии.Nicolass писал(а):Я рассчитываю органику, различные гетероциклы.
Free energy = свободная энергия = энергия Гиббса
А изменение энтропии реакции в обоих случаях одинаково?Nicolass писал(а):При этом тепловой эффект реакции мне упорно показывает то, что термодинамически более выгодно образование менее устойчивой молекулы?
Такое бывает, особо если у вас эксперимент в растворе, а расчет для газа. И от выбранного метода может многое зависеть. Попробуйте разные базисы - и увидите как меняется картина.Nicolass писал(а):это вообще нормально?
А я вот паровоз поднимал... Но не поднял.
Re: Помогите рассчитать энергию Гиббса
Nicolass
b3lyp в топку: попробуйте M062X и если железо позволяет, то базис поприличней возьмите.
b3lyp в топку: попробуйте M062X и если железо позволяет, то базис поприличней возьмите.
прозвище "Фабержé" легендарный разведчик Дроздов получил за свое уникальное умение работать с информацией, добывать ее и превращать в драгоценность высшей пробы.
Re: Помогите рассчитать энергию Гиббса
Довольно станное сочетание слов - термохимия и b3lyp/6-31G(d,p). Нужно хотя бы B3LYP/6-311G(2d,2p).Nicolass писал(а):Для термохимии я использую метод функционала плотности с b3lyp, базис 6-31 d,p Gamess.
И там дальше думать, что включать сольватацию или, может, фукционалы разные пробовать. Может, и другие методы попробовать следует.
Последний раз редактировалось o-oxhem Чт окт 20, 2011 10:09 pm, всего редактировалось 1 раз.
Re: Помогите рассчитать энергию Гиббса
Именно. Только это уже не энергия Гиббса будет, а скорее Гельмгольца (т.е., при постоянном объеме -- у нас же раствор).Nicolass писал(а): про PCM я в курсе. Только вот (простите за глупость и не знание ) что за Free energy и как мне ее применить для того, чтобы узнать какой процесс вероятнее? сумма Free energy продуктов - сумма Free energyреагентов что ли?
Вот и вся моя работа. Стеречь ребят над пропастью во ржи. (Дж. Д. Сэлинджер)
Re: Помогите рассчитать энергию Гиббса
А энергию Гиббса можно и в Гельмгольца превратить, если нужный коэффициент прибавить
Или константу на нужное умножить.
Или константу на нужное умножить.
Re: Помогите рассчитать энергию Гиббса
А зачем? Если константа равновесия при постоянном объеме (в растворе), по-любому для нее берется свободная энергия F (Гельмгольца). Это только для газов константа считается при постоянном давлении, т.е. с энергией Гиббса.
Вот и вся моя работа. Стеречь ребят над пропастью во ржи. (Дж. Д. Сэлинджер)
Re: Помогите рассчитать энергию Гиббса
Всем спасибо за помощь) С тепловым эффектом теперь все нормально, а вот по энергии Гиббса процесс запрещен((((
Возможно, просто температуру низкую взял....
Возможно, просто температуру низкую взял....
Re: Помогите рассчитать энергию Гиббса
А моим советом воспользовались?Nicolass писал(а):Всем спасибо за помощь) С тепловым эффектом теперь все нормально, а вот по энергии Гиббса процесс запрещен((((
Возможно, просто температуру низкую взял....
прозвище "Фабержé" легендарный разведчик Дроздов получил за свое уникальное умение работать с информацией, добывать ее и превращать в драгоценность высшей пробы.
Re: Помогите рассчитать энергию Гиббса
Пока нет еще, но обязательно попробую, как только свои расчеты до ума доведу. А что, вообще хорошая штука, помогает, проверено?
Кто сейчас на конференции
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 15 гостей