Articles you may be interested inContinuum solvation models are widely used to accurately estimate solvent effects on energy, structural and spectroscopic properties of complex molecular systems. The polarizable continuum model ͑PCM͒ is one of the most versatile among the continuum models because of the variety of properties that can be computed and the diversity of methods that can be used to describe the solute from molecular mechanics ͑MM͒ to sophisticated quantum mechanical ͑QM͒ post-self-consistent field methods or even hybrid QM/MM methods. In this contribution, we present a new formulation of PCM in terms of a free energy functional whose variational parameters include the continuum polarization ͑represented by the apparent surface charges͒, the solute's atomic coordinates andpossibly-its electronic density. The problem of finding the optimized geometry of the ͑polarized͒ solute, with the corresponding self-consistent reaction field, is recast as the minimization of this free energy functional, simultaneously with respect to all its variables. The numerous potential applications of this variational formulation of PCM are discussed, including simultaneous optimization of solute's geometry and polarization charges and extended Lagrangian dynamics. In particular, we describe in details the simultaneous optimization procedure and we include several numerical examples.