Hello,

I’m new to Lammps, and try to understand everything in the Langevin example script in thermal conductivity.In the thermal conductivity calculation. The example mentioned that the atoms number of the system is 8000, I think the atoms number is defined by some parameters (fcc$rho, mass, box size), unfortunately, I can’t figure it out. Why is the atoms number 8000 in the system? Can anyone tell me about it? Any help would be appreciated!

Here is the part of Langevin example input:

# sample LAMMPS input script for thermal conductivity of liquid LJ

# thermostatting 2 regions via fix Langevin

# settings

variable x equal 10

variable y equal 10

variable z equal 20

variable rho equal 0.6

variable Â Â Â Â t equal 1.35

variable rc equal 2.5

variable Â Â Â Â tlo equal 1.0

variable Â Â Â Â thi equal 1.70

#variable rho equal 0.85

#variable Â Â Â Â t equal 0.7

#variable rc equal 3.0

#variable Â Â Â Â tlo equal 0.3

#variable Â Â Â Â thi equal 1.0

# setup problem

units lj

atom_style atomic

lattice fcc ${rho}

region box block 0 $x 0 $y 0 $z

create_box 1 box

create_atoms 1 box

mass 1 1.0

velocity all create $t 87287

pair_style lj/cut ${rc}

pair_coeff 1 1 1.0 1.0

neighbor 0.3 bin

neigh_modify delay 0 every 1