File: README

package info (click to toggle)
lammps 20210122~gita77bb%2Bds1-2
  • links: PTS, VCS
  • area: main
  • in suites: bullseye
  • size: 331,996 kB
  • sloc: cpp: 802,213; python: 24,256; xml: 14,949; f90: 10,448; ansic: 8,476; perl: 4,161; sh: 3,466; fortran: 2,805; makefile: 1,250; objc: 238; lisp: 163; csh: 16; awk: 14; tcl: 6
file content (88 lines) | stat: -rw-r--r-- 2,745 bytes parent folder | download | duplicates (3)
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
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
This is a simple example of showing the computation of virtual x-ray
and electron diffraction patterns for Ni.  

In addition to the LAMMPS output, a simple visualizaiton of the electron 
diffraction data is included.  The visualization was created using the VisIt 
Open Source visualiztion software develeoped at Larrence Livermore National 
Laboraatory (https://wci.llnl.gov/simulation/computer-codes/visit/).

Below is a step by step guide to produce similar images:

a) Open the diffraction .vtk files in VisIt 

b) Click add->pseudocolor->intensity;

*Note, green name means its ready to be drawn

c) Click Draw;

*Note, you will see the full rectangular mesh grid which is what we
 are wanting.

We need to remove ghost data (those values set to -1).

d) Select intensity plot;

Click Operators->Selection->Isovolume;
Double click isovolume (under intensity plot);
Change lower bound to 0;
Click Apply;
Click Dismiss;
Click Draw;
 
*Note, now you will see only where there was valid intensity data (if you 
don't use Zone it will be a sphere).  But most intensity are low so color 
the intensity by log scale.

e) Select intensity plot

Double click Pseudocolor;
Change Scale to Log;
Check Mimimum Limit;
Input non-zero minimum limit (i.e. 1);
Click Apply;
Click Dismiss;

Note, clicking apply will change the scale automatically.  But unless
you have used Zone and a thin dR_Ewald you likely have lots of extra
data visualized that does not intersect the Ewald sphere.  Next, we
want to take a spherical slice to mimic the intersection of Ewald's
sphere with reciprocal space.

f) Select intensity plot;
Click Operators->Slicing>Spherical Slice;
Double click Spherical Slice;
Change the origin/radius to represent the center of your Ewald sphere;
      **(i.e. 39.84063 0 0) and radius (i.e. 39.84063) for lambda = 0.251 **
Click Apply;
Click Dismiss;
Click Draw;

This is the saed pattern, but we are not guaranteed to be aligned correctly 
with the zone axis.

g) On the menu bar... select Controls -> View;

Change View Normal to the zone axis of choice (i.e. -1 0 0);
Adjust Up Vector if desired (i.e. 0 1 0);
Turn off Perspective view;
Click Apply;
Click Dismiss;

*Note that for the primary x, y, and z zone axis you can use the drop 
own Align to axis.  However, it is likely that there is a bunch of extra 
annotations on the image that don't look appealing.

h) On the menu bar... select Control-> Annotation;

*Here you can change what is visualized along side the data and how
you display axis information.  I typically turn most all the display
information off.

i) Save the image

Click File-> Set save options
Click Save

** Note, further fine tuning can be achieved by adjusing the color table and 
max/min values **