Hi,

I am running a minimization using the hftn algorithm. I ran two simulations: one with a force norm convergence criterion of 1e-3 and another with a force norm convergence criterion of 1. All other parameters are the same.

However the optimization seemed to take different paths after a certain iteration step. I find this strange as the convergence criterion shouldn’t influence the path of the optimization. I was just expecting the second simulation to take the same path as the first but just stop earlier.

Please let me know if anyone knows why this may be happening.

Thanks,

Sridhar

I don’t know. Here is a simple test script that

does not do that. The min output is the same

up to the point the earliest convergence takes place.

Steve

# 2d Lennard-Jones melt and subsequent energy minimization

units lj

dimension 2

atom_style atomic

lattice sq2 0.8442

region box block 0 20 0 20 -0.1 0.1

create_box 1 box

create_atoms 1 box

mass 1 1.0

velocity all create 5.0 87287 loop geom

pair_style lj/cut 2.5

pair_coeff 1 1 1.0 1.0 2.5

pair_modify shift yes

neighbor 0.3 bin

neigh_modify delay 0 every 1 check yes

fix 1 all nve

fix 2 all enforce2d

dump 1 all atom 100 dump.min

thermo 100

run 1000

neigh_modify delay 0 every 1 check yes

dump_modify 1 every 25

thermo 10

min_style hftn

minimize 0.0 0.01 1000 10000

From the doc page, “At each (HFTN) iteration a quadratic model of the energy potential is solved by a conjugate gradient (CG) inner iteration.”

HFTN performs a convergence test for the CG iterations where a properly chose tolerance is crucial so as to ensure that CG terminates in a reasonable way. Also there is the correction to the “residual” that is based on the tolerance. Therefore depending on the ruggedness of the potential energy surface, difference tolerance may give different results.

In a nutshell, it is possible that convergence criterion (tolerance) *does* influence the path of the optimization.

Ray

Ray, thanks for your note. That explains my observations.

Sridhar