Hello all,
I am probably misunderstanding something.
I used fix langevin and its working fine and the temperature is fixed as the command states in the manual.
But I tried to use fix Brownian and the temperature (from thermo) is very large, a factor of 400!
Is there a specific way to calculate or fix the temperature when using fix Brownian?
It seems to be connected to timestep, can someone please clarify it or let me know which topic I should look into?
Best
code used:
Initialization
units lj
dimension 2
atom_style atomic
pair_style none #lj/cut 2.5
boundary p p p
System definition
region myreg block -100 100 -100 100 -0.5 0.5
create_box 2 myreg
create_atoms 1 random 2000 341341 myreg
Simulation settings
mass * 1
neigh_modify every 1 delay 1 check yes
dump mydmp all atom 1000 dump.lammpstrj
Run
thermo 10
minimize 1.0e-4 1.0e-6 1000 10000
fix mynve all nve
#fix mylgv all langevin 1.0 1.0 0.1 1530917
fix 3 all brownian 1 12908410 gamma_t 1
fix myefn all enforce2d
timestep 0.00005
thermo 10000
run 100000