Because of a non-Hermitian nature of EOM-CC, γ pq (I, I ) = γ qp (I , I), which leads to different numerical results for Ψ I |A|Ψ I and Ψ I |A|Ψ I A common resolution of this problem is to take a geometric average . However, this treatment should be modified for complex-valued tensor properties, in order to preserve phase consistency and to generalize it for a matrix. Below (section II C) we introduce a spinless density matrix, which enables generation of all phase-consistent SOCs between the selected pairs of states. We also propose a different resolution of the averaging problem.