processors 2 2 2

units lj
atom_style bond

special_bonds lj/coul 1 1 1

read_data polymer_100c100m.dat

neighbor 0.4 bin
neigh_modify every 1 delay 1
comm_modify vel yes

bond_style harmonic
bond_coeff * 1000 0.97

dump 	1 all custom 600 *_Monomer_Coordinates id mol type x y z ix iy iz
compute b1 all property/local btype batom1 batom2
compute b2 all bond/local dist engpot force
dump    2 all local 600 *_Bond_Analysis index c_b1[*] c_b2[*]

timestep 0.01
thermo 1000

pair_style dpd 1.0 1.0 122347
pair_coeff * * 25 4.5 1.0

velocity all create 1.0 17786140

fix 1 all nve
run 50000

pair_coeff * * 50.0 4.5 1.0
velocity all create 1.0 15086120
run 50
pair_coeff * * 100.0 4.5 1.0
velocity all create 1.0 15786120 
run 50
pair_coeff * * 150.0 4.5 1.0
velocity all create 1.0 15486120
run 50
pair_coeff * * 200.0 4.5 1.0
velocity all create 1.0 17986120
run 100
pair_coeff * * 250.0 4.5 1.0
velocity all create 1.0 15006120
run 100
pair_coeff * * 500.0 4.5 1.0
velocity all create 1.0 15087720
run 100
pair_coeff * * 1000.0 4.5 1.0
velocity all create 1.0 15086189
run 100

pair_style hybrid/overlay lj/cut 1.122462 dpd/tstat 1.0 1.0 1.122462 122347
pair_modify shift yes
pair_coeff * * lj/cut 1.0 1.0 1.122462
pair_coeff * * dpd/tstat 4.5 1.122462
velocity all create 1.0 1508612013
run 100
velocity all create 1.0 15086121
run 100
velocity all create 1.0 15086111
run 100
velocity all create 1.0 15086125
run 6000000
write_data polymer_100x100.data
write_restart my_restart_lj