new_parralelism

lammps.in
View Options

			# Set the units (real = A, kcal/mol, fs, g/mol, atm)
			units real

			#parameters for aluminum (E=68GPa nu=0.4) at 0K first neighbor interaction
			#define the LJ parameters
			variable epsilon equal 4.243051074
			variable sigma equal 2.54000083658
			variable rcut equal 3.89461293159
			variable latsize equal "4.0320*sqrt(2)/2"
			variable mass equal 26.98153860


			# Set style and dimension
			dimension 2
			atom_style atomic

			# Turn newton off for contact algorithm
			newton on

			# Specify the boundaries (p = periodic , s = shrink-wrapped)
			boundary s p p

			lattice hex ${latsize} origin 0.25 0.25 0
			region box block -102 102 -25 25 -2 2
			create_box 1 box
			create_atoms 1 region box

			# Assign masses
			mass * ${mass}



			# Define pair styles
			pair_style lj/cut ${rcut}
			pair_coeff * * ${epsilon} ${sigma}

			# Some computation settings for the neighbor lists
			neigh_modify delay 1

			# Do the actual integration
			fix 1 all nve