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

sanya
sanya@photonics.ru


Well, in my view, everything is wrong here.

First of all, 3-21G basis set is too small for thermochemistry. Triple or, better, quadruple-zeta basis sets augmented with diffuse and polarization functions are commonly used.

Second, RHF is also inappropriate for thermochemistry. You should use correlated methods (DFT or MP2 geometry optimization with single-point MP4 correction).

Third, and most important is the solvation model. PCM cannot account for specific solvation of the cation and functional groups of the anion and the parent ligand. You should explicitly solvate the cation, anion, and complex with sufficient amount of water molecules (5-6 for Na+, at least 1-2 molecules for each O atom of the ligand and complex, 1 water molecule for N atom) and embed this cluster in PCM. Obviously, you should take into account hydration/dehydration of Na+ during complexation (excess water molecules should also be embedded in PCM; probably, they should be combined in a cluster, because water always forms hydrogen bonds with its neighbors).

Finally, vibrational calculations (gas phase, harmonic oscillator, rigid rotor) are irrelevant for solution thermochemistry. PCM calculations directly give free energy of the solute in solution (TOTAL FREE ENERGY IN SOLVENT). To calculate the complexation free energy, you should consider the free energy of the following equilibrium (L is the ligand, m + n = p + q, (H2O)q is a cluster of q water molecules):

[Na(H2O)n]+ + [L(H2O)m]- = NaL(H2O)p + (H2O)q


[ Previous ] [ Next ] [ Index ]           Sat Jan 31 '15 7:44pm
[ Reply ] [ Edit ] [ Delete ]           This message read 402 times