add content

content/blog/my-2024.md

@@ -0,0 +1,66 @@
+---
+title: '回顾我的2024年'
+date: 2025-01-20T15:26:05+08:00
+draft: true
+summary: 弄坏了，再修好
+---
+
+2024 年对我来说是非常非常非常漫长的一年，这一年经历了很多的事情，使得我的直观感受是过去了两三年。
+回忆年初时发生的事情时，我需要反复确认才能肯定那的确是发生在仅仅一年前的事情。
+
+<br/>
+
+---
+
+我想从 2023 年的秋天或者冬天开始讲，因为我觉得“命运的齿轮”是从那天开始转动的；
+  如果要拍动画或者电影的话，或者写小说的话，从这里开始是最合适的。
+厦门的四季并不分明，所以我很难说那时是秋天还是冬天；
+  那天突然变冷了，那天之前是穿短袖感到凉爽的天气，那天开始就需要穿上外套了。
+故事就从那天开始讲。
+
+那天晚上，我的大老板的儿子给我打电话说，我小老板去世了。
+我其实不惊讶，大概半个月前我在医院看到他的时候，他的状态就已经很不好了；
+  那天突然降温，厦门也没有暖气，正是容易出事的时候。
+
+再往后没过几天，我们几个同学去参加了他的葬礼。
+这是我第一次参加熟悉的人的葬礼。
+之前有过熟悉的人去世，也参加过不熟悉的人的葬礼，这样的情况是第一次。
+葬礼的前大半段并没有什么特别的，我也没有落泪；
+  在场的大多数人我并不认识，闽南的风俗我也不熟悉，我大多处在一种等待的状态，有时帮忙摆摆板凳发瓶水。
+最后，我随着人群排队去透明棺材前看最后一眼的时候，突然被吓了一跳。
+棺材里的人脸色惨白，但嘴唇又被画成红丹丹的颜色，
+  就像是已死的人因为活人的风俗又被拉起来画上了浓妆，马上还要被驾着去参加活人的宴席那样。
+葬礼结束之后我们几个同学开车去吃饭的途中，一边聊天，我才感到悲伤，哭了出来。[^1]
+
+我和小老板身边的大多数人并不熟悉，和他的家人也只是一两面之缘。
+对我自己而言，这件事标志着我的研究生以及人生的某个阶段的彻底结束。
+我大四（大概是 2019 年）犹豫读研还是工作的时候，在校园里碰到了我后来的小老板，
+  他吹牛说组里有最先进的仪器（指八寸的碳化硅外延炉，实际上几年后才调试到基本能用），
+  夸我说我很会编程，问我要不要来他这里读研。
+我说好呀好呀，然后就来了，并且逐渐和组里的人熟悉起来。
+后来，我唯一的学姐毕业去远处工作了，唯一的学长跑去另外一个校区读博了；
+  曾经带我去万象城的、唯一的同级同学因为我搞不清的原因（绯闻，或者是抑郁之类的心理疾病）去了别的组也基本不再联系[^2]；
+  小老板得了癌症，之后经历了几轮住院和基本康复。
+在我小老板无法工作之后，（前文提到的）我的大老板和他的儿子接手了带我；
+  但对我来说，在我小老板手下工作和在他们手下工作是完全不同的。
+
+我熟悉的那个世界逐渐快速地萎缩和崩塌，最后完全消失，只留下了我一个人。
+
+<br/>
+
+---
+
+我知道上一个故事的氛围有些沉重。
+这个故事充满了理想主义的粉色泡泡，应该可以在情绪上缓冲一些。
+
+这个故事发生在 2024 年的元旦之后、春节之前。
+那时我的论文刚投稿，放了寒假，我自然
+
+
+[^1]: 葬礼给我的感觉是，我的小老板是一个符号而不是主角，葬礼是活人们以这个符号为借口聚集在一起办了一件活人想办的事情，而根本与他没有关系。
+  直到最后看到他躺在棺材里，以及和同学在车上聊天，才感到确实是他。
+[^2]: “万象城”是厦门的一个繁华的商业中心。
+  她是厦门本地人，我是外来、又很少出门逛商场的宅人，她把我带去万象城的天台时我很震撼，之后每次想到她都会想到那个场景。
+  2023 年夏天我们一起写人生的第一篇论文的过程中，不知道什么原因，她逐渐和我们（我以及我的老师们）疏远了。
+  小老板的葬礼之后我只见到过她一两次，大多时候，不论是小组里聚餐还是我单独请她吃饭，她都表示有事不能来。
+

