the example of elastic constant calculation

Dear All,
I am using LAMMPS (14 Feb 2014) and have some questions :
1.In the example of elastic constant calculation,I couldn’t find ‘atom_style atomic’ command in init.mod.However,such
a command is usually essential in the input file I used before.Why doesn’t init.mod have one?

2.This example can only calculate elastic constant at 0K,right?

  1. To get the desired configuration,some models need to be treated via some methods such as quenching and so on.Because this example can only calculate elastic constant at 0K,before I use ‘read_data’ to load the configuration above,should I need to cool down it to 0K? Cooling down to other temperatures above 0K,like 300K,seems inappropriate.

Dear All,
I am using LAMMPS (14 Feb 2014) and have some questions :
1.In the example of elastic constant calculation,I couldn't find 'atom_style
atomic' command in init.mod.However,such
   a command is usually essential in the input file I used before.Why
doesn't init.mod have one?

because atomic is the default setting for atom_style

2.This example can only calculate elastic constant at 0K,right?

the scripts use minimizations, which by definition do not consider any
dynamics, i.e. determine the 0K structures unless they get trapped in
some other local minimum. as already mentioned by steve, obtaining
elastic constants at finite temperature, is a much trickier procedure.

3. To get the desired configuration,some models need to be treated via some
methods such as quenching and so on.Because this example can only calculate
elastic constant at 0K,before I use 'read_data' to load the configuration
above,should I need to cool down it to 0K? Cooling down to other
temperatures above 0K,like 300K,seems inappropriate.

the script assumes you start from a suitable crystal structure. simply
cooling down a system is not necessarily generating that (except for
rather simple materials). better to look up proper crystallographic
data.

axel.