MaterialsStudio 论坛问答全集Word文件下载.docx

1、1，在什么系统下装，是否装了并行版本。2，计算时设置参数的地方是否选择了并行。3，程序运算时，并不是时时刻刻都要用到多个CPU 3、问：我已经成功地安装了MS3.1的Linux版本， 串行的DMol3可以成功运行。但是运行并行的时候出错。机器是双Xeon5320（四核）服务器，rsh和rlogin均开启，RHEL4.6系统。其中hosts.equiv的内容如下：localhost ibm-console machines.LINUX的内容如下：localhost:8 现在运行RunDMol3.sh时，脚本停在 $MS_INSTALL_ROOT/MPICH/bin/mpirun $nolocal

2、 -np $nproc $MS_INSTALL_ROOT/DMol3/bin/dmol3_mpi.exe $rootname $DMOL3_DATA 这一处，没法执行这一命令 并行运算时，出现以下PIxxxx（x为数字）输出 ibm-console 0 /home/www/MSI/MS3.1/DMol3/bin/dmol3_mpi.exe localhost 3 /home/www/MSI/MS3.1/DMol3/bin/dmol3_mpi.exe 请问这是什么原因？谢谢！主要是rsh中到ibm-console的没有设置 把/etc/hosts改为 127.0.0.1 localhost.lo

3、caldomain localhost ibm-console 在后面加个ibm-console 也希望对大家有帮助!4、问：在最后结果的dos图中,会显示不同电子spd的贡献,我想问的是, 假设MS考虑的原子Mg的电子组态为2p6 3s2,那么最后的dos结果中的s,p是不是就是2p,跟3s的贡献.比如更高能量的3p是否可能出现在 dos中?如果可能的话，在这种情况下,如何区分2p和3p的贡献,谢谢. A、取决于你的餍势 势里面没有3p电子，DOS怎么会有呢？自然，你的1p1s也不会出现在你的DOS中。B、我觉得原子没有d轨道但是计算出来的DOS图上有d轨道的贡献，这可能是电子跃迁所引起的，

4、虽然按照教科书上的d轨道上没有电 子，但是不排除有电子跃迁过去，但及时跃迁也应该只是很少的而且一般也都是这样。C、我考虑的是dos并非指电子出现的几率,而是能级的密度.3p构成导带能级不可能吗?没有d电子出现d的贡献可以理解是跃迁过去,那是不是也不否认3p的可能?既然不可能有3p,怎么会有1s1p?毕竟一开始计算中已经做了近似?D、第一性原理计算高于费米面的能带本来就不怎么可靠的。现在，计算的带隙小于实验值，很大的原因就来源于价带计算不准。p能量高于s不太好理解 这里有成键反键的因素 p的反键是有可能小于s的成键的 另外有时候原子电子组态没有d电子的时候计算结果也出现了d轨道 估计是做轨道投影

5、的问题。vasp就有这个问题，但是可以消除这个错误的。5、问：那位高手能够介绍一下MS中CASTEP的势能模型的适用范围？比如说LDA的（PWC,VWN），GGA的（PW91,BP,PBE,BLYP,BOP,VWN-BP,PRBE,HCTH），能够简述一下其大致差别更好 LDA局域密度近似（LDA）：局域密度近似（LDA）是第一阶梯。它仅仅采用空间点r处的电子密度n（r）来决定那点交换-相关能 密度的形式。交换-相关能密度由密度相同的均匀电子气完全确定。泛函的交换部分就准确的用均匀电子气的微分表达。各种不同的 局域密度近似（LDA）仅仅是相关部分表示方法不同，所有现代应用的局域密度泛函都基于C

6、eperly和Alders在年代对均匀电子 气总能量的Monte Carlo模拟。广义梯度近似（）：是Jacob阶梯的第二个台阶，将电子密度的梯度也作为一个独立的变量（|n（r）|），在描述交换 相关能方面，梯度引入了非定域性。泛函包含了两个主要的方向：一个称为“无参数”，泛函中新的参数通过已知形式中参数 或在其它准确理论帮助下得到。另外一个就是经验方法，未知参数来自于对实验数据的拟和或通过对原子和分子性质准确的计算。Perdew,Burke and Emzerhof（PBE）以及Perdew-Wang from 1991（PW91）是无参数的，在量子化学中广泛采用的GGA，比如 Becke,

