等身大ぬいぐるみ ラブドール 6
レス数が1000を超えています。これ以上書き込みはできません。
γ の第 j 成分 γj をそれぞれ γj ± π/4 に置き換えたものを γ±
∂ ⟨Oi⟩γ = ⟨Oi⟩γ+ − ⟨Oi⟩γ− (4.12) ∂γj
を得る. つまり, 量子状態 ρi に
アンザッツ U (γ±) を作用させて,
物理量 Oi を測定して得られた結果 ⟨Oi⟩γ± から ∂γj ⟨Oi ⟩γ を
正確に計算654 ⟨Oi ⟩γ± を有限回の物理量の測定によって推 定するので,
その真の値を得ることはできず, 統計誤差が生じることに注意しておく.
このように ∂γj ⟨Oi⟩γ を 計算する方法をパラメータシフトルール25 URPQC (γ) のような構造を持つアンザッツの構造 に関する
パラメータシフトルールを考えたが, より一般的なアンザッツの構造に対する
パラメータシフトルー ル26 勾配の情報を用いる代表的なオプティマイザーの 1 つとして,
勾配降下法が挙げられる. 勾配降下法とは, コ スト関数のパラメータを
γ(0) に適当に初期化した後
γ(t+1) ← γ(t) − α∇C(γ(t)) (4.13)
のように, 勾配方向にパラメータの更新を行うことを
何度も繰り返すことで, コスト関数の最小点を求める
ア ルゴリズムである
α ∈ R を学習率という. 上述したように,
変分量子アルゴリズムにおけるコスト関
は, γj に依らない (4.11)5241 変分量子アルゴリズム 37
数の勾配の評価では, 有限回の物理量の測定によって
∇C γ(t) を推定していることに注意しなければならな い.
このように, コスト関数の勾配の値を推定する勾配降下法を
一般に確率的勾配降下法54 学習率 α はイテレーション t に依らない定数としていたが,
Adam オプティマイザー のように, 学習率 α を
イテレーション t ごとに更新させることで,
オプティマイザーの収束性を向上54 コスト関数が, 4.2 で述べる変分量子固
有値ソルバーや Fixed input state compiling と
呼ばれる変分量子アルゴリズムのように,
CRPQC (γ) = Tr OURPQC (γ) ρURPQC (γ)† (4.14) URPQC (γ)
この設定の下, パラメータ γj に
注目してコスト関数の解析的な性質
第 t イテレーションにおけるコスト関数の値 C おいて
, γj 以外のパラメータを固定した関数
RPQC
γ(t)に (4.15)
は, γ
C(t) (γj ) := CRPQC (γ)| (t) j γj′=γj′
に依らない実数 a(t), a(t), a(t) を用いて, j 123
′
C(t) (γ ) = a(t) sin2γ + a(t) cos2γ + a(t) jj1j2j3
(4.16) 実数 a(t), a(t), a(t) は, C(t) (γ ) の独立な
3 点の値から計算できる. 例えば, 独立な
123 jj3点として,γ(t),γ(t)+π,γ(t)−π を選べばよい.
すると,C(t)(γ )=a(t)sin2γ +a(t)cos2γ +a(t) は簡単
jj4j4 jj1j2j3 98 コスト関数のあるパラメータ γj について注目してみれば,
C(t) (γj) の最小点を求めることができる.
この最小点を求めるステップを, 注目する j
パラメータ γj の添字 j を変化させて何度も繰り返すことで,
コスト関数の最小点を探索するアルゴリズムを
逐次最小化アルゴリズム54 変分量子アルゴリズ ムのコスト関数を, より一般に
C (γ) = Tr OU (γ) ρU (γ)† (4.17) とし,
第 t イテレーションにおけるコスト関数の値 C γ(tにおいて,
γj 以外のパラメータを固定した関数を
C(t) (γj) とする. 逐次最小化アルゴリズムでは,
パラメータの添字 j を変化させて C(t) (γj) の最小点を求め jj
(j ̸=j)32 |ψ (γ)⟩ = U (γ) |ψ0⟩ は, 初期状態 |ψ0⟩ に
アンザッツ U (γ) を作用させて得られる量 子状態で, 試行状態という.
試行状態は n スピン系の量子状態であるから, 2n 次元の複素内積空間の
単位ベク トルで表現される. 試行状態を愚直に古典コンピュータ上で
表現するには O(2n) ビットが必要であるが, 量子 コンピュータ上であれば
n 量子ビット65 変分量子固有値ソルバーでは, 古典コンピュータ上で は表現しきれない程に
次元の大きい複素内積空間の試行状態から基底状態を探索することができよう.
4.2.2 Fixed input state compiling
Fixed input state compiling (FISC) とは,
n 量子ビットの状態 |ψ0⟩ に作用する V と同等の計算を行う量 子ゲートを求める
変分量子アルゴリズムである [22]. すると, U (γ) をアンザッツすれば, FISC の目的は
コスト関数 Cglobal (γ) = −F 理量
V |ψ0⟩ , U (γ) |0⟩
を最小化65 関数は, 物V |ψ0⟩ = U (γ) |0⟩⊗n (4.21)
を満たす γ を求めることにある.
2 つの量子状態 V |ψ0⟩, U (γ) |0⟩⊗n 間の近さは,例えば,
忠実度 F V |ψ0⟩ , U (γ) |0⟩⊗n987 (4.21) を近似的に満たす U (γ) を求めるには,
忠実度を最大化する γ を求めればよい. これは,
⊗n2
Oglobal :=−(|0⟩⟨0|)⊗n (4.22)
を用いて, Cglobal (γ) = Tr [Oglobal |ψ (γ)⟩ ⟨ψ (γ)|] と書ける
|ψ (γ)⟩ := U (γ)† V |ψ0⟩ とした. Cglobal (γ) は, |ψ (γ)⟩ の
n 量子ビット全てを測定して |0⟩⊗n を得る確率に, −1 をかけることで得られる.
し たがって, Cglobal (γ) は −1 以上 0 以下の値をとりうる.
FISC では, Cglobal 以外のコスト関数847 物理量 Olocal を, n−1
Olocal :=−n
I ⊗|0⟩⟨0|⊗I
(4.23)
1 ⊗j ⊗n−j−1
j=0
として, FISC のコスト関数を
Clocal (γ) = Tr [Olocal |ψ (γ)⟩ ⟨ψ (γ)|] 9 Clocal (γ) は, −1
以上 0 以下の値をとり, γ が (4.21) を満たすときにのみ最小値 −1 をとる.
Clocal (γ) は, |ψ (γ)⟩ の第 j 量 子ビットのみを測定して |0⟩ が得られる
確率の, j = 0, 1, . . . , n − 1 に対する平均に −1 変分量子アルゴリズム 38
Algorithm 1 逐次最小化アルゴリズム
Require: コスト関数は, (4.14) で定義した CRPQC (γ) とする.
1: 2: 3: 4: 5: 6: 7: 8: 9:
パラメータの初期値 γ(0) を定める. while t < tmax do
for j = 1,2,...,Np do for j′ = 1,2,...,Np do
if j′ = j then
C(t) γ(t) , C(t) γ(t) + π , C(t) γ(t) − π を推定し,
a(t), a(t), a(t) を求める. jjjj4jj4 123
γ(t+1) ← arg min a(t) sinγ + a(t) cosγ + a(t) j 1j2j3
else
(t+1) γj′
10: end if 11: end for 12: t←t+1 13: end for
14: end while
← γj′654 変分量子アルゴリズム 41バレンプラトーを引き起こす量子ビット数 n の
変分量子アルゴリズムにおいて, 必要な物理量の測定回数 Ntotal
勾配の情報を用いるオプティマイザー (gradient descent) に限らず,
勾配の情報を用いないオプティ マイザー (Nelder-Mead, Powell, COBYLA) についても,
量子ビット数に対して指数的に多くの物理量の測定5487 量子ビットの数 n に対して O (log n) 程の深さの Alternating Layerd Ansatz という
クラスのア ンザッツを用いた変分量子アルゴリズムについては, (4.22) で定義した
Oglobal のように全ての量子ビットに 作用するような物理量 O を用いる場合には
バレンプラトーが起こる一方で, (4.23) で定義した Olocal ように 一部の量子ビットに
作用するような物理量 O を用いる場合にはバレンプラトーが起こらないことが示されて いる
バレンプラトーが起こらないアルゴリズムの例として, 量子畳み込みニューラルネットワー クや
アンザッツがツリーテンソルネットワーク構造を持つ量子ニューラルネットワーク
バレンプラトーの影響を軽減するためのアルゴリズムも提案され始めている
パラメータの一部をランダムに初期化し, 残りのパラメータをアンザッツが恒等演算子と
なるように選ぶパラ メータ初期化の手法である
アンザッツのパラメータを層ごとに最適化 変分量子アルゴリズムにおける
逐次最小化アルゴリズムの疑似コード.
ることを繰り返す. そのアルゴリズムの
変分量子アルゴリズムを NISQ デバイ ス上で
実装するためには, アンザッツを U (γ) = UNg UNg −1 . . . U2U1 の
ようにハードウエア上で実装可能な 基本ゲートに分解する必要があった.
そこで, NISQ デバイスへの雑音のモデルとして, (3.25) で定義した
Np
N = ⃝ (Dpi ◦ Ui) (4.18)
j=198 Dpi (pi ∈ (0, 1]) は
分極解消チャンネルとし た. この雑音のモデルに対して,
逐次最小化アルゴリズムが剛健であることを裏付ける
次の定理 4.1 定理 4.1 (4.17) で定義したコスト関数 C (γ) を考える. (4.18) で
定義した雑音のモデル N の下で計算される コスト関数を C ̃ (γ) とする.
このとき, 物理量の期待値の推定のための測定回数が無限回であならば,
任意の 自然数 t に対して, 逐次最小化アルゴリズムによって求めた
第 t イテレーション後の C (γ) の最適点と C ̃ (γ) 9 Algorithm 2 一般的な逐次最小化アルゴリズム Require:
コスト関数は (4.17) の形で表される.
1: 2: 3: 4: 5: 6: 7: 8: 9:
10: 11: 12: 13:
パラメータの初期値 γ(0) を定める. while t < tmax do
for j = 1,2,...,Np do for j′ = 1,2,...,Np do
if j′ = j then
γ(t+1) ← arg min C(t) (γj)
else
(t+1) γj′
end if end for
t←t+1 end for
end while
(t) ← γj′
jj γj 一般的な逐次最小化アルゴリズムの疑似コード.
4.2 変分量子アルゴリズムの応用 4.2.1 変分量子固有値ソルバー
変分量子固有値ソルバー (Variational Quantum Eigensolove, VQE) は,
量子系の基底状態とその固有値 (基底エネルギー) を求める
変分量子アルゴリズム9874 変分量子固有値ソルバーでは, 量子コンピュータ上で量子系を表現する必要があるので,
対象となる量子系 をスピン系にマッピングする. 例えば, ジョルダン・ウィグナー変換 や
ブラヴィ・キタエフ変換 によって, フェルミオン系をスピン系にマッピングできる.
こうして, 対象となる量子系を n スピン系にマッ ピングしたハミルトニアン H とすると,
一般にハミルトニアン H は,
(4.19)
(4.20)
n−1 3 n−1 3iiijii
H = hασα + hαβσασβ +··· i=0 α=1 i,j=0 α,β=1 定義 4.2 Γ-値確率変数 γ は一様分布に従うとする
n 量子ビットの変分量子アルゴリズムのコスト関数を
C:Γ∋γ→C(γ)∈Rとし,コスト関数CはC1 級とする
このとき,C(γ)がパラメータγj に関してバ レンプラトーで
あるとは, パラメータ γj に関するコスト関数の 1 階微分の
∂γj C (γ) の期待値が 0 で，分散が ある b > 1 を用いて
O (b−n) ∂C (γ)∂C (γ) −n Eγ ∂γ =0, Vγ ∂γ =O b
jj
定義 4.2 にチェビシェフの不等式 (補題 A.15) を
用いると, 任意の δ > 0 に対し,2ν ∂C(γ)≥δ ≤ 1E
∂C(γ) = 1V ∂C(γ) =O b−n∂γj δ2 γ ∂γj δ2 γ ∂γj
ここで, ν は γ の従う一様分布とした. (4.26) は,
Γ から一様に γ を選んだ時に, その点での γj に
関 する勾配の大きさが δ 以上である確率が,
量子ビットの数 n に対して指数625 スト関数 CRPQC (γ) がバレンプラトーか否か
定義 4.2 に基づき, γ を一様分布 ν に従う確率変数と
みなし, γj に関するコスト関数の勾配の期待値と
分散を計算する
URPQC (γ) を注 目しているパラメータ γj に
依存する部分と依存しない部分に分解して
j−1
j′=1 ∂γj C (γ) が連続であれば,
∂γj C (γ) が確率変数になる
UR (γR = (γ1,γ2,...,γj−1)) := Wj
Uj′ (γj′)Wj′ (4.27)65 V
γ
∂CRPQC (γ) ∂γj
(4n − 1)2 2 Tr[X]2
Np
UL γL = (γj+1,γj+2,...,γNp) :=
として,
URPQC (γ ) = UL (γL ) Uj (γj ) UR (γR ) 6 ∂CRPQC (γ) = Tr OUL(γL)Uj,γj (γj )UR(γR)
ρUR(γR)†Uj (γj )†UL(γL)† ∂γj †††
+Tr OUL(γL)Uj(γj)UR(γR)ρUR(γR Uj,γj(γj) UL(γL) Uj,γj (γj) := ∂γjUj(γj)
, γR, γL はそれぞれ確率変数であり,
それらに依 存する UR(γR), UL(γL) 65 j′=j+1
Uj′ (γj′)Wj′
(4.28)
∂CRPQC (γ)
Eγ ∂γ =0
j n+1 (2) ∆2n
(2) (2)
(ρ) ∆2n (O) ∆2n (Vj )
= 2
∆d (X):=Tr X − d
(2)
(4.32)0 O を (4.22) で定義した Oglobal とし, ρ を
純粋状態とし, Vj をトレースレスとすると,
(2) n (2) (2) −n
∆2n (Vj)=2 ,∆2n (O)=∆2n (ρ)=1−2 であるから,
∂CRPQC (γ)
1 −n
n 2 =O 4 (4.33)
Vγ ∂γ コスト関数 CRPQC (γ) の勾配の期待値が
0 であり, 分散が量子ビットの数 n に対して
指数的
=
2(2 +1)
j RPQC(γ)(·)URPQC(γ)† とし, p := Np j=1
N =Dp ◦U
pj とすると,
(4.35)
(4.36)
Uj (γj)
UL (γL)
Np
N=⃝ Dpj ◦Uj (4.34)
j=1 定義 4.2 に定義したバレンプラトーを, noise-induced
バレンプラトーと区別して, noise-free バレンプラトー
ということもある. noise-free バレンプラトーは,
パラメータ空間のほんとんどの領域コスト
関数の勾配が十分小さくなる現象であった
noise-induced バレンプラトーは,
コ スト関数全体が平坦になっていく NISQ デバイスの制約を超えた長時間発展シミュレーション 5.1 はじめに
量子コンピュータは, その量子ビットの数 n に対して指数的に
大きな O (2n) 次元の情報を表現できる. 一 方, 古典コンピュータは,
その古典ビットの数 n に対して O (n) 次元の情報52 NISQ デバイスの制約を超えた長時間発展シミュレーション
トロッター分解に依る時間発展シミュレーションについて述べ
る. 5.3 では, RQD に依る時間発展シミュレーションについて
5.4 では, RQD に依る格子シュウィン ガーモデルの長時間発展
シミュレーションをいかにして実現したかを述べる. 5.5 では,
サイズの小さな格子 シュウィンガーモデルに対する,
トロッター分解と RQD に依る長時間発展シミュレーション52 5.2 トロッター分解
ハミルトニアン H で表される系の実時間発展を
量子コンピュータ上で行うことを考える.
量子コンピュー タ上でハミルトニアンを表現するために,
考えたい系のハミルトニアンをスピン系のハミルトニアン
Hspin98 e−iH ∆T を O (∆t) のオーダーで近似した
量子ゲートを Utrot (∆T ) := e−iHi ∆T
e−iHspinT =(e−iHspin∆T)M =(Utrot(∆T))M +O(∆T)2(5.2)
であるから,M個のUtrot(∆T)を作用させることで,
欲しかった時間発展演算子e−iHspinT をO(∆T)87 M = T /∆T 個の量子ゲート Utrot (∆T ) が必要となる.
トロッターステップ数 M に比例して, 必要な量子ゲートの
数が増えていく. よって, 計算可 能な量子ゲートの深さが
限られている NISQ デバイスでは, トロッター分解による
長時間の時間発展シミュi(5.1)784 (Utrot (∆T))KU(γˆ1)|0⟩⊗n = U(γ2)|0⟩⊗n を
満たすように最適化されたパラメータ γ2 を求めれば よい.
この手続きを繰り返すことで, 時刻 3K∆T, 4K∆T,... における
量子状態を生成するための浅い量子 ゲート U(γˆ3), U(γˆ4), . . .
を用意することができる. そして, (S2) で, U(γˆ1), U(γˆ2), . . .
を |0⟩⊗n に作用させ ることで, 時刻 K∆T, 2K∆T NISQ デバイスの制約を超えた長時間発展シミュレーション
(a) トロッター分解による時間発展シミュレーション.
シミュレーション時間に比例して, 必要な量子回路の深さが深くなる.
(b) RQD による時間発展シミュレーション. (S1) で時間発展演算子に
対応する量子回路を浅い回路に近似したのち, (S2) で近似回路 を
用いた時間発展365 NISQ デバイスの制約を超えた長時間発展シミュレーション 50
5.4.1 格子シュウィンガーモデル
格子ウィンガーモデルとは, 1 次元空間格子上の
量子電磁力学を記述し, 高エネルギー物理のための
量子アルゴリズムのトイモデルとしてよく用いられる
[73, 74, 75, 76, 77]. n を正の偶数として, 格子間隔 a の n サイ
トの格子ウィンガーモデルのハミルトニアン Hlat は,
n−2 2 n−2 n−1
i†iθ †−iθ ga2j†
Hlat =−2a χje jχj+1 −χj+1e jχj + 2 j=0
Lj +m (−1) χjχj (5.5) j=0
j=0 χj は第 j サイトの質量 m のスタッガードフェルミオンであり,
正準反交換関係
{χ†j,χk} = δjk, {χj,χk} = 0 (5.6)
を満たす. 奇数サイトの非占有状態を電子の存在に,
偶数サイトの占有状態を陽電子の存在に対応させる.
Lj と θj は, 第 j サイトと第 j + 1 サイトのリンク上の
ゲージ場とその共役運動量95 [θj,Lk] = iδjk
を満たす. さらに, フェルミオン場とゲージ場の
相互作用の強さは結合定数 g で特徴付けらている.
(5.7)
(5.8)
ゲージ場の自由度 Lj と θj は, Hlat から取り除くことができる. まず, ガウスの法則 † 1 − (−1)j
を境界条件 L−1 = 0 の下で解くと,
Lj −Lj−1 =χjχj − 2 k
k=0
χj →
χ†kχk − 1 − (−1) 2
Lj =987 n−2 2 n−2j k2 n−1 i† † ga†1−(−1)j†
Hlat =−2a χnχn+1 −χn+1χn + 2 χkχk − 2 +m (−1) χjχj (5.11) j=0 j=0 k=0 j=0
e χj (5.10)
(5.9)
j−1
−iθk8 ジョルダン・ウィグナー変換 [61] j−1
χj → 2
(−iZk) (5.12)
Xj −iYj
k=0 1 サイトあたりの粒子の数
n−2 1
で与えられる. また, 物理量 Q
n−2 j
g2aZk+(−1)
n−1 m j
(XjXj+1 + YjYj+1) + 2 を
数密度と呼ぶことにすると, 数密度 N は,
+ 2
(−1) Zj (5.13)
を用いて,
Hspin = 4a
j=0
2
N = 1 (−1)jZj +1
n−1 1
Q=2 Zj j=0
j=0
2n
j=0
(5.14) 格子シュウィンガーモデル Hspin の時間発展演算子 e−iHspin ∆T の量子ゲート Utrot (∆T )
Hspin は,
n−1 m(−1)j ga n j HZ:= αjZZj withαjZ:= 2 +4 2−2
j=0
j=0 k=j+1
|0011>98 Xj := I⊗j ⊗ X ⊗ In−j−1, Yj := I⊗j ⊗ Y ⊗ In−j−1,
Zj := I⊗j ⊗ Z ⊗ In−j−1 とした. すると, 奇数番目の
量子ビットの |1⟩ が電子の存在に, 偶数番目の量子ビットの |0⟩ が
陽電子の存在に対応54 n−1
j=0 k=0
HXY := HZZ :=
XY XY 1
j=0 n−2
αj (XjXj+1 + YjYj+1) with αj := 4a
n−3 n−2
ZZ ZZ g2a
αjk ZjZk with αjk := 4 (n−k−1)
Hspin = HZ + HXY + HZZ
(5.17) j R 2αZ∆T Zj
(a) e−iHZ ∆T = ∏ e−iαZj ∆T Zj の
各項 e−iαZj ∆T Zj の量子ゲートによる実装. j
j H S • R 2αXY∆T • † H XjS
j+1 H S R 2αXY∆T † H ZjS
(b) e−iHXY ∆T =
∏ XY j e−iαj
XY
∆T (Xj Xj+1+Yj Yj+1) e−iHZ ∆T , e−iHX Y ∆T , e−iHZ Z ∆T の量子ゲートによる実装.
(c) e−iHZZ ∆T =
Utrot (∆T ) := e−iHZ ∆T e−iHX Y ∆T e−iHZ Z ∆T
n−1
Z
n−3 n−2 ZZ
= e−iαj ∆TZj j=0
e−iαj
∆T(XjXj+1+YjYj+1)
n−2
XY
(5.18) e−iαj ∆TZjZk (5.19)
j=0
j=0 k=j+1 [Q, Utrot (∆T )] = 0 であるから,
Utrot (∆T ) に依る時間発 展によって電荷 Q の保存則は保たれる.
5.4.2 アンザッツ: 粒子数保存アンザッツ
RQD が NISQ の雑音下でも機能するためには,
Utrot (∆T ) n 量子ビット系を記述する複素内積空間 H の
部分空間 Hn,m (m = 0,1,...,n) を n−1 n−1
Hn,m = span |ik⟩ | ik ∈ {0,1}, k=0
ik = m
(5.20)
k=09 NISQ デバイスの制約を超えた長時間発展シミュレーション
と定義し, m を粒子数と呼ぶ. すると, H = nm=0 Hn,m である.
また, Hn,m は, n 量子ビットのうち, m 量 子ビットが |1⟩,
n − m 量子ビットが |0⟩ となっている量子状態を正規直交基底として持つ
したがって, Hn,m の次元 dn,m は nCm 8 n サイト格子シュウィンガーモデルの電荷 Q の
固有値 q の固有空間は, Hn, n2 −q である.
q ≥ 0 とする. このとき, 電荷 q の系の状態は
電子が k 個, 陽電子が q + k
個 (k = 0, 1, . . . , n2 − q) 陽電子がq+k個,電子がk個の状態を量子ビット上で
表現すると,粒子数k+(n/2−(q+k))= n2 −qの量子状態
陽電子が q + k 個, 電子が k 個の状態が張る空間の
次元は, n2 Ck · n2 Cq+k であるから, 電荷 q の固有空間の
次元は,n2−q nCq+k ·nCk =nCn−q =dn/ q≥0のとき,電荷Qの固 k=02 2 2 2有値 q の固有空間は,
Hn, n2 −q であることが言えた. q < 0 の場合も同様に
示すことができる. したがって, 格 子シュウィンガーモデルの
電荷 Q の保存則は, 量子ビット上の粒子数の保存則665487 A ̃ (θ )n−2, (φ )n−2 は, 粒子数保存アンザッツ A
(5.22) (θ, φ) の第 l 層に対応
L−1n−2 n−2
An (θl,i)i=0 , (φl,i)i=0 n l,i i=0 l,i i=0
l=0
n,L9 A θl,0,φl,0 ()
A θl,1,φl,1 .
A
()
θl,n/2 , φl,n/2
. ..
()
A θl,n−2 , φl,n−2
A
()
θl,n/2−1 , φl,n/2−1 量子ビットに作用する 3 層の粒子数保存アンザッツ An=4,L=3 (θ, φ) は
粒子数保存アンザッツ An,L (θ, φ) は, 粒子数を保存する A ゲートのみから
構成されているので, 粒 子数保存アンザッツ An,L (θ, φ) もまた粒子数を
保存874 n−1
A (θ0,0, φ0,0)
A (θ1,0, φ1,0)
A (θ0,2, φ0,2)
A (θ1,2, φ1,2)
A (θ2,2, φ2,2)
A (θ0,1, φ0,1)
A (θ1,1, φ1,1)
A (θ2,1, φ2,1)
στj (5.23) τj ∈ {0, 1} は, n−1 τj = m を満たす
Xn,m |0⟩⊗n は粒子数 m の量子状態
Xn,m =
j=0 時刻 t > 0 の状態 |ψ(t)⟩ もまた固有値 q の固有ベクトルで,
|ψ (t)⟩ ∈ Hn, n2 −q を満た す. よって, |ψ(t)⟩ を, ある θˆ と φˆ を
用いて An,L(θˆ, φˆ)Xn, n2 −q |0⟩⊗n8 γ,l,i l,i
1 2
5
C(θ,φ)|θ′
l ,i
′=θ(t) l′,i′
(for(l′,i′)̸=(l,i)), φ′
l ,i
′=φ(t) l′,i′
(γ=θ) (γ=φ)
C(θ,φ)|φ′
l ,i
′=φ(t) l′,i′
(for(l′,i′)̸=(l,i)), θ′
l ,i
′=θ(t) l′,i′
C(t) (γ ) = a(t) sin2γ γ,l,i l,i 1
+ a(t) cos2γ l,i 2
+ a(t) sinγ + a(t) cosγ l,i 3 l,i 4
+ a(t) l,i 5 (t)
Cγ ,l,i (γ1 ) sin 2γ1
(t)
Cγ ,l,i (γ2 ) sin 2γ2
C(t) (γ ) = sin 2γ γ,l,i 3 3
C (t) (γ4 ) sin 2γ4 γ,l,i
cos 2γ1 cos 2γ2
γ ,l,i
sin γ1 sin γ2
cos γ1 cos γ2
(t) 1a1
(t) 1 a2
1 a(t) 3 1 a(t)
598 パラメータの初期値 (θ(0) = (θ(0))l,i, φ(0) = (φ(0))l,i) を定める. l,i l,i
while t < tmax do
for (γ,l,i) ∈ {θ,φ} × {0,1,...,L − 1} × {0,1,...,n − 2} do
for (γ′,l′,i′) ∈ {θ,φ} × {0,1,...,L − 1} × {0,1,...,n − 2} do if (γ,l,i) = (γ′,l′,i′) then
else
(t+1) γl′,i′
end if end for
t←t+1 end for
end while
← γl′,i′987 γ(t+1) ← arg min a(t) sin2γ
+ a(t) cos2γ l,i 2
+ a(t) sinγ l,i 3
+ a(t) cosγ l,i 4
l,i
γl,i (t)
1
1 2
5
+ a(t)
l,i 5 0.5 0.4 0.3
0.2 0.1 0.0
Exact Trotter
RQD (Cglobal)
0 /5 2/5 3/5 time Eθ,φ Eθ,φ
∂γ
∂Clocal(θ,φ) ∂γ
= 0, Vθ,φ
=0, Vθ,φ
∂γ
∂Clocal(θ,φ) ∂γ
= 4bγ dn−2,m−1 dn,m(dn,m+1)2
(5.32)
∂Cglobal(θ,φ) ∂Cglobal(θ,φ) 0
(m = 0, n)
(m=1, 2, ..., n−1)
0
=16bγd2n−2,m−1m(n−m) (m=1, 2, ..., n−1)
(dn,m +1)(d2n,m −1)n3
(5.33) C (γ) = Tr UL (γL)† OUL (γL)UM
(γM)UR (γR)ρUR (γR)† UM (γM)† と
表せるから, γ に関する勾配は,
∂C(γ)† † † ∂γ =Tr UL (γL) OUL (γL)UM,γ (γM)UR (γR)ρUR (γR) UM (γM)
(6.4)
(6.5) (6.6)
(6.7)
+ Tr UL (γL)† OUL (γL) UM (γM ) UR (γR) ρUR (γR)† UM,γ (γM )†
UM,γ (γM ) := ∂UM (γM ) ∂γ
UM (γM ), UM (γM )† （初コメ3号）僕は貼るだけで意見を述べないと告知しているにも関わらず
成りすましのさらに成りすましの無能なレスで多くの無能なスレ民が釣られている。
馬鹿なのですね。
今後も無能なレスは全部成りすましです。念の為。
僕は今後も無言で貼り続けますね。
勝ち負けの問題では有りません。誰とも勝負はしていません。
ぬいドーラーの居場所を永久に潰すだけです。 （今後）
後残り、5500レス分は量子力学です。
その後はフランス構造主義哲学周辺（デリダ、フーコー、ドゥルーズ、ガタリ）を7800レス予定しています。
その後はギリシャ哲学（ピタゴラス、アリストテレス、プラトン、ニーチェまで）や、ハイデガーも15000レス貼ります。
（次回スレの事前予告）
何も言わずに平均100レス連続で貼っているのが僕です。
そのあとで皆さんを怒らせる（呆れさせる）無能なレスをしているのが成りすましです。
くれぐれも今後は釣られないようにして下さい。 このスレッドは１０００を超えました。
もう書けないので、新しいスレッドを立ててくださいです。。。
life time: 7日 0時間 55分 16秒 レス数が1000を超えています。これ以上書き込みはできません。
