24小时热门版块排行榜

返回列表

戏剧

新虫 (正式写手)

应助: 0 (幼儿园)
金币: 3146.2
帖子: 562
在线: 129.6小时
虫号: 3570184
注册: 2014-12-01
专业: 理论和计算化学

[求助] gromacs用于FEP计算问题，求高手指点迷津

本人在用gromacs进行自由能扰动（FEP）模拟时候，得出的结果居然高达-400KBT（主要是Coulombic贡献的），这条多肽仅有9个氨基酸。核心参数如下，请高手指点一下。
; Free energy control stuff
free_energy             = yes
init_lambda_state       = LAMBDA_STATE
delta_lambda          = 0
calc_lambda_neighbors = 1       ; only immediate neighboring windows
couple-moltype          = Protein_chain_C
; Vectors of lambda specified here
; Each combination is an index that is retrieved from init_lambda_state for each simulation
; init_lambda_state       0    1    2    3    4    5    6    7    8    9    10    11    12    13    14    15    16    17    18    19    20    21    22    23    24    25    26    27    28    29    30    31    32    33    34    35    36    37    38    39    49    41    42    43    44    45    46    47    48    49    50    51    52    53    54    55    56    57    58    59    60    61    62    63
vdw-lambdas             = 0.00000 0.00100 0.01000 0.05000 0.10000 0.20000 0.30000 0.40000 0.50000 0.60000 0.70000 0.80000 0.90000 0.95000 0.99000 0.99900 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 ;1.00000 1.00000 1.00000 1.00000  1.00000 ; 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000  1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000  1.00000
coul-lambdas          = 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00100 0.01000 0.05000 0.10000 0.20000 0.30000 0.40000 0.50000 0.60000 0.70000 0.80000 0.90000 0.95000 0.99000 0.99900 1.00000 ; 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000  1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000  1.00000
;restraint-lambdas    = 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.000000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.000000 0.00000 ; 0.00001 0.00010 0.00100 0.01000 0.05000 0.10000 0.20000 0.30000 0.40000  0.50000 0.60000 0.70000 0.80000 0.90000 0.95000 0.99000 0.99900 0.99990 0.99999 0.999999 1.00000
couple-lambda0          = none
couple-lambda1          = vdw-q

; Options for the decoupling
sc-alpha                = 0.5
sc-coul                = no    ; linear interpolation of Coulomb (none in this case)
sc-r-power             = 6
sc-sigma                = 0.3
nstdhdl                = 100

ld_seed=-1

回复此楼

» 猜你喜欢

1楼 2019-08-09 19:57:12

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主戏剧的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 304求调剂 +4	ahbd 2026-03-14	4/200	2026-03-16 16:48 by 我的船我的海
[考研] 085600调剂 +5	漾漾123sun 2026-03-12	6/300	2026-03-16 15:58 by 漾漾123sun
[考研] 070303一志愿西北大学学硕310找调剂 +5	d如愿上岸 2026-03-12	8/400	2026-03-16 15:19 by peike
[考研] 0856求调剂 +3	刘梦微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[基金申请] NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗 20+5	NSFC2026我来了 2026-03-10	14/700	2026-03-15 23:53 by 不负韶华的虎
[考研] 求老师收留调剂 +4	jiang姜66 2026-03-14	5/250	2026-03-15 20:11 by Winj1e
[考研] 机械专硕调剂 +3	笨笨兔子 2026-03-12	3/150	2026-03-15 20:02 by 栗子粥?
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3	嘉年新程 2026-03-15	3/150	2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 0856专硕279求调剂 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015
[考研] 289求调剂 +5	步川酷紫123 2026-03-11	5/250	2026-03-15 00:45 by kruisytel
[考研] 268求调剂 +5	一定有学上- 2026-03-14	6/300	2026-03-14 22:20 by 运气yunqi
[基金申请] 面上和青基一样限30页不合理 +5	wowsunflower 2026-03-10	7/350	2026-03-14 17:21 by kingkocxr
[考研] 328求调剂 +3	5201314Lsy！ 2026-03-13	6/300	2026-03-14 15:31 by hyswxzs
[考研] 290求调剂 +4	@将就将就看 2026-03-10	8/400	2026-03-14 14:23 by 千千运气
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5	Becho359 2026-03-11	5/250	2026-03-14 11:35 by 哦哦123
[考研] 学硕285求调剂 +13	Wisjxn 2026-03-12	46/2300	2026-03-14 10:33 by JourneyLucky
[考研] 0703化学调剂 +4	快乐的香蕉 2026-03-11	4/200	2026-03-13 22:41 by JourneyLucky
[考研] 求调剂（材料与化工327） +4	爱吃香菜啦 2026-03-11	4/200	2026-03-13 22:11 by JourneyLucky
[考研] 0856化学工程280分求调剂 +4	shenzxsn 2026-03-11	4/200	2026-03-13 11:55 by ymwdoctor
[考研] 一志愿：武汉理工，材料工程，英二数二总分314 +3	2202020125 2026-03-10	4/200	2026-03-10 13:54 by xiongyaxuan