We anticipate that the present approach will enhance the accuracy of simulations of systems where molecules move between regions with very different dielectric constants, such as water and a lipid membrane, or with very different ambient electric fields, such as from water to a protein pocket with many ionized side-chains. These are settings in which large changes in electronic polarization occur, making it important to account for them in detail.