In quantum mechanics, an open system interacting with the environment is described by mixed state, which is represented by the density operator ρ(t) acting on a finite dimensional Hilbert space Hn. By the spectral theorem, ρ(t) can be decomposed as ρ(t)=∑nk=1ωk(t)|ϕk(t)⟩⟨ϕk(t)|, where ωk(t) are the eigenvalues, |ϕk(t)⟩ are the eigenvectors and ⟨ϕk(t)| are the corresponding vectors in the dual space of Hn. For simplicity, here only non-degenerate case is considered. I am now trying to assign a geometric phase for the mixed state represented by ρ(t)=∑nk=1ωk(t)|ϕk(t)⟩⟨ϕk(t)|. The geometric phase is defined as n∑k=1∫t0ωk(t′)dγk(t′)dt′dt′, where γk(t)=arg⟨ϕk(0)|ϕk(t)⟩+i∫t0⟨ϕk(t′)|˙ϕk(t′)⟩dt′ are the geometric phases for eigenvectors |ϕk(t)⟩. Here we assume that |ϕk(t)⟩ are not orthogonal to |ϕk(0)⟩ for any time t and γk(t) are smooth functions of time t. Note that γk(t) are put in differentiation of time t. There is no ambiguity caused by multiple values of γk(t). Thus the above geometric phase is well-defined. Question: How to test in experiment the geometric phase for mixed states defined above?