#C-C by REBO
# ------------------------ INITIALIZATION ----------------------------
units        metal
boundary     p p p  
atom_style   atomic
#neighbor        0.3 bin
variable    T equal 600

# ----------------------- ATOM DEFINITION ----------------------------

read_data       t2.lmp
timestep 0.001

# ------------------------ FORCE FIELDS ------------------------------
pair_style hybrid/overlay sw sw rebo  &  
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  & 
kolmogorov/crespi/z 24.0   kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0  kolmogorov/crespi/z 24.0   

# SW for layer1-WS2
pair_coeff      * *  sw 1 h-ws2.sw S1 W1 S2 S3 W2 S4 S5 W3 S6 S7 W4 S8 NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
# SW for layer2-MoS2
pair_coeff      * *  sw 2 h-mos2.sw NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S1 Mo1 S2 S3 Mo2 S4 S5 Mo3 S6 S7 Mo4 S8 NULL NULL

# SW for layer3,4-Graphene
pair_coeff  * * rebo  CH.rebo NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL C C

# Set arbitrary interactions - 0, before overlaying with KC
#pair_coeff * * lj/cut 0.0 3.4

#------------------------------------------
# S layer1-WS2: 1, 3, 4, 6, 7, 9, 10, 12
# S layer2-MoS2: 13, 15, 16, 18, 19, 21, 22, 24
# W layer1-WS2: 2, 5, 8, 11
# Mo layer2-MoS2: 14, 17, 20, 23
# C layer3-C: 25
# C layer4-C: 26
#------------------------------------------
# S-S KC interactions.
#------------------------------------------
# KC between 1 and 13,15,16,18,19,21,22,24
pair_coeff 1 13 kolmogorov/crespi/z 1 ILP.KC S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL   S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 1 15 kolmogorov/crespi/z 2 ILP.KC S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 1 16 kolmogorov/crespi/z 3 ILP.KC S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 1 18 kolmogorov/crespi/z 4  ILP.KC S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL 
pair_coeff 1 19 kolmogorov/crespi/z 5 ILP.KC S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL
pair_coeff 1 21 kolmogorov/crespi/z 6  ILP.KC S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL 
pair_coeff 1 22 kolmogorov/crespi/z 7 ILP.KC S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL
pair_coeff 1 24 kolmogorov/crespi/z 8  ILP.KC S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL

# KC between 3 and 13,15,16,18,19,21,22,24
pair_coeff 3 13 kolmogorov/crespi/z 9 ILP.KC  NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL   S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 3 15 kolmogorov/crespi/z 10 ILP.KC  NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 3 16 kolmogorov/crespi/z 11 ILP.KC  NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 3 18 kolmogorov/crespi/z 12  ILP.KC NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL 
pair_coeff 3 19 kolmogorov/crespi/z 13 ILP.KC  NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL
pair_coeff 3 21 kolmogorov/crespi/z 14  ILP.KC  NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL 
pair_coeff 3 22 kolmogorov/crespi/z 15 ILP.KC  NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL
pair_coeff 3 24 kolmogorov/crespi/z 16  ILP.KC  NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL

# KC between 4 and 13,15,16,18,19,21,22,24
pair_coeff 4 13 kolmogorov/crespi/z 17 ILP.KC  NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL   S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 4 15 kolmogorov/crespi/z 18 ILP.KC  NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 4 16 kolmogorov/crespi/z 19 ILP.KC  NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 4 18 kolmogorov/crespi/z 20  ILP.KC NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL 
pair_coeff 4 19 kolmogorov/crespi/z 21 ILP.KC  NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL
pair_coeff 4 21 kolmogorov/crespi/z 22  ILP.KC  NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL 
pair_coeff 4 22 kolmogorov/crespi/z 23 ILP.KC  NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL
pair_coeff 4 24 kolmogorov/crespi/z 24  ILP.KC  NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL

