| 查看: 4817 | 回复: 11 | |||||
| 本帖产生 1 个 模拟EPI ，点击这里进行查看 | |||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||
robertpeter2005木虫 (正式写手)
|
[交流]
关于熵变的计算，有人知道怎么计算大的体系的熵变吗？比如将近一千个残基的体系？已有6人参与
|
||||
|
关于熵变的计算，有人知道怎么计算大的体系的熵变吗？比如将近一千个残基的体系？原来在分子模拟论坛上看到有人说，可以减小大体系的尺寸来计算熵变，找到一个人问了一下，结果人家不愿意说，没空搭理我，郁闷，我有个体系急需算熵来分析构象的变化，请问大家有谁知道如何计算大体系的熵变啊？如果可以减小体系尺寸到300个残基以内，那具体该如何操作啊？我也没有见到过这样的文献，大家有见过的吗？听说可以把轨迹文件里面不相关的残基删除再保存，然后用来算熵，是用什么软件啊？vmd?这个问题困扰好长时间了，还望各位高人能帮忙想想办法，在此先谢谢了！！！ 另外有用ONIOM做QM/MM的人吗？我做了一个输入文件，但是计算的时候显示硫原子没有参数？进行不下去，但是蛋白质里面都会有硫元素的啊，蛋白两端的残基是怎么处理的呢？ |
回复此楼» 收录本帖的淘贴专辑推荐
 模型 - 模拟 - 统计 - 热力学 | 精华 |
» 猜你喜欢
连续两年医学口青年项目初审体会
已经有7人回复
-
蹲一个6月7号江苏省青基公示
已经有32人回复
-
关于2023的收入
已经有3人回复
-
求助
已经有8人回复
-
函评意见：既然你有中科院****人才项目了，就不要申请面上了。
已经有34人回复
-
给大家分享一下扎心的投稿经历
已经有11人回复
-
迷茫毕业季
已经有8人回复
-
今年香江初审出结果了吗
已经有14人回复
-
审稿交流
已经有7人回复
-
迷茫的毕业季
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- vasp 计算单原子能量失败！请教各位！
已经有8人回复
- 异氰酸酯
已经有1人回复
- 生产车间管理建议
已经有6人回复
- 有人能提供各种碳纤维的泊松比吗？
已经有6人回复
- 学习第一性原理计算的准备
已经有13人回复
- vasp计算能带，刚开始就服务器就卡了
已经有20人回复
- 难道每次计算前都咬进行收敛性测试吗？
已经有14人回复
- 两道物化题目求熵变和一个热力学证明
已经有16人回复
- 在-5度，p标准下苯的凝固过程的熵变【认真看过教材的同学来】
已经有9人回复
- 理想气体自由膨胀的熵变意味着什么？
已经有18人回复
- 【求助】用Gaussian计算分子的熵时，不同分子的熵的计算基准是一样的吗？谢谢
已经有11人回复
- 【求助】绝对熵值的计算
已经有8人回复
- 【求助】一道关于熵变计算
已经有7人回复
- 哪位大侠能告诉我如何理解熵稳定和焓稳定啊？
已经有10人回复
ChemiAndy
木虫 (正式写手)
- 模拟EPI: 40
- 应助: 52 (初中生)
- 金币: 2573.3
- 散金: 18
- 红花: 104
- 帖子: 590
- 在线: 553.6小时
- 虫号: 836907
- 注册: 2009-08-30
- 性别: GG
- 专业: 理论和计算化学
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
御剑江湖(金币+10, 模拟EPI+1): 2011-11-13 20:09:42
小木虫(金币+0.5):给个红包，谢谢回帖
御剑江湖(金币+10, 模拟EPI+1): 2011-11-13 20:09:42
|
学而不思则罔啊 1. 算熵变与体系大小有关系吗？ 2. 熵变必然包括了一个过程，你是什么变化过程导致构象数目变化的呢？ 3. 熵不是热力学平均量，是无法通过模拟直接计算的； 我没算过，但我相信必须通过计算自由能变化来倒推熵变: ΔG=ΔH-TΔS ΔG可以通过自由能微扰获取。自由能微扰计算自由能差，找一本讲分子模拟自由能的教材可以看到。英文Wiki也能告诉我们怎么做: Zwanzig equation： http://en.wikipedia.org/wiki/Free_energy_perturbation In practice, one runs a normal simulation for state A, but each time a new configuration is accepted, the energy for state B is also computed. The difference between states A and B may be in the atom types involved, in which case the ΔG obtained is for "mutating" one molecule onto another, or it may be a difference of geometry, in which case one obtains a free energy map along one or more reaction coordinates. 人家说的去掉一个残基，应该就是上面说的mutating, 是为了得到state B的能量。你要算差300个残基的熵变，得一个一个来，分300步走。自由能按Zwanzig equation积分得到。不过，前后构象差一个残基的话，误差也太大了吧。一般H原子变甲基还得分成10步走呢。搞生物模拟的真心不在乎精度吗？ 积分完自由能的同时，积分每一步前后的焓变，然后TΔS=ΔH-ΔG. 如果初末两态差300个残基，体积差肯定明显，那么每一步都要跑NPT啊。则么 ΔH=ΔU+PΔV; ΔU是前后两态的内能差。 顺便说，GROMACS直接支持热力学微扰计算自由能，不需要手动去掉残基，只要设定好哪些原子参与mutating。不要问我怎么使用GROMACS软件。我从来没算过。俺只对热力学理论感兴趣。 [ Last edited by ChemiAndy on 2011-11-11 at 20:52 ] |
3楼2011-11-12 14:09:13
ChemiAndy
木虫 (正式写手)
- 模拟EPI: 40
- 应助: 52 (初中生)
- 金币: 2573.3
- 散金: 18
- 红花: 104
- 帖子: 590
- 在线: 553.6小时
- 虫号: 836907
- 注册: 2009-08-30
- 性别: GG
- 专业: 理论和计算化学
★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
御剑江湖(金币+5): 谢谢 2011-11-13 20:10:14
小木虫(金币+0.5):给个红包，谢谢回帖
御剑江湖(金币+5): 谢谢 2011-11-13 20:10:14
|
哦，原来是MM/PBSA方法的振动熵啊。恩，可以忽略我上面的回帖，那个是岁体系有较小的变化时严格计算自由能的。不适合体系有较大的变化。 据我所知，MM/PBSA估计熵的方法很简单，仅仅是振动熵，先把蛋白取出来做优化，然后正则分析频率+如下统计公式：  nj - frequencies of normal vibrations in the molecule h - Planck's constant k - Boltzmann's constant T - absolute temperature (K) 这种估算方法要用到正则分析 normal mode analysis，需要构建所有原子对所有原子的二维矩阵，对于大的体系当然就不行了。这就是为什么对超过300个残基不能用的原因。 希望对楼主有用。 [ Last edited by ChemiAndy on 2011-11-12 at 18:33 ] |
5楼2011-11-13 14:29:13
50