Firefly and PC GAMESS-related discussion club

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

Re: Calculation of complexation constant in water using firefly

Jim Kress

Is there any particular reason you are not doing this calculation with a Molecular Dynamics (e.g. gromacs) or Quantum Molecular Dynamics (e.g. Cp2K) code?

I would think they would be better suited to your needs.


On Mon Jan 19 '15 2:12pm, Enrico wrote
>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


[ Previous ] [ Next ] [ Index ]           Mon Feb 9 '15 7:57am
[ Reply ] [ Edit ] [ Delete ]           This message read 561 times