# KC between 6 and 13,15,16,18,19,21,22,24
pair_coeff 6 13 kolmogorov/crespi/z 25 ILP.KC  NULL NULL NULL  NULL NULL S NULL NULL NULL NULL NULL NULL   S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 6 15 kolmogorov/crespi/z 26 ILP.KC  NULL NULL NULL  NULL NULL S NULL NULL NULL NULL NULL NULL   NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 6 16 kolmogorov/crespi/z 27 ILP.KC  NULL NULL NULL  NULL NULL S NULL NULL NULL NULL NULL NULL   NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 6 18 kolmogorov/crespi/z 28  ILP.KC NULL NULL NULL  NULL NULL S NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL 
pair_coeff 6 19 kolmogorov/crespi/z 29 ILP.KC  NULL NULL NULL  NULL NULL S NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL
pair_coeff 6 21 kolmogorov/crespi/z 30  ILP.KC  NULL NULL NULL  NULL NULL S NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL 
pair_coeff 6 22 kolmogorov/crespi/z 31 ILP.KC  NULL NULL NULL  NULL NULL S NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL
pair_coeff 6 24 kolmogorov/crespi/z 32  ILP.KC  NULL NULL NULL  NULL NULL S NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL

# KC between 7 and 13,15,16,18,19,21,22,24
pair_coeff 7 13 kolmogorov/crespi/z 33 ILP.KC  NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL   S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 7 15 kolmogorov/crespi/z 34 ILP.KC  NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL   NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 7 16 kolmogorov/crespi/z 35 ILP.KC  NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL   NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 7 18 kolmogorov/crespi/z 36  ILP.KC NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL 
pair_coeff 7 19 kolmogorov/crespi/z 37 ILP.KC  NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL
pair_coeff 7 21 kolmogorov/crespi/z 38  ILP.KC  NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL 
pair_coeff 7 22 kolmogorov/crespi/z 39 ILP.KC  NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL
pair_coeff 7 24 kolmogorov/crespi/z 40  ILP.KC  NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL

# KC between 9 and 13,15,16,18,19,21,22,24
pair_coeff 9 13 kolmogorov/crespi/z 41 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL   S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 9 15 kolmogorov/crespi/z 42 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL   NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 9 16 kolmogorov/crespi/z 43 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL   NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 9 18 kolmogorov/crespi/z 44  ILP.KC NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL   NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL 
pair_coeff 9 19 kolmogorov/crespi/z 45 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL   NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL
pair_coeff 9 21 kolmogorov/crespi/z 46  ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL 
pair_coeff 9 22 kolmogorov/crespi/z 47 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL
pair_coeff 9 24 kolmogorov/crespi/z 48  ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL

# KC between 10 and 13,15,16,18,19,21,22,24
pair_coeff 10 13 kolmogorov/crespi/z 49 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL   S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 10 15 kolmogorov/crespi/z 50 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL   NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 10 16 kolmogorov/crespi/z 51 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL   NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 10 18 kolmogorov/crespi/z 52  ILP.KC NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL   NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL 
pair_coeff 10 19 kolmogorov/crespi/z 53 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL   NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL
pair_coeff 10 21 kolmogorov/crespi/z 54  ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL 
pair_coeff 10 22 kolmogorov/crespi/z 55 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL
pair_coeff 10 24 kolmogorov/crespi/z 56  ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL

# KC between 12 and 13,15,16,18,19,21,22,24
pair_coeff 12 13 kolmogorov/crespi/z 57 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S   S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 12 15 kolmogorov/crespi/z 58 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S  NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 12 16 kolmogorov/crespi/z 59 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S  NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 12 18 kolmogorov/crespi/z 60  ILP.KC NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S  NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL 
pair_coeff 12 19 kolmogorov/crespi/z 61 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S  NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL
pair_coeff 12 21 kolmogorov/crespi/z 62  ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S  NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL 
pair_coeff 12 22 kolmogorov/crespi/z 63 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S  NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL
pair_coeff 12 24 kolmogorov/crespi/z 64  ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL
#------------------------------------------
# S-Mo, S-W KC interactions.
#------------------------------------------
pair_coeff 2 13 kolmogorov/crespi/z 65 ILP.KC  NULL W NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL   S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 2 15 kolmogorov/crespi/z 66 ILP.KC  NULL W NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 2 16 kolmogorov/crespi/z 67 ILP.KC  NULL W NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 2 18 kolmogorov/crespi/z 68  ILP.KC NULL W NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL 
pair_coeff 2 19 kolmogorov/crespi/z 69 ILP.KC  NULL W NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL
pair_coeff 2 21 kolmogorov/crespi/z 70  ILP.KC  NULL W NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL 
pair_coeff 2 22 kolmogorov/crespi/z 71 ILP.KC  NULL W NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL
pair_coeff 2 24 kolmogorov/crespi/z 72  ILP.KC  NULL W NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL

