For the direct analysis of the induced dipoles and the associated electrostatic potentials, fully self-consistent ground-state optimizations were performed for a randomly chosen single configuration of each system using models A-D with varying cutoff distances for models A and B. The electrostatic potentials from the permanent multipoles and induced dipoles were then calculated at the center of mass and on the van der Waals surface of the quantum subsystem for each configuration, model, and cutoff distance, as were the lowest-lying singlet excitation energies, oscillator strengths, and TPA cross-sections. Electrostatic potentials (at the center of mass), excitation energies, oscillator strengths, and TPA cross-sections were also calculated for all of the configurations of each system using models C and D. All calculations were performed using the Dalton program and PElib. The embedding potentials were created using PyFraME.