For the multi-center PAD calculations in this section, Dyson orbitals were computed using both the EOM-IP-CCSD and the EOM-EA/SF-CCSD methods described above, but with the standard aug-cc-pVTZ basis set (without augmenting it with additional highly diffuse functions). This is to avoid artifacts from the highly diffuse functions that result in electron density far from the atom centers [e.g., see the X 3   EOM-EA/SF Dyson orbital in Fig. ]. β(ε) was computed using both plane (ZS, ZO = 0) and Coulomb (ZS, ZO  0) waves to test the effect of varying the effective charges on oxygen (ZO) and sulfur (ZS), under the constraint ZO = -ZS. Since oxygen is more electronegative than sulfur, we use a negative charge on O and a positive charge on S. The results for ZS = 0, 0.05, 0.10, 0.15, and 0.20, obtained using the X 3   channel Dyson orbital (Fig. The dipole moment of SO at its X 3   equilibrium bond length is 0.610(8) a.u. (1.55 Debye). Dividing by the equilibrium bond length of SO  , 1.592 Å (from CCSD/aug-cc-pVTZ calculations), the partial charge on sulfur is estimated to be ZS = 0.20 (assuming that the dipole moment does not change from its adiabatic value). A similar estimate for the a 1  channel (D = 0.52(2) a.u.) yields ZS = 0.17. These estimates are in reasonably close to the values suggested by the above multicenter calculations. Due to the approximate treatment of the photoelectron wave function and the sensitivity of the multi-center approach to the details of the method and basis set, the better quantitative agreement for the EOM-EA/SF-CCSD calculations is likely coincidental and does not mean that the EOM-EA/SF-CCSD method is more accurate than EOM-IP-CCSD for this system.