pair_coeff 5 13 kolmogorov/crespi/z 73 ILP.KC  NULL NULL NULL NULL W NULL NULL NULL NULL NULL NULL NULL   S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 5 15 kolmogorov/crespi/z 74 ILP.KC  NULL NULL NULL NULL W NULL NULL NULL NULL NULL NULL NULL   NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 5 16 kolmogorov/crespi/z 75 ILP.KC  NULL NULL NULL NULL W NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 5 18 kolmogorov/crespi/z 76  ILP.KC NULL NULL NULL NULL W NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL 
pair_coeff 5 19 kolmogorov/crespi/z 77 ILP.KC  NULL NULL NULL NULL W NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL
pair_coeff 5 21 kolmogorov/crespi/z 78  ILP.KC  NULL NULL NULL NULL W NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL 
pair_coeff 5 22 kolmogorov/crespi/z 79 ILP.KC  NULL NULL NULL NULL W NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL
pair_coeff 5 24 kolmogorov/crespi/z 80  ILP.KC  NULL NULL NULL NULL W NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL

pair_coeff 8 13 kolmogorov/crespi/z 81 ILP.KC  NULL NULL NULL NULL NULL NULL NULL W NULL NULL NULL NULL   S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 8 15 kolmogorov/crespi/z 82 ILP.KC  NULL NULL NULL NULL NULL NULL NULL W NULL NULL NULL NULL   NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 8 16 kolmogorov/crespi/z 83 ILP.KC  NULL NULL NULL NULL NULL NULL NULL W NULL NULL NULL NULL   NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 8 18 kolmogorov/crespi/z 84  ILP.KC NULL NULL NULL NULL NULL NULL NULL W NULL NULL NULL NULL   NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL 
pair_coeff 8 19 kolmogorov/crespi/z 85 ILP.KC  NULL NULL NULL NULL NULL NULL NULL W NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL
pair_coeff 8 21 kolmogorov/crespi/z 86  ILP.KC  NULL NULL NULL NULL NULL NULL NULL W NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL 
pair_coeff 8 22 kolmogorov/crespi/z 87 ILP.KC  NULL NULL NULL NULL NULL NULL NULL W NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL
pair_coeff 8 24 kolmogorov/crespi/z 88  ILP.KC  NULL NULL NULL NULL NULL NULL NULL W NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL

pair_coeff 11 13 kolmogorov/crespi/z 89 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL W NULL   S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 11 15 kolmogorov/crespi/z 90 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL W NULL   NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 11 16 kolmogorov/crespi/z 91 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL W NULL   NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 11 18 kolmogorov/crespi/z 92  ILP.KC NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL W NULL   NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL 
pair_coeff 11 19 kolmogorov/crespi/z 93 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL W NULL   NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL
pair_coeff 11 21 kolmogorov/crespi/z 94  ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL W NULL   NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL 
pair_coeff 11 22 kolmogorov/crespi/z 95 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL W NULL   NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL
pair_coeff 11 24 kolmogorov/crespi/z 96  ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL W NULL   NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL

pair_coeff 1 14 kolmogorov/crespi/z 97 ILP.KC  S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL   NULL Mo NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 3 14 kolmogorov/crespi/z 98 ILP.KC  NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL   NULL Mo NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 4 14 kolmogorov/crespi/z 99 ILP.KC  NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL   NULL Mo NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 6 14 kolmogorov/crespi/z 100  ILP.KC NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL   NULL Mo NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL 
pair_coeff 7 14 kolmogorov/crespi/z 101 ILP.KC  NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL   NULL Mo NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 9 14 kolmogorov/crespi/z 102  ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL   NULL Mo NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL 
pair_coeff 10 14 kolmogorov/crespi/z 103 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL   NULL Mo NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 12 14 kolmogorov/crespi/z 104  ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S   NULL Mo NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL

pair_coeff 1 17 kolmogorov/crespi/z 105 ILP.KC  S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL Mo NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 3 17 kolmogorov/crespi/z 106 ILP.KC  NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL Mo NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 4 17 kolmogorov/crespi/z 107 ILP.KC  NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL Mo NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 6 17 kolmogorov/crespi/z 108  ILP.KC NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL Mo NULL NULL NULL NULL NULL NULL NULL NULL NULL 
pair_coeff 7 17 kolmogorov/crespi/z 109 ILP.KC  NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL   NULL NULL NULL NULL Mo NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 9 17 kolmogorov/crespi/z 110  ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL   NULL NULL NULL NULL Mo NULL NULL NULL NULL NULL NULL NULL NULL NULL 
pair_coeff 10 17 kolmogorov/crespi/z 111 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL   NULL NULL NULL NULL Mo NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 12 17 kolmogorov/crespi/z 112  ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S   NULL NULL NULL NULL Mo NULL NULL NULL NULL NULL NULL NULL NULL NULL

