####################################################################################################### ####################################################################################################### ############################ATOM TYPE-MASS-CHARGE ####################################################################################################### ####################################################################################################### #@ 3 atom types variable C3 equal 1 #cnt carbon variable C4 equal 2 #polymer carbon variable H equal 3 #polymer hydrogen ###########################ATOM MASSES mass ${C3} 12.010 mass ${C4} 12.010 mass ${H} 1.010 ##########################ATOM CHARGES #set type ${C3} charge 0.0000 #set type ${C4} charge -0.1060 #set type ${H} charge 0.0530 ####################################################################################################### ####################################################################################################### ############################BOND TYPE-PARAMETERS ####################################################################################################### ####################################################################################################### #@ 3 bond types bond_style class2 #C3-C3 bond_coeff 1 1.4170 470.8661 -627.6179 1327.6345 #C4-C4 bond_coeff 2 1.5300 299.6700 -501.7700 679.8100 #C4-H1 bond_coeff 3 1.1010 345.0000 -691.8900 844.6000 ####################################################################################################### ####################################################################################################### ############################ANGLE TYPE-PARAMETERS ####################################################################################################### ####################################################################################################### #@ 4 angle types angle_style class2 #C3-C3-C3 angle_coeff 1 118.9000 61.0226 -34.9931 0.0000 angle_coeff 1 bb 68.2856 1.4170 1.4170 angle_coeff 1 ba 28.8708 28.8708 1.4170 1.4170 #C4-C4-H1 angle_coeff 2 110.7700 41.4530 -10.6040 5.1290 angle_coeff 2 bb 3.3872 1.5300 1.1010 angle_coeff 2 ba 20.7540 11.4210 1.5300 1.1010 #C4-C4-C4 angle_coeff 3 112.6700 39.5160 -7.4430 -9.5583 angle_coeff 3 bb 0.0000 1.5300 1.5300 angle_coeff 3 ba 8.0160 8.0160 1.5300 1.5300 #H1-C4-H1 angle_coeff 4 107.6600 39.6410 -12.9210 -2.4318 angle_coeff 4 bb 5.3316 1.1010 1.1010 angle_coeff 4 ba 18.1030 18.1030 1.1010 1.1010 ####################################################################################################### ####################################################################################################### ############################DIHEDRAL TYPE-PARAMETERS ####################################################################################################### ####################################################################################################### #@4 dihedral type dihedral_style class2 #C3-C3-C3-C3 dihedral_coeff 1 8.3667 0.0000 1.2000 0.0000 0.0000 0.0000 dihedral_coeff 1 mbt 27.5989 -2.3120 0.0000 1.4170 dihedral_coeff 1 ebt -0.1185 6.3204 0.0000 -0.1185 6.3204 0.0000 1.4170 1.4170 dihedral_coeff 1 at 1.9767 1.0239 0.0000 1.9767 1.0239 0.0000 118.9000 118.9000 dihedral_coeff 1 aat 0.0000 118.9000 118.9000 dihedral_coeff 1 bb13 53.000 1.4170 1.4170 #C4-C4-C4-H1 dihedral_coeff 2 0.0000 0.0000 0.0316 0.0000 -0.1681 0.0000 dihedral_coeff 2 mbt -14.8790 -3.6581 -0.3138 1.5300 dihedral_coeff 2 ebt 0.2486 0.2422 -0.0925 0.0814 0.0591 0.2219 1.5300 1.1010 dihedral_coeff 2 at -0.2454 0.0000 -0.1136 0.3113 0.4516 -0.1988 112.6700 110.7700 dihedral_coeff 2 aat -16.1640 112.6700 110.7700 dihedral_coeff 2 bb13 0.0000 1.5300 1.1010 #H1-C4-C4-H1 dihedral_coeff 3 -0.1432 0.0000 0.0617 0.0000 -0.1530 0.0000 dihedral_coeff 3 mbt -14.2610 -0.5322 -0.4864 1.5300 dihedral_coeff 3 ebt 0.2130 0.3120 0.0777 0.2130 0.3120 0.0777 1.1010 1.1010 dihedral_coeff 3 at -0.8085 0.5569 -0.2466 -0.8085 0.5569 -0.2466 110.7700 110.7700 dihedral_coeff 3 aat -12.5640 110.7700 110.7700 dihedral_coeff 3 bb13 0.0000 1.1010 1.1010 #C4-C4-C4-C4 dihedral_coeff 4 0.0000 0.0000 0.0514 0.0000 -0.1430 0.0000 dihedral_coeff 4 mbt -17.7870 -7.1877 0.0000 1.5300 dihedral_coeff 4 ebt -0.0732 0.0000 0.0000 -0.0732 0.0000 0.0000 1.5300 1.5300 dihedral_coeff 4 at 0.3886 -0.3139 0.1389 0.3886 -0.3139 0.1389 112.6700 112.6700 dihedral_coeff 4 aat -22.0450 112.6700 112.6700 dihedral_coeff 4 bb13 0.0000 1.5300 1.5300 ####################################################################################################### ####################################################################################################### ############################IMPROPER TYPE-PARAMETERS ####################################################################################################### ####################################################################################################### #@4 improper types improper_style class2 #C3-C3-C3-C3 improper_coeff 1 7.1794 0.0000 improper_coeff 1 aa 0.0000 0.0000 0.0000 118.9000 118.9000 118.9000 #C4-C4-C4-H1 improper_coeff 2 0.0000 0.0000 improper_coeff 2 aa 0.2738 -0.4825 0.2738 110.7700 107.6600 110.7700 #H1-C4-C4-H1 improper_coeff 3 0.0000 0.0000 improper_coeff 3 aa -1.3199 -1.3199 0.1184 112.6700 110.7700 110.7700 #C4-C4-C4-C4 improper_coeff 4 0.0000 0.0000 improper_coeff 4 aa -0.3157 -0.3157 -0.3157 107.6600 107.6600 107.6600 ####################################################################################################### ####################################################################################################### ############################PAIR POTENTIALS ####################################################################################################### ####################################################################################################### pair_style hybrid/overlay lj/class2 12.0 coul/msm 12.0 #####CNT INTERACTION pair_coeff * * coul/msm #####POLYMER CNT INTERACTION pair_coeff ${C3} ${C3} lj/class2 0.06400 4.0100 pair_coeff ${C3} ${C4} lj/class2 0.05878 4.0100 pair_coeff ${C3} ${H} lj/class2 0.03577 3.5025 #####POLYMER INTERACTION pair_coeff ${C4} ${C4} lj/class2 0.054 4.0100 pair_coeff ${H} ${H} lj/class2 0.0200 2.9950 pair_coeff ${C4} ${H} lj/class2 0.03286 3.5025 kspace_style ewald 1.0e-4