张量网络编程学习笔记(5-1):量子多体含时演化,Trotter-Suzuki近似

本文仅从本人复习需要出发:作为对首都师范大学冉仕举老师课程内容的整理、归纳以及补充 

课程详情:张量网络PyThon编程:3.4 量子线路模块化编程(a) ADGate类_哔哩哔哩_bilibili

 冉仕举老师本人空间:StringCNU的个人空间_哔哩哔哩_bilibili

详见上一讲:最小代数本征值问题,以及其对应的本征态。有明确的统计意义。

当Srodinger-Equation 的 空间部分可以分离出来,|\phi> = \psi(x) |\phi(t)>

将算符的积分写成黎曼和的形式,固定每个切片\hat{H}(t_k),布置一系列门,该层门当中H 视作不变

 事实上,算符的指数求和并不能严格等价成简单的连乘形式,会引入Trotter 误差。

 以下摘录了 P73-P76,对 Trotter 误差的 推导:

(1)编时乘积传播子:\hat{U}(t) \approx \prod_{k=0}^{K-1} \hat{U}(k\tau,k\tau+\tau)= \prod_{k=0} ^{K-1} \hat{U_{tk}}(\tau)=\prod_{k=0}^{K-1} e^{-i\tau \hat{H}(t_k)}

(2)对于每层传播子: \underset{\tau \rightarrow 0}{lim}\hat{U_{t_k}}(\tau) = \underset{\tau \rightarrow 0}{lim} e^{-i\tau\sum_{ij}^{}\hat{H_{t_k}}_{ ij}} = \underset{\tau\rightarrow 0}{lim} \prod_{ij}^{} e^{-i\tau \hat{H_{t_k}}_{ ij}} + O(f(\tau))

(3)连乘 -演化算符的差距:也类似 e^{-i\tau \hat{H}(t_k)} : \tau * H 的求和 - 连乘,累加为:K \times O(f(\tau))

(4)Trotter 误差大概就是举AB两个简单的例子描述这种连乘-指数累加的差距:

e^{\tau(\hat{A}+\hat{B})} - e^{\tau \hat{A}}e^{\tau \hat{B}} \approx - \frac{\tau^2}{2} [\hat{A},\hat{B}] = O(\tau^2),故整套含时演化下来,误差控制可以控制在 t/\tau \times O(\tau^2) = t \times O(\tau)

同理,我们可以用有限步骤,将演化算符邻接交换,stairs - bricks 结果相差O(\tau^2),brick更紧凑

值得一提的是,令 t = i \beta ,演化算符写成 e^{-\beta \hat{H}} :

 * 对着参数演化到 \beta \to \infty ,e^{-\beta \hat{H}} \to |\phi_0><\phi_0|

 * \underset{\beta \to \infty}{lim} e^{-\beta\hat{H}}|\phi> \to \frac{1}{Z} |\phi_0>系统会落到基态上,

但毕竟没有可能实现 \beta \to \infty ,找出出现不动点的地方: \varsigma = ||\phi>-|\phi'>| < \varepsilon

同时要观察到演化算符| e^{-\beta \hat{H}}| < 1,属于压缩映射:  

 * 我们称之为虚时演化。 直接(整块地)计算 e^{-\beta \hat{H}}的虚时演化,每层复杂度~ O(2^{2N}) ,

 * 因此:每层还是考虑和原有含时演化一样布置门,

只是  \hat{U}^{i,j}(\tau) = e^{-i\tau \hat{H}_{ij}} \to \hat{U}^{i,j}(\tau) = e^{-\tau \hat{H}_{ij}} ,就可以用来计算基态了!

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值