git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@13739 f3b2605a-c512-4ea7-a41b-209d697bcdaa

This commit is contained in:
sjplimp 2015-07-27 18:12:57 +00:00
parent 9ed741e904
commit 317744626c
8 changed files with 69 additions and 78 deletions

View File

@ -69,11 +69,11 @@ is computed from the structure factor F using the equations:
position of each atom, fj are atomic scattering factors.
</P>
<P>Diffraction intensities are calculated on a three-dimensional mesh of
reciprocal lattice nodes. The mesh spacing is defined either (I) by
the entire simulation domain or (II) manually using selected values as
reciprocal lattice nodes. The mesh spacing is defined either (a) by
the entire simulation domain or (b) manually using selected values as
shown in the 2D diagram below.
</P>
<CENTER><A HREF = "saed_mesh.jpg"><IMG SRC = "JPG/saed_mesh_small.jpg"></A>
<CENTER><A HREF = "JPG/saed_mesh.jpg"><IMG SRC = "JPG/saed_mesh_small.jpg"></A>
</CENTER>
<P>For a mesh defined by the simulation domain, a rectilinear grid is
constructed with spacing <I>c</I>*inv(A) along each reciprocal lattice
@ -89,7 +89,7 @@ will defined using the <I>c</I> values for the spacing along each reciprocal
lattice axis. Note that manual mapping of the reciprocal space mesh is
good for comparing diffraction results from multiple simulations; however
it can reduce the likelihood that Bragg reflections will be satisfied
unless small spacing parameters <B><0.05 Angstrom^(-1)</B> are implemented.
unless small spacing parameters <0.05 Angstrom^(-1) are implemented.
Meshes with manual spacing do not require a periodic boundary.
</P>
<P>The limits of the reciprocal lattice mesh are determined by the use of
@ -103,9 +103,9 @@ denote the z1=h, z2=k, and z3=l (in a global since) zone axis of an
intersecting Ewald sphere. Diffraction intensities will only be
computed at the intersection of the reciprocal lattice mesh and a
<I>dR_Ewald</I> thick surface of the Ewald sphere. See the example 3D
intestiety data and the intersection of a <B>010</B> zone axis in the below image.
intestiety data and the intersection of a [010] zone axis in the below image.
</P>
<CENTER><A HREF = "saed_ewald_intersect.jpg"><IMG SRC = "JPG/saed_ewald_intersect_small.jpg"></A>
<CENTER><A HREF = "JPG/saed_ewald_intersect.jpg"><IMG SRC = "JPG/saed_ewald_intersect_small.jpg"></A>
</CENTER>
<P>The atomic scattering factors, fj, accounts for the reduction in
diffraction intensity due to Compton scattering. Compute saed uses
@ -139,7 +139,7 @@ type. Valid chemical symbols for compute saed are:
Tl: Pb: Bi: Po: At:
Rn: Fr: Ra: Ac: Th:
Pa: U: Np: Pu: Am:
Cm: Bk: Cf:
Cm: Bk: Cf:tb(c=5,s=:)
</P>
<P>If the <I>echo</I> keyword is specified, compute saed will provide extra
reporting information to the screen.
@ -160,11 +160,7 @@ options.
</P>
<P><B>Restrictions:</B>
</P>
<P>This command is part of the USER-DIFFRACTION package. It is only
enabled if LAMMPS was built with that package. See the <A HREF = "Section_start.html#start_3">Making
LAMMPS</A> section for more info.
</P>
<P>This command does not work for triclinic cells.
<P>The compute_saed command does not work for triclinic cells.
</P>
<P><B>Related commands:</B>
</P>

View File