7、Lee,Parr and Yang（BLYP）是经验性。LYP校正采用了密度的二阶Laplace算符，因此严格上讲属于Jacob阶梯的第三阶，但通常 仍然归类为GGA. 6、求助Linux上MS的安装问题 我在RHES4.6上安装MS3.0 但是在castep和dmol下的RunDMol3.sh大小为零，没有东西， 同时在share下也没有bin目录，也找不到ms_setuo.sh脚本 请问是什么原因？相应的出现这样的错误信息：/home/sript/MSI/3.0/3rdParty/regxpcom: error while loading shared libraries: libcxa.

8、so.3: cannot open shared object file: No such file or directory 假如你有Materials Studio 3.1的ISO版本以及licene.dat。在linux机器上的用户名为firefox。下面是安装Materials Studio3.1 linux版本的步骤 1、假定把Materials Studio3.1的iso文件mount到/home/firefox/ms目录，切换到root帐号下，然后用下面的命令如下：mount -o loop /home/firefox/MS_MODELING.iso /home/firefox/

9、ms 2、文件mount上去后，退出root帐号，回到firefox帐号。然后到/home/firefox/ms/UNIX目录，运行下面的：./Install 3、开始按提示安装了 1）、指定你要Materials Studio安装到的目录，比如/home/firefox/MaterialsStudio 2）、下面就是指定你所需要安装的模块了，根据自己的需要选择，当然也可以全部选上 3）、然后指定你License_Pack所在的目录，我们假定也是放到/home/firefox/MaterialsStudio/License_Pack目录中，那么先建 立这样的目录：mkdir /home/fir

10、efox/MaterialsStudio/License_Pack 接着按提示Enter。开始安装Licens_Pack了。4）、下面把你所得到的licene文件msilic.lic拷贝到/home/firefox/MaterialsStudio/License_Pack/licenses目录中，并改名为 msilicense.dat。5）、把/home/firefox/MaterialsStudio/License_Pack/licenses/msilicense.dat完整的文件名按提示输入。这些完成后。就基本 安装好了。下面就是把gateway运行起来。4、到/home/firefox/

11、MaterialsStudio/Gateway目录，运行下面的命令：./msgateway_control_18888 start 提示说Gateway start succeeded - running as process xxxx。后，就启动成功了。5、在ie的地址栏输入，http:/127.0.0.1:18888/ ; 则页面会显示该机器上已经安装好的计算模块 基本就安装成功了。（注意之处：不要用root安装它）。7、问：我想对NaYF4的电子云分布进行模拟计算,请各位高手指点怎么计算?如果有NaYF4的结构的话,请分享一下. 如果是晶体的话，可以找相关的文献看晶格常数和原子坐标，自己

12、构建模型后计算 如果没有文献，可以上无机晶体结构数据库ICSD查询，如果查到，直接导出来 就可以用ms算了 电子云分布算一下电子密度就可以了 可以用dmol和castep 如果要某个轨道的电子云，dmol算homolumo上下若干个轨道 castep算能带 最后那visulizer看电子云密度图就ok了 8、问：各位高手请教一下用demolJ计算分子轨道怎么看HUMO,LUMO图，其中的黄色蓝色代表什么啊，谢谢了 dmol计算的时候选择properties里面的orbitals 把homo和lumo选上 然后计算，算完以后打开analysis对话框 选择orbitals，就看到你算的homo和

13、lumo了，然后import就在模型中可以看见homo和lumo了 9、求助有MS并行计算经验的朋友解答一下 我虽然在LINUX CLUSTER并行成功了，但是并行的时候CPU计算利用率很不均衡也不稳定 有时候100%有时候0% 而且计算到一定程度的时候会变换运行程序在LINUX中的进程名 换了一次之后CPU使用就开始出现不均衡的情况，有一个CPU 90%+，别的都是2-5%的样子进程处与sleeping状态只有那一个CPU在RUN 另外说一下我每个节点的硬件配置相同，用的是CASTEP模块，K点数是CPU数的整数倍，但是算能带的时候不知道如何调节K点数目， 不均衡不知道是不是算能带的时候导致

14、K点在每个CPU上分布不均造成的 还想问一下 我各个节点都是同一个交换机下的机器，NFS的效率如何评估 A、你应该看看内存的使用率,一般如果内存不足,会造成你说的现象 B、所有节点空闲内存都在50%以上啊，发现只有一个节点的CPU占用的厉害 别的不超过10% C、你选用的几个cpu进行计算,如果是一个,那就是你发生的问题原因,如果不是一个,而是选择几个cpu进行计算,那么应该是节点之间 的通讯问题了,这需要进行机群的的调节了 10、问：CASTEP能算pdos的自旋情况吗 可以计算体系的自旋态密度 alpha和beta分别为自旋向上和自旋向下 total spin是总自旋 计算的时候必须把自旋

