Abstract Calculating the wave function and the properties of a molecule in solution is a common procedure today that can be performed with a variety of methods. However, the study of solvated systems did not begin with the establishment of the Schrodinger equation in 1926 and the development of the first methods in computational chemistry. In fact, it took until 1973 to see the first self-consistent field equations that included the electrostatic potential of the solvent. It took another two decades for computational methods and solvent models to become consolidated, with an explosion of applications in the 1990s. This article briefly describes that fascinating story.