软件学习
  • gromacs-mdp-Temperature coupling

    tau-t 是每个分子随机化之间的平均时间。目前仅使用 velocity Verlet 实现。对于 velocity Verlet 积分器,nsttcouple 设置为 1。nh-chain-length : 对于 velocity Verlet 积分器,integrator=md 仅支持 1。默认情况下,NH链变量的数据不会输入到 edr 文件,但可以使用 print-nose-hoover-chain-variables 选项打开。print-nose-hoover-chain-variables:no: 不在能量文件中存储Nose-Hoover链变量。

    2023-03-14 0 0 0 0
  • gromacs-mdp-Ewald

    fourier-nx/fourier-ny/fourier-nz: 使用 Ewald 时 reciprocal space 中波矢量的最大幅值,或使用 PME 或 P3M 时的网格大小。请注意,通过 gmx mdrun 中的 PME tuning,这些网格大小可以随着 rcoulomb 的扩展而减小。ewald-rtol : ewald-rtol 给出了 rcoulomb 处 Ewald-shifted direct potential 的相对强度。ewald-geometry:3d: 在所有三个维度中执行Ewald求和。epsilon-surface : 这控制了对三维 Ewald 求和的偶极子校正。

    2023-03-14 0 0 0 0
  • gromacs-mdp-Tables

    table-extension [nm]: 将非键合势能查找表扩展到最大截止距离之外。当 user table 用于 静电 和/或 VdW 时,这里可以给出应使用单独用户表的成对能量组。这两个能量组将按照它们在 energygrps 中的定义顺序附加到表文件名,并用下划线分隔。例如,如果 energygrps = Na Cl Sol,energygrp-table = Na Na Na Cl,gmx mdrun 将读取table_Na_Na.xvg 和 table_Na_Cl.xvg 以及将用于所有其他能量组对的 table.xvg。

    2023-03-14 0 0 0 0
  • gromacs-mdp-VdW

    vdwtype:Cut-off: 具有成对列表半径 rlist 和 VdW 截止 rvdw 的平面截止,其中 rlist=rvdw。网格尺寸以与静电相同的方式用 fourierspacing 控制,插值顺序用 pme-order 控制。Shift: 2022.2中已被弃用,并被使用 vdw modifier=Force-switch 的 vdwtype=Cut-off 替代。LJ电势在全范围内被降低,并且力在 rvdw-switch 与 rvdw 之间平滑衰减至零。LJ电势正常输出至 rvdw-switch,之后在 rvdw 处切换至零。Potential-shift-Verlet:仅对能量应用长程色散校正。

    2023-03-14 0 0 0 0
  • gromacs-mdp-Electrostatics

    实空间截止 rcoulomb 应等于 rlist。在 reciprocal space 中使用的波矢量的最高幅度由 fourierspacing 控制。direct/reciprocal space 的相对精度由 ewald rtol 控制。Direct space 类似于Ewald sum,而 reciprocal space 使用 FFT 执行。库仑截止为 rcoulomb 的广义反应场,其中 rcoulumb≥rlist。Reaction-Field-zero 通过使电势零点超过截止值来解决这个问题。切换库仑电势可能会导致严重的伪影,建议使用 Reaction-Field-zero。coulomb-modifier:Potential-shift: 将库仑电势增减一个常数,使其在截止点处为零。

    2023-03-14 0 0 0 0
  • gromacs-mdp-Neighbor searching

    当设置 dynamics 和 verlet-buffer-tolerance 时,nstlist 实际上是一个最小值,gmx mdrun 可能自动增加,除非将其设置为1。这仅适用于截止值小于 rlist 的 switched, shifted or user potentials。若要模拟无截止,请将所有 cut-offs 相关选项和 nstlist 设置为 0。periodic-molecules:no: 分子是有限的,可以使用 fast molecular PBC。应当就是指 cutoff-scheme=Verlet 的情况,这一点在 5.0.4 中指出)。请注意,PP-PME 负载平衡可能会自动启用双程力评估。

    2023-03-14 0 0 0 0
  • gromacs-mdp-Output control

    nstvout [steps]: 将速度写入输出轨迹文件的间隔步数,其余同上。该值应为 nstcalcenergy 的倍数。请注意,每一步 MD step 的精确和以及波动模 nstcalcenergy 的值都存储在能量文件中,因此当 nstenergy1 时,gmx energy 也可以报告精确的能量平均值和波动。compressed-x-precision [real]: 写入压缩轨迹文件的精度。compressed-x-grps 写入压缩轨迹文件的索引组,默认情况下写入整个系统。

    2023-03-14 0 0 0 0
  • gromacs-mdp-Test particle insertion

    rtpi [nm]:测试粒子插入半径,请参见 integrator=tpi 及 integrator=tpic。

    2023-03-14 0 0 0 0
  • gromacs-mdp-Shell molecular dynamics

    emtol [kJ mol-1 nm-1]: 当最大力小于 emtol 时,最小化收敛。fcstep [ps2]: 用于优化柔性约束的步长。设置时应当使 fcstep=mu/,其中mu是柔性约束中两个粒子的缩减质量,d2V/dq2 是约束方向上电势的二阶导数。幸运的是这一数字在不同柔性约束之间不会发生太大变动,尽管迭代次数和运行时长对 fcstep 非常敏感。

    2023-03-14 0 0 0 0
  • gromacs-mdp-Energy minimization

    emtol [kJ mol-1 nm-1]: 当最大力小于该值时,表明最小化收敛。nstcgsteep [steps]: 使用共轭梯度能量最小化时执行1个最陡下降步骤的频率。nbfgscorr : 用于 l-bfgs 最小化的校正步骤数。更高的数字更准确,但速度较慢。

    2023-03-14 0 0 0 0
技术支持: CLOUD | 管理登录
seo seo