15、极化spin polarized 选上才能画自旋态密度 11、问：大家帮忙看看这是什么问题 我再用demol计算动力学的时候,不管怎么设置参数它老出现这个提示 Error: SCF iterations not converged in 50 iterations Message: SCF not converging. Choose Use Smearing on DMol3 SCF panel or set Occupation Thermal in the input file You may also need to change spin or use symmetry Resubm

16、it DMol3 DMol3 job failed DMol3 exiting 请教高手这是什么原因啊?SCF not converging. Choose 无法收敛，加大计算步数。你错误的原因是体系在50步内不能收敛,你可以增加收敛计算步数,比如增加到100 系统提示你改变和 的设置参数（这两个参数在设置对话框中）,还提出建议让你改变自旋和对 称设置来试一试。不收敛，解决方法一加大圈数，DMol3 calculationselectronic标签下More选项SCF选项，中Max. SCF Cycles加大，最多 1000圈，如果1000圈收敛不了 就用方法二，加大smearing 在SCF

17、标签下把orbital occypancy下面的勾打上，然后设置smearing 设置重小到大逐渐增大，smearing如果很大，虽然能收敛但是最后的构型有问题。做MS计算，我的感觉是没有必要一步算到位，就是说你的计算精度（能量，力等等）可以一步一步来，或者是smearing慢慢变小，这 样容易收敛些。12、新手请教几个不专业的问题 1.（1x1）CO中（1x1）是什么意思？2.态密度、能态密度和电子态密度的区别？3. K-resolved是什么意思？4.如何在半导体上建金属表面？5.如何翻译？supercell geometry the Vienna ab-initio simulation

18、 package （VASP） the generalized gradient approximation the plane wave expansion structural relaxations a tight-binding linear muffintin orbital （TB-LMTO） method the TB-LMTO method the principal-layer Greens function technique the K-resolved transmission 6K-resolved DOS at the Fermi energy这句话怎么理解啊？7用

19、MS能不能计算一个体系的电导？如何去算？1. （11）表示的是表面的构型，22是（11）的四倍，也就是x和y上分别用两倍单位平移定义表面的晶格常数a和b，c轴与表面 方向平行，c的大小取决于原子层和真空层厚度 2.应该是没有区别，可能是称呼不同而已。3.不清楚什么意思 4. 用ms里面的layer builder构建界面，基底为半导体，半导体上长金属层是通过layer builder构建的 5. supercell geometry超胞结构 the Vienna ab-initio simulation package （VASP） 维也纳从头算模拟软件包 the generalized gr

20、adient approximation广义梯度近似 the plane wave expansion平面波展开 structural relaxations结构弛豫 a tight-binding linear muffintin orbital （TB-LMTO） method紧束缚线性丸盒轨道 the TB-LMTO method紧束缚线性丸盒轨道方法 the principal-layer Greens function technique？the K-resolved transmission？6.不是很清楚，没有上下文 7 先计算出态密度。然后通过由玻尔兹曼方程推导出的电导率公式来

21、算,见固体物理。提供两篇参考文献是做输运系数的计算的：Georg K.H.Madsen, David J.Singh. Computer Physics Communications 175 （2006 6771. Philip B. Aleen ,Warren E.Pickett, Henry Karkauer. Physcical Review B 37 （1988） 27212726. 13、请问： 在用CASTEP计算材料空位缺陷时，如何才能实现只取一个空位呢？降低对称性，再DELETE一个原子就可以了 在build下的symtrey有make P1 不是降低对称性，是构建超胞 mak

22、e supper cell 然后就可以删去一个原子了 正常的话如果晶体是原胞，直接删除一个原子可能造成很大影响 可以建大一点的超晶胞（比如222超晶胞）然后删掉一个原子，这样比较合理， 毕竟缺陷浓度是很低的。14、求助 MS是否可以做原子数在10-20个左右的簇合物的结构模拟 我想得到一些原子数在10-20个左右的一元或二元原子簇合物的低能构型。MS里的分子模拟功能是否可以模拟出一些能量局部最小化的构型？如果不行的话，麻烦熟悉这方面的虫子给推荐个软件，最好是免费的。我的要求很简单，就只要找到低能量构型即可。手动直接构建实在是太麻烦了，而且说服力很低。谢谢指教！A、可以，dmol模块可以达到你的

