## Firefly and PC GAMESS-related discussion club

Learn how to ask questions correctly

** Re^2: Calculation of complexation constant in water using firefly **
Enrico

e.chiavazza@gmail.com

Thx very much for your suggestions Sanya, i knew that the basis set i mentioned where very small, but i have have constrains on the computation power so i chose such low level to make some initial screening...anyway i'm performing my task with B3LYP/6-31+G(p) for geometry optimization and B3LYP/6-311+G(2p,d) for single point dpcm and all the molecules have explicit 4-6 molecules of water pre-optimized by molecular dynamics.

I'll let you know if i'll end up with any useful information for my work!

have a nice day

EC

On Sat Jan 31 '15 7:44pm, sanya wrote

-------------------------------------

>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, (H_{2}O)_{q} is a cluster of q water molecules):

>[Na(H_{2}O)_{n}]^{+} + [L(H_{2}O)_{m}]^{-} = NaL(H_{2}O)_{p} + (H_{2}O)_{q}

Mon Feb 9 '15 2:36pm

This message read **423** times