On the other hand, the models that include the SO dipole into the theoretical framework yield a much-improved agreement with the experiment. These The computed photoelectron angular distributions are sensitive to the long-range behavior of wave functions and, therefore, the exact method for calculating the Dyson orbitals. At a relatively low-level treatments (up to double excitations), the EOM-IP approach using the anionic reference yields more diffuse Dyson orbitals and quantitatively different () curves compared to the EOM-EA/SF method using the neutral reference. For most of the calculations in this work, the EOM-IP Dyson orbital calculations yield the more accurate () curves. This is consistent with a known behavior of these methods in calculations using closed-shell references, e.g., the effect of triple excitations is more important for EOM-EA than for EOM-IP. Hence, we expect the differences between the two approaches to be reduced when triple excitations are included.