Sub-Geometric Phases in Density Matrices

This study presents the generalization of geometric phases in density matrices. We show that the extended sub-geometric phase has an unified expression during the adiabatic or nonadiabatic process and establish the relations between them and the usual Berry or Aharonov-Anandan phases. We also demonstrate the influence of sub-geometric phases on the physical observables. Finally, the above treatment is used to investigate the geometric phase in a mixed state.

with instantaneous eigenstate ψ R t ( ( )) n and eigen energy E R t ( ( )) n , and the Aharonov-Anandan (AA) phase 7 in the nonadiabatic evolution of wave function with Hamiltonian H R t ( ( )), are canceled by their conjugate and therefore do not contribute to their density matrices. This motivates us to find a geometric phase that can appear in the density matrix and allows us to address the effects of geometric phase on the physical observables.
The geometric phase in the wave function is usually observed via interference experiments. Does the geometric phase produce other observable effects? Several works had been done toward this goal. In 1983, Thouless et al. demonstrated the influence of geometric phase on the adiabatic particle transport process 8 . They also elaborated the quantized Hall conductance via geometric phase 9 . Based on the path integral formalism, Huratsuji et al. suggested that the geometric phase may affect the quantization rule 5 and hence the quantized energy. To illuminate the influence of geometric phase on the acoustoelectric current induced by surface acoustoelectric waves, we developed a position dependent geometric phase, and exploited it to study this quantized acoustoelectric current 10 . Further, can we expect this geometric phase to have effect on other physical observables? To our knowledge, this is an open question and worth to explore. The usual geometric phases do not appear in the density matrices; however, the physical observables are always calculated via the statistical average by the density matrix. Thus, the geometric phase that can exist in the density matrix must be determined.

the Sub-Geometric phases
Initially, we studied the adiabatic evolution of a system H R t ( ( )) with parameter R t ( ), whose eigenequation is To obtain the Berry phase, the instantaneous eigenstate ψ R t ( ( )) n should be calculated. However, it is seem not trivial because the instantaneous eigenequation is difficult to solve except in some special cases, which is discussed below.
If we denote the Hamiltonian at initial time as H R ( (0)) with its eigenequation ψ ψ = H R E ( (0)) n n n 0 0 0 , the time-dependent Hamiltonian can be divided as follows: where ∆H R t ( ( )) describes the difference between the Hamiltonian H R t ( ( )) and its initial value H R ( (0)), which is usually small when the system varies adiabatically. We can now perform the expansion of the instantaneous eigenfunction by the complete basis { } R t c R t e ( ( )) ( ( )) , (2) n k In the case of adiabatic evolution, ∆H R t ( ( )) may be regarded as perturbation, therefore, c nk may take the form of the first-order approximation as follows: c R t e HR t dt It should be ponied out that Eqs (1-3) are still suitable for the description of nonadiabatic case, in which we can not solve Eq. (3) by the perturbative theory again.
Considering Eq. (2), if we introduce a notation such that R t satisfy the Schrödinger equation, f(t) can be solved as follows: constitutes a complete basis. The instantaneous eigenfunction can also expanded as follows: nk , we will not write them out in the following. The phase in the above expression is just the geometric phase associated with the basis | 〉 and is name d as t he sub -ge ome t r i c phas e, w he re as t he phas e is the dynamical phase. Following Eq. (5), its corresponding density matrix can be represented as follows: is just the relative sub-geometric phase between states ψ | 〉 k 0 and l 0 ψ | 〉 in the density matrix. Using Eq. (2), the sub-geometric phase can be connected with the usual Berry phase, which can be expressed as , and the sub-geometric phase is , we find that the Berry phase can be regarded as the summation of sub-geometric phases with the following probability: which relates the Berry phase and the sub-geometric phase proposed herein.
In fact, although the Berry phase will cancel with its conjugate in the expression of density operator when the wavefunction evolves as ψ ≅ ∫ in the adiabatic case, it should be noted that this wavefunction is an approximate form, its exact expression should be a e m R t (R(t)) e ( ( )) ,  . Employing the exact expression (8), one can obtain its density operator as:  (8) and (9) are very similar to Eqs (5) and (6) for the sub-geometric phase, they are based on the instantaneous eigenstate representation, while the sub-geometric phases are based on the initial eigenstate representation, but they are the same in the physical nature. It can be seen that the Berry phase shown in expression (9) is the relative Berry phase between the instantaneous eigenstates R t n ( ( ) ) ′ and m R t ( ( ) ) ′ , they appear only in the off-diagonal elements of the density matrix like the sub-geometric phase.
Compared with the adiabatic evolution, the above sub-geometric phase is suitable for the description of the nonadiabatic case, in which the Hamiltonian H(t) can also be written as H Ht H(t) ( ) 0 = + ∆ , and its wave function may still be expanded similar to Eq. (5) as follows: is the sub-geometric phase. We can now see that the expression has the same form as that in the adiabatic case. It should be emphasized that ∆H t ( ) is not small in the nonadiabatic procedure; therefore, the coefficients c t { ( )} k in the above expression cannot be obtained via perturbative theory as in the adiabatic case. However, they still satisfy Eq. (3), and we can solve them from these coupled differential equations. Beyond this, the sub-geometric phase in the density matrix is also analogous to the adiabatic case (6), in which the relative geometric phase indicates the coherence between states ψ | 〉 k 0 and ψ | 〉 l 0 . Therefore, the sub-geometric phase reported herein has a unified expression during adiabatic and nonadiabatic processes.
Generally, the nonadiabatic AA phase is not easy to calculate owing to the unknown function t ( ) ψ ∼ . However, it can be explicitly expressed by the coefficient c t ( ) k in Eq. (2), where the total phase is The AA phase β can be finally obtained as φ β τ = − α τ ( ) ( ). Therefore, we can provide an alternative way to calculate the AA phase.
Since the sub-geometric phase may appear in the density matrix during adiabatic and nonadiabatic processes, it will affect the physical observables through the statistical average by the density matrix, 〈 〉 = ρA A tr( ) . As stated above, the relative sub-geometric phases in the density matrix will not vanish. Instead, they will produce observable effects on these physical quantities, which answers the question raised at the second paragraph.
Employing the sub-geometric phases in the density matrix, we proceed to discuss the geometric phase in a mixed state. For a mixed state p k k k ρ ρ = ∑ , where k ρ is the density matrix of the pure state k ( 1, 2 ) k ψ | 〉 = … and p k is the probability for k ψ | 〉 appearing in the mixed state, there is no fixed relative phase including the relative geometric phase between different pure states. Despite some of the treatments proposed [11][12][13][14][15][16][17] for the geometric phase in a mixed state, defining the geometric phase in a mixed state is difficult. Due to the density matrix k ρ for a pure state ψ | 〉 k comprising the sub-geometric phases, we can naturally define the sub-geometric phases in the density matrix of a mixed state using the proposed method. Certainly, there is still no relative geometric phases between different pure states { } k ψ | 〉 , the sub-geometric phases only exist within each pure state.

