I am a physics undergraduate. I just read through Sakurai's section 3.4 (3rd ed) on Density Operators and Pure vs. Mixed states.
Are pure states as kets generalized to mixed states as operators in order to retain all the properties we know and love of states as kets while affording the ability to introduce relative population weights of such states that do not interfere with such ket properties (like interference)?
For example, if we tried to stick with the ket state formalism, then we wouldn't be able to distinguish relative population weights from complex amplitudes and thus between a more "classical" probability from the "quantum" probability?
Additionally, would we write a general density operator as
$$\rho = \sum_j w_j\left(\sum_{\textbf{i}_n} c_{\textbf{i}_n}|\textbf{i}_n\rangle \langle \textbf{i}_n |\right)$$
where $| \textbf{i}_n \rangle = |i_1\rangle \otimes |i_2 \rangle \otimes ...|i_n\rangle$ and $\sum_i w_i = 1$?
In other words, would we write an arbitrary general density matrix as a convex combination of arbitrary pure state density operators?