Помогите рассчитать энергию Гиббса

[Yurii]

Пока нет еще, но обязательно попробую, как только свои расчеты до ума доведу. А что, вообще хорошая штука, помогает, проверено?

Проверено: B3LYP бред выдавал даже в хорошем базисе. А M062X кореллировал с MP2. Забудьте о B3LYP - это уже маргинальщиной попахивает.

[Nicolass]

про PCM я в курсе. Только вот (простите за глупость и не знание ) что за Free energy и как мне ее применить для того, чтобы узнать какой процесс вероятнее? сумма Free energy продуктов - сумма Free energyреагентов что ли?

Именно. Только это уже не энергия Гиббса будет, а скорее Гельмгольца (т.е., при постоянном объеме -- у нас же раствор).

а не подскажите, могу ли я сравнить энергии Гельмгольца двух продуктов реакции и выявить, какой процесс термодинамически более вероятен? А то мне долго считать реагенты с учетом растворителя, а результат нужен поскорее. Может это сгодится для начального приближения или, все же нет?

[sanya1024]

Если реагенты одни и те же, то достаточно сравнить энергии продуктов. А Вы точно уверены, что реакция контролируется термодинамически, а не кинетически?

[Nicolass]

Если реагенты одни и те же, то достаточно сравнить энергии продуктов. А Вы точно уверены, что реакция контролируется термодинамически, а не кинетически?

Я это подозреваю. Просто вычислю разность между тепловыми эффектами, при различных температурах. если разность не равна константе, то все хорошо. Кинетический анализ хотелось бы сделать, но сами понимаете, там переходное состояние) вообще, 87% на то, что контроль термодинамический.

[sanya1024]

Просто если там здоровенный барьер, то никакого толку от разности тепловых эффектов (или даже свободных энергий) не будет, все будет лимитироваться барьером. И тогда от поисков переходного состояния никуда Вам не деться.

[Nicolass]

[
А моим советом воспользовались?

Yurii, подскажите, а этот функционал в Firefly есть или он только в Gamess Us ? в FF я его не нашел(

[avelon]

gamess us, nwchem

[Yurii]

А моим советом воспользовались?

Yurii, подскажите, а этот функционал в Firefly есть или он только в Gamess Us ? в FF я его не нашел(

Посчитайте предварительно в g09, и если все будет тип-топ, для очистки совести пересчитайте в Gamess Us. Учтите, что g09 под Linux x64 считает раза в 2 быстрей, чем под виндами, памяти можно немерянно использовать, что для Вашей задачи актуально, т.к. для нее необходимо использовать хороший или лучше очень хороший базис, да и под Linux версия посвежей (B01).

[asp]

У меня тоже проблема с расчетом изменения энергии Гиббса и энтальпии реакции. Если я все правильно понял, то считать надо так. Кусок out файла из Природы:
E = -325.854161 E0 = -325.611863

ZPVE = 0.242297 a.u. = 152.0440 kcal/mol

----------------------------------------------------------------------
T = 298.150 K___ln(Q)______S________U________H ______G
______________________cal/mol*K__kcal/mol__kcal/mol __kcal/mol
----------------------------------------------------------------------
translational____18.1846__41.1049 __0.8887___1.4812___-10.7742
rotational _____14.7975 __32.3867___0.8887 __ 0.8887__-8.7674
vibrational _____ 19.7807 __83.6468__165.2634 __165.2634 __140.3241
total___________52.7629__157.1384 _167.0409__167.6334__120.7826
----------------------------------------------------------------------

ΔG=Gпродукты-Gреагенты
G=H-TS
S=Stotal
H=E+ ZPVE+Htotal
или
H=E0+Htotal
Энтальпия реакции ΔH=Hпродукты-Hреагенты
Вопрос, что делать с Gtotal?
Amge, почему Gtotal суммируется с E?

[EvgeniX]

Может быть вам программу для ПРИРОДы выслать для самопроверки? Как раз считает термодинамику по out-Файлам, только сырая пре-альфа версия.

ΔG=Gпродукты-Gреагенты

mA+nB => xC+yD
ΔG = x*G(C) + y*G(D) - m*G(A) - n*G(B)

G=H-TS

В принципе, для обработки данных эту формулу можно и не использовать.

S=Stotal

Правильно.

H=E+ ZPVE+Htotal
или
H=E0+Htotal

Не правильна ни первая, ни вторая формула. Htotal уже включает ZPVE.
Т.е. H=E+Htotal

Вопрос, что делать с Gtotal?
Amge, почему Gtotal суммируется с E?

Потому, что практически во всех термодинамических функциях явно или неявно присутствует внутренняя энергия (U0 = E +ZPE).
G=E+Gtotal

[asp]

Всё, разобрался, спасибо.