Extra:Intersecting planes

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 $\scriptstyle xy$ plane, and thus has normal vector:

$\mathbf {n} _{1}={\begin{bmatrix}0&0&1\end{bmatrix}}$

The first plane intersects the origin. The second plane also intersects the origin, but is tilted about the $\scriptstyle x$-axis by $\scriptstyle \alpha$, such that its normal is:

$\mathbf {n} _{2}={\begin{bmatrix}0&-\sin \alpha \ &\cos \alpha \end{bmatrix}}$

We are interested in quantities that are a particular distance ($\scriptstyle q$) from the origin. Imagine a vector of length $\scriptstyle q$ lying in plane 1, rotated about the Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle \scriptstyle z} axis by $\scriptstyle \phi$ (i.e. the angular distance from the Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle \scriptstyle y} -axis is $\scriptstyle \phi$):

Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle \mathbf {v} _{1}={\begin{bmatrix}q\sin \phi &q\cos \phi \ &0\end{bmatrix}}}

The second vector (Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle \scriptstyle \mathbf {v} _{2}} ) is lying in plane 2. We call Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle \scriptstyle \alpha _{r}} the angle between Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle \scriptstyle \mathbf {v} _{1}} and Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle \scriptstyle \mathbf {v} _{2}} . The specified geometry uniquely defines Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle \scriptstyle \alpha _{r}} in terms of the angle between the planes ($\scriptstyle \alpha$) and the amount of rotation of the vectors ($\scriptstyle \phi$) within their respective planes. In particular, Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle \scriptstyle \mathbf {v} _{2}} can be thought of as Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle \scriptstyle \mathbf {v} _{2b}} rotated about $\scriptstyle \mathbf {n} _{2}$ by $\scriptstyle \phi$, where Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle \scriptstyle \mathbf {v} _{2}b} is the vector in plane 2 without any $\scriptstyle \phi$ rotation (i.e. lying in the Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle \scriptstyle yz} plane):

Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle \mathbf {v} _{2b}={\begin{bmatrix}0&q\cos \alpha \ &q\sin \alpha \end{bmatrix}}}

In general, rotation of a vector Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle \scriptstyle \mathbf {v} _{\mathrm {start} }={\begin{bmatrix}x&y&z\end{bmatrix}}} about an arbitrary unit-vector Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle \scriptstyle \mathbf {n} ={\begin{bmatrix}u&v&w\end{bmatrix}}} gives (1, 2):

$\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}}$

In this particular case, we thus expect:

${\begin{alignedat}{2}\mathbf {v} _{2}&={\begin{bmatrix}(-wy+vz)\sin \phi \\v(vy+wz)(1-\cos \phi )+y\cos \phi \\w(vy+wz)(1-\cos \phi )+z\cos \phi \end{bmatrix}}\\&={\begin{bmatrix}(-\cos \alpha q\cos \alpha +-\sin \alpha q\sin \alpha )\sin \phi \\-\sin \alpha (-\sin \alpha q\cos \alpha +\cos \alpha q\sin \alpha )(1-\cos \phi )+q\cos \alpha \cos \phi \\\cos \alpha (-\sin \alpha q\cos \alpha +\cos \alpha q\sin \alpha )(1-\cos \phi )+q\sin \alpha \cos \phi \end{bmatrix}}\\&=q{\begin{bmatrix}-(\cos ^{2}\alpha +\sin ^{2}\alpha )\sin \phi \\\sin ^{2}\alpha (\cos \alpha -\cos \alpha )(1-\cos \phi )+\cos \alpha \cos \phi \\\cos ^{2}\alpha (-\sin \alpha +\sin \alpha )(1-\cos \phi )+\sin \alpha \cos \phi \end{bmatrix}}\\&=q{\begin{bmatrix}-\sin \phi \\\cos \alpha \cos \phi \\\sin \alpha \cos \phi \end{bmatrix}}\end{alignedat}}$

Note that we replace $\scriptstyle \phi$ by $\scriptstyle -\phi$ to force the same orientation convention in the definition of rotating Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle \scriptstyle \mathbf {v} _{1}} and Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle \scriptstyle \mathbf {v} _{2}} :

${\begin{alignedat}{2}\mathbf {v} _{2}&=q{\begin{bmatrix}\sin \phi \\\cos \alpha \cos \phi \\\sin \alpha \cos \phi \end{bmatrix}}\end{alignedat}}$

The angle between Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle \scriptstyle \mathbf {v} _{1}} and $\scriptstyle \mathbf {v} _{2}$ is Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle \scriptstyle \alpha _{r}} :

