应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » 其他 » 【其他】征集本版成功求助贴
251/3返回列表123下一页
查看: 1805  |  回复: 89
查看全部回帖@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

wuli8

荣誉版主 (知名作家)

…………

优秀版主优秀版主

[交流] 【其他】征集本版成功求助贴

【活动目的】
经常光临“第一性原理版”的虫子都很清楚，本版是以求助贴为主的。有的求助贴，求助成功了，但是有的没有。这样就很混乱，时常有2次求助的现象发生。这样很不利于大家交流。为了解决这一问题特向各位征集本版成功求助贴。
【跟帖方式】
[求助贴链接]（必填项）
[求助的问题]（必填项）
[求助问题的归类]（必填项）
[问题的正确答案]（必填项）
[你认为正确的答案]（可选项）
【注意事项】
1、发帖前请查看有没有别人已经发过您的内容，重复贴无效。
2、问题的答案要正确无误，如果原帖没有正确答案，你可以进行补充。明显错误的答案，为无效贴。
【相关奖励】
每一个有效的征集奖励5个金币，特别优秀者奖励Vip0.1个。



PS:帖子前缀为【求助成功】的表示别人已经在此处发过帖子，请不要重复发帖。谢谢大家的参与！

[ Last edited by wuli8 on 2009-6-25 at 17:24 ]
回复此楼

» 猜你喜欢
…………
1楼2009-06-02 10:04:47
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

mozhui

★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
wuli8(金币+5,VIP+0):感谢参与 6-4 11:52
http://muchong.com/bbs/viewthread.php?tid=1329674
其他
静电相互作用能和电子的交换-关联能有什么区别？本人初学，两者我觉的都是库伦相互作用能，不知道是从哪个角度来区分的。谢谢

2楼：
密度函数理论中将实际的n个电子之间的库仑相互作用分解为两个部分：
1.电荷分布为￡的空间连续电荷的相互作用能
2.真实的电子之间的库仑作用与1中的能量差
其中第二部分即为通常所说的交换关联能，在用独立电子体系近似时交换关联能通常还包括了用独立电子系的总动能替换实际总动能所导致的偏差~
这么做的原因是一般情况下真实电子之间的库仑相互作用并不能写成电子密度的显函数
3楼：
交换能和关联能分开似乎是在hartree-fock单电子近似时才用的，在密度泛函理论中都是合在一起定义的
在hartree-fock近似下：
交换能是指在由电子的交换反对称性（此时的体系波函数为单电子波函数的slater行列式）引起的相对不考虑交换反对称性（此时波函数为单电子波函数的简单乘积）时的能量差
关联能则指由于一般情况下实际多体系统波函数不能分解为slater行列式，而我们采用slater行列式时所引入的偏差

都是本人回的，感觉差不多吧~
赞一下(11人) 回复此楼
38楼2009-06-04 11:48:08
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

mozhui