pair_coeff 1 20 kolmogorov/crespi/z 113 ILP.KC  S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL Mo NULL NULL NULL NULL NULL NULL
pair_coeff 3 20 kolmogorov/crespi/z 114 ILP.KC  NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL Mo NULL NULL NULL NULL NULL NULL
pair_coeff 4 20 kolmogorov/crespi/z 115 ILP.KC  NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL Mo NULL NULL NULL NULL NULL NULL
pair_coeff 6 20 kolmogorov/crespi/z 116  ILP.KC NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL Mo NULL NULL NULL NULL NULL NULL 
pair_coeff 7 20 kolmogorov/crespi/z 117 ILP.KC  NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL Mo NULL NULL NULL NULL NULL NULL
pair_coeff 9 20 kolmogorov/crespi/z 118  ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL Mo NULL NULL NULL NULL NULL NULL 
pair_coeff 10 20 kolmogorov/crespi/z 119 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL   NULL NULL NULL NULL NULL NULL NULL Mo NULL NULL NULL NULL NULL NULL
pair_coeff 12 20 kolmogorov/crespi/z 120  ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S   NULL NULL NULL NULL NULL NULL NULL Mo NULL NULL NULL NULL NULL NULL

pair_coeff 1 23 kolmogorov/crespi/z 121 ILP.KC  S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL Mo NULL NULL NULL
pair_coeff 3 23 kolmogorov/crespi/z 122 ILP.KC  NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL Mo NULL NULL NULL
pair_coeff 4 23 kolmogorov/crespi/z 123 ILP.KC  NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL Mo NULL NULL NULL
pair_coeff 6 23 kolmogorov/crespi/z 124  ILP.KC NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL Mo NULL NULL NULL 
pair_coeff 7 23 kolmogorov/crespi/z 125 ILP.KC  NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL Mo NULL NULL NULL
pair_coeff 9 23 kolmogorov/crespi/z 126  ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL Mo NULL NULL NULL 
pair_coeff 10 23 kolmogorov/crespi/z 127 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL Mo NULL NULL NULL
pair_coeff 12 23 kolmogorov/crespi/z 128  ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S   NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL Mo NULL NULL NULL

#------------------------------------------
# W-Mo KC interactions.
#------------------------------------------
# W layer1: 2, 5, 8, 11
# Mo layer2: 14, 17, 20, 23

# KC between 2 and 14,17,20,23
pair_coeff 2 14 kolmogorov/crespi/z 129  ILP.KC NULL W NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL   NULL Mo NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 2 17 kolmogorov/crespi/z 130  ILP.KC NULL W NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL Mo NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 2 20 kolmogorov/crespi/z 131  ILP.KC NULL W NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL Mo NULL NULL NULL NULL NULL NULL
pair_coeff 2 23 kolmogorov/crespi/z 132  ILP.KC NULL W NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL Mo NULL NULL NULL
#
#
## KC between 5 and 14,17,20,23
pair_coeff 5 14 kolmogorov/crespi/z 133  ILP.KC NULL NULL NULL NULL W NULL NULL NULL NULL NULL NULL NULL   NULL Mo NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 5 17 kolmogorov/crespi/z 134  ILP.KC NULL NULL NULL NULL W NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL Mo NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 5 20 kolmogorov/crespi/z 135  ILP.KC NULL NULL NULL NULL W NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL Mo NULL NULL NULL NULL NULL NULL
pair_coeff 5 23 kolmogorov/crespi/z 136  ILP.KC NULL NULL NULL NULL W NULL NULL NULL NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL Mo NULL NULL NULL
#
## KC between 8 and 14,17,20,23
pair_coeff 8 14 kolmogorov/crespi/z 137  ILP.KC NULL NULL NULL NULL NULL NULL NULL W NULL NULL NULL NULL   NULL Mo NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 8 17 kolmogorov/crespi/z 138 ILP.KC NULL NULL NULL NULL NULL NULL NULL W NULL NULL NULL NULL   NULL NULL NULL NULL Mo NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 8 20 kolmogorov/crespi/z 139  ILP.KC NULL NULL NULL NULL NULL NULL NULL W NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL Mo NULL NULL NULL NULL NULL NULL
pair_coeff 8 23 kolmogorov/crespi/z 140  ILP.KC NULL NULL NULL NULL NULL NULL NULL W NULL NULL NULL NULL   NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL Mo NULL NULL NULL
#
## KC between 11 and 14,17,20,23
pair_coeff 11 14 kolmogorov/crespi/z 141  ILP.KC NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL W NULL   NULL Mo NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 11 17 kolmogorov/crespi/z 142  ILP.KC NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL W NULL   NULL NULL NULL NULL Mo NULL NULL NULL NULL NULL NULL NULL NULL NULL
pair_coeff 11 20 kolmogorov/crespi/z 143  ILP.KC NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL W NULL   NULL NULL NULL NULL NULL NULL NULL Mo NULL NULL NULL NULL NULL NULL
pair_coeff 11 23 kolmogorov/crespi/z 144  ILP.KC NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL W NULL   NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL Mo NULL NULL NULL

