units              	lj
atom_style         	full
neighbor           	0.36 bin
neigh_modify       	delay 2

pair_style         	lj/cut 2.5

bond_style         	fene
special_bonds    	lj  0.0 1.0 1.0
read_data 			data_file.data

pair_coeff         	* * 1.0 1.0 2.5
bond_coeff        	* 30.0 1.5 1.0 1.0
pair_modify        	shift yes

timestep        	0.010
reset_timestep 		0

fix 				1 all npt temp 1.0 1.0 2.0 iso 0.0 0.0 10.0
thermo     			1000
thermo_style 		custom step temp pe etotal epair pxx pyy pzz lx ly lz
run 				100000

unfix 				1

# Store final simulation box length for strain calculation
variable 			tmp equal "lx"
variable 			L0 equal ${tmp}
variable 			srate_load equal 0.00001
print 				"Initial Length, L0: ${L0}"

# Perform Deformation
reset_timestep 		0

fix 	  			1 all nvt temp 1.0 1.0 2.0
fix 				2 all deform 1 x erate ${srate_load} y volume z volume
compute 			p all pressure thermo_temp
fix 				3 all ave/time 4 25 100 c_p[1] c_p[2] c_p[3] c_p[4] c_p[5] c_p[6] file press.1

run    				10000000