#2D box simulation # Variable definition variable cutf equal 0.0109 variable neigh equal 0.1*${cutf} #Initialization units si atom_style atomic dimension 2 boundary p p p #Atom definition region box block 0 0.2675 0 0.2675 -1.5e-4 1.5e-4 create_box 3 box create_atoms 1 random 245000 3432435 NULL # Force field pair_style yukawa 1.1432e3 ${cutf} pair_coeff * * 3.2847e-22 #Settings neighbor ${neigh} bin mass * 6.99e-13 compute new3d all temp # equilibrate velocity all create 5.6275e2 5812775 temp new3d units box fix eq1 all nvt temp 5.6275e2 5.6275e2 0.02 thermo 10 timestep 0.0002 thermo_modify temp new3d # Radial distribution function compute myRDF all rdf 1000 fix 7 all ave/time 1 100 100 c_myRDF file tmp.rdf mode vector run 100