python 遍历矩阵_多体动力学演化python入门——quantum many-body scars 和稀疏矩阵后续...

好久没更新了，肚子里也没什么货，就算python稀疏矩阵的最后一篇吧。之前的

路飞的哥哥：多体物理python入门——Ising模型和稀疏矩阵​zhuanlan.zhihu.com

计算了本征值，也就是平衡态问题。这次是动力学演化的演示代码。量子多体物理的动力学演化是个很大的方向（本人不是特别了解，可能有误），大致有多体局域化（MBL），热化（ETH），scar。多体局域化算大方向了，一般研究无序系统。热化比较小众，不是很了解。scar是又新又小众。大致上我理解，对于非无序系统（有序系统），大部分态演化都会热化，特殊情况会有scar。总体来讲非平衡态问题还是有很多未开垦的荒野。

物理我也讲不清楚，不献丑了，详细的scar可以参考文献。

Quantum many-body scars​arxiv.org

大致上，就是某些系统，在某些初态下，经历的时间演化比较特殊。经典上，如题目图中台球所示，一般的入射角度（蓝色）会使台球（不考虑损耗）遍历空间所有点，而存在特殊的入射角度（黑色）台球只会在一定空间里来回往复。

文献主要的内容是，对于PXP模型

特殊初态

（实心表示0，对应

；空心表示1，对应

）的演化，呈周期性，部分图如下

为演示动力学演化，只重复最后一个图，代码如下，权做参考

实现哈密顿量

实现态，和态演化，基本是一两分钟能算完，如果L为16的话，等最后的一行代码的时间能下一把云顶之弈。基本逻辑是实现

最后作图

文章中用的iTEBD，其实用精确解算12个点也能取得很好的效果，基本一样，在我看来。一般来讲，一个好的工作，数值解只是用来辅助说明。真正有意义的在于你对其中物理机制的理解，据此建立有效理论甚至经典理论解释。如果你的工作都是数值结果说明的话，建议用比较现代的DMRG，iTEBD等算法及软件。稀疏矩阵，我感觉20个点到顶了。16个点我还心疼我的笔记本。

我不知道为什么最近好多人关注我，有关注我的评论说下原因，我挺懵逼的。好像也不是文章的原因吧，都没几个赞。其实个人水平不咋地，我平时也用不到稀疏矩阵，只是上次无聊现学的。如果代码对你有帮助，或者里面的物理引起了你的兴趣，实乃我的荣幸。看情况更新吧，如果有很多赞的话。下次更新也是其他内容了，这个稀疏矩阵感觉就那么多内容。多看看文档就行。quantum many-body scar我翻译为量子多体隐痕，感觉有点味道，不知道对不对。祝大家科研顺利。