真实情况似乎是,我们不愿意承认在世界上还有其他生物具有与人类同等的智慧或者道德,不愿意承认其他生物与我们具有相同的重要性,并且认为其他生物命中注定在世界上扮演某些额外角色…… 我们不打算承认它们或许与我们类似,我们宁可相信除人之外的其他生物都是愚蠢的,只能靠本能生存和活动。 —— 梅特林克《花的智慧》

# 背景设定和定性分析

考虑薛定谔方程

itψ(t)=H(t)ψ(t),i\frac{\partial}{\partial t}|\psi(t)\rangle = H(t)|\psi(t)\rangle,

其中H(t)CN×NH(t)\in \mathbb{C}^{N\times N} 是一个 Hermitian matrix,其中N=2nN=2^n
我们希望找到一个 gate complexity poly(n,t,1/ϵ)poly(n,t,1/\epsilon) 的量子算法,来模拟ψ(0)ψ(t)|\psi(0)\rangle\mapsto |\psi(t)\rangle
(其中我们允许算法输出和ψ(t)|\psi(t)\rangle 的 norm 差ϵ\leq \epsilon
如果这个算法能找到,我们说H(t)H(t) 是可以被 efficiently simulate 的。

下面我们考虑定时 Hamiltonian H(t)H(t),即H(t)=HH(t)=H 是一个常数矩阵。

  1. (Local) 如果HH acts on O(1)O(1) qubits,那么HH 是可以被 efficiently simulate 的。
  2. (Scalability) 如果HH 是可以被 efficiently simulate 的,那么c=poly(n)c=poly(n)cHc\cdot H 也可以被 efficiently simulate 的。
  3. (Unitary) 如果HH 是可以被 efficiently simulate 的,那么UHUUHU^\dagger 也可以被 efficiently simulate 的,其中UU 是 unitary。
  4. (Addition) 如果H1H_1H2H_2 是可以被 efficiently simulate 的,那么H1+H2H_1+H_2 也可以被 efficiently simulate 的。(注意H1,H2H_1,H_2 可能不交换)
  5. (Commutator) 如果H1H_1H2H_2 是可以被 efficiently simulate 的,那么i[H1,H2]i[H_1,H_2] 也可以被 efficiently simulate 的。
  6. (Diagonal) 如果HH 是对角的,那么HH 是可以被 efficiently simulate 的。
  7. (Sparse) 如果HHO(1)O(1)-sparse 的,那么HH 是可以被 efficiently simulate 的。这里会有一些 matrix access oracle,关于 oracle 的 query complexity 也是 polynomial 的。

Example: 对于一个粒子的 Hamiltonian H=p22m+V(x)H=\frac{p^2}{2m}+V(x),其中p2=d2dx2p^2=-\frac{d^2}{dx^2} 是 square of the momentum operator,p22m\frac{p^2}{2m} 是动能部分,V(x)V(x) 是势能部分。势能部分是对角的,P2P^2 可以在 Fourier transformation(unitary)下是对角的,因此HH 是可以被 efficiently simulate 的。

# 定量分析

(一阶 Lie-Trotter-Suzuki Formula)
考虑H=A+BH=A+B,其中A,BA,B 都是 Hermitian。令

U(t)=eiHt,U1(t)=eiAteiBt,U(t)=e^{-iHt},\qquad U_1(t)=e^{-iAt}e^{-iBt},

考虑泰勒展开带余项,

f(t)=f(0)+f(0)t+0tf(s)(ts)ds,f(t)=f(0)+f'(0)t+\int_0^t f''(s)(t-s)ds,

那么是 Hermitian 时,矩阵也可以套入:

U(t)=Ii(A+B)t+0t(A2+AB+BA+B2)ei(A+B)s(ts)ds,U1(t)=Ii(A+B)t+0teiAs(A22ABB2)eiBs(ts)ds.\begin{aligned} & U(t)=I-i(A+B)t+\int_0^t -\left (A^2+AB+BA+B^2\right )e^{-i(A+B)s}(t-s)ds,\\ & U_1(t)=I-i(A+B)t+\int_0^t e^{-iAs}\left (-A^2-2AB-B^2\right )e^{-iBs}(t-s)ds. \end{aligned}

注意到eiAt,eiBte^{-iAt},e^{-iBt} 是 unitary,放在 norm 里可以丢掉(这里不是完全丢掉,需要放缩 2 倍UABV2AB\|UA-BV\|\leq 2\|A-B\| 其中U,VU,V 是 unitary,需要仔细分析),所以

U(t)U1(t)0t2ABBA(ts)ds20t(A+B)2(ts)ds=t2(A+B)2.\|U(t)-U_1(t)\|\leq\int_0^t2\|AB-BA\|(t-s)ds\leq 2\int_0^t (\|A\|+\|B\|)^2(t-s)ds=t^2(\|A\|+\|B\|)^2.

因此,对于分割ll,如果希望通过迭代逼近:

ei(A+B)t(eiAt/leiBt/l)llt2l2(A+B)2ϵ,\left \|e^{-i(A+B)t}-\left (e^{-iAt/l}e^{-iBt/l}\right )^l \right \|\leq l\cdot\frac{t^2}{l^2}(\|A\|+\|B\|)^2\leq\epsilon,

那么l=O(t2(A+B)2ϵ)l=O\left (\frac{t^2(\|A\|+\|B\|)^2}{\epsilon}\right )

(二阶 Lie-Trotter-Suzuki Formula / Gilbert Strang Formula)
考虑H=A+BH=A+B,其中A,BA,B 都是 Hermitian。令

U2(t)=eiAt/2eiBteiAt/2,U_2(t)=e^{-iAt/2}e^{-iBt}e^{-iAt/2},

U2(t)U_2(t)U(t)U(t) 的二阶展开是一样的:U(0)=U2(0),U(0)=U2(0),U(0)=U2(0)U(0)=U_2(0), U'(0)=U_2'(0), U''(0)=U_2''(0),所以要考虑三阶展开。最后可以得到

U(t)U2(t)t3(A+B)3.\|U(t)-U_2(t)\|\leq t^3(\|A\|+\|B\|)^3.

类似地,我们知道

U(t)U2l(t/l)t3(A+B)3l2.\|U(t)-U_2^l(t/l)\|\leq \frac{t^3(\|A\|+\|B\|)^3}{l^2}.

这就比之前的一阶 Trotter formula 要好。

高阶 Suzuki 公式的话

U2k(t)=U2k2(skt)2U2k2((14sk)t)U2k2(skt)2,U_{2k}(t)=U_{2k-2}(s_kt)^2\cdot U_{2k-2}((1-4s_k)t)\cdot U_{2k-2}(s_kt)^2,

其中

sk=14412k1.s_k=\frac{1}{4-4^{\frac{1}{2k-1}}}.

误差是

U(t)U2k(t)t2k+1(A+B)2k+1(K2k+1+1)(2k+1)!,\|U(t)-U_{2k}(t)\|\leq \frac{t^{2k+1}(\|A\|+\|B\|)^{2k+1}(K^{2k+1}+1)}{(2k+1)!},

其中K=25k1K=2\cdot 5^{k-1}。那么类似前面分析,

l=O((tK(A+B))1+12kϵ12k).l=O\left (\frac{(tK(\|A\|+\|B\|))^{1+\frac{1}{2k}}}{\epsilon^{\frac{1}{2k}}}\right ).

实际上已知A,BA,B 和时间tt 后,可以取到一个最优的kk,使得ll 是最小的:

k=max1,log(t(A+B)/ϵ)log5.k=\max\left \lceil 1,\sqrt{\frac{\log (t(\|A\|+\|B\|)/\epsilon)}{\log 5}}\right \rceil.

所以并不是越高阶越好,高阶可能常系数指数爆炸。[Paper:A Theory of Trotter Error, 用 commutator 去 bound]

Interacting picture
考虑H=A+BH=A+B,而且H1>>H2\|H_1\|>>\|H_2\| 时。考虑

itψ(t)=Hψ(t),i\frac{\partial}{\partial t}|\psi(t)\rangle=H|\psi(t)\rangle,

定义

ψ1(t)eiH1tψ(t),|\psi_1(t)\rangle\triangleq e^{iH_1t}|\psi(t)\rangle,

那么

itψ1(t)=i(iH1)eiH1tψ(t)+eiH1tHψ(t)=(H1eiH1t+eiH1t(H1+H2))ψ(t)=eiH1tH2ψ(t)=eiH1tH2eiH1tψ1(t).\begin{aligned} i\frac{\partial}{\partial t}|\psi_1(t)\rangle&=i(iH_1)e^{iH_1t}|\psi(t)\rangle+e^{iH_1t}H|\psi(t)\rangle\\ &=(-H_1e^{iH_1t}+e^{iH_1t}(H_1+H_2))|\psi(t)\rangle\\ &=e^{iH_1t}H_2|\psi(t)\rangle\\ &=e^{iH_1t}H_2e^{-iH_1t}|\psi_1(t)\rangle. \end{aligned}

这就定义了一个新的 Hamiltonian 的薛定谔方程。此时模拟的复杂性可以把H1\|H_1\| 从复杂度里剔除,相当于对H2H_2 做了个 unitary 变换。