python实现: VMC做一维谐振子
摘要：
(未经允许,不得转载) 参考文献：J. M. Thijssen, 《Computational Physics》 1. Formalism 为了简化，设一维谐振子的哈密顿量为 ˆH=−d2dx2+x2 引入拟设波函数 $$ \psi(x, \a 阅读全文
posted @ 2021-11-29 20:09 luyi07 阅读(482) 评论(0) 推荐(0) 编辑
实反对称矩阵正则化
摘要：
这是上次一个小文献笔记(https://www.cnblogs.com/luyi07/p/15442971.html)里一个定理的实践。 1. 实数反对称矩阵 M 所有矩阵元为实数，并且有反对称性 M^\top = - M。 2. 反对称阵的正则形式 如果反对称矩阵 M 是块 阅读全文
posted @ 2021-11-19 16:57 luyi07 阅读(505) 评论(0) 推荐(0) 编辑
PVPC kb3g pn/upn Ti44 LAP 脚本
摘要：
随便记在这吧 LAPkb3g.sh #!/bin/bash # ./LAPkb3g.sh 1 1 example/kb3g.jsp example/kb3g.pn y 10 0.001 0.001 projection_pf0202_p.dat projection_pf0202_n.dat nuc 阅读全文
posted @ 2021-11-17 20:23 luyi07 阅读(58) 评论(0) 推荐(0) 编辑
小矩阵相乘效率对比：lapack, cblas, 手写函数
摘要：
我们需要做很多很多小矩阵相乘（维数只有几十），但是次数很多，所以用哪个矩阵库的函数对我们很重要。这里写一个很小的测试代码，测试lapack(包含着朴素的blas)，cblas，还有手写函数，对比它们做小矩阵相乘的效率。 对于给定的维数，这三种办法，每种都做１０００次方阵相乘AB = C，每次相乘 阅读全文
posted @ 2021-11-16 18:22 luyi07 阅读(869) 评论(0) 推荐(0) 编辑
Ising模型的Metropolis模拟中的组态阻塞和振荡
摘要：
Ising模型哈密顿量为 \hat{H} = - J \sum_{i,j相邻} s_i s_j - B \sum_i s_i 若取 J=1, B=0，做 4 x 4 一共 16 个 spin 的 Ising 模型模拟，可以得到很平缓的相变曲线，如下图 曲线最左侧有些点低于低于1. 阅读全文
posted @ 2021-11-07 11:15 luyi07 阅读(106) 评论(0) 推荐(0) 编辑
python 画能级图
摘要：
核结构经常要画能级图，所以记一下这个代码吧。 # plt 绘制 energy 中的能级，按照 0+,2+,4+,6+,8+,... 的顺序，画在 xleft 处， # 如果 flaglabel == True, 则把 0+, 2+, ..., 画在 xlabel 处。 def EnergyLevel 阅读全文
posted @ 2021-10-31 11:25 luyi07 阅读(721) 评论(0) 推荐(0) 编辑
python 调用 shell 命令，制作用户界面
摘要：
1. python调用shell指令: os.system 格式为 os.system('command') import os os.system('make') flag0 = flag[0]; flag1 = flag[1]; flag2 = flag[2]; flag3 = flag[3]; 阅读全文
posted @ 2021-10-29 10:02 luyi07 阅读(422) 评论(0) 推荐(0) 编辑
打包 python 程序，变成一个可执行文件
摘要：
1. 安装 pyinstaller pip3 install pyinstaller 2. 改装 matplotlib 3.1.1 现在的 matplotlib 版本都高于 3.2，而打包工具 pyinstaller需要低于 3.2 的版本，所以需要卸载 matplotlib，然后再安装低版本 pi 阅读全文
posted @ 2021-10-29 09:23 luyi07 阅读(331) 评论(0) 推荐(0) 编辑
宇称投影
摘要：
1. 宇称投影算符 宇称投影算符\hat{P}作用在任意函数f(x)得到 \hat{P} f(x) = f(-x), 三维空间中则为 \hat{P} f(\vec{r}) = f(-\vec{r}). 如果有 $\hat{P} f(\vec{r}) = f(-\vec{r 阅读全文
posted @ 2021-10-26 12:06 luyi07 阅读(371) 评论(0) 推荐(0) 编辑
Bruno Zumino, "Normal Forms of Complex Matrices", Journal of Mathematical Physics 3, 1055 (1962)
摘要：
有了 Gantmacher 的书（前两篇随笔摘录的内容），Bruno Zumino 的这篇文章很容易 follow，矩阵部分的数学可能个把两个小时就跟着推完了。还有一点关于量子场论的内容，我看不懂，暂时也不需要。 网上查了一下，Bruno Zumino 是超弦理论的提出者之一！没想到和这么重要的人面 阅读全文
posted @ 2021-10-23 13:06 luyi07 阅读(133) 评论(0) 推荐(0) 编辑