@ -59,11 +59,11 @@ Here, K is the location of the reciprocal lattice node, rj is the
position of each atom, fj are atomic scattering factors.
Diffraction intensities are calculated on a three-dimensional mesh of
reciprocal lattice nodes. The mesh spacing is defined either (I) by
the entire simulation domain or (II) manually using selected values as
reciprocal lattice nodes. The mesh spacing is defined either (a) by
the entire simulation domain or (b) manually using selected values as
shown in the 2D diagram below.
:c,image(JPG/saed_mesh_small.jpg,saed_mesh.jpg)
:c,image(JPG/saed_mesh_small.jpg,JPG/saed_mesh.jpg)
For a mesh defined by the simulation domain, a rectilinear grid is
constructed with spacing {c}*inv(A) along each reciprocal lattice
@ -79,7 +79,7 @@ will defined using the {c} values for the spacing along each reciprocal
lattice axis. Note that manual mapping of the reciprocal space mesh is
good for comparing diffraction results from multiple simulations; however
it can reduce the likelihood that Bragg reflections will be satisfied
unless small spacing parameters [<0.05 Angstrom^(-1)] are implemented.
unless small spacing parameters <0.05 Angstrom^(-1) are implemented.
Meshes with manual spacing do not require a periodic boundary.
The limits of the reciprocal lattice mesh are determined by the use of
@ -93,9 +93,9 @@ denote the z1=h, z2=k, and z3=l (in a global since) zone axis of an
intersecting Ewald sphere. Diffraction intensities will only be
computed at the intersection of the reciprocal lattice mesh and a
{dR_Ewald} thick surface of the Ewald sphere. See the example 3D
intestiety data and the intersection of a [010] zone axis in the below image.
intestiety data and the intersection of a \[010\] zone axis in the below image.
:c,image(JPG/saed_ewald_intersect_small.jpg,saed_ewald_intersect.jpg)
:c,image(JPG/saed_ewald_intersect_small.jpg,JPG/saed_ewald_intersect.jpg)
The atomic scattering factors, fj, accounts for the reduction in
diffraction intensity due to Compton scattering. Compute saed uses
@ -129,7 +129,7 @@ type. Valid chemical symbols for compute saed are:
Tl: Pb: Bi: Po: At:
Rn: Fr: Ra: Ac: Th:
Pa: U: Np: Pu: Am:
Cm: Bk: Cf:
Cm: Bk: Cf:tb(c=5,s=:)
If the {echo} keyword is specified, compute saed will provide extra
@ -151,11 +151,7 @@ All array values calculated by this compute are "intensive".
[Restrictions:]
This command is part of the USER-DIFFRACTION package. It is only
enabled if LAMMPS was built with that package. See the "Making
LAMMPS"_Section_start.html#start_3 section for more info.
This command does not work for triclinic cells.
The compute_saed command does not work for triclinic cells.
[Related commands:]

View File

