# simulacion de estructura INAS pura comprimida (con paredes y sin paredes y luego sin comprimir) dimension 3 boundary s s s newton on neighbor 2 nsq processors 2 2 2 #neigh_modify delay 5 atom_style molecular pair_style sw units metal #lectura de estructura read_data /home/lsi/bernardo/data.prueba group bloque id <> 1 388772 #esto se coloca para que el suelo y el techo sean fijos z, se ha comprobado quitandolo del grupo #bloque y la simulacion es la misma no hay cambios. group suelo id <> 1 5372 #group techo id 6 12 18 24 30 36 42 48 54 60 66 72 78 84 90 96 102 108 114 120 126 132 138 144 150 156 162 168 174 180 186 192 198 204 210 216 522 528 534 540 546 552 558 564 570 576 582 588 594 600 606 612 618 624 630 636 642 648 654 660 666 #bloques para fijar ejes x e y #group xhigh id 1009 1010 1011 1012 1053 1054 1055 1056 3259 3260 3261 3262 3303 3304 3305 3306 #group xlow id 1 2 3 4 13 14 15 16 33 34 35 36 61 62 63 64 97 98 99 100 141 142 143 144 193 194 195 196 253 254 255 256 1157 1158 1159 1160 1173 1174 1175 1176 1197 1198 1199 1200 1229 1230 1231 1232 1269 1270 1271 1272 1317 1318 1319 1320 1373 1374 1375 1376 2276 2277 2278 2279 2288 2289 2290 2291 2308 2309 2310 2311 2336 2337 2338 2339 2372 2373 2374 2375 2416 2417 2418 2419 2468 2469 2470 2471 2528 2529 2530 2531 3406 3407 3408 3418 3419 3420 3436 3437 3438 3460 3461 3462 3490 3491 3492 3526 3527 3528 3568 3569 3570 #group ylow id 1153 1154 1155 1156 1161 1162 1163 1164 1177 1178 1179 1180 1201 1202 1203 1204 1233 1234 1235 1236 1273 1274 1275 1276 1321 1322 1323 1324 1377 1378 1379 1380 1437 1438 1439 1440 1497 1498 1499 1500 1557 1558 1559 1560 1632 1633 1634 1635 1707 1708 1709 1710 1782 1783 1784 1785 1857 1858 1859 1860 1932 1933 1934 1935 1992 1993 1994 1995 3403 3404 3405 3409 3410 3411 3421 3422 3423 3439 3440 3441 3463 3464 3465 3493 3494 3495 3529 3530 3531 3571 3572 3573 3616 3617 3618 3661 3662 3663 3706 3707 3708 3771 3772 3773 3836 3837 3838 3901 3902 3903 3966 3967 3968 4031 4032 4033 4076 4077 4078 #group yhigh id 253 254 255 256 313 314 315 316 373 374 375 376 453 454 455 456 533 534 535 536 613 614 615 616 693 694 695 696 773 774 775 776 833 834 835 836 893 894 895 896 953 954 955 956 1005 1006 1007 1008 1049 1050 1051 1052 1085 1086 1087 1088 1113 1114 1115 1116 1133 1134 1135 1136 1145 1146 1147 1148 1149 1150 1151 1152 2528 2529 2530 2531 2588 2589 2590 2591 2648 2649 2650 2651 2723 2724 2725 2726 2798 2799 2800 2801 2873 2874 2875 2876 2948 2949 2950 2951 3023 3024 3025 3026 3083 3084 3085 3086 3143 3144 3145 3146 3203 3204 3205 3206 3255 3256 3257 3258 3299 3300 3301 3302 3335 3336 3337 3338 3363 3364 3365 3366 3383 3384 3385 3386 3395 3396 3397 3398 3399 3400 3401 3402 pair_coeff * * /home/lsi/bernardo/GANINN.sw Ga N In # initial velocities compute new bloque temp #compute 2 all stress/atom velocity bloque create 0.1 888823 temp new # fixes #los setforce 0 se colocan para que el suelo y techo en z sean fijos. fix 2 suelo setforce 0.0 0.0 0.0 #fix 3 techo setforce NULL NULL 0.0 #si quitamos los de arriba dejamos libre la simulacion sin paredes fijas #fix 4 xhigh setforce 0.0 NULL NULL #fix 5 xlow setforce 0.0 NULL NULL #fix 6 ylow setforce NULL 0.0 NULL #fix 7 yhigh setforce NULL 0.0 NULL #a continuacion se ponen los fix para dejar fijas las paredes de x y z son 4 fix. fix 8 all nve # run thermo 1000 thermo_style custom step temp c_new[1] c_new[2] c_new[3] c_new[4] c_new[5] c_new[6] vol press lx ly lz pxx pyy pzz pxy pxz pyz dump 1 all atom 10000 dump.ganinn #dump 2 all custom 250 tmp.stress tag c_2[1] c_2[2] c_2[3] c_2[4] c_2[5] c_2[6] #run 1000 #dump_modify 1 every 5 #thermo 5 #min_style cg/fr #min_modify lineiter 5 #min_modify linestyle scan dmin 0.0001 dmax 0.1 minimize 0 1000 1000000