zoukankan      html  css  js  c++  java
  • charmm学习历程2

    最近几天setup已经成功了。说说mini、solvation和equilibrium

    minimization是能量最小化,运行之后相比于之前H的坐标发生改变,达到能量最小的状态。但是我的经历总是那么不顺。mini之后我的首尾出问题了,原子之间夹在一起。

    导师叫我改了一下inp:define heavy sele all .and. .not. ( bynu 1:28 .or. bynu 12547:12569 .or. resname TIP3 .or. hydrogen) end,bynu就是将残基不再受限,因此最后首尾两个残基就是好的了。

    不知道我的问题出在哪里,但是解决了就好

    solvation这一段比较机械,一般不会出问题,把自己底物上的一个原子设置成中心坐标(0,0,0)然后水就在其边缘形成一个圆球,droplet这个英文单词太形象了。
    coordinate translate xdirection -108.573 ydirection 38.918 zdirection -155.551  

    设置成相反数就好了。

    生成的文件检查一遍,出来的PDB实在好看,艺术家的魅力。

    solvation完成后改脚本inp,跑一遍equilibrium,我又出问题了,heat300K跑了一点就直接终止,我也不知道为啥,写的时候正在改inp,估计首位又得不固定。

    charmm论坛上的一个博主rmv实在太棒,我问一个问题他基本全部都摸清楚了。我以后也要成为这样的人吧。加油。希望这次equilibrium成功。

  • 相关阅读:
    k8s-[排查记录]解决节点无法查看pod日志
    k8s kube-proxy模式
    容器网络
    k8s-使用kubeadm安装集群
    k8s-Deployment重启方案
    k8s-NetworkPolicy-网络策略
    nodejs 解析终端特殊字符
    fluentd 日志自定义字段解析
    题目笔记 CF 1494b
    CF1225D Power Products(分解质因子 哈希)
  • 原文地址:https://www.cnblogs.com/dengwenhao/p/13130042.html
Copyright © 2011-2022 走看看