@ -46,8 +46,8 @@
<PRE>compute 1 all xrd 1.541838 Al O 2Theta 0.087 0.87 c 1 1 1 LP 1 echo
compute 2 all xrd 1.541838 Al O 2Theta 10 100 c 0.05 0.05 0.05 LP 1 manual
</PRE>
<PRE>fix 1 all ave/histo 1 1 1 0.087 0.87 250 c_1[1] mode vector weights c_1[2] file Rad2Theta.xrd
fix 2 all ave/histo 1 1 1 10 100 250 c_2[1] mode vector weights c_2[2] file Deg2Theta.xrd
<PRE>fix 1 all ave/histo/weights 1 1 1 0.087 0.87 250 c_1[1] c_1[2] mode vector file Rad2Theta.xrd
fix 2 all ave/histo/weights 1 1 1 10 100 250 c_2[1] c_2[2] mode vector file Deg2Theta.xrd
</PRE>
<PRE>
</PRE>
@ -55,11 +55,11 @@ fix 2 all ave/histo 1 1 1 10 100 250 c_2[1] mode vector weights c_2[2] file Deg2
</P>
<P>Define a computation that calculates x-ray diffraction intensity as described
in <A HREF = "#Coleman">(Coleman)</A> on a mesh of reciprocal lattice nodes defined
by the entire simulation domain (or manually) using simulated radiation
by the entire simulation domain (or manually) using a simulated radiation
of wavelength lambda.
</P>
<P>The x-ray diffraction intensity I at each reciprocal lattice point
is computed from the structure factor F using the equations:
<P>The x-ray diffraction intensity, I, at each reciprocal lattice point, k,
is computed from the structure factor, F, using the equations:
</P>
<CENTER><IMG SRC = "Eqs/compute_xrd1.jpg">
</CENTER>
@ -73,14 +73,14 @@ is computed from the structure factor F using the equations:
position of each atom, fj are atomic scattering factors, LP is the
Lorentz-polarization factor, and theta is the scattering angle of
diffraction. The Lorentz-polarization factor can be turned off using
the <I>LP</I> switch.
the optional <I>LP</I> keyword.
</P>
<P>Diffraction intensities are calculated on a three-dimensional mesh of
reciprocal lattice nodes. The mesh spacing is defined either (I)
by the entire simulation domain or (II) manually using selected values as
reciprocal lattice nodes. The mesh spacing is defined either (a)
by the entire simulation domain or (b) manually using selected values as
shown in the 2D diagram below.
</P>
<CENTER><A HREF = "xrd_mesh.jpg"><IMG SRC = "JPG/xrd_mesh_small.jpg"></A>
<CENTER><A HREF = "JPG/xrd_mesh.jpg"><IMG SRC = "JPG/xrd_mesh_small.jpg"></A>
</CENTER>
<P>For a mesh defined by the simulation domain, a rectilinear grid is
constructed with spacing <I>c</I>*inv(A) along each reciprocal lattice
@ -159,7 +159,7 @@ type. Valid chemical symbols for compute xrd are:
Ac3+: Th: Th4+: Pa: U:
U3+: U4+: U6+: Np: Np3+:
Np4+: Np6+: Pu: Pu3+: Pu4+:
Pu6+: Am: Cm: Bk: Cf:
Pu6+: Am: Cm: Bk: Cf:tb(c=5,s=:)
</P>
<P>If the <I>echo</I> keyword is specified, compute xrd will provide extra
reporting information to the screen.
@ -182,15 +182,11 @@ options.
</P>
<P><B>Restrictions:</B>
</P>
<P>This command is part of the USER-DIFFRACTION package. It is only
enabled if LAMMPS was built with that package. See the <A HREF = "Section_start.html#start_3">Making
LAMMPS</A> section for more info.
</P>
<P>This command does not work for triclinic cells.
<P>The compute_xrd command does not work for triclinic cells.
</P>
<P><B>Related commands:</B>
</P>
<P><A HREF = "compute_ave_histo.html">fix ave/histo</A>,
<P><A HREF = "fix_ave_histo.html">fix ave/histo</A>,
<A HREF = "compute_saed.html">compute saed</A>
</P>
<P><B>Default:</B>

View File

@ -37,8 +37,8 @@ keyword = {2Theta} or {c} or {LP} or {manual} or {echo} :l
compute 1 all xrd 1.541838 Al O 2Theta 0.087 0.87 c 1 1 1 LP 1 echo
compute 2 all xrd 1.541838 Al O 2Theta 10 100 c 0.05 0.05 0.05 LP 1 manual :pre
fix 1 all ave/histo 1 1 1 0.087 0.87 250 c_1\[1\] mode vector weights c_1\[2\] file Rad2Theta.xrd
fix 2 all ave/histo 1 1 1 10 100 250 c_2\[1\] mode vector weights c_2\[2\] file Deg2Theta.xrd :pre
fix 1 all ave/histo/weights 1 1 1 0.087 0.87 250 c_1\[1\] c_1\[2\] mode vector file Rad2Theta.xrd
fix 2 all ave/histo/weights 1 1 1 10 100 250 c_2\[1\] c_2\[2\] mode vector file Deg2Theta.xrd :pre
:pre
@ -46,11 +46,11 @@ fix 2 all ave/histo 1 1 1 10 100 250 c_2\[1\] mode vector weights c_2\[2\] file
Define a computation that calculates x-ray diffraction intensity as described
in "(Coleman)"_#Coleman on a mesh of reciprocal lattice nodes defined
by the entire simulation domain (or manually) using simulated radiation
by the entire simulation domain (or manually) using a simulated radiation
of wavelength lambda.
The x-ray diffraction intensity I at each reciprocal lattice point
is computed from the structure factor F using the equations:
The x-ray diffraction intensity, I, at each reciprocal lattice point, k,
is computed from the structure factor, F, using the equations:
:c,image(Eqs/compute_xrd1.jpg)
:c,image(Eqs/compute_xrd2.jpg)
@ -61,14 +61,14 @@ Here, K is the location of the reciprocal lattice node, rj is the
position of each atom, fj are atomic scattering factors, LP is the
Lorentz-polarization factor, and theta is the scattering angle of
diffraction. The Lorentz-polarization factor can be turned off using
the {LP} switch.
the optional {LP} keyword.
Diffraction intensities are calculated on a three-dimensional mesh of
reciprocal lattice nodes. The mesh spacing is defined either (I)
by the entire simulation domain or (II) manually using selected values as
reciprocal lattice nodes. The mesh spacing is defined either (a)
by the entire simulation domain or (b) manually using selected values as
shown in the 2D diagram below.
:c,image(JPG/xrd_mesh_small.jpg,xrd_mesh.jpg)
:c,image(JPG/xrd_mesh_small.jpg,JPG/xrd_mesh.jpg)
For a mesh defined by the simulation domain, a rectilinear grid is
@ -148,8 +148,7 @@ type. Valid chemical symbols for compute xrd are:
Ac3+: Th: Th4+: Pa: U:
U3+: U4+: U6+: Np: Np3+:
Np4+: Np6+: Pu: Pu3+: Pu4+:
Pu6+: Am: Cm: Bk: Cf:
Pu6+: Am: Cm: Bk: Cf:tb(c=5,s=:)
If the {echo} keyword is specified, compute xrd will provide extra
reporting information to the screen.
@ -172,15 +171,11 @@ All array values calculated by this compute are "intensive".
[Restrictions:]
This command is part of the USER-DIFFRACTION package. It is only
enabled if LAMMPS was built with that package. See the "Making
LAMMPS"_Section_start.html#start_3 section for more info.
This command does not work for triclinic cells.
The compute_xrd command does not work for triclinic cells.
[Related commands:]
"fix ave/histo"_compute_ave_histo.html,
"fix ave/histo"_fix_ave_histo.html,
"compute saed"_compute_saed.html
[Default:]

