3. Density operator