Phoenix in fire
从上周三开完组会一直到现在都还没怎么缓过神来，不能再这样下去了😟，得主动出击。早上起来忧心忡忡的，高分子物理没心情上，吃完早饭直接去鼎新馆吧。昨晚大概M$入门了，建模后提交了任务，但是可能INCAR里面参数设置得不够合理，今天下午提交的任务才收敛。尝试想了一些方案但是被师兄否定了，提高电解质中离子电导率的机理也许是清楚的，但是没有一篇明晰的文献手把手教我应该挑选那些元素掺杂来提高离子电导率，也许永远不会吧🥺。Keep on going，再试试吧，多读读文献，也许我现在的困境和英语只有60分的同学所面临的困境相似吧？——都是词汇量不够？ Anyway, Just keep on going! 如果走不出迷堆明天就去请教一下魏老师。 合抱之木，生于毫末；九层之台，起于累土；千里之行，始于足下
阳光总在风雨后
¶歌词 123456789101112131415161718192021222324252627282930313233343536373839404142434445[00:00.60]许美静 - 阳光总在风雨后[00:01.60]词：陈佳明[00:02.60]曲：陈佳明[00:24.34]人生路上甜苦和喜忧[00:29.24]愿与你分担所有[00:34.15]难免曾经跌倒和等候[00:39.30]要勇敢的抬头[00:43.95]谁愿藏躲在避风的港口[00:48.85]宁有波涛汹涌的自由[00:53.91]愿是你心中灯塔的守候[00:59.22]在迷雾中让你看透[01:06.93]阳光总在风雨后[01:11.89]乌云上有晴空[01:16.84]珍惜所有的感动[01:21.89]每一份希望在你手中[01:26.89]阳光总在风雨后[01:32.00]请相信有彩虹[01:37.05]风风雨雨都接受[01:41.81]我一直会在你的左右[02:11.74]人生路上甜苦和喜忧[02:16.98]愿与你分担所有[02:21.35]难免曾经跌倒和等候[02:26.71]要勇敢地抬头[0...
快乐无罪
¶歌词 123456789101112131415161718192021222324252627282930313233343536373839404142[00:00.43]快乐无罪 - 许美静[00:01.62]词：许常德/陈佳明[00:03.13]曲：陈建宁/陈政卿[00:19.82]你要我快乐[00:23.76]我非常快乐[00:27.22]在不想太多的时候[00:34.32]吻你是烟火 吻后是寂寞[00:41.32]爱是我犯的差错[00:49.38]我要你什么 我贪你什么[00:56.73]在你想消失的时候[01:04.12]想你是黑洞 它没有尽头[01:10.72]就算流星也该坠落[01:17.55]快乐有泪 怎会有泪[01:21.23]我哭到心碎[01:25.52]不解幸福会被永恒连累[01:29.16]不解刹那是你要的美[01:32.60]快乐有罪 怎会有罪[01:36.29]我怎么拒绝[01:40.55]也许到最后你我都反悔[01:44.21]留下空荡的城市黑夜[01:47.60]让我撒野[02:22.36]我要你什么 我贪你什么[02:29.56]在你想消失...
截断能ENCUT的选取
ENCUT是用来指定平面波基组能量截止值的参数 以上是wiki对ENCUT的描述，POTCAR文件中已经包含了元素的ENMAX和ENMIN，这也是互联网上对ENCUT的选择方案对通常是1.3倍的ENMAX 侯柱峰老师提供了另外一种方法并指出通过该方法选择的ENCUT通常能满足1.3倍ENMAX，供参考使用 截断能指定了用于波函数展开的平面波基组的截断能量，此能量越大则用来描述波函数的平面波基组越多，精度越高，但计算也越耗时。 可以通过计算测试来选择合适的ENCUT。 以下是用于测试的脚本： 123456789101112131415#!/bin/shrm WAVECARfor i in 150 200 250 300 350 400docat > INCAR <<!SYSTEM = Si-DiamondENCUT = $iISTART = 0 ; ICHARG = 2ISMEAR = -5PREC = Accurate!echo "ENCUT = $i eV" ; time vaspE=`grep "TOTEN" OU...
表面能的计算
¶关于表面能 根据VASP官方手册Hands on Session III，表面能的定义如下，大师兄的文章也提到了这个定义( Nanoscale, 2017, 9, 13089–13094) #mjx-bb96b1c8{ display:contents; mjx-assistive-mml { user-select: text !important; clip: auto !important; color: rgba(0,0,0,0); } mjx-container[jax="SVG"] { direction: ltr; } mjx-container[jax="SVG"] > svg { overflow: visible; min-height: 1px; min-width: 1px; } mjx-container[jax="SVG"] > svg a { fill: blue...
POTCAR赝势的选取
本文为一个查阅性的记录文章，也可以说是搬运，内容取自侯柱锋老师的使用VASP的个人经验手册 元素 普通赝势 (eV) 高硬赝势 (eV) 软赝势 (eV) d 电子作为半芯态来处理赝势 (eV) B 318 700 (B_h) 250 (B_s) - C 400 700 (C_h) 273 (C_s) - N 400 700 (N_h) 250 (N_s) - O 400 700 (O_h) 250 (O_s) - F 400 700 (F_h) 250 (F_s) - Al 240 295 (Al_h) - - Si 245 380 (Si_h) - - P 270 390 (P_h) - - S 280 402 (S_h) - - Cl 280 409 (Cl_h) - - Ga 134 404 (Ga_h) - 282 (Ga_d) Ge 173 410 (Ge_h) - 287 (Ge_d) As 208 - - - Se 211 - - - Br 216 - - - In 95 - - 239 (...
以p(1x1)Cu(111)为例分析表面弛豫的结果
本文其实没啥新的内容，只是对上两节的一个总结分析 在表面弛豫前（slab模型优化计算前） 12345678910111213SYSTEM = p(1x1)Cu(111)1.0 2.5701000690 0.0000000000 0.0000000000 -1.2850500345 2.2257719501 0.0000000000 0.0000000000 0.0000000000 21.2954006195 Cu 4Selective DynamicsDirect 0.000000000 0.000000000 0.000000000 F F F 0.666670000 0.333330000 0.098540000 F F F 0.333330000 0.666670000 0.197080000 T T T ...
p(1x1)Cu(111)的优化计算
¶POSCAR的修改 与单点计算不同，优化计算可能涉及到离子步迭代、晶格常数的改变。但Cu(111)晶面的优化计算并不只是在INCAR里添加参数IBRION，优化计算确实允许原子弛豫，但在slab模型中，只允许表面的两层原子允许弛豫，而下面的那些原子则直接固定住 这样做的物理意义在于：在真实的环境下，催化剂体相看做是不变的，只有表面的原子参与催化反应 前面在乙醇分子的振动频率分析里面已经学过了，添加以下选择性动力学标记Slective Dynamics即可 PS：看分数坐标即可知道哪些是表层哪些是体相 12345678910111213CONTCAR\(1\1\1)1.0 2.5701000690 0.0000000000 0.0000000000 -1.2850500345 2.2257719501 0.0000000000 0.0000000000 0.0000000000 21.2954006195 Cu 4Selective ...
Cu(111) slab模型的单点计算
首先对Cu单胞进行结构优化（bulk计算），然后使用优化的结构进行切割，再进行计算（由于晶胞已经优化过，所以本次计算其实是单点计算）。大师兄提到： Cu(111) slab 模型（POSCAR）的制作流程（复习一下）。**注意：**前面我们用的是Conventional cell，下面用的是Primitive Cell。 一般来说用Conventional cell， FCC的金属可以用Primitive Cell 但对于其他体系，通过Primitive Cell切出来表面模型有问题 我准备两种都验证一下。 ¶使用Conventional Cell建立slab模型 ¶切割结构 将bulk计算完毕的CONTCAR文件用Vesta打开，再转换为.cif文件使用Materials Studio打开 Step1：Build-Surface-Cleave Surface Angstrom是单位埃，需要切割的晶面是(111)，相应的平面方程为 #mjx-19b68b5{ display:contents; mjx-a...
Cu单胞的bulk计算
bulk计算，直译为“体相计算”，指的是模拟理想的三维无限大晶体（无表面、界面或缺陷，不模拟则体系太大算不动），利用 周期性边界条件 描述材料的体相性质，主要用于研究材料的本征物理化学性质（如晶格结构、电子能带、热力学稳定性等） ¶准备POSCAR 打开materials studio-File-Import-Structures-Metals-pure metals-Cu.xsd，然后export为cif文件 用vesta打开Cu.cif文件，再export data，导出保存为POSCAR文件 ¶准备KPOINTS 12345KPOINTS for Cu-concentional-cell-bulk-calc0Gamma13 13 130 0 0 ¶准备INCAR 123456789101112SYSTEM = Cu-conventional-cell-bulk-calculationISMEAR = 0SIGMA = 0.1ENCUT = 700ISIF = 3NSW = 200NELM = 200IBRION = 2POTIM = 0.1EDIFFG = -0.01ALG...