View File

@ -76,7 +76,7 @@
<PRE>fix 1 all ave/histo 100 5 1000 0.5 1.5 50 c_myTemp file temp.histo ave running
fix 1 all ave/histo 100 5 1000 -5 5 100 c_thermo_press[2] c_thermo_press[3] title1 "My output values"
fix 1 all ave/histo 1 100 1000 -2.0 2.0 18 vx vy vz mode vector ave running beyond extra
fix 1 all ave/histo/weight 1 1 1 10 100 2000 c_XRD<B>1</B> c_XRD<B>2</B>
fix 1 all ave/histo/weight 1 1 1 10 100 2000 c_XRD[1] c_XRD[2]
</PRE>
<P><B>Description:</B>
</P>

View File

@ -60,7 +60,7 @@ keyword = {mode} or {file} or {ave} or {start} or {beyond} or {overwrite} or {ti
fix 1 all ave/histo 100 5 1000 0.5 1.5 50 c_myTemp file temp.histo ave running
fix 1 all ave/histo 100 5 1000 -5 5 100 c_thermo_press\[2\] c_thermo_press\[3\] title1 "My output values"
fix 1 all ave/histo 1 100 1000 -2.0 2.0 18 vx vy vz mode vector ave running beyond extra
fix 1 all ave/histo/weight 1 1 1 10 100 2000 c_XRD[1] c_XRD[2] :pre
fix 1 all ave/histo/weight 1 1 1 10 100 2000 c_XRD\[1\] c_XRD\[2\] :pre
[Description:]

View File

@ -17,7 +17,7 @@
</PRE>
<UL><LI>ID, group-ID are documented in <A HREF = "fix.html">fix</A> command
<LI>ave/time/saed = style name of this fix command
<LI>saed/vtk = style name of this fix command
<LI>Nevery = use input values every this many timesteps
@ -62,7 +62,7 @@ outside the <I>Kmax</I> range assigned in the compute saed. The ghost data is
assigned a value of -1 and can be removed setting a minimum isovolume
of 0 within the vizualiziton software. SAED images can be created by
visualizing a spherical slice of the data that is centered at
R_Ewald*<B>h k l</B>/norm(<B>h k l</B>), where R_Ewald=1/lambda.
R_Ewald*[h k l]/norm([h k l]), where R_Ewald=1/lambda.
</P>
<P>The group specified within this command is ignored. However, note that
specified values may represent calculations performed by saed computes
@ -102,16 +102,24 @@ diffraction intensity outputs.
</P>
<P>By default the header contains the following information (with example data):
</P>
<P># vtk DataFile Version 3.0 c_SAED
<PRE># vtk DataFile Version 3.0 c_SAED
Image data set
ASCII
DATASET STRUCTURED_POINTS
DIMENSIONS 337 219 209
ASPECT_RATIO 0.00507953 0.00785161 0.00821458
ORIGIN -0.853361 -0.855826 -0.854316 \n POINT_DATA 15424827
ORIGIN -0.853361 -0.855826 -0.854316
POINT_DATA 15424827
SCALARS intensity float
LOOKUP_TABLE default
...data
...data
</PRE>
<P>In this example, kspace is sampled across a 337 x 219 x 209 point mesh
where the mesh spacing is approximately 0.005, 0.007, and 0.008
inv(length) units in the k1, k2, and k3 directions, respectively.
The data is shifted by -0.85, -0.85, -0.85 inv(length) units so that
the origin will lie at 0, 0, 0. Here, 15,424,827 kspace points are
sampled in total.
</P>
<HR>
@ -170,12 +178,8 @@ minimization</A>.
</P>
<P><B>Restrictions:</B>
</P>
<P>This command is part of the USER-DIFFRACTION package. It is only
enabled if LAMMPS was built with that package. See the <A HREF = "Section_start.html#start_3">Making
LAMMPS</A> section for more info.
</P>
<P>The attributes for fix_saed_vtk must match the values assigned in the
associated <A HREF = "compute_saed.txt">compute_saed</A> command.
associated <A HREF = "compute_saed.html">compute_saed</A> command.
</P>
<P><B>Related commands:</B>
</P>

View File

@ -13,7 +13,7 @@ fix saed/vtk command :h3
fix ID group-ID saed/vtk Nevery Nrepeat Nfreak c_ID attribute args ... keyword args ... :pre
ID, group-ID are documented in "fix"_fix.html command :ulb,l
ave/time/saed = style name of this fix command :l
saed/vtk = style name of this fix command :l
Nevery = use input values every this many timesteps :l
Nrepeat = # of times to use input values for calculating averages :l
Nfreq = calculate averages every this many timesteps :l
@ -54,7 +54,7 @@ outside the {Kmax} range assigned in the compute saed. The ghost data is
assigned a value of -1 and can be removed setting a minimum isovolume
of 0 within the vizualiziton software. SAED images can be created by
visualizing a spherical slice of the data that is centered at
R_Ewald*[h k l]/norm([h k l]), where R_Ewald=1/lambda.
R_Ewald*\[h k l\]/norm(\[h k l\]), where R_Ewald=1/lambda.
The group specified within this command is ignored. However, note that
specified values may represent calculations performed by saed computes
@ -100,10 +100,18 @@ ASCII
DATASET STRUCTURED_POINTS
DIMENSIONS 337 219 209
ASPECT_RATIO 0.00507953 0.00785161 0.00821458
ORIGIN -0.853361 -0.855826 -0.854316 \n POINT_DATA 15424827
ORIGIN -0.853361 -0.855826 -0.854316
POINT_DATA 15424827
SCALARS intensity float
LOOKUP_TABLE default
...data
...data :pre
In this example, kspace is sampled across a 337 x 219 x 209 point mesh
where the mesh spacing is approximately 0.005, 0.007, and 0.008
inv(length) units in the k1, k2, and k3 directions, respectively.
The data is shifted by -0.85, -0.85, -0.85 inv(length) units so that
the origin will lie at 0, 0, 0. Here, 15,424,827 kspace points are
sampled in total.
:line
@ -162,12 +170,8 @@ minimization"_minimize.html.
[Restrictions:]
This command is part of the USER-DIFFRACTION package. It is only
enabled if LAMMPS was built with that package. See the "Making
LAMMPS"_Section_start.html#start_3 section for more info.
The attributes for fix_saed_vtk must match the values assigned in the
associated "compute_saed"_compute_saed.txt command.
associated "compute_saed"_compute_saed.html command.
[Related commands:]