Firefly and PC GAMESS-related discussion club

Learn how to ask questions correctly  
We are NATO-free zone

Calculation of complexation constant in water using firefly


Hello to everyone,
I'm trying to estimate a trend in the complexation constant of sodium ion with several polycarboxylic ligand in water, i know that i should calculate the dG value of the complex from the G value of the reagents, i read several paper on which they explain the need of corrections that account for the solvatation contribution and for the thermal factor (in order to bring the system from 0 K to 298 K and 1 bar). In order to do that I've found the optimized structure of the participant molecule under PCM solvent modeling and than i run an hessian calculation to obtain the thermochemistry at RT but if i use the gibbs energy that i obtained from the thermochemistry i obtain a positive dG for the complex that do not reflect at all the experimental behavior so now i have some question to which i couldn't find an answer elsewhere:

- what is the correct energy value that i have to consider from the firefly output?
- what are the correct correction factor (Ezpe,Hcorr,Scorr,solvent) that i need to get a formally correct dG for the complexation reaction?

solved these basic issues, in your opinion is it better to take into account explicitly the first solvatation sphere of the sodium ion (6 water molecules)?

you can find attached a set of output file for the smaller complex of my serie,
thank you


This message contains the 1298 kb attachment
[ ] Output files for complexation reaction

[ Previous ] [ Next ] [ Index ]           Mon Jan 19 '15 2:12pm
[ Reply ] [ Edit ] [ Delete ]           This message read 731 times