fire_minimization
Differences
This shows you the differences between two versions of the page.
| 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 technique, which 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)
Page Tools
Except where otherwise noted, content on this wiki is licensed under the following license: CC Attribution-Noncommercial-Share Alike 4.0 International
