# LAMMPS Input File

# Licensed under CC BY 4.0

# A Set of Tutorials for the LAMMPS Simulation Package (LiveCoMS, 2025)

# By Simon Gravelle, Cecilia M. S. Alvares, Jacob R. Gissinger, and Axel Kohlmeyer

# Please cite doi.org/10.33011/livecoms.6.1.3037

# Find more on GitHub: https://github.com/lammpstutorials

variable sigma equal 3.405

variable epsilon equal 0.238

variable U0 equal 1.5*${epsilon}

variable dlt equal 1.0

variable x0 equal 10.0

units real

atom_style atomic

pair_style lj/cut $(2^(1/6)*v_sigma)

pair_modify shift yes

boundary p p p

region myreg block -50 50 -15 15 -50 50

create_box 1 myreg

create_atoms 1 random 200 34134 myreg overlap 3 maxtry 50

mass * 39.95

pair_coeff * * ${epsilon} ${sigma}

variable U atom ${U0}*atan((x+${x0})/${dlt})-${U0}*atan((x-${x0})/${dlt})

variable F atom ${U0}/((x-${x0})^2/${dlt}^2+1)/${dlt}-${U0}/((x+${x0})^2/${dlt}^2+1)/${dlt}

fix myadf all addforce v_F 0.0 0.0 energy v_U

fix mynve all nve

fix mylgv all langevin 119.8 119.8 500 30917

region mymes block -${x0} ${x0} INF INF INF INF

variable n_center equal count(all,mymes)

thermo_style custom step temp etotal v_n_center

thermo 2000

dump viz all image 5000 myimage-*.ppm type type &

shiny 0.1 box yes 0.01 view 180 90 zoom 6 size 1600 500 fsaa yes

dump_modify viz backcolor white acolor 1 cyan adiam 1 3 boxcolor black

timestep 2.0

run 50000

reset_timestep 0

thermo 200000

compute cc1 all chunk/atom bin/1d x 0.0 2.0

fix myac all ave/chunk 100 20000 2000000 &

cc1 density/number file free-sampling.dat

run 2000000