[gmx-users] How to estimate the reduced unit temperature?

[DeChang Li]

Dear gmx-users,

   I have simulated a Lennard-Jones system that used reduced units. In the
system it only contains one atom type B, the non-bonded parameter I wrote
below:

[ nonbond_params ]
  ; i    j func          c6           c12             sigma
epsilon
    B    B    1 0.400000000   0.400000000   ; --->     1.0000000
0.1000000




; Berendsen temperature coupling is on in one groups
Tcoupl              =  berendsen
tau_t               =   0.5
tc-grps             =  B_atoms
ref_t               =  10.0

   In the simulation I used Berendsen temperature coupling, and set the
ref_t = 10. Here I got a reduced temperature T'=0.083, is it right? And if
the real interaction strength of atom B is epsilon = 5 KJ/mol, how I can
estimate the actual temperature of my simulation system?

   Any suggestions will be help, thank you very much!



Dechang Li
Li.DC06 at gmail.com
2006-11-20
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20061120/498162dc/attachment.html>