# 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 #neigh_modify delay 5 atom_style molecular pair_style sw units metal #lectura de estructura read_data /home/lsi/bernardo/pruebas/MINIMICENUEVO/INASstress/data.INAS555 group bloque id <> 1 1166 #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 7 13 19 25 31 37 43 49 55 61 67 73 79 85 91 97 103 109 115 121 127 133 139 145 151 157 163 169 175 181 187 193 199 205 211 517 523 529 535 541 547 553 559 565 571 577 583 589 595 601 607 613 619 625 631 637 643 649 655 661 #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 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 group xlow id 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 group ylow id 1 2 3 4 5 6 37 38 39 40 41 42 73 74 75 76 77 78 109 110 111 112 113 114 145 146 147 148 149 150 181 182 183 184 185 186 367 368 369 370 371 397 398 399 400 401 427 428 429 430 431 457 458 459 460 461 487 488 489 490 491 group yhigh id 31 32 33 34 35 36 67 68 69 70 71 72 103 104 105 106 107 108 139 140 141 142 143 144 175 176 177 178 179 180 211 212 213 214 215 216 392 393 394 395 396 422 423 424 425 426 452 453 454 455 456 482 483 484 485 486 512 513 514 515 516 pair_coeff * * /home/lsi/bernardo/pruebas/MINIMICENUEVO/INASstress/InGaAs.sw In As # 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] vol press lx ly lz pxx pyy pzz pxy pxz pyz dump 1 all atom 250 dump.inas555 dump 2 all custom 100 tmp.stress tag c_2[1] c_2[2] c_2[3] #run 10000 #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