reaxff forcefield

Dear LAMMPS users and developers

I simulate a reaction of CO and H2 on NI(111) surface using reaxff forcefiled.

However, the CO cannot absorb on the NI(111) surface using reaxff provide by "Phys. Rev. B 87, 104114 (2013)"

This reaxff forcefiled contain C H O and Ni four atom type

Now, I do not know the what wrong in my simulations
Forcefield
This is my input file
units real
dimension 3
atom_style charge
boundary p p f

read_data data.498200

pair_style reax/c lmp_control
pair_coeff * * ffield.reax.FC 1 2 9 3

fix 2 all qeq/reax 1 0.0 10.0 1e-6 reax/c

neighbor 5.0 bin
neigh_modify delay 10 check yes

#region mybottom block -2.000 800.000 -2.000 800.000 -1.000 5.5000 units box
#group bot region mybottom

group bot id 1 7 9 15 17 23 25 31 33 39 41 47 49 55 57 63 65 71 73 79 81 87 89 95 97 103 105 111 113 119 121 127 129 135 137 143 145 151 153 159 161 167 169 175 177 183 185 191 193 199 201 207 209 215 217 223 225 231 233 239 241 247 249 255 257 263 265 271 273 279 281 287 289 295 297 303 305 311 313 319 321 327 329 335 337 343 345 351 353 359 361 367 369 375 377 383 385 391 393 399 401 407 409 415 417 423 425 431 433 439 441 447 449 455 457 463 465 471 473 479 481 487 489 495 497 503 505 511 513 519 521 527 529 535 537 543 545 551 553 559 561 567 569 575 577 583 585 591 593 599 601 607 609 615 617 623 625 631 633 639 641 647 649 655 657 663 665 671 673 679 681 687 689 695 697 703 705 711 713 719 721 727 729 735 737 743 745 751 753 759 761 767 769 775 777 783 785 791 793 799 801 807 809 815 817 823 825 831 833 839 841 847 849 855 857 863 865 871 873 879 881 887 889 895 897 903 905 911 913 919 921 927 929 935 937 943 945 951 953 959 961 967 969 975 977 983 985 991 993 999 1001 1007 1009 1015 1017 1023 1025 1031 1033 1039 1041 1047 1049 1055 1057 1063 1065 1071 1073 1079 1081 1087 1089 1095 1097 1103 1105 1111 1113 1119 1121 1127 1129 1135 1137 1143 1145 1151 1153 1159 1161 1167 1169 1175 1177 1183 1185 1191 1193 1199 1201 1207 1209 1215 1217 1223 1225 1231 1233 1239 1241 1247 1249 1255 1257 1263 1265 1271 1273 1279 1281 1287 1289 1295 1297 1303 1305 1311 1313 1319 1321 1327 1329 1335 1337 1343 1345 1351 1353 1359 1361 1367 1369 1375 1377 1383 1385 1391 1393 1399 1401 1407 1409 1415 1417 1423 1425 1431 1433 1439 1441 1447 1449 1455 1457 1463 1465 1471 1473 1479 1481 1487 1489 1495 1497 1503 1505 1511 1513 1519 1521 1527 1529 1535 1537 1543 1545 1551 1553 1559 1561 1567 1569 1575 1577 1583
group uper id 2 3 4 5 6 8 10 11 12 13 14 16 18 19 20 21 22 24 26 27 28 29 30 32 34 35 36 37 38 40 42 43 44 45 46 48 50 51 52 53 54 56 58 59 60 61 62 64 66 67 68 69 70 72 74 75 76 77 78 80 82 83 84 85 86 88 90 91 92 93 94 96 98 99 100 101 102 104 106 107 108 109 110 112 114 115 116 117 118 120 122 123 124 125 126 128 130 131 132 133 134 136 138 139 140 141 142 144 146 147 148 149 150 152 154 155 156 157 158 160 162 163 164 165 166 168 170 171 172 173 174 176 178 179 180 181 182 184 186 187 188 189 190 192 194 195 196 197 198 200 202 203 204 205 206 208 210 211 212 213 214 216 218 219 220 221 222 224 226 227 228 229 230 232 234 235 236 237 238 240 242 243 244 245 246 248 250 251 252 253 254 256 258 259 260 261 262 264 266 267 268 269 270 272 274 275 276 277 278 280 282 283 284 285 286 288 290 291 292 293 294 296 298 299 300 301 302 304 306 307 308 309 310 312 314 315 316 317 318 320 322 323 324 325 326 328 330 331 332 333 334 336 338 339 340 341 342 344 346 347 348 349 350 352 354 355 356 357 358 360 362 363 364 365 366 368 370 371 372 373 374 376 378 379 380 381 382 384 386 387 388 389 390 392 394 395 396 397 398 400 402 403 404 405 406 408 410 411 412 413 414 416 418 419 420 421 422 424 426 427 428 429 430 432 434 435 436 437 438 440 442 443 444 445 446 448 450 451 452 453 454 456 458 459 460 461 462 464 466 467 468 469 470 472 474 475 476 477 478 480 482 483 484 485 486 488 490 491 492 493 494 496 498 499 500 501 502 504 506 507 508 509 510 512 514 515 516 517 518 520 522 523 524 525 526 528 530 531 532 533 534 536 538 539 540 541 542 544 546 547 548 549 550 552 554 555 556 557 558 560 562 563 564 565 566 568 570 571 572 573 574 576 578 579 580 581 582 584 586 587 588 589 590 592 594 595 596 597 598 600 602 603 604 605 606 608 610 611 612 613 614 616 618 619 620 621 622 624 626 627 628 629 630 632 634 635 636 637 638 640 642 643 644 645 646 648 650 651 652 653 654 656 658 659 660 661 662 664 666 667 668 669 670 672 674 675 676 677 678 680 682 683 684 685 686 688 690 691 692 693 694 696 698 699 700 701 702 704 706 707 708 709 710 712 714 715 716 717 718 720 722 723 724 725 726 728 730 731 732 733 734 736 738 739 740 741 742 744 746 747 748 749 750 752 754 755 756 757 758 760 762 763 764 765 766 768 770 771 772 773 774 776 778 779 780 781 782 784 786 787 788 789 790 792 794 795 796 797 798 800 802 803 804 805 806 808 810 811 812 813 814 816 818 819 820 821 822 824 826 827 828 829 830 832 834 835 836 837 838 840 842 843 844 845 846 848 850 851 852 853 854 856 858 859 860 861 862 864 866 867 868 869 870 872 874 875 876 877 878 880 882 883 884 885 886 888 890 891 892 893 894 896 898 899 900 901 902 904 906 907 908 909 910 912 914 915 916 917 918 920 922 923 924 925 926 928 930 931 932 933 934 936 938 939 940 941 942 944 946 947 948 949 950 952 954 955 956 957 958 960 962 963 964 965 966 968 970 971 972 973 974 976 978 979 980 981 982 984 986 987 988 989 990 992 994 995 996 997 998 1000 1002 1003 1004 1005 1006 1008 1010 1011 1012 1013 1014 1016 1018 1019 1020 1021 1022 1024 1026 1027 1028 1029 1030 1032 1034 1035 1036 1037 1038 1040 1042 1043 1044 1045 1046 1048 1050 1051 1052 1053 1054 1056 1058 1059 1060 1061 1062 1064 1066 1067 1068 1069 1070 1072 1074 1075 1076 1077 1078 1080 1082 1083 1084 1085 1086 1088 1090 1091 1092 1093 1094 1096 1098 1099 1100 1101 1102 1104 1106 1107 1108 1109 1110 1112 1114 1115 1116 1117 1118 1120 1122 1123 1124 1125 1126 1128 1130 1131 1132 1133 1134 1136 1138 1139 1140 1141 1142 1144 1146 1147 1148 1149 1150 1152 1154 1155 1156 1157 1158 1160 1162 1163 1164 1165 1166 1168 1170 1171 1172 1173 1174 1176 1178 1179 1180 1181 1182 1184 1186 1187 1188 1189 1190 1192 1194 1195 1196 1197 1198 1200 1202 1203 1204 1205 1206 1208 1210 1211 1212 1213 1214 1216 1218 1219 1220 1221 1222 1224 1226 1227 1228 1229 1230 1232 1234 1235 1236 1237 1238 1240 1242 1243 1244 1245 1246 1248 1250 1251 1252 1253 1254 1256 1258 1259 1260 1261 1262 1264 1266 1267 1268 1269 1270 1272 1274 1275 1276 1277 1278 1280 1282 1283 1284 1285 1286 1288 1290 1291 1292 1293 1294 1296 1298 1299 1300 1301 1302 1304 1306 1307 1308 1309 1310 1312 1314 1315 1316 1317 1318 1320 1322 1323 1324 1325 1326 1328 1330 1331 1332 1333 1334 1336 1338 1339 1340 1341 1342 1344 1346 1347 1348 1349 1350 1352 1354 1355 1356 1357 1358 1360 1362 1363 1364 1365 1366 1368 1370 1371 1372 1373 1374 1376 1378 1379 1380 1381 1382 1384 1386 1387 1388 1389 1390 1392 1394 1395 1396 1397 1398 1400 1402 1403 1404 1405 1406 1408 1410 1411 1412 1413 1414 1416 1418 1419 1420 1421 1422 1424 1426 1427 1428 1429 1430 1432 1434 1435 1436 1437 1438 1440 1442 1443 1444 1445 1446 1448 1450 1451 1452 1453 1454 1456 1458 1459 1460 1461 1462 1464 1466 1467 1468 1469 1470 1472 1474 1475 1476 1477 1478 1480 1482 1483 1484 1485 1486 1488 1490 1491 1492 1493 1494 1496 1498 1499 1500 1501 1502 1504 1506 1507 1508 1509 1510 1512 1514 1515 1516 1517 1518 1520 1522 1523 1524 1525 1526 1528 1530 1531 1532 1533 1534 1536 1538 1539 1540 1541 1542 1544 1546 1547 1548 1549 1550 1552 1554 1555 1556 1557 1558 1560 1562 1563 1564 1565 1566 1568 1570 1571 1572 1573 1574 1576 1578 1579 1580 1581 1582 1584 1585 1586 1587 1588 1589 1590 1591 1592 1593 1594 1595 1596 1597 1598 1599 1600 1601 1602 1603 1604 1605 1606 1607 1608 1609 1610 1611 1612 1613 1614 1615 1616 1617 1618 1619 1620 1621 1622 1623 1624 1625 1626 1627 1628 1629 1630 1631 1632 1633 1634 1635 1636 1637 1638 1639 1640 1641 1642 1643 1644 1645 1646 1647 1648 1649 1650 1651 1652 1653 1654 1655 1656 1657 1658 1659 1660 1661 1662 1663 1664 1665 1666 1667 1668 1669 1670 1671 1672 1673 1674 1675 1676 1677 1678 1679 1680 1681 1682 1683 1684 1685 1686 1687 1688 1689 1690 1691 1692 1693 1694 1695 1696 1697 1698 1699 1700 1701 1702 1703 1704 1705 1706 1707 1708 1709 1710 1711 1712 1713 1714 1715 1716 1717 1718 1719 1720 1721 1722 1723 1724 1725 1726 1727 1728 1729 1730 1731 1732 1733 1734 1735 1736 1737 1738 1739 1740 1741 1742 1743 1744 1745 1746 1747 1748 1749 1750 1751 1752 1753 1754 1755 1756 1757 1758 1759 1760 1761 1762 1763 1764

reset_timestep 0
#min_style cg
#minimize 1e-25 1e-25 5000 10000
velocity all create 1000 4559 rot yes dist gaussian

dump myrec all atom 25000 4981000f.lammpstrj
thermo_style custom step etotal temp press
thermo_modify flush yes

fix 3 all wall/reflect zlo EDGE zhi EDGE
fix 1 uper nvt temp 1000.0 1000.0 5.0
fix feeze bot setforce 0.0 0.0 0.0

thermo 10
timestep 0.25

      restart 20000 restart.*.equil

      run 100000000

Dear Wang

You mentioned that ReaxFF contain C H O Ni but in the input file you have mentioned 9 in the pair_coeff.
pair_coeff * * ffield.reax.FC 1 2 9 3

Best

vasu

Dear Wang

You mentioned that ReaxFF contain C H O Ni but in the input file you have mentioned 9 in the pair_coeff.
pair_coeff * * ffield.reax.FC 1 2 9 3

Not necessarily since we do not know how these four elements are being assigned in the actual data file…

Best

vasu

This particular paper and parameterization was developed for another purpose (which is titled: Thermal properties of fluorinated graphene). Are you sure this force field description is suitable for CO/H2 adsorption on Ni surfaces?