为什么我用dmol3优化晶面模型后再复制分子过去上层原子会消失

y1p1p1

有没有大佬知道，我用没优化的晶面试过，复制过去不会消失，原图是bmp的，不会转jpg

卡开发发

模型的display style尽可能从一开始就用original模式试试。暂时的结局方案可以用redefine lattice直接填个单位矩阵进去，这样试试看。

y1p1p1

卡开发发 发表于 2021-11-5 22:05
模型的display style尽可能从一开始就用original模式试试。暂时的结局方案可以用redefine lattice直接填个 ...

大佬你好，我刚刚都切换成了original，果然原型是一模一样的，可是单位矩阵是啥意思，我不懂呀

y1p1p1

卡开发发 发表于 2021-11-5 22:05
模型的display style尽可能从一开始就用original模式试试。暂时的结局方案可以用redefine lattice直接填个 ...

我是顺着师兄的论文开始先对晶胞进行优化，切了晶面加了真空层扩胞之后再次优化不收敛，没办法只好加了smearing，一开始勾选optimize cell的时候原子层融合到一起去了，复制过去分子原子层倒是不消失，可是也没法儿进行吸附，后来我就没勾选optimize cell，然后复制过去上层原子就消失了

y1p1p1

卡开发发 发表于 2021-11-5 22:05
模型的display style尽可能从一开始就用original模式试试。暂时的结局方案可以用redefine lattice直接填个 ...

我主要是想复现一下师兄的论文内容，师兄那上面倒是没有出现原子层消失的现象

y1p1p1

卡开发发 发表于 2021-11-5 22:05
模型的display style尽可能从一开始就用original模式试试。暂时的结局方案可以用redefine lattice直接填个 ...

大佬这个redifine lattice我刚刚对原晶面和优化后的晶面都试过了，还是不行呀，用了之后上面的原子层就消失了

卡开发发

y1p1p1 发表于 2021-11-6 16:47
大佬这个redifine lattice我刚刚对原晶面和优化后的晶面都试过了，还是不行呀，用了之后上面的原子层就消 ...

应该是你复制之后有些原子可能因为显示模式的原因所以不显示了，所以才建议你一开始就用original的模式。redefine lattice里面的选项默认填的就是单位矩阵，也就是((1,0,0),(0,1,0),(0,0,1))。

其他的问题我看不到你前后操作的图片，所以给不了更多帮助。

y1p1p1

卡开发发 发表于 2021-11-6 17:10
应该是你复制之后有些原子可能因为显示模式的原因所以不显示了，所以才建议你一开始就用original的模式。 ...

这样的嘛，我看看调整显示风格能不能变回来

卡开发发

y1p1p1 发表于 2021-11-6 17:18
这样的嘛，我看看调整显示风格能不能变回来

redefine lattice的操作简单说下原理，道理上一般每次做redefine lattice最后一步会将全部的原子的分数坐标对1取余，保证全部的原子全都落在晶格内，当填写单位矩阵进去的时候晶格矢量其实还是保持原样，但是原子就全部按照前面说的方式分数坐标被调整成0～1的范围。

y1p1p1

卡开发发 发表于 2021-11-6 17:22
redefine lattice的操作简单说下原理，道理上一般每次做redefine lattice最后一步会将全部的原子的分数坐 ...

incell和默认的都试过了，还没优化的晶面复制分子没事，优化之后的原子消失了一层，本来我以为可能是因为没有选上optimize cell，不过我上周真空层加大再优化把分子复制过去准备做吸附外层原子也消失了，并且删除氢气分子也没变回来，对于这点我猜可能是氢气复制过去显示默认变成了Incell。我刚刚先对初始晶面进行了redifine lattice，再对优化后的晶面模型进行redifine lattice的时候外层原子就消失了；又只对优化后的晶面模型进行redifine lattice，结果用完外层原子还是直接消失了，就不明白是为啥了

y1p1p1

卡开发发 发表于 2021-11-6 17:22
redefine lattice的操作简单说下原理，道理上一般每次做redefine lattice最后一步会将全部的原子的分数坐 ...

就像这样

卡开发发

y1p1p1 发表于 2021-11-6 18:21
就像这样

少一层的原因是一开始你用default显示，那么边界上的原子就在两侧重复显示了，但是实际上并没有这一层，这也是为什么一开始强调建议用oirigin模式的原因。

其他的方面， 你的层高太薄了，真空层也不够大，dmol3的算法计算量几乎与真空层尺寸无关，所以对表面来说放心地加。

y1p1p1

卡开发发 发表于 2021-11-6 18:30
少一层的原因是一开始你用default显示，那么边界上的原子就在两侧重复显示了，但是实际上并没有这一层， ...

师兄论文里我看是有这一层的，这是怎么做到的呀

y1p1p1

卡开发发 发表于 2021-11-6 18:30
少一层的原因是一开始你用default显示，那么边界上的原子就在两侧重复显示了，但是实际上并没有这一层， ...

还有就是他这个吸附模拟是怎么进行的，看起来好像不是直接复制分子过去，我之前按照教程复制分子然后用forcite进行动力学模拟也是只有一个分子的吸附，和我想象中的大量分子的吸附不太一样，还有asorption来进行动力学模拟好像可以通过选项直接把分子添加到晶面上，可是接下来的蒙特卡罗模拟就不知道怎么做了

卡开发发

y1p1p1 发表于 2021-11-6 19:17
师兄论文里我看是有这一层的，这是怎么做到的呀

切面的时候thickness=1.5试试，但是不建议这样，最好按整数来。

如果是一般的第一原理计算或者力场模拟按道理是把分子复制过去就行。adsorpition应该是各自存在一个xsd文件，然后adsorption似乎可以指定，但我不做这块不大了解，你可以自己摸索或者问更了解这块的人。