计算化学公社

标题: 使用VMD构建上下两层板的channel结构 [打印本页]
作者
Author: tjuptz    时间: 2019-8-21 12:52
标题: 使用VMD构建上下两层板的channel结构
本帖最后由 tjuptz 于 2019-10-17 21:59 编辑

前面我整理了论坛中关于固体表面cluster/slab模型的问题讨论以及界面建模考虑表面弛豫或重构时考虑的层数问题讨论。MD中建立表面模型，可以用MS切面扩胞加真空管后直接使用，这样固体表面在盒子最底部；也可以在GROMACS中用editconf命令对盒子尺寸调整，默认情况下会把表面挪到盒子中间，这样适用于某些体系。有时候我们想建立类似三明治结构，两层表面模型中间真空用来填充一些分子，比如研究流体输运或材料剪切问题时。http://bbs.keinsci.com/thread-13806-1-1.html中有人问这种模型怎么建，我尝试写了个VMD脚本，如下：
  1. display rendermode GLSL
  2. display depthcue off
  3. display projection Orthographic

  5. mol new 1.gro
  6. pbc box
  7. rotate x by -90
  8. mol modstyle 0 top CPK 1.000000 0.300000 12.000000 12.000000

  10. puts "Use like 'transform 30 x 2'. The first value is the distance bewteen two slabs and the last value is for mode selection. 'x' and 'y' representing X and Y axis are equivalent in isotropic condition!"
  11. puts "The built model would be saved automatically into 2.gro. You should copy the coordinates from 10.gro to 2.gro manually!"
  12. puts "CAUTION:cut mode without manual operation after cut!"
  13. puts "CAUTION:split mode without manual operation after split!"

  15. proc transform {d {axis x} {isoXY 2}} {
  16. set dis {0 0}
  17. lappend dis $d
  18. set sel [atomselect top all]
  19. $sel writegro 10.gro
  20. set vec [measure center $sel]
  21. if {$isoXY==1} {
  22. #inversion mode:isotropic for X/Y direction
  23.         $sel moveby [vecinvert $vec]
  24.         $sel move [transaxis $axis 180]
  25.         $sel moveby $vec
  26.         $sel moveby $dis
  27.         $sel writegro 2.gro
  28. } elseif {$isoXY==0} {
  29. #inversion mode:anisotropic for X/Y direction
  30.         $sel moveby [vecinvert $vec]
  31.         $sel move [transaxis $axis 180]
  32.         $sel move [transaxis z 180]
  33.         $sel moveby $vec
  34.         $sel moveby $dis
  35.         $sel writegro 2.gro
  36. } else {
  37. #translation mode
  38.         $sel moveby $dis
  39.         $sel writegro 2.gro
  40. }
  41. }

  43. proc cut {d} {
  44. #cut mode
  45. set sel [atomselect top all]
  46. lassign [lindex [measure minmax $sel] 1] a b c
  47. set range "z<=($c-$d)/2 or z>=($c+$d)/2"
  48. mol modselect 0 top $range
  49. set model [atomselect top $range]
  50. $model writegro 2.gro
  51. }

  53. proc split {d} {
  54. #split mode
  55. set sel [atomselect top all]
  56. lassign [lindex [measure minmax $sel] 1] a b c
  57. set dis {0 0}
  58. lappend dis $d
  59. set range "z>=$c/2"
  60. set upper [atomselect top $range]
  61. $upper moveby $dis
  62. $sel writegro 2.gro
  63. }
复制代码

把脚本和改名为1.gro的原始模型放到VMD目录下使用。切完除了窗口显示以外，会自动保存到VMD自动目录下为2.gro。10.gro是原结构的副本，用于使用transfrom时与2.gro合并。

固体有不同情况，1.点阵点是单个原子的晶体如金属、单原子离子盐，2.点阵点是原子但共价连接的晶体如金刚石、硅，3.点阵点包含多个分子片段的晶体如碳酸盐、分子晶体，4.无规堆垛、无定形或短程有序长程无序的片层。该脚本给出了三个命令，功能上存在冗余，我是想到了多种情况就都写进来了，可认为适用于不同目的。split命令是在MS建好的带真空模型上把底层slab劈开并将上半部分挪到上层；cut命令是把建好的不带真空的厚厚slab切掉中间层；transform命令则是把底部slab变换成顶部slab，又分了几种情况：直接移动对应isoXY不等于0或1，倒置后再移动又分为XY向同性isoXY=1与异性isoXY=0。
注意：对于晶体而言，晶层经过重复总会有一个周期，所以不同命令可能建的模一样；对2.情况，使用cut和split需要恰当处理切除原子造成的悬挂键问题；对3.情况，使用cut或split，适当添加“same resid as……”语句保证片段完整性。MS建模时，Z尺寸应为d+2c，d是真空厚度，c是slab厚度。

用什么命令还是要看你想模拟什么情景。比如我对建好的底部slab最表面进行了嫁接改性，为了上表面也是如此而且不用再做一遍，用transform就好了，至于XY异性时到底让上下板呈镜面还是倒反另考虑。

PS：tcl语言对空格太敏感了，格式困扰了好久





欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3