24小时热门版块排行榜

>论坛更新日志 (2751)
>虫友互识 (259)
>文献求助 (180)
>导师招生 (157)
>考博 (111)
>休闲灌水 (71)
>论文投稿 (52)
>论文道贺祈福 (45)
>硕博家园 (37)
>博后之家 (35)
>基金申请 (33)
>教师之家 (30)
>公派出国 (23)
>招聘信息布告栏 (22)
>考研 (21)
>找工作 (20)

返回列表

戏剧

新虫 (正式写手)

应助: 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

回复此楼

» 猜你喜欢

存款400万可以在学校里躺平吗已经有26人回复
拟解决的关键科学问题还要不要写已经有8人回复
最失望的一年已经有7人回复
推荐一本书已经有16人回复
国自然申请面上模板最新2026版出了吗？已经有20人回复
26申博已经有3人回复
请教限项目规定已经有3人回复
基金委咋了？2026年的指南还没有出来？已经有10人回复
基金申报已经有6人回复
疑惑？已经有5人回复

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

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

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

返回列表