[讨论交流] 求助关于minimize 命令的问题 相关搜索: minimize, 命令, 求助

sunysb

I just found one strange problem about the minimize command.

The minimize command considers energy due to fixes which can act as
constraints or apply force to atoms. The Lennard-Jones potential is
used between this two balls.

Now I created one problem. There are two balls, and one ball is fixed
the other one is applied an external force. Then I use the minimize
command to find the lowest potential energy position. I think at the
lowest potential point this system must be at balance state. So the
external force must be equal to the attraction( or repulsion) force to
make the potential energy minimum. But the lammps does not follow this
rule. Please run the below code.

I used the command fix_modify but it does not give the right result.

# minimization test

echo                screen
dimension           3
units               real
atom_style          atomic
boundary            f f f

pair_style          lj/cut   100
neighbor            100  bin
neigh_modify        every   1

#read in coordinate information
read_data           data.try
group upper         id    2
group lower         id    1


velocity            lower set 0 0 0 units box
velocity            upper set 0 0 0 units box
fix                 1   lower planeforce 0 1 0

#displace_atoms      upper move 0 1.107 0 units box


fix                 3 upper addforce 0 2.3 0
fix_modify          3 energy yes
minimize            1.0e-19 1.0e-19 100000000 100000000


timestep            1
thermo              1
dump                2   upper   custom   1   dump.try2  y fy vy
dump                1   all xyz 1 dump.try1
run                 100




Below is the file data.try

2 atoms  
1 atom types

-300 300 xlo xhi
-300 300 ylo yhi
-300 300 zlo zhi


Masses
  
1 1

Pair Coeffs

1  1.0  1.0   100

Atoms

1   1    0.0    0.0      0
2   1    0.0    1.1225      0

Thank you very much!

redream

找steve吧，呵呵，这个是算法问题。

minimize只能用于周期性边界条件，而且fix rigid的原子都会给你打乱的。

alan8410

请问楼主：Pair Coeffs  中

1  1.0  1.0   100

各代表什么含义？

stonyhit

minimize可以用于非周期边界条件，我试过纳米压痕（压头方向非周期边界），没问题的。