content/blog/smear-in-dft.md

@@ -0,0 +1,85 @@
+---
+title: 'DFT中的涂抹'
+date: 2025-01-17T12:36:30+08:00
+draft: true
+summary: ISMEAR应该如何设置？
+---
+
+VASP 里有两个神秘的参数：`ISMEAR` 和 `SIGMA`。
+初次接触 VASP 的时候我曾经注意到过这个参数，但因为看不懂所以没有深究，然后就把这个事情忘记了。
+前几天又有人问我这个问题，我就注意到了这个盲点，于是来研究了一下。
+
+这两个参数控制的是 KS 轨道上在费米能级附近电子如何分布：
+  尽管按照严格的 DFT 理论，应该是费米能级以下的轨道都被占满、费米能级以上的轨道都空，
+  但实际计算的时候，有时会做一些平滑过渡，
+  比如说费米能级下面的第一个轨道放 $1.9$ 个电子、费米能级上面的第一个轨道放 $0.1$ 个电子。
+为什么要这样做？乃至于，如果你不了解 DFT 的原理，你可能会问：KS 轨道是什么玩意儿？
+这个文章就是要解释这个问题。
+
+需要先说明的一点是：这两个参数是不“物理”的，
+  也就是说，按照物理图景来思考的话，就不应该有这个平滑过渡，`SIGMA` 应该永远是 $0$ 才对[^1]。
+设计这个参数只是为了缓解计算细节中的问题，所以之后的解释也都是从计算的角度来解释的。
+
+## 一点科普：VASP 计算的原理
+
+用 DFT 做计算的人并不见得理解 DFT 的原理（我最开始也是这样，这并不罕见）。
+这里简单介绍一下原理，之后才能进一步解释为何要涂抹，如何涂抹更合适。
+这一节的内容需要基础的量子力学知识，基础的就够了。
+
+VASP（乃至所有基于 DFT 的第一性原理计算软件）核心目的都是求解关于电子的薛定谔方程：
+
+$$
+  \hat{H} \Psi(\mathbf{r_1}, \mathbf{r_2}, \ldots, \mathbf{r_N})
+    = E \Psi(\mathbf{r_1}, \mathbf{r_2}, \ldots, \mathbf{r_N})
+$$
+
+这里 $E$ 是总能量，$\hat{H}$ 是哈密顿算符。
+$\Psi$ 是所有电子（不是单个电子！）的波函数，$\mathbf{r_i}$ 是第 $i$ 个电子的位置。
+也就是说，如果体系有一千个电子的话（这在基于平面波基矢的 DFT 计算中是常见的），这个波函数就是一个三千个维度的复变函数。
+这么复杂的函数现在的计算机是没法处理的，稍等会需要想办法来绕过这个问题（实际上求解的是另外一个与它等价的方程）；
+  但形式上来说，这就是我们想要解的方程（解不解得出来另说），所有的计算结果都是求解这个方程得到波函数之后，再得到的副产品。
+例如，所谓“弛豫得到体系最稳定的结构”就是知道波函数之后，再计算一下各个原子核的受力，然后把原子核顺着受力方向移动一点，
+  再重新求解这个方程，如此循环往复，直到所有原子核的受力都小到可以忽略不计为止。
+
+上述哈密顿算符 $\hat{H}$ 至少得包含三项：电子动能、电子和原子核之间的库伦吸引、电子之间的库伦排斥。也就是：
+$$
+  \hat{H} =
+    \sum_{i} \frac{-\hbar^2}{2m} \nabla_i^2
+    - \sum_{i, I} \frac{Z_I e^2}{4\pi\varepsilon_0\left|\mathbf{r}_i - \mathbf{R}_I\right|} \
+    + \sum\_{i\neq j} \frac{e^2}{8\pi\varepsilon_0\left|\mathbf{r}_i - \mathbf{r}_j\right|}
+$$
+当然你也可以把更多的东西包含进来（比如原子核之间的斥力，磁场与自旋的耦合，等等），我们对此不做过多讨论，
+  因为这和之后要说的事情没有关系，我们就假定讨论的体系是这样简单的体系。
+这里写出这个表达式来只是为了告诉读者：
+  到现在为止，哈密顿量和薛定谔方程与量子力学课堂上学的并没有太大不同；但很快我们就要把它改造得面目全非了，不这样做的话没法写程序计算。
+
+这个改造的起点是这样一个结论：
+  **基态电荷密度蕴含了体系所有的信息**，包括但不限于基态能量、各个激发态的能量，
+  各种状态下的波函数、原子核的位置和所带电荷，以及你今天早上吃的豆浆还是油条。
+总之，它就是字面意思上的“所有”。
+如果你要搜索这个结论的证明的话，大概可以看到一个非常古老、非常简洁但又非常“耐人寻味”的“证明”[^2]。
+但这个证明并不能说服我，我总觉得那就是个循环论证；我感到信服的证明是在某本 DFT 教学的书上看到的[^3]，它的证明思路是这样的：
+* 薛定谔方程（哈密顿算符）蕴含了体系的一切信息，这是量子力学的基本假设；
+* 对于上述体系，只需要三个信息就可以写出哈密顿算符的表达式：各个原子核的位置、各个原子核所带电荷、电子的个数；
+* 从基态电荷密度出发，这三个信息都是求出来的：
+  * 电子个数：在全空间中积分一下就得到了。
+  * 原子核位置：解析的表达式在书中有给出，没记错的话，大概就是，原子核存在于电荷密度的倒数不连续的地方。
+  * 原子核所带电荷：解析的表达式在书中也有给出。
+* 因此，基态电荷密度已经蕴含了体系的所有信息。
+
+这个结论并没有马上指出我们要如何写代码做计算，而是给出了一个可能：
+  虽然波函数复杂得可怕，但我们可能并不需要真的求解出它；
+  我们只要找到绕过波函数求解出基态电荷密度（这只是一个三维的实数函数）的方法，同样可以得到所有想要的信息（基态能量，能带结构，等），
+    整个过程或许不需要考虑那个三千维的波函数。
+事实上也确实是这样做的。这就是“密度泛函理论”（density functional theory，DFT）名字的由来。
+
+这个方法是这样做的：将上述哈密顿量“略微”修改，使得它描述的不再是实际的体系，而是另外一个假象的体系；
+  其中的粒子也不是电子，而是另外一种假想的粒子。
+这个假象的体系有这样两个特点：其一是，它的**基态电荷密度与实际体系相同**，这样求解出它的基态电荷密度就等于求解出了实际体系的基态电荷密度；
+  其二是，**假想粒子之间没有相互作用**，这样假想体系中的
+
+
+[^1]: 尽管电子也应该服从费米-狄拉克分布，但一来 DFT 本身是在基态（绝对零度）下才严格成立的理论，有限温度下的计算结果是没有什么理论来支撑的；
+  二来实际生活中的温度（最多几千开尔文）对电子状态影响其实不大，它主要影响的是原子核的振动状态。
+[^2]: 好像是叫什么什么第一定理第二定理之类的，我记不清了。
+[^3]: 英文的，书名我忘记了，有空再找找补充进来。