23、要求 B、感谢楼上各位。不过我琢磨了半天，好像无论是Gaussian还是MS,都是要给定基本定型的构型，然后进行优化后获得更精炼一些的构型。我想做这样一件事： 已知原子数为 10 的构型，第 11 个原子任意放置，然后程序自动变化第11个原子的位置，甚至变化 11 个原 子的位置 （既第11个原子对其他10原子的构型进行微扰）得到不同能量的局部的低能构型。自己琢磨的，肯定很多不足，不是否有相关的参考资料。还请赐教， 15、ms discover停电后如何处理 ms discover算了一个星期，还有三天就算完了，结果停电全没了 重开机，启动ms，提示有任务非正常关掉，载入auto recove

24、r，job explorer里面还显示未算完的任务，server console打不开，无 法查看是否正在算 查任务管理器，没有找到discover.exe进程，只有msstudio的进程，而且cpu占用高 手动停止job，然后从discover的dynamics模块下选中use restart result，但是结果显示完全从头开始算 没钱买ups，很不走运，两个星期内，这种停电出现两次了 请教大家有什么好的方法，从断电时的数据开始接着算 谢谢 如果你在dynamic中trajectory的选项里save选择了full，那么在重启动时，在restart选项里，选中use restart da

25、te选项，就 可以了。16、求助：黄色和蓝色的电子云代表什么问题 请教各位一个问题，在MS的Dmol3下计算能量，在分析里查看HOMO和LUMO轨道的电子云的计算结果时，有两种颜色，默认的是黄色和 蓝色。这种颜色的不同代表什么意思？还有一般用HOMO和LUMO轨道的电子云说明什么问题啊？我看很多文献说什么（-）派-派跃迁就是3.1415926所代表的那个希腊字母 ，这里不知道怎么输入。还有（-）西格玛-西格玛等等，学习化学的人应该能够了解，可以从这上面看出来吗？A、其实两种颜色的电子云不是代表homo和lumo的，只是说这个地方波函数是正号还是负号。因为对于电子云来说，其实是波函 数的平方，这

26、个地方为了区别波函数的正负号，所以引进了两种颜色的电子云 B、在分子进行反应时，最高占据轨道（HOMO）和最低空轨道（LUMO）其有特别重要的作用。分子轨道对称性守恒原理认为：在反应过程中，当有反应物分子轨道转变成产物的分子轨道时，其对称性始终保持不变，即始终都属 于同一个不可约表示。如果反应物基态的分子轨道变成的是产物的基态的分子轨道，那么反应就是热允许的。（意即温和加热就可进行） 如果变成的是激发态的分子轨道，就是热禁阻的（意即不可能用温和加热的方法促其进行）。如果反应物受光激发后由电子占据的分子轨道（激发态轨道）转变成产物的基态分子轨道，这一过程就是光允许的（意即用适当的光 照就可使其进

27、行）。若转变成产物的更高激发态分子轨道，就是光禁阻的（意即无法用光照促其进行）。17、请问MS计算能带的问题 由CASTEP或DMol计算得到的能带，要看能带的带系宽度，是怎样看？是默认费米能级为0吗？那就是只看0附近的宽度吗？A、0两边的宽度，最高占据和最低未占据之间的宽度就是带隙 B、我没有仔细的看过castep的能带，但是看过课题组的人做过报告。我印象中，对于非金属，castep的费米能级是在价带顶，找出价带顶最大的那个能量和导带底对应的最小的能量之差，就是带隙吧。至于有的说dos看带隙，这个不大可靠，要看你取了多大的展宽因子，展宽因子大，非金属也变成金属了，没带隙。另外对于dmol来说

28、，何来什么能带 全部分子能级，不过可能有gap，要从homo-lumo的差值去算gap C、横轴是布里渊区的高对称点，与晶体的对称性有关 纵轴表示能带的能量。带隙可以分成直接带隙和间接带隙 直接带隙就是0eV上下最高和最低两个能带在某个高对称点（如G，X、R点等）之间的能量差。间接带隙就是0eV上能量最低的能带上的最低点的能量值减去0eV下能量最高的能带的最高点的能量值。18、什么是Overlap populations？在一篇文章上看见 Table 2 Bond distances （A ） and Overlap populations （/E/in brackets） TiO2 TiO2/Sn6c TiO2/Sn5c TiO2/Sn6-in TiO2/Sn5-in ObMe6c 1.831 （172） 1.963 （220） 1.83 （165） 1.831 （172） 1.832 （171） OpMe6c 2.066 （86） 2.107 （32） 1.98