User Tools

Site Tools


fire_minimization

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
fire_minimization [2014/07/08 12:31]
prokop
fire_minimization [2014/07/08 14:20] (current)
Line 2: Line 2:
 ====== FIRE minimization ====== ====== FIRE minimization ======
  
-FIRE ( ''​iquench = -6''​ ) is molecular dynamics based minimization technique which use only gradient of potential ( force ), but it often provides performance comparable or even superior ​tor quasi-second-order Quasi-Newton algorithms such as BFGS. It differs from quenched molecular dynamics in the way how velocity is updated (see. pseudocode or referecnes ​for more details). An other difference from quenched-MD (iquench=-1) is that mass of all atoms is set to 4.0 atomic units ( can be modified by FIRE_mass parameter) so Gold and Hydrogen atoms have the same inertia. In comparison to BFGS, FIRE algorithm is more robust ( and often also faster ) if the potential is far from being quadratic form, and in case of long narrow bend valley. For this reason it is well suited in case of soft degrees of freedom (such as dihedral angles, intermolecular electrostatic forces, hydrogen bond etc.)+FIRE ( ''​iquench = -6''​ ) is molecular dynamics based minimization techniquewhich use only gradient of the energy ​potential ( force ), but it often provides performance comparable or even superior ​to quasi-second-order Quasi-Newton algorithms such as BFGS. It differs from quenched molecular dynamics in the way how velocity is updated (see. pseudocode or references ​for more details). An other difference from standard ​quenched-MD (iquench=-1) is that mass of all atoms is set to 4.0 atomic units ( can be modified by FIRE_mass parameter) so Gold and Hydrogen atoms have the same inertia. In comparison to BFGS, FIRE algorithm is more robust ( and often also faster ) if the potential is far from being quadratic form, and in case of long narrow bend valley. For this reason it is well suited in case of soft degrees of freedom (such as dihedral angles, intermolecular electrostatic forces, hydrogen bond etc.)
  
 === Settings === === Settings ===
  
-Default ​paramenters ​of FIRE algorithm are good for most of the cases and it is not necessary to set them manually. The only parameter which should be set in ''​fireball.in''​ is time step ''​dt = 0.5 - 1.0''​ femtosecond.+Default ​parameters ​of FIRE algorithm are good for most of the cases and it is not necessary to set them manually. The only parameter which should be set in ''​fireball.in''​ is time step ''​dt = 0.5 - 1.0''​ femtosecond.
  
-In case you want to try your luck and play with the parameters of algorithm it is optionally possible to provide file 'FIRE.opt' of the folowing format (with default ​paramenters ​as examples):+In case you want to try your luck and play with the parameters of algorithm it is optionally possible to provide file 'FIRE.optional' of the folowing format (with default ​parameters ​as examples):
  
  1.1  ! FIRE_finc ​  ... increment time step if dot(f,v) is positive  1.1  ! FIRE_finc ​  ... increment time step if dot(f,v) is positive
Line 43: Line 43:
  
 === Results === === Results ===
 +
 +{{:​fire_hcooh.png| }}
 +
 +{{:​fire_ch-kink.png|}}
 +
 +{{:​fire_h2o2_fixed.png|}}
fire_minimization.1404815508.txt.gz ยท Last modified: 2014/07/08 12:31 (external edit)