The stationary Schrödinger equation can be cast in the form Hρ = Eρ, where H is the system's Hamiltonian and ρ is the system's density matrix. We explore the merits of this form of the stationary Schrödinger equation, which we refer to as SSEρ, applied to many-body systems with symmetries. For a nondegenerate energy level, the solution ρ of the SSEρ is merely a projection on the corresponding eigenvector. However, in the case of degeneracy ρ in nonunique and not necessarily pure. In fact, it can be an arbitrary mixture of the degenerate pure eigenstates. Importantly, ρ can always be chosen to respect all symmetries of the Hamiltonian, even if each pure eigenstate in the corresponding degenerate multiplet spontaneously breaks the symmetries. This and other features of the solutions of the SSEρ can prove helpful by easing the notations and providing an unobscured insight into the structure of the eigenstates. We work out the SSEρ for the system of spins 1/2 with Heisenberg interactions. Eigenvalue problem for quantum observables other than Hamiltonian can also be formulated in terms of density matrices. We provide an analytical solution to one of them, S 2 ρ = S(S + 1)ρ, where S is the total spin of N spins 1/2, and ρ is chosen to be invariant under permutations of spins. This way we find an explicit form of projections to the invariant subspaces of S 2 . Finally, we note that the anti-Hermitian part of the SSEρ can be used to construct sum rules for temperature correlation functions, and provide an example of such sum rule.
03J 12 3J 23 3J 13 J 12 −2J 12 J 23 J 13 J 23 J 12 −2J 23 J 13 J 13 J 12 J 23 −2J 13