计算化学公社

 找回密码 Forget password
 注册 Register
Views: 139|回复 Reply: 0
打印 Print 上一主题 Last thread 下一主题 Next thread

[Lammps] 求助 在lammps中使用eABF方法计算自由能如何改善采样?

[复制链接 Copy URL]

1

帖子

0

威望

47

eV
积分
48

Level 2 能力者

各位老师们好,我想要获得纳米流体中两个纳米颗粒之间自由能与纳米颗粒距离之间的关系,盒子大小是60*60*120,我为了使纳米颗粒加速收敛使用fix spring/self 限制了他在xz方向上的运动,步长为1fs时算了2ns的时间后自由能还没收敛,在两个纳米颗粒靠的比较近时无法取样而变成平台一样的形状,想问下有什么方法可以改善取样,加速自由能的收敛?
自由能图片和colvars文件在下方:
#eABF方法对两个纳米氧化铝颗粒进行自由能采集

Colvarstrajfrequency 1000
Colvarsrestartfrequency 500000
indexFile non1.ndx
indexFile non2.ndx

colvar {
    name d
    width 0.2

    extendedlagrangian on
    extendedFluctuation 0.1
    extendedTemp 1200
    extendedTimeConstant 200

    subtractAppliedForce on

    distance {
        group1 {
            indexGroup non1
        }
        group2 {
            indexGroup non2
        }
        forceNoPBC yes
    }
    lowerBoundary 1
    upperBoundary 21
}
harmonicWalls {
    name mywall
    colvars d
    lowerWalls 1
    upperWalls 21
    forceConstant 20
}

abf {
    name abf1
    colvars d
    applyBias yes
    fullSamples 500
    outputEnergy on
    outputFreq 100000
    historyFreq 100000
    writeCZARwindowFile on
    CZARestimator on
}



74ba1de9be7a26955cee89814371d8a.png (6.76 KB, 下载次数 Times of downloads: 6)

74ba1de9be7a26955cee89814371d8a.png

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-14 04:13 , Processed in 0.153205 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list