1、使用IRI函数同时考察反应过程的共价和非共价相互作用,相应的cub文件已经使用Multiwfn产生好了。 2、在VMD目录下写个文本文件isoall.tcl,准备用VMD程序运行isoall.tcl来产生每个cube文件对应的等值面图像文件,但运行后只能产生一帧。
3、VMD错误提示见图片,isoall.tcl文件如下:
color Display Background white
color scale method BGR
set isoval 0.5
axes location Off
for {set i 1} {$i<=114} {incr i} {
set name IRC[format %04d $i]
puts "Processing f1_$name.cub and f2_$name.cub..."
mol default style CPK
mol new f1_$name.cub
mol addfile f2_$name.cub
mol representation CPK 1.0 0.3 18.0 16.0
mol addrep top
mol representation Isosurface 1.0 1 0 0 1 1
mol color Volume 0
mol addrep top
mol scaleminmax top 1 -0.04 0.02
mol modstyle 0 0 CPK 0.700000 0.300000 18.000000 16.000000
color scale midpoint 0.666
color scale method BGR
color Display Background white
axes location Off
display depthcue off
display rendermode GLSL
light 3 on
color Element N iceblue
mol modcolor 0 0 Element