An example for two States System
Next, the sub-geometric phase is illustrated using two states system with the Hamiltonian www.nature.com/scientificreports www.nature.com/scientificreports/

experimental observation of Sub-Geometric phase
Consider a spin-1/2 particle under a rotating magnetic field Fig. 1 is the magnetic moment of the particle. When the particle passes through the solenoid, it will experiences a periodic rotating magnetic field → B 1 (t). The Hamiltonian for this system can be written as to study. According to the notation in Eqs (10) and (13) In order to observe the sub-geometric phase γ 2 in Eq. (15), we try to devise an experimental setup which is depicted in Fig. 2 The spin-1/2 particles from one source are split into two beams which pass through two solenoids, respectively, then they meet together on the screen and interfere with each other. The magnetic field subject to the two solenoids rotating with the opposite phase, that It is known that the wavefunction on the screen is superposed as has the opposite sub-geometric phase than ψ − t ( ) 1 , then the intensity of interference can be obtained by use of the expression (15)  we can see that the interference is determined by the sub-geometric phases γ 1 and γ 2 , so the sub-geometric phase can be observed by measuring the intensity of interference, which indicates that our sub-geometric phase is physical observable.

Summary and Discussion
As we know that the usual geometric phase, no matter what Berry phase or AA phase, can not appear in the density matrix because it is canceled by its conjugate. To explore the influence of geometric phase on the physical observables, in this paper, we have shown that a generalized sub-geometric phase can exist in the density matrix, therefore have effect on the physical observables averaged by the density matrix. This sub-geometric phase has an unified expression during the adiabatic and nonadiabatic process, it can also be extended to study the geometric phase in a mixed state. It should be pointed out that the above sub-geometric phase in a mixed state is only investigated at zero temperature. In fact, it can be extended to study a mixed state at finite temperature, i.e., with a density matrix as follows: Figure 2. The experimental setup for the interference of two beams of particles with sub-geometric phases which is induced by the rotating magnetic fields produced by the two solenoids with the opposite phase.