断熱定理 (英語 : adiabatic theorem )は、ハミルトニアン がゆるやかに時間変化する状況では、ある時刻に系を(その時刻での)ハミルトニアンの一つの固有状態として用意した場合に、縮退がない場合には、時間発展した後の状態は当初の固有状態に対応する固有状態であるという定理である[ 1] 
 
  
   
ハミルトニアンが時間に依存しない場合、エネルギーが
       
        
         
          
          
           E
            
           
           
            n
             
            
           
          
         
       
        {\displaystyle E_{n}}
         
        
       
       
        
         
          
          
           |
            
           
          
          
           ψ
            
           
           
            n
             
            
           
         
          ⟩
           
          
         
       
        {\displaystyle |\psi _{n}\rangle }
         
        
       
       
        
         
          
          
           |
            
           
          
          
           Ψ
            
           
           
            n
             
            
           
         
          (
           
         
          t
           
         
          )
           
         
          ⟩
           
          
         
       
        {\displaystyle |\Psi _{n}(t)\rangle }
         
        
       
 
  
   
    
        
         
          
           
           
            |
             
            
           
           
            Ψ
             
            
            
             n
              
             
            
          
           (
            
          
           t
            
          
           )
            
          
           ⟩
            
          
           =
            
          
           exp
            
          
           
            
          
           (
            
          
           −
            
           
           
            i
             
            
           
           
            E
             
            
            
             n
              
             
            
          
           t
            
           
           
            /
             
            
          
           ℏ
            
          
           )
            
           
           
            |
             
            
           
           
            ψ
             
            
            
             n
              
             
            
          
           ⟩
            
           
          
        
         {\displaystyle |\Psi _{n}(t)\rangle =\exp(-\mathrm {i} E_{n}t/\hbar )|\psi _{n}\rangle }
          
         
         
    
  と記述できる。ハミルトニアンが時間に依存する場合、時刻ごとに対応するハミルトニアンが異なるわけだが、ハミルトニアンがエルミート演算子であることは変わらないので、いずれの時刻においてもスペクトル分解が可能であり、そのような、各時刻tでのハミルトニアンと対応する固有値と固有状態について、 
 
  
   
    
        
         
          
           
            
             
             
              H
               
             
              ^
               
              
             
            
          
           (
            
          
           t
            
          
           )
            
           
           
            |
             
            
          
           ψ
            
          
           (
            
          
           t
            
          
           )
            
          
           ⟩
            
          
           =
            
           
           
            E
             
            
            
             n
              
             
            
          
           (
            
          
           t
            
          
           )
            
           
           
            |
             
            
           
           
            ψ
             
            
            
             n
              
             
            
          
           (
            
          
           t
            
          
           )
            
          
           ⟩
            
           
          
        
         {\displaystyle {\hat {H}}(t)|\psi (t)\rangle =E_{n}(t)|\psi _{n}(t)\rangle }
          
         
         
    
  と表しておくことにする。エルミート演算子の固有ベクトルに対しては、正規直交性 
 
  
   
    
        
         
          
          
           ⟨
            
           
           
            ψ
             
            
            
             n
              
             
            
          
           (
            
          
           t
            
          
           )
            
           
           
            |
             
            
           
           
            ψ
             
            
            
             m
              
             
            
          
           (
            
          
           t
            
          
           )
            
          
           ⟩
            
          
           =
            
           
           
            δ
             
            
            
             n
              
            
             m
              
             
            
           
          
        
         {\displaystyle \langle \psi _{n}(t)|\psi _{m}(t)\rangle =\delta _{nm}}
          
         
         
    
  を課すことが可能である。 
  さて、ここでシュレーディンガー方程式  
 
  
   
    
        
         
          
           
           
            i
             
            
          
           ℏ
            
           
            
             
             
              d
               
              
             
              
              
               d
                
               
             
              t
               
              
             
            
           
           
            |
             
            
          
           ψ
            
          
           (
            
          
           t
            
          
           )
            
          
           ⟩
            
          
           =
            
           
            
             
             
              H
               
             
              ^
               
              
             
            
           
           
            |
             
            
          
           ψ
            
          
           (
            
          
           t
            
          
           )
            
          
           ⟩
            
           
          
        
         {\displaystyle \mathrm {i} \hbar {\frac {\mathrm {d} }{\mathrm {d} t}}|\psi (t)\rangle ={\hat {H}}|\psi (t)\rangle }
          
         
         
    
  について考察し、断熱定理を満たすような形で状態が時間発展することを確かめたい。 
  まず、エルミート演算子の固有状態はヒルベルト空間の完全系を張っていることを利用して、解となる状態は、それぞれの時刻にそれぞれの時刻での固有状態で展開できることを用いて、 
 
  
   
    
        
         
          
           
           
            |
             
            
          
           Ψ
            
          
           (
            
          
           t
            
          
           )
            
          
           ⟩
            
          
           =
            
           
           
            ∑
             
            
            
             n
              
             
            
           
           
            c
             
            
            
             n
              
             
            
          
           (
            
          
           t
            
          
           )
            
          
           exp
            
          
           
            
          
           (
            
           
           
            i
             
            
           
           
            θ
             
            
            
             n
              
             
            
          
           (
            
          
           t
            
          
           )
            
          
           )
            
           
           
            |
             
            
           
           
            ψ
             
            
            
             n
              
             
            
          
           (
            
          
           t
            
          
           )
            
          
           ⟩
            
           
          
        
         {\displaystyle |\Psi (t)\rangle =\sum _{n}c_{n}(t)\exp(\mathrm {i} \theta _{n}(t))|\psi _{n}(t)\rangle }
          
         
         
    
  と記述しておくことにする。ただし、以後の計算の都合、位相について 
 
  
   
    
        
         
          
           
           
            θ
             
            
            
             n
              
             
            
          
           (
            
          
           t
            
          
           )
            
          
           =
            
          
           −
            
          
           1
            
           
           
            /
             
            
          
           ℏ
            
           
           
            ∫
             
            
            
             0
              
             
            
            
             t
              
             
            
           
           
            E
             
            
            
             n
              
             
            
          
           (
            
           
           
            t
             
           
            ′
             
            
          
           )
            
           
           
            d
             
            
           
           
            t
             
           
            ′
             
            
           
          
        
         {\displaystyle \theta _{n}(t)=-1/\hbar \int _{0}^{t}E_{n}(t')\mathrm {d} t'}
          
         
         
    
  の関数を導入して調整することにしておく。 
  次に、この解をシュレディンガー方程式に代入し、
       
        
         
         
          ⟨
           
          
          
           ψ
            
           
           
            m
             
            
           
          
          
           |
            
           
          
         
       
        {\displaystyle \langle \psi _{m}|}
         
        
       
 
  
   
    
        
         
          
           
           
            ∑
             
            
            
             n
              
             
            
           
           
            c
             
            
            
             n
              
             
            
           
           
            δ
             
            
            
             n
              
            
             m
              
             
            
          
           exp
            
          
           
            
          
           (
            
           
           
            i
             
            
           
           
            θ
             
            
            
             n
              
             
            
          
           )
            
          
           =
            
          
           −
            
           
           
            ∑
             
            
            
             n
              
             
            
           
           
            c
             
            
            
             n
              
             
            
          
           ⟨
            
           
           
            ψ
             
            
            
             m
              
             
            
           
           
            |
             
            
           
            
             
              
              
               ψ
                
               
               
                n
                 
                
               
             
              ˙
               
              
             
            
          
           ⟩
            
          
           exp
            
          
           
            
          
           (
            
           
           
            i
             
            
           
           
            θ
             
            
            
             n
              
             
            
          
           )
            
           
          
        
         {\displaystyle \sum _{n}c_{n}\delta _{nm}\exp(\mathrm {i} \theta _{n})=-\sum _{n}c_{n}\langle \psi _{m}|{\dot {\psi _{n}}}\rangle \exp(\mathrm {i} \theta _{n})}
          
         
         
    
  
   
    
        
         
          
           
            
             
              
              
               c
                
              
               ˙
                
               
              
             
            
            
             m
              
             
            
          
           (
            
          
           t
            
          
           )
            
          
           =
            
          
           −
            
           
           
            ∑
             
            
            
             n
              
             
            
          
           ⟨
            
           
           
            ψ
             
            
            
             m
              
             
            
           
           
            |
             
            
           
            
             
              
              
               ψ
                
               
               
                m
                 
                
               
             
              ˙
               
              
             
            
          
           ⟩
            
          
           exp
            
          
           
            
          
           (
            
           
           
            i
             
            
          
           (
            
           
           
            θ
             
            
            
             n
              
             
            
          
           −
            
           
           
            θ
             
            
            
             m
              
             
            
          
           )
            
          
           )
            
           
          
        
         {\displaystyle {\dot {c}}_{m}(t)=-\sum _{n}\langle \psi _{m}|{\dot {\psi _{m}}}\rangle \exp(\mathrm {i} (\theta _{n}-\theta _{m}))}
          
         
         
    
  の条件が得られる。ここで、固有方程式(i)の両辺をtで微分する。 
 
  
   
    
        
         
          
           
            
             
             
              H
               
             
              ˙
               
              
             
            
           
           
            |
             
            
           
           
            ψ
             
            
            
             n
              
             
            
          
           ⟩
            
          
           +
            
          
           H
            
           
           
            |
             
            
           
            
             
              
              
               ψ
                
              
               ˙
                
               
              
             
            
            
             n
              
             
            
          
           ⟩
            
          
           =
            
           
            
             
              
              
               E
                
               
               
                n
                 
                
               
             
              ˙
               
              
             
            
           
           
            |
             
            
           
           
            ψ
             
            
            
             n
              
             
            
          
           ⟩
            
          
           +
            
           
           
            E
             
            
            
             n
              
             
            
           
           
            |
             
            
           
            
             
              
              
               ψ
                
              
               ˙
                
               
              
             
            
            
             n
              
             
            
          
           ⟩
            
           
          
        
         {\displaystyle {\dot {H}}|\psi _{n}\rangle +H|{\dot {\psi }}_{n}\rangle ={\dot {E_{n}}}|\psi _{n}\rangle +E_{n}|{\dot {\psi }}_{n}\rangle }
          
         
         
    
  
   
    
        
         
          
          
           n
            
          
           ≠
            
          
           m
            
           
          
        
         {\displaystyle n\neq m}
          
         
        
        
         
          
          
           ⟨
            
           
           
            ψ
             
            
            
             m
              
             
            
           
           
            |
             
            
           
          
        
         {\displaystyle \langle \psi _{m}|}
          
         
         
    
  
   
    
        
         
          
          
           ⟨
            
           
           
            ψ
             
            
            
             m
              
             
            
           
           
            |
             
            
           
            
             
             
              H
               
             
              ˙
               
              
             
            
           
           
            |
             
            
           
           
            ψ
             
            
            
             n
              
             
            
          
           ⟩
            
          
           =
            
          
           (
            
           
           
            E
             
            
            
             n
              
             
            
          
           −
            
           
           
            E
             
            
            
             m
              
             
            
          
           )
            
          
           ⟨
            
           
           
            ψ
             
            
            
             m
              
             
            
           
           
            |
             
            
           
            
             
              
              
               ψ
                
               
               
                m
                 
                
               
             
              ˙
               
              
             
            
          
           ⟩
            
           
          
        
         {\displaystyle \langle \psi _{m}|{\dot {H}}|\psi _{n}\rangle =(E_{n}-E_{m})\langle \psi _{m}|{\dot {\psi _{m}}}\rangle }
          
         
         
    
  を得られる。また、(vii)式に(iv)式の結果を用いることによって、 
 
  
   
    
        
         
          
           
            
             
              
              
               c
                
              
               ˙
                
               
              
             
            
            
             m
              
             
            
          
           (
            
          
           t
            
          
           )
            
          
           =
            
          
           −
            
           
           
            c
             
            
            
             m
              
             
            
          
           ⟨
            
           
           
            ψ
             
            
            
             m
              
             
            
           
           
            |
             
            
           
            
             
              
              
               ψ
                
               
               
                m
                 
                
               
             
              ˙
               
              
             
            
          
           ⟩
            
          
           −
            
           
           
            ∑
             
            
            
             n
              
            
             ≠
              
            
             m
              
             
            
           
           
            c
             
            
            
             n
              
             
            
           
            
             
             
              ⟨
               
              
              
               ψ
                
               
               
                m
                 
                
               
              
              
               |
                
               
              
               
                
                
                 H
                  
                
                 ˙
                  
                 
                
               
              
              
               |
                
               
              
              
               ψ
                
               
               
                n
                 
                
               
             
              ⟩
               
              
             
              
              
               E
                
               
               
                n
                 
                
               
             
              −
               
              
              
               E
                
               
               
                m
                 
                
               
              
             
            
          
           exp
            
          
           
            
          
           (
            
          
           −
            
           
           
            i
             
            
           
           
            /
             
            
          
           ℏ
            
           
           
            ∫
             
            
            
             0
              
             
            
            
             t
              
             
            
          
           [
            
           
           
            E
             
            
            
             n
              
             
            
          
           (
            
           
           
            t
             
           
            ′
             
            
          
           )
            
          
           −
            
           
           
            E
             
            
            
             m
              
             
            
          
           (
            
           
           
            t
             
           
            ′
             
            
          
           )
            
          
           ]
            
           
           
            d
             
            
           
           
            t
             
           
            ′
             
            
          
           )
            
           
          
        
         {\displaystyle {\dot {c}}_{m}(t)=-c_{m}\langle \psi _{m}|{\dot {\psi _{m}}}\rangle -\sum _{n\neq m}c_{n}{\langle \psi _{m}|{\dot {H}}|\psi _{n}\rangle \over {E_{n}-E_{m}}}\exp(-\mathrm {i} /\hbar \int _{0}^{t}[E_{n}(t')-E_{m}(t')]\mathrm {d} t')}
          
         
         
    
  という関係を得ることができる。 
  さて、ここで断熱近似、すなわち、ハミルトニアンの時間変化が十分に小さいということで、
       
        
         
         
          ⟨
           
          
          
           ψ
            
           
           
            m
             
            
           
          
          
           |
            
           
          
           
            
            
             H
              
            
             ˙
              
             
            
           
          
          
           |
            
           
          
          
           ψ
            
           
           
            n
             
            
           
         
          ⟩
           
         
          =
           
         
          0
           
          
         
       
        {\displaystyle \langle \psi _{m}|{\dot {H}}|\psi _{n}\rangle =0}
         
        
       
 
  
   
    
        
         
          
           
           
            c
             
            
            
             m
              
             
            
          
           (
            
          
           t
            
          
           )
            
          
           =
            
           
           
            c
             
            
            
             m
              
             
            
          
           (
            
          
           0
            
          
           )
            
          
           exp
            
          
           
            
          
           (
            
           
           
            i
             
            
           
           
            γ
             
            
            
             m
              
             
            
          
           (
            
          
           t
            
          
           )
            
          
           )
            
           
          
        
         {\displaystyle c_{m}(t)=c_{m}(0)\exp(\mathrm {i} \gamma _{m}(t))}
          
         
         
    
  という解が得られる。ただし、 
 
  
   
    
        
         
          
           
           
            γ
             
            
            
             m
              
             
            
          
           (
            
          
           t
            
          
           )
            
          
           ≡
            
           
           
            i
             
            
           
           
            ∫
             
            
            
             0
              
             
            
            
             t
              
             
            
          
           ⟨
            
           
           
            ψ
             
            
            
             m
              
             
            
          
           (
            
          
           t
            
          
           )
            
           
           
            |
             
            
           
            
             
              
              
               ψ
                
              
               ˙
                
               
              
             
            
            
             m
              
             
            
          
           (
            
           
           
            t
             
           
            ′
             
            
          
           )
            
          
           ⟩
            
           
           
            d
             
            
           
           
            t
             
           
            ′
             
            
           
          
        
         {\displaystyle \gamma _{m}(t)\equiv \mathrm {i} \int _{0}^{t}\langle \psi _{m}(t)|{\dot {\psi }}_{m}(t')\rangle \mathrm {d} t'}
          
         
         
    
  と記述されるもので、断熱発展が周期的なものであれば、γはBerry位相 となる。(x)式の結果を(iv)式に適用することで、 
 
  
   
    
        
         
          
           
           
            |
             
            
           
           
            Ψ
             
            
            
             n
              
             
            
          
           (
            
          
           t
            
          
           )
            
          
           ⟩
            
          
           =
            
          
           exp
            
          
           
            
          
           (
            
          
           −
            
           
           
            i
             
            
           
           
            θ
             
            
            
             n
              
             
            
          
           (
            
          
           t
            
          
           )
            
          
           t
            
           
           
            /
             
            
          
           ℏ
            
          
           )
            
          
           exp
            
          
           
            
          
           (
            
          
           −
            
           
           
            i
             
            
           
           
            γ
             
            
            
             n
              
             
            
          
           (
            
          
           t
            
          
           )
            
           
           
            /
             
            
          
           ℏ
            
          
           )
            
           
           
            |
             
            
           
           
            ψ
             
            
            
             n
              
             
            
          
           (
            
          
           t
            
          
           )
            
          
           ⟩
            
           
          
        
         {\displaystyle |\Psi _{n}(t)\rangle =\exp(-\mathrm {i} \theta _{n}(t)t/\hbar )\exp(-\mathrm {i} \gamma _{n}(t)/\hbar )|\psi _{n}(t)\rangle }
          
         
         
    
  という結論を得ることができる。