|
98 | 98 | "cell_type": "markdown", |
99 | 99 | "metadata": {}, |
100 | 100 | "source": [ |
| 101 | + "From looking at the path tree, we can get a sense of how well our sampling is going. Some questions to consider:\n", |
| 102 | + "\n", |
101 | 103 | "* Do you see acceptance of both forward and backward shots?\n", |
102 | 104 | "* How many decorrelated paths do you see? \n", |
103 | | - "* How many accepted moves between decorrelated paths? How many total moves?" |
| 105 | + "* How many accepted moves between decorrelated paths? How many total moves?\n", |
| 106 | + "\n", |
| 107 | + "If you don't see acceptance of both forward and backward shots, that might indicate a problem in either the state definitions or in the sampling (which might be remedied by using a more advanced shooting point selection algorithm). The number MC steps between decorrelated paths gives a sense of how quickly your simulation is exploring path space. If you see very few decorrelated paths (or none at all), that may indicate that you should switch to a more advanced algorithm such as spring shooting." |
104 | 108 | ] |
105 | 109 | }, |
106 | 110 | { |
|
152 | 156 | "```" |
153 | 157 | ] |
154 | 158 | }, |
| 159 | + { |
| 160 | + "cell_type": "markdown", |
| 161 | + "metadata": {}, |
| 162 | + "source": [ |
| 163 | + "Next let's visualize the shooting point from one of your MC steps. There are two reasons we do this: First, this gives you more practice visualizing OPS data with NGLView. Second, this shows you how to extract the shooting point snapshot from a MC step. Analyzing the shooting points can give ideas on how to troubleshoot problems with sampling (e.g., if you had very few accepted backward shots)." |
| 164 | + ] |
| 165 | + }, |
155 | 166 | { |
156 | 167 | "cell_type": "code", |
157 | 168 | "execution_count": null, |
|
326 | 337 | "cell_type": "markdown", |
327 | 338 | "metadata": {}, |
328 | 339 | "source": [ |
329 | | - "* Estimate (by eye) what the path length should be if you ensured that fixed-length TPS would include 99% of the paths. How much longer is that than the average path length here?" |
| 340 | + "* Estimate (by eye) what the path length should be if you ensured that fixed-length TPS would include 99% of the paths. How much longer is that than the average path length here? What does that mean about the simulation time required for fixed-length TPS as compared the flexible-length TPS?" |
330 | 341 | ] |
331 | 342 | }, |
332 | 343 | { |
|
407 | 418 | "* Which one appears to be more frequent? \n", |
408 | 419 | "* Do they look like completely separate channels, or does it look like some paths involves parts of both?\n", |
409 | 420 | "\n", |
410 | | - "Next, you'll find a representative trajectory for each channel. This will just involve trial and error until you find one." |
| 421 | + "Next, you'll find a representative trajectory for each channel. This will just involve trial and error until you find one. You're likely to also find trajectories that at least partly cross into the other channel -- these two channels are not very strictly separated." |
411 | 422 | ] |
412 | 423 | }, |
413 | 424 | { |
|
438 | 449 | "cell_type": "markdown", |
439 | 450 | "metadata": {}, |
440 | 451 | "source": [ |
441 | | - "If you have the full simulation data, you could visualize these trajectories with NGLView. (Full simulation data not provided online, due to size), using the same approach as above." |
| 452 | + "If you have the full simulation data, you could visualize these trajectories with NGLView, using the same approach as above. (Full simulation data not provided online, due to size.)" |
442 | 453 | ] |
443 | 454 | }, |
444 | 455 | { |
|
465 | 476 | "name": "python", |
466 | 477 | "nbconvert_exporter": "python", |
467 | 478 | "pygments_lexer": "ipython3", |
468 | | - "version": "3.7.3" |
| 479 | + "version": "3.7.10" |
469 | 480 | }, |
470 | 481 | "toc": { |
471 | 482 | "base_numbering": 1, |
|
0 commit comments