#
## KC betwen layer1(WS2) and C layers(3,4)
pair_coeff 1  25   kolmogorov/crespi/z 145 ILP.KC  S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL C NULL
pair_coeff 3  25   kolmogorov/crespi/z 146 ILP.KC  NULL NULL  S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL C NULL
pair_coeff 4  25   kolmogorov/crespi/z 147 ILP.KC  NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL C NULL
pair_coeff 6  25   kolmogorov/crespi/z 148 ILP.KC  NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL C NULL
pair_coeff 7  25   kolmogorov/crespi/z 149 ILP.KC  NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL C NULL
pair_coeff 9  25   kolmogorov/crespi/z 150 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL C NULL
pair_coeff 10 25   kolmogorov/crespi/z 151 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL C NULL
pair_coeff 12 25   kolmogorov/crespi/z 152 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL C NULL
pair_coeff 2  25   kolmogorov/crespi/z 153 ILP.KC  NULL W NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL C NULL
pair_coeff 5  25   kolmogorov/crespi/z 154 ILP.KC  NULL NULL NULL NULL W NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL C NULL
pair_coeff 8  25   kolmogorov/crespi/z 155 ILP.KC  NULL NULL NULL NULL NULL NULL NULL W NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL C NULL
pair_coeff 11 25   kolmogorov/crespi/z 156 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL W NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL C NULL
pair_coeff 1  26   kolmogorov/crespi/z 157 ILP.KC  S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL  NULL C
pair_coeff 3  26   kolmogorov/crespi/z 158 ILP.KC  NULL NULL  S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL C
pair_coeff 4  26   kolmogorov/crespi/z 159 ILP.KC  NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL  NULL C
pair_coeff 6  26   kolmogorov/crespi/z 160 ILP.KC  NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL  NULL C
pair_coeff 7  26   kolmogorov/crespi/z 161 ILP.KC  NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL  NULL C
pair_coeff 9  26   kolmogorov/crespi/z 162 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL  NULL C
pair_coeff 10 26   kolmogorov/crespi/z 163 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL  NULL C
pair_coeff 12 26   kolmogorov/crespi/z 164 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL  NULL C
pair_coeff 2  26   kolmogorov/crespi/z 165 ILP.KC  NULL W NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL C
pair_coeff 5  26   kolmogorov/crespi/z 166 ILP.KC  NULL NULL NULL NULL W NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL C
pair_coeff 8  26   kolmogorov/crespi/z 167 ILP.KC  NULL NULL NULL NULL NULL NULL NULL W NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL C
pair_coeff 11 26   kolmogorov/crespi/z 168 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL W NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL C
#
## KC betwen layer2(MoS2) and C layers(3,4)
pair_coeff 13  25   kolmogorov/crespi/z 169 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL C NULL
pair_coeff 15  25   kolmogorov/crespi/z 170 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL  S NULL NULL NULL NULL NULL NULL NULL NULL NULL C NULL
pair_coeff 16  25   kolmogorov/crespi/z 171 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL C NULL
pair_coeff 18  25   kolmogorov/crespi/z 172 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL C NULL
pair_coeff 19  25   kolmogorov/crespi/z 173 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL C NULL
pair_coeff 21  25   kolmogorov/crespi/z 174 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL C NULL
pair_coeff 22  25   kolmogorov/crespi/z 175 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL C NULL
pair_coeff 24  25   kolmogorov/crespi/z 176 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S C NULL
pair_coeff 14  25   kolmogorov/crespi/z 177 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL Mo NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL C NULL
pair_coeff 17  25   kolmogorov/crespi/z 178 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL Mo NULL NULL NULL NULL NULL NULL NULL C NULL
pair_coeff 20  25   kolmogorov/crespi/z 179 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL Mo NULL NULL NULL NULL C NULL
pair_coeff 23  25   kolmogorov/crespi/z 180 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL Mo NULL C NULL
pair_coeff 13  26   kolmogorov/crespi/z 181 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL C
pair_coeff 15  26   kolmogorov/crespi/z 182 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL  S NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL C
pair_coeff 16  26   kolmogorov/crespi/z 183 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL NULL NULL C
pair_coeff 18  26   kolmogorov/crespi/z 184 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL NULL C
pair_coeff 19  26   kolmogorov/crespi/z 185 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL NULL NULL C
pair_coeff 21  26   kolmogorov/crespi/z 186 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL NULL C
pair_coeff 22  26   kolmogorov/crespi/z 187 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL NULL NULL C
pair_coeff 24  26   kolmogorov/crespi/z 188 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL S NULL C
pair_coeff 14  26   kolmogorov/crespi/z 189 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL Mo NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL C
pair_coeff 17  26   kolmogorov/crespi/z 190 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL Mo NULL NULL NULL NULL NULL NULL NULL NULL C
pair_coeff 20  26   kolmogorov/crespi/z 191 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL Mo NULL NULL NULL NULL NULL C
pair_coeff 23  26   kolmogorov/crespi/z 192 ILP.KC  NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL NULL Mo NULL NULL C

