摘要： (未经允许,不得转载) 参考文献：J. M. Thijssen, 《Computational Physics》 1. Formalism 为了简化，设一维谐振子的哈密顿量为 $$ \hat{H} = - \frac{d^2}{dx^2} + x^2 $$ 引入拟设波函数 $$ \psi(x, \a 阅读全文

摘要： 这是上次一个小文献笔记(https://www.cnblogs.com/luyi07/p/15442971.html)里一个定理的实践。 1. 实数反对称矩阵 \(M\) 所有矩阵元为实数，并且有反对称性 \(M^\top = - M\)。 2. 反对称阵的正则形式 如果反对称矩阵 \(M\) 是块 阅读全文

摘要： 随便记在这吧 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 阅读全文

摘要： 我们需要做很多很多小矩阵相乘（维数只有几十），但是次数很多，所以用哪个矩阵库的函数对我们很重要。这里写一个很小的测试代码，测试lapack(包含着朴素的blas)，cblas，还有手写函数，对比它们做小矩阵相乘的效率。 对于给定的维数，这三种办法，每种都做１０００次方阵相乘$AB = C$，每次相乘 阅读全文

摘要： 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. 阅读全文

摘要： 核结构经常要画能级图，所以记一下这个代码吧。 # plt 绘制 energy 中的能级，按照 0+,2+,4+,6+,8+,... 的顺序，画在 xleft 处， # 如果 flaglabel == True, 则把 0+, 2+, ..., 画在 xlabel 处。 def EnergyLevel 阅读全文

摘要： 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]; 阅读全文

摘要： 1. 安装 pyinstaller pip3 install pyinstaller 2. 改装 matplotlib 3.1.1 现在的 matplotlib 版本都高于 3.2，而打包工具 pyinstaller需要低于 3.2 的版本，所以需要卸载 matplotlib，然后再安装低版本 pi 阅读全文

摘要： 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 阅读全文

摘要： 有了 Gantmacher 的书（前两篇随笔摘录的内容），Bruno Zumino 的这篇文章很容易 follow，矩阵部分的数学可能个把两个小时就跟着推完了。还有一点关于量子场论的内容，我看不懂，暂时也不需要。 网上查了一下，Bruno Zumino 是超弦理论的提出者之一！没想到和这么重要的人面 阅读全文