★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
wuli8(金币+5,VIP+0):谢谢回帖 6-4 21:30
http://muchong.com/bbs/viewthread.php?tid=1365608&fpage=1
其他
【求助】第一性是啥东东？
2楼：DFT密度泛函
3楼：上帝先创造男人，所以第一性是男人:-(

不过wuli8讲的不很全，HF也不是DFT，但也算是第一原理。第一原理就是说从基本的物理原理比如薛定谔方程出发，不加任何经验参数而算出体系的性质。
换句话说，分子动力学里面的势能一般用经验势，就不是第一原理了。
但是因为第一原理能做的事情有限，比如强关联体系就很难处理，所以会加一些经验的东西来拓展处理能力和应用领域。打个比方说，第一原理是男人，那些经验的东西就是伟哥。
5楼：
据我所知，第一性原理好像包括两个基本方程，一个是薛定谔方程，另一个是波松方程。
按照量子化学，传统第一性原理包括HF,POST-HF,MPn,CCSD,CASSCF等等，这是传统的量子化学方法，其中有些精度甚至高于DFT，比如CCSD(T)。
7楼：第一原理方法，就是不依赖于任何经验常数的方法。第一原理又叫第一性原理，指的是一些公理性的原理，以及一些硬性规定的法则或假设等。
量子化学的第一原理就是Schrödinger方程，以及量子力学的几个基本假设，除此之外，只需要几个基本的物理常数， 如电子质量、核电荷数（原子序号）、普郎克常数、玻耳兹曼常数、真空光速等，即可开展计算，演绎出物质的结构和各种物理性质。而具体的计算方法则有很多，如五楼所说。当然为了实现这些方法，也不可避免地有一些近似和假设。

大家说得都很好，综合起来就差不多了~
赞一下(11人) 回复此楼
39楼2009-06-04 13:24:26
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

mozhui

★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
wuli8(金币+5,VIP+0): 6-4 21:29
http://muchong.com/bbs/viewthread.php?tid=1359976&fpage=2

Vasp&MedeA

【求助】vasp作动力学用哪个系综？我想研究团簇在室温下的稳定性用哪个系综？要设那些关键词？ 谢谢 我用的SMASS=0 TEBEG=300 TEEND=300，可发现OSZICAR里面的温度是变化的，是什么原因阿？ 请教  
3楼：NVT
5楼：IBRION=0  
POTIM
SMASS
TEBEG=
TEEND=
6楼：SMASS 的值决定你的系综
8楼：温度波动是正常的，那是nose hoover热浴法的特点， 不过你要保证你的温度波动幅度较小，
赞一下(11人) 回复此楼
40楼2009-06-04 13:30:42
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

mozhui

★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
wuli8(金币+5,VIP+0): 6-4 21:29
http://muchong.com/bbs/viewthread.php?tid=1338535&fpage=1

Vasp&MedeA

【求助】VASP里面的费米能级表示什么？

OUTCAR里面有一个参数，用grep “fermi" OUTCAR可以得到，比如：

E-fermi :   1.0372     XC(G=0):  -6.7744     alpha+bet : -8.3695

求教：这个Fermi Energy表示的是
1.价带顶的位置 还是
2.Fermi-Dirac分布中占据能级的几率是1/2的位置？

谢谢～

3楼： 电子正好填充到该能量尚
4楼： 我觉得是前者吧
电子的最高占据轨道
6楼： 同意前者。
7楼： 电子正好填充到该能量尚
8楼： 0温下还有什么1/2占据？把T=0代近fermi-dirac公式，整个就一个阶跃函数好吧！自己看公式去。。

另外，9楼介绍了相互作用费米子的分布情况，赞一个，不过与本问题关系不大~
VASP总是取电子的最高占据态为费米能级，因而对于半导体，VASP所取的费米能级为价带顶，与通常的定义不同，这点要小心~
赞一下(11人) 回复此楼
41楼2009-06-04 15:29:01
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

mozhui

★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
wuli8(金币+5,VIP+0):感谢参与 6-4 21:29
http://muchong.com/bbs/viewthread.php?tid=1347934&fpage=2

Vasp&MedeA

【求助】求助，CHGCAR的单位?
2楼：e/A3
3楼：这个问题，我自己已经搞定了。 应该是e／A3 在乘以fft的点数的。
赞一下(11人) 回复此楼
42楼2009-06-04 15:38:16
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

mozhui

★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
wuli8(金币+5,VIP+0):感谢参与 6-4 21:28
http://muchong.com/bbs/viewthread.php?tid=1335547&fpage=3

Vasp&MedeA

【求助】求教如何画态密度图？

搜了一下，看到一般的步骤为：
第一步，就是对研究体系进行驰豫，得到稳定的结构和体积。
第二步，就是进行静态计算，获得电子密度等
第三步，非自恰计算获得DOSCAR

我已经对体系弛豫完毕。但是不知道第2、3步中INCAR里具体的那些参数需要调整。

2楼：
second step:
SYSTEM =XXX
PREC = Accurate
ISTART=0
ENCUT = 1000.0    ＃在作结构优化时，需要收敛测试。用那个就行。
EDIFF = 1e-5
ISMEAR =1          # 1 metal ; 0 non-metal
SIGMA = 0.2
LCHARG = T
LWAVE = T

POSCAR KPOINTS POTCAR 都用结构优化时的。POSCAR 用你想做的那个压力点的CONTCAR。
third step：
SYSTEM = XXX
PREC = Accurate
ISTART=1
ENCUT = 1000.0
EDIFF = 1e-5
ICHARG=11
ISMEAR=-5
SIGMA = 0.2
LORBIT=10    # 计算s p d 时用10 想算用px py 11。
第二步计算是在结构优化的结果上进行的，所以开始第二步的时候，将第一步中的输入文件INCAR, POTCAR, POSCAR, KPOINT 以及 C* 文件放入静态自洽计算中去，并且将CONTCAR 拷贝到 POSCAR中，然后运行VASP。计算结果中的Fermi能是准确的，需要记下。
得到的DOSCAR ， 就是你想要的了。
4楼：
1.结构优化，即普通计算
    2.将优化之后产生的CONTCAR文件改成POSCAR, 做静态计算.输入文件包括KPOINT,POTCAR,POSCAR,CHGCAR,INCAR;INCAR中ICHARG=2,ISTART=0,ISMEAR=-5,LCHARG = .TURE.;
    3.将静态计算产生的CHARGCAR文件拷贝到输入目录中,INCAR,POSCAR,POTCAR,KPOINTS和CHARGCAR, INCAR需做如下改动, ICHARG=11,ISMEAR=-5,LORBIT=11;并增加K点.
赞一下(11人) 回复此楼
43楼2009-06-04 15:45:02
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

mozhui

★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
wuli8(金币+5,VIP+0):感谢，参与 6-5 20:51
http://muchong.com/bbs/viewthread.php?tid=1341003&fpage=2

Vasp&MedeA

【求助】一个VASP安装文件有多大？

弱弱地问一下！
一个vasp安装盘中的文件有多大？
2楼：弱弱的说一下,安装文件不到1.5M
4楼：加上势文件不超过100M
5楼：人家卖给你的是源代码，可不是打包好的商业软件。虽然只有1.5M，但软件售价又不是和文件大小相关的。再说VASP售价也算厚道了
6楼：还是1.5M
9楼：势函数是单独下载的.
安装文件就1.5M
赞一下(11人) 回复此楼
45楼2009-06-04 21:46:55
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

mozhui

★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
wuli8(金币+5,VIP+0): 6-5 20:48
http://muchong.com/bbs/viewthread.php?tid=1311400&fpage=3

Vasp&MedeA

求助】关于VASP测试的命令输入

例如在做截断能测试时 写好脚本 按照下面提示在终端输入命令
==================
（二）计算任务执行方法：
输入：dos2unix optencut.sh
bash optencut.sh
==================
结果就出错了
没办法继续算下去

这是为什么呢？？？那该怎么写入呢？？？
对FeISMEAR能设置成-5吗？

问题脚本：
========================
（一）输入文件：
1、用脚本程序optencut.sh代替INCAR文件：
rm WAVECAR
for i in 100 150 200 250 300 350 400 450 500 550 600
do
cat > INCAR < SYSTEM = bcc Fe
ENCUT = $i
ISTART = 0 ; ICHARG = 2
ISMEAR = -5
PREC = Accurate
!
echo "ENCUT = $i eV"; time vasp
E=$(grep "TOTEN" OUTCAR | tail -1| awk '{printf "%12.6f \n", $5}')
echo $i $E >> comment
done
===========================

6楼：要将 optencut.sh 文件属性改为可执行文件，用命令：  chmod +x optencut.sh

下一步输入命令：  ./optencut.sh   就可以运行了.
8楼：金属的话一般用１和２　但做静态计算的时候设成-５
赞一下(11人) 回复此楼
46楼2009-06-04 21:52:03
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

mozhui

★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
wuli8(金币+2,VIP+0): 6-5 20:49
http://muchong.com/bbs/viewthread.php?tid=1293203&fpage=4

Vasp&MedeA

【讨论】paw势过于陈旧，会导致计算会出现错误么？
现在我使用的paw势都是98年到01年的。尤其是用的比较多的O 的paw-GGA。结果计算的体系在自洽过程中，很难收敛。已经调整过mixing的方法等促进收敛的方法。依然没有作用。后来按照手册上测试单原子的赝势，发现能量在几百毫电子伏，与期望的1~10m eV有很大出入、
请问大家用的paw势都是哪一年出版的呢？有没有遇到过同样的问题

2楼：VASP里面都是01年造的吧，我觉得和时间关系不大，都是成熟的理论，赝势主要靠测试，如果不行就换吧，
3楼：我觉得，老的、很多人用过的，才更保险啊。除非你用的是很早就造出来但从来没有人用过的。
赞一下(11人) 回复此楼
47楼2009-06-04 21:57:08
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

mozhui

★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
wuli8(金币+2,VIP+0): 6-5 20:49
http://muchong.com/bbs/viewthread.php?tid=1306532&fpage=4

Vasp&MedeA

【求助】谁知道新版VASP5.2的价钱？
2楼：1000欧元从4.6升级。全新购买4000欧元
3楼：一个licence，升级1000；
两个licence，升级2000.
赞一下(11人) 回复此楼
48楼2009-06-04 22:00:01
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wuli8 的主题更新
251/3返回列表123下一页
信息提示
关闭
请填处理意见
关闭