1、关于 K 点1. 应当使用多少个 k 网格？很难一般地回答，只能给出一般建议。注意：一定要检查 k 网格，首先用较粗糙的网格计算，接下来用精细的网格计算。通过比较两次的结果，决定选用较粗糙的网格，或是继续进行更精细网格的计算，直到达到收敛。金属体系需要精细的网格，绝缘体使用很少的 k点通常就可以。小单胞需要精细格点，大单胞很可能不需要。因此：单位晶胞内原子数很多（比如 40-60 个）的绝缘体，可能仅需要一个（移动后的）k 点。另一方面，面心立方的铝可能需要上万个 k 点以获得好的 DOS。对于孤立原子或分子的超晶胞，仅需要在 Gamma点计算。对于表面（层面）的超晶胞计算，仅需要（垂直于表面
2、）z 方向上有 1 个 k 点。甚至可以增加晶格参数 c，这样即使对精细格点，沿 z 方向上也只产生一个 k 点（产生 k点后，不要忘记再把 c 改回） 。2. 当体系没有出现时间反演对称操作时，是否加入？大多数情况下的回答是“是” ，只有包含自旋-轨道耦合的自旋极化（磁性）计算除外。这时，时间反演对称性被破坏（+k 和-k 的本征值可能不同） ，因此决不能加入时间反演对称性。3. 是否移动 k 网格？（只对某些格子类型有效）“移动”k 网格意味着把所有产生的 k 点增加(x,x,x)，把那些位于高对称点（或线）上的k 点移动到权重更大的一般点上。通过这种方法（也即众所周知的“特殊 k 点方法
3、” ）可以产生等密度的，k 点较少的网格。通常建议移动。只有一点注意：当对半导体的带隙感兴趣时（通常位于 Gamma，X，或 BZ 边界上的其它点） ，使用移动的网格将不会得到这些高对称性的点，因此得到的带隙和预期结果相比或大或小。这个问题的解决：用移动的网格做SCF 循环，但对 DOS 计算，改用精细的未移动网格。关于 k 空间布点的问题，建议参阅以下文献 Phys.Rev.B 49,16223 (1994) 如何构建缺陷晶体结构晶体结构改成 P1，然后去掉想抹去的原子就可以了在 ms 中如何做空穴对于金属缺陷，是直接剪切一个原子？个人经验：就是直接把原子去掉就；如果不是正版软件，有可能出现
4、同时去掉其他同位置的原子，如果这种情况，就重新定义，问题就不会出现了还有，一般考虑孔穴的时候，都要标明哪些原子的迟豫，具体为什么不知道，国外的文献有提到希望有做空位的一起多讨论我：PDOS 选项计算 DOS 时，选择 PDOS，可以画出 s,p,d 轨道的 DOS，但无法画出某一个原子的 s,p,d 图关于 PDOS 的 Chart 中求积分的问题在用 Castep 计算出 PDOS 后,如何在 Chart 中对曲线局部进行积分?将 Chart 输出为 cav 格式,然后在 excel 中求和 ?简单，把数据导出，在 Origin 里作图，程序里有积分微分卷积功能，在数据分析下面。作图时选取积
5、分范围。优化结构算能带一般需要优化结构。如果选择实验的参数，全部固定的话就不需要了如何做二维电子密度图MS 结果文件夹中*.grd 文件内存储的是三维空间各点的电荷密度值，利用这个数据就可以得到二维的电荷密度等值线图,应该有专门的软件能画,不过用 matlab 编自己编程序也不难，其中关键命令是 contourslice，实现在某一平面内绘制等高线对 DOS 图的分析根据 DOS 的积分曲线可以计算出，对于表面 Nb和 C 原子，大约有 16.8%和 14.8%的电子态集中在 4.02.0eV 的区域，而对体相原子则分别为 6.4%和 6.0%表面吸附我做 H 在 ZnO 上吸附。刚开始时后我
6、构建的吸附构型忘记 impose symmetry 了，Groupname 是 P1。在第二次计算的时候 我加上了 symmetry。两次计算差别出来了： （1） 首先是 impose symmetry 后，supercell 中的原子位置由原来的现面跑到了上面，也就是和真空层换了一下位置！ 而且吸附原子竟然 超出了supercell 的 白色线框！ （2）比较两次计算的 DOS，一模一样；但是两次计算的 Band Structure 却有很大差异！回答：1）只是显示问题（2） Band Structure 有很大差异是指那种差异？使用不同对称性计算能带时，默认计算的 K 点是不同的，所以图像
7、肯定不同。如果你确定是选择计算了同样的 K 点，能带结构仍然不同，那可能是采用对称后结构变化导致的。如何计算结构中某一元素的分波态密度1. 我最近看了关于氧气锌的论文，上面有锌原子的分波态密度图，可是我怎么都没算出来，一直得到氧和锌原子的分波态密度，希望哪位大侠指教2. 在计算性质选择了计算 density of state 是，对话框下面有一个口 calculate PDOS,把它勾选上，计算成功后，在分析就能看到总态密度和分态密度，随你选择。3. 按照楼上的指导计算成功后，需要哪个原子的 PDOS,就选中哪个原子，然后在 analysis里 DOS 项前打勾，partial 项打勾，再 v
8、iew 就 ok 了。 有关能带分析能带图分析能带图的横坐标是在模型对称性基础上取的 K 点。为什么要取 K 点呢？因为晶体的周期性使得薛定谔方程的解也具有了周期性。按照对称性取 K 点，可以保证以最小的计算量获得最全的能量特征解。能带图横坐标是 K 点，其实就是倒格空间中的几何点。其中最重要也最简单的就是 gamma 那个点，因为这个点在任何几何结构中都具有对称性，所以在castep 里，有个最简单的 K 点选择，就是那个 gamma 选项。纵坐标是能量。那么能带图应该就是表示了研究体系中，各个具有对称性位置的点的能量。我们所得到的体系总能量，应该就是整个体系各个点能量的加和。 记得氢原子的
9、能量线吧？能带图中的能量带就像是氢原子中的每条能量线都拉宽为一个带。通过能带图，能把价带和导带看出来。在 castep 里，分析能带结构的时候给定 scissors 这个选项某个值，就可以加大价带和导带之间的空隙，把绝缘体的价带和导带清楚地区分出来。 DOS 叫态密度，也就是体系各个状态的密度，各个能量状态的密度。从 DOS 图也可以清晰地看出带隙、价带、导带的位置。要理解 DOS，需要将能带图和 DOS 结合起来。分析的时候，如果选择了 full，就会把体系的总态密度显示出来，如果选择了 PDOS，就可以分别把体系的 s、 p、d、f 状态的态密度分别显示出来。还有一点要注意的是，如果在分析
10、的时候你选择了单个原子，那么显示出来的就是这个原子的态密度。否则显示的就是整个体系原子的态密度。要把周期性结构能量由于微扰裂分成各个能带这个概念印在脑袋里。 最后还有一点，这里所有的能带图和 DOS 的讨论都是针对体系中的所有电子展开的。研究的是体系中所有电子的能量状态。根据量子力学假设，由于原子核的质量远远大于电子，因此奥本海默假设原子核是静止不动的，电子围绕原子核以某一概率在某个时刻出现。我们经常提到的总能量，就是体系电子的总能量。 这些是我看书的体会，不一定准确，大家多多批评啊！如何分析第一原理的计算结果转自：http:/ 转自量化网 http:/ 摘要：本文总结了对于第一原理计算工作的
11、结果分析的三个重要方面，以及各自的若干要点用第一原理计算软件开展的工作，分析结果主要是从以下三个方面进行定性/定量的讨论：1、电荷密度图（charge density） ； 2、能带结构（Energy Band Structure） ； 3、态密度（Density of States，简称 DOS） 。 电荷密度图是以图的形式出现在文章中，非常直观，因此对于一般的入门级研究人员来讲不会有任何的疑问。唯一需要注意的就是这种分析的种种衍生形式，比如差分电荷密图（def-ormation charge density）和二次差分图（difference charge density）等等，加自旋极化
12、的工作还可能有自旋极化电荷密度图（spin-polarized charge density） 。所谓“差分”是指原子组成体系（团簇）之后电荷的重新分布， “二次” 是指同一个体系化学成分或者几何构型改变之后电荷的重新分布，因此通过这种差分图可以很直观地看出体系中个原子的成键情况。通过电荷聚集（accumulation）/ 损失（depletion）的具体空间分布，看成键的极性强弱；通过某格点附近的电荷分布形状判断成键的轨道（这个主要是对 d 轨道的分析，对于 s 或者 p 轨道的形状分析我还没有见过） 。分析总电荷密度图的方法类似，不过相对而言，这种图所携带的信息量较小。 能带结构分析现在在
13、各个领域的第一原理计算工作中用得非常普遍了。但是因为能带这个概念本身的抽象性，对于能带的分析是让初学者最感头痛的地方。关于能带理论本身，我在这篇文章中不想涉及，这里只考虑已得到的能带，如何能从里面看出有用的信息。首先当然可以看出这个体系是金属、半导体还是绝缘体。判断的标准是看费米能级和导带（也即在高对称点附近近似成开口向上的抛物线形状的能带）是否相交，若相交，则为金属，否则为半导体或者绝缘体。对于本征半导体，还可以看出是直接能隙还是间接能隙：如果导带的最低点和价带的最高点在同一个 k 点处，则为直接能隙，否则为间接能隙。在具体工作中，情况要复杂得多，而且各种领域中感兴趣的方面彼此相差很大，分析
14、不可能像上述分析一样直观和普适。不过仍然可以总结出一些经验性的规律来。主要有以下几点： 1） 因为目前的计算大多采用超单胞（supercell）的形式，在一个单胞里有几十个原子以及上百个电子，所以得到的能带图往往在远低于费米能级处非常平坦，也非常密集。原则上讲，这个区域的能带并不具备多大的解说/阅读价值。因此，不要被这种现象吓住，一般的工作中，我们主要关心的还是费米能级附近的能带形状。 2） 能带的宽窄在能带的分析中占据很重要的位置。能带越宽，也即在能带图中的起伏越大，说明处于这个带中的电子有效质量越小、非局域（non-local）的程度越大、组成这条能带的原子轨道扩展性越强。如果形状近似于抛
15、物线形状，一般而言会被冠以类 sp 带（sp-like band）之名。反之，一条比较窄的能带表明对应于这条能带的本征态主要是由局域于某个格点的原子轨道组成，这条带上的电子局域性非常强，有效质量相对较大。 3） 如果体系为掺杂的非本征半导体，注意与本征半导体的能带结构图进行对比，一般而言在能隙处会出现一条新的、比较窄的能带。这就是通常所谓的杂质态(doping state)，或者按照掺杂半导体的类型称为受主态或者施主态。 4） 关于自旋极化的能带，一般是画出两幅图：majority spin 和 minority spin。经典的说，分别代表自旋向上和自旋向下的轨道所组成的能带结构。注意它们在
16、费米能级处的差异。如果费米能级与 majority spin 的能带图相交而处于 minority spin 的能隙中，则此体系具有明显的自旋极化现象，而该体系也可称之为半金属（half metal） 。因为 majority spin 与费米能级相交的能带主要由杂质原子轨道组成，所以也可以此为出发点讨论杂质的磁性特征。5） 做界面问题时，衬底材料的能带图显得非常重要，各高对称点之间有可能出现不同的情况。具体地说，在某两点之间，费米能级与能带相交；而在另外的 k 的区间上，费米能级正好处在导带和价带之间。这样，衬底材料就呈现出各项异性：对于前者，呈现金属性，而对于后者，呈现绝缘性。因此，有的工
17、作是通过某种材料的能带图而选择不同的面作为生长面。具体的分析应该结合试验结果给出。 （如果我没记错的话，物理所薛其坤研究员曾经分析过$beta$-Fe 的(100)和(111)面对应的能带。有兴趣的读者可进一步查阅资料。 ） 原则上讲，态密度可以作为能带结构的一个可视化结果。很多分析和能带的分析结果可以一一对应，很多术语也和能带分析相通。但是因为它更直观，因此在结果讨论中用得比能带分析更广泛一些。简要总结分析要点如下： 1） 在整个能量区间之内分布较为平均、没有局域尖峰的 DOS，对应的是类 sp 带，表明电子的非局域化性质很强。相反，对于一般的过渡金属而言，d 轨道的 DOS 一般是一个很大
18、的尖峰，说明 d 电子相对比较局域，相应的能带也比较窄。 2） 从 DOS 图也可分析能隙特性：若费米能级处于 DOS 值为零的区间中，说明该体系是半导体或绝缘体；若有分波 DOS 跨过费米能级，则该体系是金属。此外，可以画出分波（PDOS）和局域（LDOS）两种态密度，更加细致的研究在各点处的分波成键情况。 3） 从 DOS 图中还可引入“赝能隙” （pseudogap ）的概念。也即在费米能级两侧分别有两个尖峰。而两个尖峰之间的 DOS 并不为零。赝能隙直接反映了该体系成键的共价性的强弱：越宽，说明共价性越强。如果分析的是局域态密度（LDOS） ，那么赝能隙反映的则是相邻两个原子成键的强弱
19、：赝能隙越宽，说明两个原子成键越强。上述分析的理论基础可从紧束缚理论出发得到解释：实际上，可以认为赝能隙的宽度直接和 Hamiltonian 矩阵的非对角元相关，彼此间成单调递增的函数关系。 4） 对于自旋极化的体系，与能带分析类似，也应该将 majority spin 和 minority spin 分别画出，若费米能级与 majority 的 DOS 相交而处于 minority 的 DOS 的能隙之中，可以说明该体系的自旋极化。 5） 考虑 LDOS，如果相邻原子的 LDOS 在同一个能量上同时出现了尖峰，则我们将其称之为杂化峰（hybridized peak） ，这个概念直观地向我们展
20、示了相邻原子之间的作用强弱。 以上是本人基于文献调研所总结的一些关于第一原理工作的结果分析要点。期冀能对刚进入这个领域内的科研工作者有所启发。受本人的水平所限，文章的内容可能会有理论上的不足甚至错误之处，希望大家指出，共同发展第一原理计算物理的方法和研究内容。 希望大家多读几遍。 http:/cmt.dur.ac.uk/sjc/Castep_Publications/pubblications.html看看人家用 castep 怎样选择赝势和 K 点的关于出错信息的解读（对自己解决问题有帮助）Geomrtry Optimization 时，到 77 步（共 100 步）就没有反应，我等了一天终
21、于把它停止了，文件中有下列出错信息，请大侠帮忙解读一下。精度是 Fine，之前在 coarse 下可以24 步就完成。 1、Warning: There are no empty bands for at least one kpoint and spin; this may slow the convergence and/or lead to an inaccurate groundstate. If this warning persists, you should consider increasing nextra_bands and/or reducing smearing_wid
22、th in the param file. Recommend using nextra_bands of 6 to 12. 2、BFGS: Warning - Repeated consecutive reset of inverse Hessian BFGS: without satisfying convergence criteria which BFGS: looks like BFGS has run out of search directions. BFGS: Warning - Lets try allowing some uphill steps and see if BF
23、GS: we can get around this barrier. BFGS: Warning - It is possible that the system may now converge to BFGS: a stationary point OTHER than the desired minimum. BFGS: Hint - this may be an indication that either: BFGS: a) you are using a poor guess at geom_frequency_est BFGS: and/or geom_modulus_est,
24、 or BFGS: b) you are using unrealistic convergence criteria. BFGS: Suggest therefore that you consider changing them! 3、energies not corrected for finite basis set 什么意思？要修改那些参数？回答：你可以在 property 中的 SCF 增加 empty band 和减小 smearing 来试试关于结合能 形成能的问题我不太明白形成能 结合能和总能之间的关系?是否可以指教回答：举例来说吧 。对于 Si 块体材料，通过 CASTEP
25、 计算得到的 能量就叫总能。如果你想把它切开，得到（100）面，通过添加真空层 构建一个 spercell，计算得到的也叫总能。 刚才对于 Si 块体材料的总能乘上系数再减去 Si（100）面总能 这就是 Si（100）面的形成能。 乘上的系数 保证 他们的原子总数一致。结合能： 如果在 Si(100)面上吸附了一个 Mn 原子（假如你关心磁性） ， 和原来一样计算出总能后， 减去 Si(100)面的 总能 再减去 Mn 原子的总能 ，得到的就是 binding energy 那位高手能够介绍一下 MS 中 CASTEP 的势能模型的适用范围？比如说 LDA 的（PWC,VWN） ，GGA 的
26、（PW91,BP,PBE,BLYP,BOP,VWN-BP,PRBE,HCTH） ，能够简述一下其大致差别更好LDA 局域密度近似（LDA）：局域密度近似（LDA）是第一阶梯。它仅仅采用空间点 r 处的电子密度 n（r）来决定那点交换-相关能密度的形式。交换-相关能密度由密度相同的均匀电子气完全确定。泛函的交换部分就准确的用均匀电子气的微分表达。各种不同的局域密度近似（LDA）仅仅是相关部分表示方法不同，所有现代应用的局域密度泛函都基于Ceperly 和 Alders 在年代对均匀电子气总能量的 Monte Carlo 模拟。广义梯度近似（）：是 Jacob 阶梯的第二个台阶，将电子密度的梯度也
27、作为一个独立的变量（|n(r)|） ，在描述交换相关能方面，梯度引入了非定域性。泛函包含了两个主要的方向：一个称为“无参数” ，泛函中新的参数通过已知形式中参数或在其它准确理论帮助下得到。另外一个就是经验方法，未知参数来自于对实验数据的拟和或通过对原子和分子性质准确的计算。Perdew,Burke and Emzerhof(PBE)以及 Perdew-Wang from 1991(PW91)是无参数的，在量子化学中广泛采用的 GGA，比如 Becke,Lee,Parr and Yang(BLYP)是经验性。 LYP 校正采用了密度的二阶 Laplace 算符，因此严格上讲属于Jacob 阶梯的
28、第三阶，但通常仍然归类为 GGA.（见后面英文介绍）请问： 在用 CASTEP 计算材料空位缺陷时，如何才能实现只取一个空位呢*在 build 下的 symtrey 有 make P1*不是降低对称性，是构建超胞make supper cell然后就可以删去一个原子了正常的话如果晶体是原胞，直接删除一个原子可能造成很大影响可以建大一点的超晶胞(比如 222 超晶胞)然后删掉一个原子，这样比较合理，毕竟缺陷浓度是很低的。*都可以啊！先看一下 show symmetry 是 P1 就可以直接删一个，如果不是就换成 P1，一般 make supercell 后自动变成 P1 了，你可以根据自己所需的
29、浓度选择其中一种方式MS 无法收敛请大家帮忙看看了,掺杂 Eu 的体系,一开始提示增加 empty band,减小 smear,都做了,把精度设置到 coarse,可最后还是在 100 轮后提示无法收敛.- - SCF*Warning* max. SCF cycles performed but system has not reached the groundstate.Current total energy, E = -10141.84970726 eVCurrent free energy (E-TS) = -10141.84970726 eV(energies not correct
30、ed for finite basis set)NB est. 0K energy (E-0.5TS) = -10141.84970726 eV*Warning: electronic minimisation did not converge duringfinite basis set correction.*Checkpoint file cannot be written.Error calculate_finite_basis : Convergence failed when doing finite basis set correction.Error calculate_fin
31、ite_basis : Convergence failed when doing finite basis set correction.*Warning: electronic minimisation did not converge duringfinite basis set correction.*Checkpoint file cannot be written.1 MPI Abort by user Aborting program !1 Aborting program!0 MPI Abort by user Aborting program !0 Aborting prog
32、ram!回答：1.增加 smearing 容易收敛，但是太大可能结果会有些问题。2.增加循环的圈数，3.增加空带数目，4.体系有自旋单电子的，要设置自旋极化，如果设置非自旋极化，可能不收敛5. 如果是金属体系，比较难收敛，可以试试修改 k 点和 empty band 试试看。1. 几何优化可以在多大范围内改变原子的位置啊？感觉每次优化后的位置调整不是很多啊？2. 大家在算能带和态密度的时候的时候是用单点能计算给出能带还是几何优化给出能带啊？回答：问题 1几何优化的目的是寻找压力最小的几何结构，原子位置改变不多是因为你建的构形比较合理，比如直接从软件数据库中导入的结构在 0 压力下改变很小，如果
33、加个压力就会变化大一些。问题 2一般都是优化的时候算，如果你想算特定构形的能带或其它性质，就可以用单点能；【表面模拟】castep 计算中遇到的几个问题PW91 Perdew-Wang generalized-gradient approximation Perdew and Wang (1992) BP Becke exchange plus Perdew correlation Becke (1988), Perdew and Wang (1992) PBE Perdew-Burke-Ernzerhof correlation Perdew et al. (1996) RPBE Revis
34、ed PBE functional by Hammer et al. Hammer et al. (1999) HCTH Hamprecht, Cohen, Tozer and Handy functional Boese and Handy (2001) BLYP Becke exchange plus Lee-Yang-Parr correlation Becke (1988), Lee et al. (1988) BOP Becke One Parameter functional Tsuneda et al. (1999) VWN-BP BP functional with the l
35、ocal correlation replaced by the VWN functional. Vosko et al. (1980), Becke (1988), Perdew and Wang (1992) VWN: The Vosko-Wilk-Nusair (VWN) functional is the most popular LSD correlation potential. It uses a fit to accurate numerical results (by Ceperly and Alder) of a uniform electron gas. Ceperley
36、 and Alder performed quantum Monte Carlo calculations on a uniform electron gas at low and high spin limits for several electron densities. VWN uses the Pade interpolation procedure to fit the CA results for both the para and ferro states and for low and high densities. DMol3 uses the best VWN (so c
37、alled “Fit“) parameters.PWC: The Perdew-Wang (PWC) functional is a recent parameterization of the Ceperley and Alder data, which corrects some VWN problems with fitting. PWC is the default functional for DMol3 calculations.The local spin-density (LSD) approximation accurately predicts structures, vi
38、brations, and relative energies of covalent systems; however, bond energies are seriously overestimated. The local DFT should not be used for systems with weak bonds, such as hydrogen bonds. These problems with the LSD method can be corrected to a large extent by using the so-called gradient-correct
39、ed (or nonlocal) functionals.P91, BP, BLYP, BOP: DMol3 supports several nonlocal exchange and correlation functionals. The most popular, the Becke exchange functional (B88) is used in conjunction with the Perdew-Wang correlation functional (BP) or the Lee-Yang-Parr correlation functional (BLYP). The
40、 so-called generalized gradient corrected (GGA) functional, by Perdew and Wang (P91) was derived by considering low and high density regimes and by enforcing various summation rules.PBE: The PBE (Perdew, Burke and Enzerhof) functional (1996, 1997) is another GGA functional in which all the parameter
41、s (other then those in its LDA component) are fundamental constants. The exchange part of this functional is similar to the Becke formula (1986), and the correlation part is close to the Perdew-Wang functional (1986).This functional has a strong physical background, reliable numerical performance an
42、d it is frequently used in DFT calculations.RPBE: More recently, Hammer, Hansen and Norskov (1999) proposed a modified version of the PBE formula that improves considerably thermochemical results. So far this functional, called RPBE, has been mainly used for solids. HCTH: The HCTH functional, named
43、for the authors initials (Hamprecht et al., 1998), represents a “pragmatic“ philosophy in designing a DFT functional. Assuming that the exact functional will never be found, they propose a flexible form of gradient corrected functional that is fitted to the training set of molecules. The default for
44、 the current implementation, the so called HCTH/407 functional, that was obtained by fitting to the set of the 407 atomic and molecular systems (Boese and Handy, 2001). This functional was found to predict a much improved thermochemistry for inorganic and hydrogen bonded systems. The standalone keyw
45、ord for this functional is hcth407. The original HCTH-type functionals can be activated by keyword hcth93 (Hamprecht et al, 1998) and hcth147 (Boese et al. 2000). Although the NLSD methods are significantly better than the LSD method, particularly in studying chemical reactions, the NLSD methods may
46、 still lead to reaction barriers that are too low.VWN-BP: The VWN-BP functional is recommended for COSMO and COSMO-RS studies.Materials-Studio 论坛问答全集（精选众多论坛讨论贴）1、问：用 MS 构造晶体时要先确立空间群,可是那些空间群的代码是啥意思啊,看不懂,我想做的是聚乙烯醇的晶体,嘿嘿,也不知道去哪可以查到它的空间群 答：A、要做晶体，首先要查询晶体数据，然后利用晶体数据再建立模型。晶体数据来源主要是文献，或者一些数据库，比如 CCDC。你都不知道
47、这个晶体是怎么样的，怎么指定空间群呢？要反过来做事情哦：）B、我不知道你指示的代码是数字代码还是字母代码,数字代码它对应了字母的代码，而字母的代码它含盖了一些群论的知识(晶系,对称操作等),如果要具体了解你的物质或者材料属于那一个群，你可以查阅一下相关的手册,当然你要了解一些基本的群论知识.MS 自带了一些材料的晶体结构,你可以查询一下. 2、问：各位高手，我用 ms 中的 castep 进行运算。无论 cpu 是几个核心，它只有一个核心在工作。这个怎么解决呢？答：请先确认以下几个问题：1，在什么系统下装，是否装了并行版本。2，计算时设置参数的地方是否选择了并行。3，程序运算时，并不是时时刻刻
48、都要用到多个 CPU 3、问：我已经成功地安装了 MS3.1 的 Linux 版本，串行的 DMol3 可以成功运行。但是运行并行的时候出错。机器是双 Xeon5320(四核)服务器，rsh 和 rlogin 均开启，RHEL4.6 系统。其中 hosts.equiv 的内容如下：localhostibm-consolemachines.LINUX 的内容如下：localhost:8现在运行 RunDMol3.sh 时，脚本停在$MS_INSTALL_ROOT/MPICH/bin/mpirun $nolocal -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.exelocalhost 3 /home/www/MSI/MS3
