There is no option in G09 to get the solvation energy directly from output. Run a gas phase calculation separately and see the difference.  

There is no option in G09 to get the solvation energy directly from output. Run a gas phase calculation separately and see the difference.  

Mr. Ilic,

Or by G98W.

Try this: SCRF=(PCM,SMD,SC,DoVacuum,Solvent=Acetonitrile)

It works for me

Just run two calculations.. One in solvent/PCM and another in gas phase.. Then take the difference in energy. This is what I usually do. 

I think M. Manathunga is right.

Let us focus attention to an example of computations of a single H2O molecule in water using "direct-output" computations and the manual approach recommended by Mr. Manathunga:

Example: H2O (HF/STO-3G)

"Direct-output" computations:

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

Self-consistent PCM results

===========================

(a.u.) = -74.960089

(a.u.) = -74.964414

(a.u.) = -74.964568

(a.u.) = -74.959929

(a.u.) = -74.964569

(Unpolarized solute)-Solvent (kcal/mol) = -2.71

(Polarized solute)-Solvent (kcal/mol) = -2.91

Solute polarization (kcal/mol) = 0.10

Total electrostatic (kcal/mol) = -2.81

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

DeltaG (solv) (kcal/mol) = -2.81

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

And the 'manual approach':

E-GP = -74.9600885 a.u.

E-PCM-H2O = -74.9645688 a.u.

D = -0.0044803

Dx627.51 = -0.0044803x625.51 = -2.811433053 kcal.mol-1

Thank you so much for your help and tips, everyone!