#

group ws2 type 1:12
group mos2 type 13:24

# ------------------------- SETTINGS ---------------------------------
neighbor     2 bin
neigh_modify delay 0 one 10000

shell mkdir dump_gBx_600K
shell cd dump_gBx_600K

#----------- Energy Minimization---------------------------
thermo_style custom step temp etotal press pxx pyy pzz lx ly lz
thermo 10000

dump 1 all custom 100 dump.min.defo.* id type x y z 
min_style cg
minimize 1e-8 1e-8 5000 10000
undump 1

#---------------Energy Equilibriation--------------------------
reset_timestep 0
velocity        all create $T 87287 mom yes rot yes dist gaussian            

fix 1 all nve
fix 2 all temp/rescale 10 $T $T 1.0 0.5
dump 1 all custom 50000 dump.nve.defo.* id type x y z 

run 50000        
unfix 1
unfix 2
undump 1

#-----------------Pressure Equilibriation-----------------------
reset_timestep 0
fix 1 all npt temp $T $T 50.0 x 0.0 0.0 10 y 0.0 0.0 10 drag 0.0
dump 1 all custom 50000 dump.npt.defo.* id type x y z              
run 100000 
unfix 1
undump 1

#--------------- Equilibriation--------------------------
reset_timestep 0
fix 1 all nvt temp $T $T 50
dump 1 all custom 50000 dump.nvt.defo.* id type x y z 
run 50000      
unfix 1
undump 1

#--------------------Storing Initial length---------------------
variable tmp equal "lx"
variable L0 equal ${tmp}
print "Initial Length, L0: ${L0}"
#variable strain equal "v_srate/1e12"

shell cd ..
shell mkdir deform_gBx_600K
shell cd deform_gBx_600K

#### Computes Required
compute csym all centro/atom 12
compute 2 all stress/atom NULL virial
compute mytemp all temp
compute 11 all reduce sum c_2[1]
compute 12 all reduce sum c_2[2]
compute 13 all reduce sum c_2[3]
compute 14 all reduce sum c_2[4]
compute 15 all reduce sum c_2[5]
compute 16 all reduce sum c_2[6]

#------------Deform------------------------------

reset_timestep 0
fix 1 all deform 1 x erate 0.001 units box 
fix 2 all nvt temp $T $T 50
#change direction along stress direction
#fix 2 all npt temp $T $T 50 y 0.0 0.0 10
fix 5 ws2 spring couple mos2 100 NULL NULL -6.137 0
fix 3 all ave/atom 1 100 100 c_2[*]

variable strain equal "(lx - v_L0)/v_L0"
variable p1 equal "v_strain"
variable p2 equal "-pxx/10000"
variable sf equal "lz/16.54"
variable m2 equal "v_p2*v_sf"

dump 2 all custom 1000 dump.defo.* id type x y z c_csym c_2[1] c_2[2] c_2[3] c_2[4] c_2[5] c_2[6]
thermo_modify lost ignore flush yes
fix def_print2 all print 1000 "${p1} ${m2}" file 600K_GWMG_B_x.txt screen no
run 500000 
print "DONE"