Home > Community > Thermochemistry calculations by using Gaussian?
Upvote

28

Downvote
+ Gaussian
+ Temperature
+ Thermochemistry
+ Pressure
Posted by
Maqsood Ahmed

Thermochemistry calculations by using Gaussian?

Danilo Renzi  Follow

Yes, it's possible.

More

Upvote

VOTE

Downvote
Andreas Nikolaou  Follow

Thank you Bojidarka B. Ivanova and Debasish Mandal

More

Upvote

VOTE

Downvote
Cross-dressing  Follow

Dear Syamlal Sk ,
You did not specify exactly what you want to calculate. However, you should be warned about quality of thermochemistry calculated by Gaussian.
1) In the way described above, the frequencies will be calculated as harmonic and using them in statistical thermodynamics calculations without proper scaling will introduce significant offsets.
2) Gaussian only uses RRHO model (or 1DHR correction when employing a right keyword) which can be unsuitable for more complex molecules than water.
3) Since the calculations are performed for ideal gas, I'm not sure about effect of pressure you want to evaluate. That one will be only simple dependence in entropic term, so you can evaluate it in excel...


More

Upvote

VOTE

Downvote
Dave Hunter  Follow

Mr. Debasish Mandal ,

Here i couldn't find any keywords works for temperature and pressure calculations?
Can you suggest any keywords?

More

Upvote

VOTE

Downvote
António José Rodrigues Rebelo  Follow
Upvote

VOTE

Downvote
David Bangley  Follow

This works. Perform the freq calculation at default temperature and pressure. KEEP the checkpoint file! Then recalculate at the desired temperature and pressure with the freq=(readfc,readiso) command.

For example, if you have calculated:

%chk=water.chk
# hf/6-31g** opt freq
water

0 1
H
O 1 0.99
H 2 0.99 1 105.


Then you get thermochemistry at 298.15 K and 1.0 atm. To recalculate at 500 K and 2.0 atm, submit the following commands:

%chk=water.chk
#hf/6-31g** geom=allcheck freq=(readfc,readiso)
water recalc'd

500.0 2.0
1
16
1

The last few lines are the isotope numbers you want to use--Gaussian converts these into exact masses.






More

Upvote

VOTE

Downvote
Claudio Giomini  Follow

following

More

Upvote

VOTE

Downvote
Bobby Baucom  Follow

This works. Perform the freq calculation at default temperature and pressure. KEEP the checkpoint file! Then recalculate at the desired temperature and pressure with the freq=(readfc,readiso) command.

For example, if you have calculated:

%chk=water.chk
# hf/6-31g** opt freq
water

0 1
H
O 1 0.99
H 2 0.99 1 105.


Then you get thermochemistry at 298.15 K and 1.0 atm. To recalculate at 500 K and 2.0 atm, submit the following commands:

%chk=water.chk
#hf/6-31g** geom=allcheck freq=(readfc,readiso)
water recalc'd

500.0 2.0
1
16
1

The last few lines are the isotope numbers you want to use--Gaussian converts these into exact masses.






More

Upvote

VOTE

Downvote
Cyndi Root  Follow

Yes it is possible. Please go through the following link
http://gaussian.com/freq/


More

Upvote

VOTE

Downvote