${\begin{alignedat}{2}\cos \alpha _{r}&={\frac {\mathbf {v} _{1}\cdot \mathbf {v} _{2}}{|\mathbf {v} _{1}||\mathbf {v} _{2}|}}\\&={\frac {(q^{2}\sin \phi \sin \phi )+(q^{2}\cos \phi \cos \alpha \cos \phi )+(0)}{(q)(q)}}\\&=\sin ^{2}\phi +\cos ^{2}\phi \cos \alpha \\\alpha _{r}&=\cos ^{-1}\left[\cos ^{2}\phi \cos \alpha +\sin ^{2}\phi \right]\end{alignedat}}$

Distance between two planes

The distance between Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle \scriptstyle \mathbf {v} _{1}} and $\scriptstyle \mathbf {v} _{2}$ is $\scriptstyle d$:

${\begin{alignedat}{2}\sin \left({\frac {\alpha _{r}}{2}}\right)&={\frac {d/2}{q}}\\d&=2q\sin \left({\frac {\alpha _{r}}{2}}\right)\\&=2q\sin \left({\frac {1}{2}}\cos ^{-1}\left[\cos ^{2}\phi \cos \alpha +\sin ^{2}\phi \right]\right)\end{alignedat}}$

Alternatively:

${\begin{alignedat}{2}d^{2}&=2q^{2}-2q^{2}\cos \alpha _{r}\\d&=q{\sqrt {2\left(1-\cos \alpha _{r}\right)}}\\d&=q{\sqrt {2\left(1-\cos \left(\cos ^{-1}\left[\cos ^{2}\phi \cos \alpha +\sin ^{2}\phi \right]\right)\right)}}\\d&=q{\sqrt {2\left(1-\cos ^{2}\phi \cos \alpha +\sin ^{2}\phi \right)}}\end{alignedat}}$

If the two vectors do not have equal length:

${\begin{alignedat}{2}d^{2}&=q_{1}^{2}+q_{2}^{2}-2q_{1}q_{2}\cos \alpha _{r}\\d&={\sqrt {q_{1}^{2}+q_{2}^{2}-2q_{1}q_{2}\left(\cos ^{2}\phi \cos \alpha +\sin ^{2}\phi \right)}}\\\end{alignedat}}$

Generalized intersection

Imagine reciprocal-space scattering that is a ring; more specifically a pseudo-toroid with Gaussian-like decay. The intensity overall is:

$I=\exp \left[-(q_{rr}-q_{0})^{2}/(2\sigma _{q}^{2})\right]\exp \left[-q_{rz}^{2}/(2\sigma _{q}^{2})\right]$

Where we use the subscript r to denote the reciprocal-space coordinate system, and $\scriptstyle q_{rr}={\sqrt {q_{rx}^{2}+q_{ry}^{2}}}$. The plane of the detector (i.e. the Ewald plane) is denoted by d:

$\mathbf {v} _{d}={\begin{bmatrix}q_{dx}&q_{dy}&0\end{bmatrix}}$

We set the symmetry axis in realspace (detector coordinate system) to be the $\scriptstyle x$-axis. The reciprocal-space is tilted by Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \scriptstyle \chi_0} (about the Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \scriptstyle y} -axis), before the 'powder' rotation about the $\scriptstyle x$-axis (where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \scriptstyle \alpha} goes from Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \scriptstyle -\pi} to Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \scriptstyle +\pi} ). Consider an initial vector:

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{v}_{d0} = \begin{bmatrix} q \sin \phi & q \cos \phi & 0 \end{bmatrix} }

The 1st rotation (about Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \scriptstyle y} -axis by Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \scriptstyle \chi_0} ) involves:

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{alignat}{2} \mathbf{v}_{d1} & = R_y(\chi_0) \mathbf{v}_{d0} \\ & = \begin{bmatrix} \cos \chi_0 & 0 & \sin \chi_0 \\ 0 & 1 & 0 \\ -\sin \chi_0 & 0 & \cos \chi_0 \\ \end{bmatrix} \begin{bmatrix} q \sin \phi \\ q \cos \phi \\ 0 \\ \end{bmatrix} \\ & = \begin{bmatrix} q \sin \phi \cos \chi_0 \\ q \cos \phi \\ - q \sin \phi \sin \chi_0 \\ \end{bmatrix} \end{alignat} }

The 2nd rotation (about the Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \scriptstyle q_{dx}} -axis by Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \scriptstyle \alpha} ) occurs with respect to the vector:

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \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} }

The second rotation can again be thought of in general terms as Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \scriptstyle \mathbf{v}_{\mathrm{start}} = \begin{bmatrix} x & y & z \end{bmatrix}} rotated by Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \scriptstyle \theta} about an unit-vector Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \scriptstyle \mathbf{u} = \begin{bmatrix} u & v & w \end{bmatrix}} .

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \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} }

In this case:

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \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} }

This reduces the problem to:

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \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} }

More specifically:

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \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} }

This vector describes the possible positions of the intersecting detector-plane, in the sample's reciprocal-space.