|
|
(21 intermediate revisions by the same user not shown) |
Line 97: |
Line 97: |
| | | |
| | | |
− | ==Generalized distance between two vectors== | + | ==Generalized intersection== |
| Imagine [[reciprocal-space]] [[scattering]] that is a ring; more specifically a pseudo-toroid with Gaussian-like decay. The intensity overall is: | | Imagine [[reciprocal-space]] [[scattering]] that is a ring; more specifically a pseudo-toroid with Gaussian-like decay. The intensity overall is: |
| ::<math> | | ::<math> |
− | I = \exp \left [ -(q_{rx}^2 + q_{ry}^2)/(2 \sigma_q^2) \right ] \exp \left [ -q_{rz}^2/(2 \sigma_q^2) \right ] | + | I = \exp \left [ -(q_{rr}-q_0)^2/(2 \sigma_q^2) \right ] \exp \left [ -q_{rz}^2/(2 \sigma_q^2) \right ] |
| </math> | | </math> |
− | Where we use the subscript ''r'' to denote the reciprocal-space coordinate system. The plane of the [[detector]] (i.e. the [[Ewald sphere|Ewald plane]]) is denoted by ''d'': | + | Where we use the subscript ''r'' to denote the reciprocal-space coordinate system, and <math>\scriptstyle q_{rr} = \sqrt{q_{rx}^2+q_{ry}^2}</math>. The plane of the [[detector]] (i.e. the [[Ewald sphere|Ewald plane]]) is denoted by ''d'': |
| ::<math> | | ::<math> |
| \mathbf{v}_{d} = \begin{bmatrix} q_{dx} & q_{dy} & 0 \end{bmatrix} | | \mathbf{v}_{d} = \begin{bmatrix} q_{dx} & q_{dy} & 0 \end{bmatrix} |
Line 131: |
Line 131: |
| \end{alignat} | | \end{alignat} |
| </math> | | </math> |
| + | The 2nd rotation (about the <math>\scriptstyle q_{dx}</math>-axis by <math>\scriptstyle \alpha</math>) occurs with respect to the vector: |
| + | ::<math> |
| + | \begin{alignat}{2} |
| + | \mathbf{u}_{d1} & = R_y(\chi_0) \mathbf{u}_{d0} \\ |
| + | & = \begin{bmatrix} |
| + | \cos \chi_0 & 0 & \sin \chi_0 \\ |
| + | 0 & 1 & 0 \\ |
| + | -\sin \chi_0 & 0 & \cos \chi_0 \\ |
| + | \end{bmatrix} |
| + | \begin{bmatrix} |
| + | 1 \\ |
| + | 0 \\ |
| + | 0 \\ |
| + | \end{bmatrix} \\ |
| + | & = \begin{bmatrix} |
| + | \cos \chi_0 \\ |
| + | 0 \\ |
| + | - \sin \chi_0 \\ |
| + | \end{bmatrix} |
| + | \end{alignat} |
| + | </math> |
| + | |
| + | The second rotation can again be thought of in general terms as <math>\scriptstyle \mathbf{v}_{\mathrm{start}} = \begin{bmatrix} x & y & z \end{bmatrix}</math> rotated by <math>\scriptstyle \theta</math> about an unit-vector <math>\scriptstyle \mathbf{u} = \begin{bmatrix} u & v & w \end{bmatrix}</math>. |
| + | ::<math> |
| + | \mathbf{v}_{\mathrm{end}} = \begin{bmatrix} u(ux +vy+wz)(1-\cos \theta) + x \cos \theta + (-wy+vz)\sin \theta \\ |
| + | v(ux+vy+wz)(1 - \cos \theta) + y \cos \theta + (wx-uz) \sin \theta\\ |
| + | w(ux+vy+wz)(1-\cos \theta) + z \cos \theta + (-vx+uy) \sin \theta \end{bmatrix} |
| + | </math> |
| + | In this case: |
| + | ::<math> |
| + | \begin{alignat}{2} |
| + | \mathbf{v}_{\mathrm{start}} & |
| + | = \begin{bmatrix} x \\ y \\ z \end{bmatrix} |
| + | = \mathbf{v}_{d1} |
| + | = \begin{bmatrix} q \sin \phi \cos \chi_0 \\ q \cos \phi \\ -q \sin \phi \sin \chi_0 \end{bmatrix} \\ |
| + | \mathbf{u} & |
| + | = \begin{bmatrix} u \\ v \\ w \end{bmatrix} |
| + | = \mathbf{u}_{d1} |
| + | = \begin{bmatrix} \cos \chi_0 \\ 0 \\ -\sin \chi_0 \end{bmatrix} |
| + | \end{alignat} |
| + | </math> |
| + | This reduces the problem to: |
| + | ::<math> |
| + | \mathbf{v}_{\mathrm{end}} = |
| + | \begin{bmatrix} u(ux +wz)(1-\cos \alpha) + x \cos \alpha + -wy \sin \alpha \\ |
| + | y \cos \alpha + (wx-uz) \sin \alpha\\ |
| + | w(ux+wz)(1-\cos \alpha) + z \cos \alpha + uy \sin \alpha \end{bmatrix} |
| + | </math> |
| + | |
| + | More specifically: |
| + | ::<math> |
| + | \begin{alignat}{2} |
| + | \mathbf{v}_{d2} & = |
| + | \begin{bmatrix} |
| + | \cos \chi_0 (\cos \chi_0 q \sin \phi \cos \chi_0 +\sin \chi_0 q \sin \phi \sin \chi_0)(1-\cos \alpha) + q \sin \phi \cos \chi_0 \cos \alpha +\sin \chi_0 q \cos \phi \sin \alpha \\ |
| + | q \cos \phi \cos \alpha + (-\sin \chi_0 q \sin \phi \cos \chi_0 + \cos \chi_0 q \sin \phi \sin \chi_0 ) \sin \alpha\\ |
| + | -\sin \chi_0 ( \cos \chi_0 q \sin \phi \cos \chi_0 +\sin \chi_0 q \sin \phi \sin \chi_0 )(1-\cos \alpha) -q \sin \phi \sin \chi_0 \cos \alpha + \cos \chi_0 q \cos \phi \sin \alpha \end{bmatrix} \\ |
| + | & = |
| + | q \begin{bmatrix} |
| + | \cos \chi_0 \sin \phi(\cos^2 \chi_0 +\sin^2 \chi_0 )(1-\cos \alpha) + \sin \phi \cos \chi_0 \cos \alpha +\sin \chi_0 \cos \phi \sin \alpha \\ |
| + | \cos \phi \cos \alpha + \sin \phi (-\sin \chi_0 \cos \chi_0 + \cos \chi_0 \sin \chi_0 ) \sin \alpha\\ |
| + | -\sin \chi_0 \sin \phi ( \cos^2 \chi_0 +\sin^2 \chi_0 )(1-\cos \alpha) - \sin \phi \sin \chi_0 \cos \alpha + \cos \chi_0 \cos \phi \sin \alpha \end{bmatrix} \\ |
| + | & = |
| + | q \begin{bmatrix} |
| + | \cos \chi_0 \sin \phi(1-\cos \alpha) + \sin \phi \cos \chi_0 \cos \alpha +\sin \chi_0 \cos \phi \sin \alpha \\ |
| + | \cos \phi \cos \alpha\\ |
| + | -\sin \chi_0 \sin \phi (1-\cos \alpha) - \sin \phi \sin \chi_0 \cos \alpha + \cos \chi_0 \cos \phi \sin \alpha \end{bmatrix} \\ |
| + | |
| + | \end{alignat} |
| + | </math> |
| + | This vector describes the possible positions of the intersecting detector-plane, in the sample's reciprocal-space. |
Latest revision as of 17:16, 23 December 2015
A common problem in scattering is to consider the intersection of various planes (representing the Ewald sphere, reciprocal space, etc.).
Angle between two planes
The general case for the angle between two planes is well known. Consider a particular case where we want to know how the angle between two planes depends on the direction/orientation of a third plane/vector that intersects the first two. I.e. what is the minimal angle between two planes along a 'certain direction' (what is the angle between two vectors that both lie on the third plane, and which lie on planes 1 and 2, respectively).
One of the planes represents reciprocal-space scattering (e.g. mostly localized to a plane); the other represents the detector. We are interested in the angle between them so that we can calculate the distance between them, so that we can compute 'how much' scattering is seen on the detector. To make this concrete, plane 1 lies in the plane, and thus has normal vector:
The first plane intersects the origin. The second plane also intersects the origin, but is tilted about the -axis by , such that its normal is:
We are interested in quantities that are a particular distance () from the origin. Imagine a vector of length lying in plane 1, rotated about the axis by (i.e. the angular distance from the -axis is ):
The second vector () is lying in plane 2. We call the angle between and . The specified geometry uniquely defines in terms of the angle between the planes () and the amount of rotation of the vectors () within their respective planes. In particular, can be thought of as rotated about by , where is the vector in plane 2 without any rotation (i.e. lying in the plane):
In general, rotation of a vector about an arbitrary unit-vector gives (1, 2):
In this particular case, we thus expect:
Note that we replace by to force the same orientation convention in the definition of rotating and :
The angle between and is :
Distance between two planes
The distance between and is :
Alternatively:
If the two vectors do not have equal length:
Generalized intersection
Imagine reciprocal-space scattering that is a ring; more specifically a pseudo-toroid with Gaussian-like decay. The intensity overall is:
Where we use the subscript r to denote the reciprocal-space coordinate system, and . The plane of the detector (i.e. the Ewald plane) is denoted by d:
We set the symmetry axis in realspace (detector coordinate system) to be the -axis. The reciprocal-space is tilted by (about the -axis), before the 'powder' rotation about the -axis (where goes from to ). Consider an initial vector:
The 1st rotation (about -axis by ) involves:
The 2nd rotation (about the -axis by ) occurs with respect to the vector:
The second rotation can again be thought of in general terms as rotated by about an unit-vector .
In this case:
This reduces the problem to:
More specifically:
This vector describes the possible positions of the intersecting detector-plane, in the sample's reciprocal-space.