鈴木=トロッター分解

演算子の和の指数演算子を各々の演算子の指数演算子の積に分解する公式

物理学において、鈴木=トロッター分解(すずき=トロッターぶんかい、: Suzuki–Trotter decompostion)、または鈴木=トロッター公式: Suzuki-Trotter formula)とは演算子の和の指数演算子を各々の演算子の指数演算子の積に分解する公式[1][2][3]量子系時間発展分配関数の計算に応用される。数学者ヘイル・トロッター英語版バナッハ空間半群の研究において、トロッター公式を与えた[4]。日本の物理学者 鈴木増雄トロッター公式を拡張した指数積の高次分解公式を研究し、系統的に求める手法を提案した。

概要

編集

量子力学における量子状態の時間発展や量子統計力学における分配関数の計算ではハミルトニアンの指数演算子が現れる。これを数値計算するには二つの非可換な演算子 A,B の指数演算子 ex(A+B) の計算が必要になる。ここでx は実数または複素数のパラメータである。トロッター公式によれば、指数演算子は

 

を満たし、

 

が成り立つ。但し、O( )ランダウの記号である。 より高次の m 次分解公式

 

が得られると、パラメータ x が小さくない場合にも、

 

であり、十分 n を大きくとれば、ex(A+B) の高精度の近似が得られる。 こうした指数演算子の積の形での近似計算では、量子力学での波動関数の時間発展が満たすユニタリ性や古典力学でのハミルトン力学系英語版の時間発展が満たすシンプレクティック性が保証されるなどの利点がある。このような指数演算子の分解公式を鈴木=トロッター分解という[1]。 例えば、

 

は2次の公式の1つの例である また、Ruthの公式として知られている

 

は6次の公式の例の1つである。 一般に m 次分解公式 Fm(xA, xB) を与えるには、条件を満たすパラメータの組 {ti}i=1,…,M を定める必要がある。こうしたパラメータの組の求め方には、AB交換子積のなす自由リー代数の理論に基づく方法や漸化式による方法、量子解析による方法がある[2][3]

漸化式による構成

編集

Qm-1(x)=Qm-1(xA,xB)ex(A+B)m-1 次分解公式であるとき、漸化式

 

で定義される Qm(x)m 次分解公式となるためにpm,1,…pm,s が条件

 

を満たせば良い。 m が奇数 2l-1 である場合、この条件を満たすpm,1,…pm,s の実数の組が存在するが、m が偶数 2l である場合、この上記の条件を満たす実数の組は存在しない。 しかしながら、一般に奇数 2l-1 次の分解公式 S2l-1(x) が対称性の条件

 

を満たすとき、これは自動的に 2l 次の分解公式となる。よって、奇数の場合にこの漸化式を適用していくことでパラメータを実数とする高次の公式が導かれる。但し、Sm(x) が対称性の条件を満たすようにパラメータには pm,s+1-i =pm,i の条件を課すものとする。

例えば、s=3 の場合、奇数 m における上記の条件を満たすパラメータの組として、

 
 

をとることができる。2次の公式

 

は対称性の条件 S2(x)S2(−x)=1 を満すことから

 
 

S3(x) を構成できる。

脚注

編集

参考文献

編集

論文

編集
  • Trotter, H. F. (1959). “On the product of semi-groups of operators”. Proceedings of the American Mathematical Society 10 (4): 545–551. doi:10.2307/2033649. ISSN 0002-9939. MR0108732. 
  • Suzuki, Masuo (1976). “Generalized Trotter's formula and systematic approximants of exponential operators and inner derivations with applications to many-body problems”. Comm. Math. Phys. 51: 183–190. doi:10.1007/bf01609348. 

書籍

編集

関連項目

編集