Processing math: 100%

随笔分类 -  原子核结构

原子核四极形变的 python 可视化
摘要:把原子核当作液滴,讨论较小的形变,经过一系列处理,可以用公式处理: \[R = R_0 \left\{ 1 + \beta \sqrt{\frac{5}{16\pi}} (\cos\gamma (3\cos^2\theta -1) + \sqrt{3}\sin\gamma\sin^2\theta\c 阅读全文

posted @ 2024-04-08 19:35 luyi07 阅读(72) 评论(0) 推荐(0) 编辑

Bloch-Messiah定理的证明
摘要:参考书:Peter Ring, Peter Shuck, "The nuclear many-body problem", Appendix E. Appendix E 中有 Bloch-Messiah 定理的证明,非常简要,我发现难以读懂,就花了些功夫推明白以后,整理了笔记。录在这里,免得搞掉了。 阅读全文

posted @ 2024-03-11 18:33 luyi07 阅读(52) 评论(0) 推荐(0) 编辑

python sympy模块计算 Clebsch-Gorden 系数
摘要:(未经允许,不得转载) 角动量代数还是挺常用的,今天本想自己写一个 python 的 CG 系数函数,又觉得可能已经有现成的包了。 google 搜了一下,第一条就是 sympy 官网的相关文档,快准爽。学习了一下官网文档,简单记在这。 1. 安装 sympy 我用的 pycharm,所以按照常规的 阅读全文

posted @ 2022-09-29 11:26 luyi07 阅读(407) 评论(0) 推荐(0) 编辑

贝塔衰变 phase-space factor
摘要:All these are probably well known, I just take notes from these references below, and figure out some derivations. References: Christian Illiadis,《Nuc 阅读全文

posted @ 2022-08-13 18:38 luyi07 阅读(204) 评论(0) 推荐(0) 编辑

单极对力下 BCS 方程求解
摘要:BCS方程可用于描述单满壳原子核的基态,这个笔记总结BCS方程的构成和特征,并展示BCS方程数值解的算法,讨论算法中的收敛性问题。 最近实验学家测量了136Xe的价质子在不同轨道上的占据数,所以我借此做了实例演示,在唯象的单极对力哈密顿量下解BCS方程,与最新的实验数据相比较。 阅读全文

posted @ 2022-08-12 16:24 luyi07 阅读(155) 评论(0) 推荐(0) 编辑

Lanczos sum rule
摘要:假设跃迁算符是 ˆO,如果是 E2,ˆO 的第三分量可以是 2,1,0,1,2。 现在对于给定的初态 |i,想要求它的某种 sum rule: $$ I = \int^\infty_0 d \omega S(\omega) g(\omega). 阅读全文

posted @ 2022-07-29 09:33 luyi07 阅读(42) 评论(0) 推荐(0) 编辑

General Seniority 学习笔记(1): BCS
摘要:General Seniority 学习笔记(1): BCS 我学习了参考文献 [1,2],把里面的核心公式推导核对整理了,因为觉得有点意思。做完笔记我顺手写了个代码,还没来得及核对。 下一步可以考虑投影 broken-pair 的优化和投影,即 seniority 取次极小的内秉态做投影。 参考文 阅读全文

posted @ 2022-06-16 18:30 luyi07 阅读(79) 评论(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) 编辑

宇称投影
摘要:1. 宇称投影算符 宇称投影算符ˆP作用在任意函数f(x)得到 ˆPf(x)=f(x), 三维空间中则为 ˆPf(r)=f(r). 如果有 $\hat{P} f(\vec{r}) = f(-\vec{r 阅读全文

posted @ 2021-10-26 12:06 luyi07 阅读(353) 评论(0) 推荐(0) 编辑

原子核壳模型程序 BigStick 脚本:sd壳能谱+B(E2)
摘要:以下脚本使用可执行程序: bigstick.x, tropic1b.x, genstrength.x 使用 sd 壳相关文件: sd.sps, usdb.int 结果文件: pxnx.res, pxnx.e2.str 使用了 awk 语句:awk 'BEGIN{print "'$A'"^0.1666 阅读全文

posted @ 2021-10-14 20:28 luyi07 阅读(105) 评论(0) 推荐(0) 编辑

格式转换:有效相互作用 -> P+Q相互作用
摘要:这是以前写的笔记,散在不同的地方,现在收集在这个随笔里,供以后自己和合作者翻阅。写了这个笔记以后,我还会在相应代码里加上注释,因为没有人能记住大串公式的每一个细节,所以只能把笔记、博客用作外部的记忆存储器。 1. 壳模型相互作用: pn格式 在 pn 格式下,单体 + 2 体相互作用写作 \[ \b 阅读全文

posted @ 2021-07-01 15:25 luyi07 阅读(380) 评论(0) 推荐(0) 编辑

原子核壳模型程序 BigStick 的用法
摘要:BigStick 公开版本可在 https://github.com/cwjsdsu/BigstickPublick.git 下载。其中 docs/ 与 util/docs/ 下都有详细的文档说明,这里只是记下简单运用的步骤。 下面的脚本是根据 util/example/ 里面的脚本修改的,可以自动 阅读全文

posted @ 2021-07-01 14:45 luyi07 阅读(172) 评论(0) 推荐(1) 编辑

c++ 中的下三角阵矩阵元标记
摘要:若有一个下三角阵,不包括对角元,其矩阵元为 Aij,    i=1,,n1,    j=0,,i1. 那么,可以用一维数组储存这些矩阵元: a[k]=Aij,   k=i(i1)/2+j. 这样可以将 阅读全文

posted @ 2021-06-29 19:32 luyi07 阅读(186) 评论(0) 推荐(0) 编辑

BCS方程和Bogoliubov变换
摘要:BCS 方程和 Bogoliubov 是处理对关联的经典内容。一半是学习,一半是为了教学,这里我写下了较为完整的推导过程。我欣赏 Bogoliubov 变换,它把 BCS 波函数变为准粒子真空,在此基础上,构造的准粒子激发可以看做是准粒子的自由气体,这很漂亮,虽然是在粒子数不守恒的 BCS 框架下, 阅读全文

posted @ 2021-06-28 17:36 luyi07 阅读(1916) 评论(0) 推荐(0) 编辑

核结构单体跃迁算符
摘要:这都是教科书上都有的内容,我只是整理整理,把一些约定也统一一下,方便以后写代码的时候参照。 1. 约定 1.1 约化矩阵约定 \[ \langle J'M' | s \sigma | J M \rangle = \langle J' || s || J \rangle (JM,s \sigma| J 阅读全文

posted @ 2021-06-23 12:26 luyi07 阅读(473) 评论(0) 推荐(1) 编辑

投影模型中的单体跃迁
摘要:1. 投影框架下的本征态 如果变分得到的能量最低波函数为 |PC,这里 PC 是 pair condensate 的缩写,表示我们最近做的对凝聚组态,但下面的讨论并不局限于这一种组态。 如果要得到角动量为 (J,M) 的近似本征态,我们先从 |PC 阅读全文

posted @ 2021-06-10 15:45 luyi07 阅读(258) 评论(0) 推荐(0) 编辑

投影技术之线性代数投影:LAP
摘要:之前写了投影技术的理论基础,即转动算符、D矩阵、投影算符等理论知识:https://www.cnblogs.com/luyi07/p/14586631.html 又写了两种投影技术的思路与比较:https://www.cnblogs.com/luyi07/p/14823838.html 在这个随笔中 阅读全文

posted @ 2021-06-07 20:32 luyi07 阅读(555) 评论(0) 推荐(0) 编辑

核结构中的投影技术
摘要:1. 回顾:角动量投影 之前写过一篇随笔,整理了角动量投影技术的推导:https://www.cnblogs.com/luyi07/p/14586631.html 投影算符定义为 \[ \hat{P}^J_{LK} | \Phi \rangle = \frac{ 2 J + 1 }{ 8 \pi^2 阅读全文

posted @ 2021-05-28 21:14 luyi07 阅读(251) 评论(0) 推荐(0) 编辑

gnuplot 绘制球谐函数图
摘要:因为对形变感兴趣,最近在看 P. Ring & P. Schuck 的书。在形变很小(β<<1)时,椭球体的表面可以用低阶球谐函数来近似描述。所以,不妨用gnuplot画两个图出来看看,比较直观,也许以后展示的时候也用得着。 1. 球谐函数 使用 Griffiths 的量子力学教 阅读全文

posted @ 2021-04-28 11:03 luyi07 阅读(1088) 评论(0) 推荐(0) 编辑

非正交归一基:Hill-Wheeler 方程
摘要:这个很简单,但是还是简单记一下,供以后翻阅。 1. 正交化 现有线性无关的非正交归一基矢 |i,i=1,2,,n.把这些基矢做内积,得到矩阵 Aoveij=i|j,然后对角化,得到 \[ 阅读全文

posted @ 2021-04-26 15:21 luyi07 阅读(819) 评论(0) 推荐(0) 编辑

导航

< 2025年2月 >
26 27 28 29 30 31 1
2 3 4 5 6 7 8
9 10 11 12 13 14 15
16 17 18 19 20 21 22
23 24 25 26 27 28 1
2 3 4 5 6 7 8
点击右上角即可分享
微信分享提示