本帖最后由 Lacrimosa 于 2023-6-15 09:50 编辑 fightingllll 发表于 2023-6-14 10:54 用这个贴子里的代码肯定实现不了,而且VMD做这个有可能会很慢。建议自己编程实现或者找一些有类似分析功能python包。自编代码的话,首先你得明确水团簇的判定条件是什么,然后找出每一帧中的所有水团簇,输出该团簇的水分子数量,最后按数量和频率做直方图就行了。编程方面有困难可以求助chatGPT |
Lacrimosa 发表于 2022-11-7 18:26 老师你好 我想通过VMD获得在水溶液中形成水团簇的数量 以及不同水团簇尺寸的百分比 可以实现吗 |
202306141054475700..png (35.76 KB, 下载次数 Times of downloads: 52)
Lacrimosa 发表于 2022-11-7 18:26 好的,十分感谢 |
YPL 发表于 2022-11-6 17:58 会直接在tkconsole里输出结果,如果你需要输出到文件用下面这三个命令: set outfile [open ./out.dat w] puts $outfile "Output data" close $outfile |
本帖最后由 YPL 于 2022-11-6 20:04 编辑 Lacrimosa 发表于 2020-11-4 18:22 您好,老师。我刚学VMD,这个在TK console中运行完上述语句之后,输出文件在哪儿呢? |
dongdong 发表于 2021-4-24 17:33 measure $sel center weight mass是求质心 measure $sel center是几何中心 |
Lacrimosa 发表于 2020-11-4 18:22 老师您好: 如果想求水团簇的几何中心应该如何改这个代码啊? |
Lacrimosa 发表于 2020-12-10 12:49 感谢感谢 |
elpa 发表于 2020-12-9 23:30 不是文件名,而是vmd中的selection 参考:http://sobereva.com/504 |
dongdong 发表于 2020-11-4 18:59 问一个愚蠢的问题。请问“将water替换为你自己的选择语句”,这里water是指自己的文件名还是? |
Lacrimosa 发表于 2020-11-4 18:22 哇,非常实用,谢谢 ![]() |
|
本帖最后由 Lacrimosa 于 2021-4-24 19:41 编辑 打开VMD Main-Extensions-Tk Console输入以下内容(set select "water"此处将water替换为你自己的选择语句) #--------------------------------------------------- set outfile [open center.dat w] set select "water" #--------------------------------------------------- set n [molinfo top get numframes] set sel [atomselect top "$select"] for { set i 1 } { $i <= $n } { incr i } { set sel [atomselect top "$select" frame $i] set p [measure center "$sel" weight mass] puts $outfile $p } close $outfile |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-25 16:44 , Processed in 0.203184 second(s), 25 queries , Gzip On.