Abstract
Starting from the screening in conductors, an algorithm for the accurate calculation of dielectric screening effects in solvents is presented, which leads to rather simple explicit expressions for the screening energy and its analytic gradient with respect to the solute coordinates. Thus geometry optimization of a solute within a realistic dielectric continuum model becomes practicable for the first time. The algorithm is suited for molecular mechanics as well as for any molecular orbital algorithm. The implementation into MOPAC and some example applications are reported.