• 分子动力学中步长的选取


    一般来说,时间步长应该设为分子运动的最小振动周期的1/10 左右为宜。在化学分子中振动周期最短的是碳氢键的伸缩振动,其特征振动频率大约是3000cm-1,其振动周期大约是10fs。所以通常的分子动力学时间步长最高不能高于1fs。

    http://blogold.chinaunix.net/u2/85704/showart_1884250.html


    The choice of timestep
    A key parameter in the integration algorithms is the integration timestep. To make the best use of the computer time, a large timestep should be used. However, too large a timestep causes instability and inaccuracy in the integration process.


    Relation of timestep to molecular vibration

    The timestep used depends on the model as well as the integrators. The main limitation imposed by the model is the highest-frequency motion that must be considered. A vibrational period must be split into at least 8-10 segments for models to satisfy the Verlet assumption that the velocities and accelerations are constant over the timestep used.

    In most organic models, the highest vibrational frequency is that of C-H bond stretching, whose period is on the order of 10-14 s (10 fs). The integration timestep should therefore be about 0.5-1 fs. If you use the SHAKE or RATTLE constraint algorithm (Constraints during dynamics simulations), a longer timestep is possible.

    If you are studying simple model liquids or solids and are not interested in internal modes, much longer timesteps may be used, e.g., up to 20 fs. A timestep of about 5 fs should be adequate for ionic material models.

    Appropriate for the integrator

    The timestep must also be appropriate to the integrator. For the ABM4 method, the timestep should be about half that needed for the Verlet algorithm. The Runge-Kutta-4 method seems to require a much smaller timestep than the other methods.

    http://www.chem.cmu.edu/courses/09-560/docs/msi/ffbsim/5_Dynamics.html#348259

    Choosing the time step
    All the dynamics objects need a time step. Choosing it too small will waste computer time, choosing it too large will make the dynamics unstable, typically the energy increases dramatically (the system “blows up”). If the time step is only a little to large, the lack of energy conservation is most obvious in Velocity Verlet dynamics, where energy should otherwise be conserved.

    Experience has shown that 5 femtoseconds is a good choice for most metallic systems. Systems with light atoms (e.g. hydrogen) and/or with strong bonds (carbon) will need a smaller time step.

    All the dynamics objects documented here are sufficiently related to have the same optimal time step.
    https://wiki.fysik.dtu.dk/ase/ase/md.html

  • 相关阅读:
    MS SQL数据类型
    ASP中调用存储过程、语法、写法-sql server数据库
    系統存儲過程System stored procedures
    ASP调用带参数存储过程的几种方式
    解密SQLSERVER2000存储过程,函数,视图,触发器
    SQL Script 執行的時間
    各種正則表达式示例
    12款操作系统安装全程图解(菜鸟必备宝典)
    内存 : DDR2与DDR
    數據庫分頁
  • 原文地址:https://www.cnblogs.com/Simulation-Campus/p/8663427.html
Copyright © 2020-2023  润新知