|
395 | 395 | },
|
396 | 396 | {
|
397 | 397 | "cell_type": "code",
|
398 |
| - "execution_count": 11, |
| 398 | + "execution_count": 13, |
399 | 399 | "id": "4898aaa3-926a-4c15-9c4f-29acd0a3425a",
|
400 | 400 | "metadata": {},
|
401 | 401 | "outputs": [],
|
402 | 402 | "source": [
|
403 | 403 | "# Do the same with water\n",
|
404 |
| - "wat = u.select_atoms(\"resname SOL\")\n", |
| 404 | + "wat = u.select_atoms(\"resname SOL and name OW\")\n", |
405 | 405 | "site_occu_wat = []\n",
|
406 | 406 | "for ts in u.trajectory:\n",
|
407 | 407 | " wat_array = sf.state_detect(wat, s5_z_cutoff, r_cutoff, s0_r_cutoff) \n",
|
|
410 | 410 | },
|
411 | 411 | {
|
412 | 412 | "cell_type": "code",
|
413 |
| - "execution_count": 12, |
| 413 | + "execution_count": 14, |
414 | 414 | "id": "421c3e83-05f9-4b31-84c9-be3950713483",
|
415 | 415 | "metadata": {},
|
416 | 416 | "outputs": [
|
|
446 | 446 | " <tbody>\n",
|
447 | 447 | " <tr>\n",
|
448 | 448 | " <th>0</th>\n",
|
449 |
| - " <td>[56425, 60299, 60300, 60301]</td>\n", |
| 449 | + " <td>[60299]</td>\n", |
450 | 450 | " <td>[]</td>\n",
|
451 | 451 | " <td>[]</td>\n",
|
452 | 452 | " <td>[]</td>\n",
|
453 | 453 | " <td>[]</td>\n",
|
454 |
| - " <td>[31255, 32579, 32580, 32581, 40145, 40146, 401...</td>\n", |
| 454 | + " <td>[40145, 44588]</td>\n", |
455 | 455 | " </tr>\n",
|
456 | 456 | " <tr>\n",
|
457 | 457 | " <th>1</th>\n",
|
|
460 | 460 | " <td>[]</td>\n",
|
461 | 461 | " <td>[]</td>\n",
|
462 | 462 | " <td>[]</td>\n",
|
463 |
| - " <td>[36476, 36477, 39977, 39978, 39979, 41862, 522...</td>\n", |
| 463 | + " <td>[61961]</td>\n", |
464 | 464 | " </tr>\n",
|
465 | 465 | " <tr>\n",
|
466 | 466 | " <th>2</th>\n",
|
|
469 | 469 | " <td>[]</td>\n",
|
470 | 470 | " <td>[]</td>\n",
|
471 | 471 | " <td>[]</td>\n",
|
472 |
| - " <td>[36476, 36477, 36478, 37471, 41862, 47115, 490...</td>\n", |
| 472 | + " <td>[49097, 51185, 51449]</td>\n", |
473 | 473 | " </tr>\n",
|
474 | 474 | " <tr>\n",
|
475 | 475 | " <th>3</th>\n",
|
476 |
| - " <td>[35385, 35386, 48137, 48138]</td>\n", |
| 476 | + " <td>[48137]</td>\n", |
477 | 477 | " <td>[]</td>\n",
|
478 | 478 | " <td>[]</td>\n",
|
479 | 479 | " <td>[]</td>\n",
|
480 | 480 | " <td>[]</td>\n",
|
481 |
| - " <td>[31295, 31296, 31297, 37469, 37471, 40427, 404...</td>\n", |
| 481 | + " <td>[31295, 43874, 55148, 58355]</td>\n", |
482 | 482 | " </tr>\n",
|
483 | 483 | " <tr>\n",
|
484 | 484 | " <th>4</th>\n",
|
485 |
| - " <td>[35384, 35385, 35386, 42595, 50423, 50424]</td>\n", |
| 485 | + " <td>[35384, 50423]</td>\n", |
486 | 486 | " <td>[]</td>\n",
|
487 | 487 | " <td>[]</td>\n",
|
488 | 488 | " <td>[]</td>\n",
|
489 | 489 | " <td>[]</td>\n",
|
490 |
| - " <td>[37016, 37017, 37018, 38117, 38118, 38119, 399...</td>\n", |
| 490 | + " <td>[37016, 49463]</td>\n", |
491 | 491 | " </tr>\n",
|
492 | 492 | " <tr>\n",
|
493 | 493 | " <th>5</th>\n",
|
494 |
| - " <td>[42764, 42765, 42766, 47197, 53675, 53676, 53677]</td>\n", |
| 494 | + " <td>[42764, 53675]</td>\n", |
495 | 495 | " <td>[]</td>\n",
|
496 | 496 | " <td>[]</td>\n",
|
497 | 497 | " <td>[]</td>\n",
|
498 | 498 | " <td>[]</td>\n",
|
499 |
| - " <td>[30694, 36476, 36477, 36478, 36844, 37470, 384...</td>\n", |
| 499 | + " <td>[40154, 48152, 58601]</td>\n", |
500 | 500 | " </tr>\n",
|
501 | 501 | " <tr>\n",
|
502 | 502 | " <th>6</th>\n",
|
503 |
| - " <td>[55301, 55302]</td>\n", |
| 503 | + " <td>[55301]</td>\n", |
504 | 504 | " <td>[]</td>\n",
|
505 | 505 | " <td>[]</td>\n",
|
506 | 506 | " <td>[]</td>\n",
|
507 | 507 | " <td>[]</td>\n",
|
508 |
| - " <td>[30692, 30693, 30694, 38414, 38415, 38416, 390...</td>\n", |
| 508 | + " <td>[38414, 39023, 59456]</td>\n", |
509 | 509 | " </tr>\n",
|
510 | 510 | " <tr>\n",
|
511 | 511 | " <th>7</th>\n",
|
512 |
| - " <td>[41922, 41923, 43610, 43611, 43612]</td>\n", |
| 512 | + " <td>[43610]</td>\n", |
513 | 513 | " <td>[]</td>\n",
|
514 | 514 | " <td>[]</td>\n",
|
515 | 515 | " <td>[]</td>\n",
|
516 | 516 | " <td>[]</td>\n",
|
517 |
| - " <td>[30692, 30693, 30694, 32579, 32580, 32581, 364...</td>\n", |
| 517 | + " <td>[37358, 39977, 40088]</td>\n", |
518 | 518 | " </tr>\n",
|
519 | 519 | " <tr>\n",
|
520 | 520 | " <th>8</th>\n",
|
521 |
| - " <td>[37128, 41121, 41122, 54627]</td>\n", |
522 | 521 | " <td>[]</td>\n",
|
523 | 522 | " <td>[]</td>\n",
|
524 | 523 | " <td>[]</td>\n",
|
525 | 524 | " <td>[]</td>\n",
|
526 |
| - " <td>[32580, 39023, 39025, 40088, 40089, 40090, 547...</td>\n", |
| 525 | + " <td>[]</td>\n", |
| 526 | + " <td>[40088, 55916]</td>\n", |
527 | 527 | " </tr>\n",
|
528 | 528 | " <tr>\n",
|
529 | 529 | " <th>9</th>\n",
|
530 |
| - " <td>[33339, 37620, 41318, 41319, 41320]</td>\n", |
| 530 | + " <td>[41318]</td>\n", |
531 | 531 | " <td>[]</td>\n",
|
532 | 532 | " <td>[]</td>\n",
|
533 | 533 | " <td>[]</td>\n",
|
534 | 534 | " <td>[]</td>\n",
|
535 |
| - " <td>[30775, 33674, 33675, 33676, 37121, 37122, 371...</td>\n", |
| 535 | + " <td>[33674, 37121]</td>\n", |
536 | 536 | " </tr>\n",
|
537 | 537 | " <tr>\n",
|
538 | 538 | " <th>10</th>\n",
|
539 |
| - " <td>[34520, 38763]</td>\n", |
| 539 | + " <td>[34520]</td>\n", |
540 | 540 | " <td>[]</td>\n",
|
541 | 541 | " <td>[]</td>\n",
|
542 | 542 | " <td>[]</td>\n",
|
543 | 543 | " <td>[]</td>\n",
|
544 |
| - " <td>[32392, 33140, 33141, 33142, 33674, 33675, 336...</td>\n", |
| 544 | + " <td>[33140, 33674, 55745]</td>\n", |
545 | 545 | " </tr>\n",
|
546 | 546 | " </tbody>\n",
|
547 | 547 | "</table>\n",
|
548 | 548 | "</div>"
|
549 | 549 | ],
|
550 | 550 | "text/plain": [
|
551 |
| - " S0 S1 S2 S3 S4 \\\n", |
552 |
| - "0 [56425, 60299, 60300, 60301] [] [] [] [] \n", |
553 |
| - "1 [] [] [] [] [] \n", |
554 |
| - "2 [] [] [] [] [] \n", |
555 |
| - "3 [35385, 35386, 48137, 48138] [] [] [] [] \n", |
556 |
| - "4 [35384, 35385, 35386, 42595, 50423, 50424] [] [] [] [] \n", |
557 |
| - "5 [42764, 42765, 42766, 47197, 53675, 53676, 53677] [] [] [] [] \n", |
558 |
| - "6 [55301, 55302] [] [] [] [] \n", |
559 |
| - "7 [41922, 41923, 43610, 43611, 43612] [] [] [] [] \n", |
560 |
| - "8 [37128, 41121, 41122, 54627] [] [] [] [] \n", |
561 |
| - "9 [33339, 37620, 41318, 41319, 41320] [] [] [] [] \n", |
562 |
| - "10 [34520, 38763] [] [] [] [] \n", |
563 |
| - "\n", |
564 |
| - " SCav \n", |
565 |
| - "0 [31255, 32579, 32580, 32581, 40145, 40146, 401... \n", |
566 |
| - "1 [36476, 36477, 39977, 39978, 39979, 41862, 522... \n", |
567 |
| - "2 [36476, 36477, 36478, 37471, 41862, 47115, 490... \n", |
568 |
| - "3 [31295, 31296, 31297, 37469, 37471, 40427, 404... \n", |
569 |
| - "4 [37016, 37017, 37018, 38117, 38118, 38119, 399... \n", |
570 |
| - "5 [30694, 36476, 36477, 36478, 36844, 37470, 384... \n", |
571 |
| - "6 [30692, 30693, 30694, 38414, 38415, 38416, 390... \n", |
572 |
| - "7 [30692, 30693, 30694, 32579, 32580, 32581, 364... \n", |
573 |
| - "8 [32580, 39023, 39025, 40088, 40089, 40090, 547... \n", |
574 |
| - "9 [30775, 33674, 33675, 33676, 37121, 37122, 371... \n", |
575 |
| - "10 [32392, 33140, 33141, 33142, 33674, 33675, 336... " |
| 551 | + " S0 S1 S2 S3 S4 SCav\n", |
| 552 | + "0 [60299] [] [] [] [] [40145, 44588]\n", |
| 553 | + "1 [] [] [] [] [] [61961]\n", |
| 554 | + "2 [] [] [] [] [] [49097, 51185, 51449]\n", |
| 555 | + "3 [48137] [] [] [] [] [31295, 43874, 55148, 58355]\n", |
| 556 | + "4 [35384, 50423] [] [] [] [] [37016, 49463]\n", |
| 557 | + "5 [42764, 53675] [] [] [] [] [40154, 48152, 58601]\n", |
| 558 | + "6 [55301] [] [] [] [] [38414, 39023, 59456]\n", |
| 559 | + "7 [43610] [] [] [] [] [37358, 39977, 40088]\n", |
| 560 | + "8 [] [] [] [] [] [40088, 55916]\n", |
| 561 | + "9 [41318] [] [] [] [] [33674, 37121]\n", |
| 562 | + "10 [34520] [] [] [] [] [33140, 33674, 55745]" |
576 | 563 | ]
|
577 | 564 | },
|
578 |
| - "execution_count": 12, |
| 565 | + "execution_count": 14, |
579 | 566 | "metadata": {},
|
580 | 567 | "output_type": "execute_result"
|
581 | 568 | }
|
|
0 commit comments