求助，tcl脚本加入重新判断连接关系语句后getbonds仍然不变

哒哒哒哒w

大家好，我在处理一段分子动力学数据，想使用tcl脚本得到每一帧键的连接情况：

1. set selnow [atomselect top "name O"]
   
2. for {set i 0} {$i <= 50000} {incr i} {
   
3. $selnow frame $i
   
4. mol bondsrecalc all
   
5. topo retypebonds
   
6. set bonds [$selnow getbonds]
   
7. puts "$i $bonds"
   
8. }

复制代码

但输出的结果每一帧的连接情况都相同：

1. Info) Determining bond structure from distance search ...
   
2. Info) Bond count: 113
   
3. 0 0 0 0 0 7 7 7 7 14 14 14 14 21 21 21 21 28 28 28 28 {7 35} 35 35 35 42 42 42 42 49 49 49 49 56 56 56 56 63 63 63 63 70 70 70 70 77 77 77 77 84 84 84 84 91 91 91 91 98 98 98 98 105 105 105 110 110 110 110 117 117 117 117 124 124 124 124 131 131 131 131 138 138 138 138 145 145 145 145 152 152 152 152 159 159 159 159 166 166 166 166 173 173 173 173 180 180 180 {105 180} 187 187 187 187
   
4. Info) Determining bond structure from distance search ...
   
5. Info) Bond count: 113
   
6. 1 0 0 0 0 7 7 7 7 14 14 14 14 21 21 21 21 28 28 28 28 {7 35} 35 35 35 42 42 42 42 49 49 49 49 56 56 56 56 63 63 63 63 70 70 70 70 77 77 77 77 84 84 84 84 91 91 91 91 98 98 98 98 105 105 105 110 110 110 110 117 117 117 117 124 124 124 124 131 131 131 131 138 138 138 138 145 145 145 145 152 152 152 152 159 159 159 159 166 166 166 166 173 173 173 173 180 180 180 {105 180} 187 187 187 187
   
7. Info) Determining bond structure from distance search ...
   
8. Info) Bond count: 113
   
9. 2 0 0 0 0 7 7 7 7 14 14 14 14 21 21 21 21 28 28 28 28 {7 35} 35 35 35 42 42 42 42 49 49 49 49 56 56 56 56 63 63 63 63 70 70 70 70 77 77 77 77 84 84 84 84 91 91 91 91 98 98 98 98 105 105 105 110 110 110 110 117 117 117 117 124 124 124 124 131 131 131 131 138 138 138 138 145 145 145 145 152 152 152 152 159 159 159 159 166 166 166 166 173 173 173 173 180 180 180 {105 180} 187 187 187 187
   
10. Info) Determining bond structure from distance search ...
    
11. Info) Bond count: 113
    
12. 3 0 0 0 0 7 7 7 7 14 14 14 14 21 21 21 21 28 28 28 28 {7 35} 35 35 35 42 42 42 42 49 49 49 49 56 56 56 56 63 63 63 63 70 70 70 70 77 77 77 77 84 84 84 84 91 91 91 91 98 98 98 98 105 105 105 110 110 110 110 117 117 117 117 124 124 124 124 131 131 131 131 138 138 138 138 145 145 145 145 152 152 152 152 159 159 159 159 166 166 166 166 173 173 173 173 180 180 180 {105 180} 187 187 187 187
    
13. Info) Determining bond structure from distance search ...
    
14. Info) Bond count: 113

复制代码

我肉眼观察时，dynamic bonds绝对是有变化的，请大家告诉我问题出在哪谢谢！

sobereva

加上mol reanalyze top试试

哒哒哒哒w

sobereva 发表于 2024-3-12 09:25
加上mol reanalyze top试试

感谢sob老师！成功解决了