LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Mailing List Archives
Re: [lammps-users] Temperature suddendly drop down to zero
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [lammps-users] Temperature suddendly drop down to zero


From: Axel Kohlmeyer <akohlmey@...24...>
Date: Wed, 13 Sep 2017 12:41:13 -0400



On Wed, Sep 13, 2017 at 2:18 AM, kia krn <kiakrn@...24...> wrote:
Hello all
 
I am doing simulation for strain effect on Radiation damage in Aluminum and i used NVT ensemble. the simulations run well for 1KeV an 5KeV Primary knock-on atom but for 10 KeV Primary Knock-on atom the temperature drop down to zero and atoms freeze (in dump file i see zero velocities). i tried to reduce timestep from 0.0002 to 0.0001 but problem remains same.I used a box of size (70 * 70 *70).

​have you visualized your trajectory and noted what is happening around step 170 / 180?
you possibly should do thermo output more frequently and also include the potential energy.

axel.​


 

fix           relax all nvt temp 300.0 300.0 1.0

fix         2 all deform 1 z erate 0.0005 units box remap x

*******Log OUTPUT*******

 Step Atoms Temp f_ThermoAverad[1] f_ThermoAverad[2] Press TotEng Volume Lx Ly Lz 
   130   665500    354.03748    354.05728    -7251.937   -7251.8198   -3.3235097     11359930    224.05334    224.05334    226.29388 
     140   665500    353.77219    353.84129   -7245.9865   -7243.0179   -3.3235096     11359930    224.05334    224.05334    226.29388 
     150   665500     353.8682    353.79938   -7243.9829   -7247.0272   -3.3235097     11359930    224.05334    224.05334    226.29388 
     160   665500    353.88971    353.95054   -7251.0697   -7249.6294   -3.3235099     11359930    224.05334    224.05334    226.29388 
     170   665500    344.44038    349.29312   -6815.9655   -6286.5693   -3.3235184     11359930    224.05334    224.05334    226.29388 
     180   665500            0            0   -10200.298   -10231.186   -3.3696401     11359930    224.05334    224.05334    226.29388 
     190   665500            0            0   -10370.032   -10404.169   -3.3702055     11359930    224.05334    224.05334    226.29388 
     200   665500            0            0   -10538.619    -10571.88   -3.3707539     11359930    224.05334    224.05334    226.29388 
     210   665500            0            0   -10703.282   -10735.828   -3.3712886     11359930    224.05334    224.05334    226.29388 
     220   665500            0            0    -10864.52    -10896.41   -3.3718109     11359930    224.05334    224.05334    226.29388 
     230   665500            0            0   -11022.572   -11053.846   -3.3723218     11359930    224.05334    224.05334    226.29388 
     240   665500            0            0   -11177.606   -11208.289   -3.3728218     11359930    224.05334    224.05334    226.29388 
     250   665500            0            0    -11329.72   -11359.827   -3.3733114     11359930    224.05334    224.05334    226.29388 
     260   665500            0            0   -11478.997   -11508.544    -3.373791     11359930    224.05334    224.05334    226.29388 
     270   665500            0            0   -11625.501   -11654.505   -3.3742608     11359930    224.05334    224.05334    226.29388 
     280   665500            0            0   -11769.309   -11797.774   -3.3747212     11359930    224.05334    224.05334    226.29388 
     290   665500            0            0   -11910.468   -11938.413   -3.3751724     11359930    224.05334    224.05334    226.29388 
     300   665500            0            0   -12049.042   -12076.474   -3.3756147     11359930    224.05334    224.05334    226.29388 
     310   665500            0            0   -12185.071   -12212.002   -3.3760483     11359930    224.05334    224.05334    226.29388 
     320   665500            0            0   -12318.617   -12345.056   -3.3764735     11359930    224.05334    224.05334    226.29388 
     330   665500            0            0   -12449.719   -12475.676   -3.3768905     11359930    224.05334    224.05334    226.29388 
     340   665500            0            0   -12578.436   -12603.921   -3.3772994     11359930    224.05334    224.05334    226.29388 
     350   665500            0            0   -12704.806   -12729.825   -3.3777005     11359930    224.05334    224.05334    226.29388 
     360   665500            0            0   -12828.884   -12853.449    -3.378094     11359930    224.05334    224.05334    226.29388 
     370   665500            0            0   -12950.702    -12974.82     -3.37848     11359930    224.05334    224.05334    226.29388 
     380   665500            0            0   -13070.324   -13094.007   -3.3788588     11359930    224.05334    224.05334    226.29388 
     390   665500            0            0   -13187.776   -13211.033   -3.3792304     11359930    224.05334    224.05334    226.29388 
     400   665500            0            0   -13303.109   -13325.947   -3.3795952     11359930    224.05334    224.

------------------------------------------------------------------------------
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot
_______________________________________________
lammps-users mailing list
lammps-users@...396...sourceforge.net
https://lists.sourceforge.net/lists/listinfo/lammps-users




--
Dr. Axel Kohlmeyer  akohlmey@...24...  http://goo.gl/1wk0
College of Science & Technology, Temple University, Philadelphia PA, USA
International Centre for Theoretical Physics, Trieste. Italy.