Difference between revisions of "Form Factor:Pyramid"

From GISAXS
Jump to: navigation, search
(Created page with "==Equations== For pyramid of base edge-length 2''R'', and height ''H''. The angle of the pyramid walls is <math>\alpha</math>. If <math>H < R/ \tan\alpha</math> then the pyram...")
 
Line 430: Line 430:
 
\right|^2 \sin\theta\mathrm{d}\theta\mathrm{d}\phi \\
 
\right|^2 \sin\theta\mathrm{d}\theta\mathrm{d}\phi \\
  
 +
\end{alignat}
 +
</math>
 +
 +
==Regular Pyramid==
 +
A regular pyramid (half of an octahedron) has faces that are equilateral triangles (each vertex is 60°). The 'corner-to-edge' distance along each triangular face is then:
 +
:<math> d_{face,c-e} = R \tan(60^{\circ}) = \sqrt{3} R</math>
 +
This makes the height:
 +
:<math>
 +
\begin{alignat}{2}
 +
(d_{face,c-e})^2 & = (H)^2 + (R)^2 \\
 +
H^2  & = (d_{face,c-e})^2 - (R)^2\\
 +
H  & = \sqrt{ (\sqrt{3} R)^2  - (R)^2 }\\
 +
  & = \sqrt{ 3 R^2  - R^2 }\\
 +
  & = \sqrt{ 2 } R \\
 +
\end{alignat}
 +
</math>
 +
 +
So that the pyramid face angle, <math>\alpha</math> is:
 +
:<math>
 +
\begin{alignat}{2}
 +
\tan(\alpha) & = \frac{ H }{ R } \\
 +
  \alpha & = \arctan \left( \frac{\sqrt{ 2 } R}{R} \right) \\
 +
  & = \arctan( \sqrt{2} ) \\
 +
  & \approx 0.9553 \\
 +
  & \approx 54.75^{\circ}
 +
\end{alignat}
 +
</math>
 +
 +
The square base of the pyramid has edges of length 2''R''. The distance from the center of the square to any corner is ''H'', such that:
 +
:<math>
 +
\begin{alignat}{2}
 +
\cos(45^{\circ}) & = \frac{R}{H} \\
 +
H & = \frac{R}{ 1/\sqrt{2} } \\
 +
  & = \sqrt{2} R
 +
\end{alignat}
 +
</math>
 +
 +
 +
===Surface Area===
 +
For a non-truncated, regular pyramid, each face is an equilateral triangle (each vertex is 60°). So each face:
 +
:<math>
 +
\begin{alignat}{2}
 +
S_{face}
 +
  & = 2 \times \left( \frac{ R R \tan(60^{\circ}) }{2} \right) \\
 +
  & = R^2 \sqrt{3}
 +
\end{alignat}
 +
</math>
 +
The base is simply:
 +
:<math>
 +
\begin{alignat}{2}
 +
S_{base}
 +
  & = 2 R \times 2 R \\
 +
  & = 4 R^2
 +
\end{alignat}
 +
</math>
 +
Total:
 +
 +
:<math>
 +
\begin{alignat}{2}
 +
S_{pyr}
 +
  & = 4 \times R^2 \sqrt{3} + 4 R^2 \\
 +
  & = 4(1 + \sqrt{3}) R^2
 +
\end{alignat}
 +
</math>
 +
===Volume===
 +
For a regular pyramid, the height <math>H=\sqrt{2}R</math> and <math>\tan(\alpha)=H/R = \sqrt{2}</math>:
 +
:<math>
 +
\begin{alignat}{2}
 +
V_{pyr}
 +
  & = \frac{4}{3} \tan (\alpha) \left[ R^3 - \left( R - \frac{H}{ \tan (\alpha)} \right)^3 \right] \\
 +
  & = \frac{4}{3} \sqrt{2} \left[ R^3 - \left( R - \frac{ \sqrt{2} R }{ \sqrt{2}} \right)^3 \right] \\
 +
  & = \frac{4\sqrt{2}}{3} R^3 \\
 
\end{alignat}
 
\end{alignat}
 
</math>
 
</math>

Revision as of 17:00, 13 June 2014

Equations

For pyramid of base edge-length 2R, and height H. The angle of the pyramid walls is . If then the pyramid is truncated (flat top).

  • Volume
  • Projected (xy) surface area

Form Factor Amplitude

where

Isotropic Form Factor Intensity

Derivations

Form Factor

For a pyramid of base-edge-length 2R, side-angle , truncated at H (along z axis), we note that the in-plane size of the pyramid at height z is:

Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle R_{z}=R-{\frac {z}{\tan \alpha }}}

Integrating with Cartesian coordinates:

A recurring integral is (c.f. cube form factor):

Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle {\begin{alignedat}{2}f_{x}(q_{x})&=\int _{-R_{z}}^{R_{z}}e^{iq_{x}x}\mathrm {d} x\\&=\int _{-R_{z}}^{R_{z}}\left[\cos(q_{x}x)+i\sin(q_{x}x)\right]\mathrm {d} x\\&=-{\frac {2}{q_{x}}}\sin(q_{x}R_{z})\\&=-2R_{z}\mathrm {sinc} (q_{x}R_{z})\\\end{alignedat}}}

Which gives:

This can be simplified automated solving. For a regular pyramid, we obtain:

Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle {\begin{alignedat}{2}F_{pyr}(\mathbf {q} )&={\frac {4{\sqrt {2}}}{q_{x}q_{y}}}{\frac {\left({\begin{array}{l}-q_{y}\left(-q_{x}^{2}+q_{y}^{2}-2q_{z}^{2}\right)\cos(q_{y}R)\sin(q_{x}R)\\\,\,\,\,-q_{x}\cos(q_{x}R)\left(2i{\sqrt {2}}q_{y}q_{z}\cos(q_{y}R)+\left(q_{x}^{2}-q_{y}^{2}-2q_{z}^{2}\right)\sin(q_{y}R)\right)\\\,\,\,\,+i{\sqrt {2}}q_{z}\left(2e^{i{\sqrt {2}}q_{z}R}q_{x}q_{y}-\left(q_{x}^{2}+q_{y}^{2}-2q_{z}^{2}\right)\sin(q_{x}R)\sin(q_{y}R)\right)\end{array}}\right)}{q_{x}^{4}+(q_{y}^{2}-2q_{z}^{2})^{2}-2q_{x}^{2}(q_{y}^{2}+2q_{z}^{2})}}\end{alignedat}}}

Form Factor near q=0

qy

When :

Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle {\begin{alignedat}{2}q_{1}&=q_{3}\\q_{2}&=q_{4}\\K_{1}&=K_{3}\\K_{2}&=K_{4}\\\end{alignedat}}}

So:

qx

When :

Since sinc is an even function:

Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle {\begin{alignedat}{2}K_{1}&=\,\,+{\text{sinc}}(q_{1}H)e^{+iq_{1}H}+\,\,{\text{sinc}}(q_{2}H)e^{-iq_{2}H}=K_{3}\\K_{2}&=-i{\text{sinc}}(q_{1}H)e^{+iq_{1}H}+i{\text{sinc}}(q_{2}H)e^{-iq_{2}H}=K_{4}\\K_{3}&=\,\,+{\text{sinc}}(q_{2}H)e^{-iq_{2}H}+\,\,{\text{sinc}}(q_{1}H)e^{+iq_{1}H}=K_{1}\\K_{4}&=-i{\text{sinc}}(q_{2}H)e^{-iq_{2}H}+i{\text{sinc}}(q_{1}H)e^{+iq_{1}H}=K_{2}\end{alignedat}}}

And:

qz

When Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle q_{z}=0} :

So:

q

When :

Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle {\begin{alignedat}{2}q_{1}&=q_{2}=q_{3}=q_{4}=0\\\end{alignedat}}}

So:

Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle {\begin{alignedat}{3}K_{1}&=+1+1&=2\\K_{2}&=-i+i&=0\\K_{3}&=+1+1&=2\\K_{4}&=-i+i&=0\end{alignedat}}}

And:

qx and qy

When Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle q_{x}=q_{y}=0} :

Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle {\begin{alignedat}{3}q_{1}&=q_{3}&=+{\frac {q_{z}}{2}}\\q_{2}&=q_{4}&=-{\frac {q_{z}}{2}}\\\end{alignedat}}}

So:

Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle {\begin{alignedat}{2}F_{pyr}(\mathbf {q} )&={\frac {H}{q_{x}q_{y}}}\left({\begin{array}{l}\cos \left[(q_{x}-q_{y})R\right]K_{1}\\\,\,\,\,+\sin \left[(q_{x}-q_{y})R\right]K_{2}\\\,\,\,\,-\cos \left[(q_{x}+q_{y})R\right]K_{1}\\\,\,\,\,-\sin \left[(q_{x}+q_{y})R\right]K_{2}\end{array}}\right)\\\end{alignedat}}}

To analyze the behavior in the limit of small and , we consider the limit of Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle q\to 0} where Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle q_{x}=q_{y}=q} . We replace the trigonometric functions by their expansions near zero (keeping only the first two terms):

Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle {\begin{alignedat}{2}\lim _{q\to 0}F_{pyr}(\mathbf {q} )&={\frac {H}{qq}}\left({\begin{array}{l}\cos \left[(q-q)R\right]K_{1}\\\,\,\,\,+\sin \left[(q-q)R\right]K_{2}\\\,\,\,\,-\cos \left[(q+q)R\right]K_{1}\\\,\,\,\,-\sin \left[(q+q)R\right]K_{2}\end{array}}\right)\\&={\frac {H}{q^{2}}}\left({\begin{array}{l}\left[1-{\frac {((q-q)R)^{2}}{2!}}+\cdots \right]K_{1}\\\,\,\,\,+\left[(q-q)R-{\frac {((q-q)R)^{3}}{3!}}+\cdots \right]K_{2}\\\,\,\,\,-\left[1-{\frac {((q+q)R)^{2}}{2!}}+\cdots \right]K_{1}\\\,\,\,\,-\left[(q+q)R-{\frac {((q-q)R)^{3}}{3!}}+\cdots \right]K_{2}\end{array}}\right)\\&={\frac {H}{q^{2}}}\left({\begin{array}{l}\left[1-{\frac {((q-q)R)^{2}}{2!}}-1+{\frac {((q+q)R)^{2}}{2!}}\right]K_{1}\\\,\,\,\,+\left[(q-q)R-{\frac {((q-q)R)^{3}}{3!}}-(q+q)R+{\frac {((q-q)R)^{3}}{3!}}\right]K_{2}\\\end{array}}\right)\\&={\frac {H}{q^{2}}}\left({\begin{array}{l}\left[{\frac {((2q)R)^{2}}{2!}}-{\frac {((q-q)R)^{2}}{2!}}\right]K_{1}\\\,\,\,\,+\left[(q-q)R-(2q)R\right]K_{2}\\\end{array}}\right)\\&={\frac {(2qR)^{2}}{2!}}{\frac {HK_{1}}{q^{2}}}-{\frac {((q-q)R)^{2}}{2!}}{\frac {HK_{1}}{q^{2}}}+(q-q)R{\frac {HK_{2}}{q^{2}}}-2qR{\frac {HK_{2}}{q^{2}}}\\&={\frac {4R^{2}HK_{1}}{2}}-{\frac {R^{2}HK_{1}}{2}}{\frac {(q-q)^{2}}{q^{2}}}+RHK_{2}{\frac {(q-q)}{q^{2}}}-{\frac {2RHK_{2}}{q}}\\&=2R^{2}HK_{1}\end{alignedat}}}

Note that since is symmetric . When and are small (but not zero and not necessarily equal), many of the above arguments still apply. It remains that Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle K_{2}\approx K_{4}\approx 0} , and:

Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle {\begin{alignedat}{2}\lim _{(q_{x},q_{y})\to 0}F_{pyr}(\mathbf {q} )&={\frac {HK_{1}}{q_{x}q_{y}}}\left(\cos \left[(q_{x}-q_{y})R\right]-\cos \left[(q_{x}+q_{y})R\right]\right)\\&={\frac {HK_{1}}{q_{x}q_{y}}}\left(\left[1-{\frac {((q_{x}-q_{y})R)^{2}}{2!}}+\cdots \right]-\left[1-{\frac {((q_{x}+q_{y})R)^{2}}{2!}}+\cdots \right]\right)\\&={\frac {HK_{1}}{q_{x}q_{y}}}\left({\frac {(q_{x}+q_{y})^{2}R^{2}}{2!}}-{\frac {(q_{x}-q_{y})^{2}R^{2}}{2!}}\right)\\&={\frac {HR^{2}K_{1}}{2q_{x}q_{y}}}\left((q_{x}+q_{y})^{2}-(q_{x}-q_{y})^{2}\right)\\\end{alignedat}}}

Isotropic Form Factor Intensity

To average over all possible orientations, we note:

and use:

Regular Pyramid

A regular pyramid (half of an octahedron) has faces that are equilateral triangles (each vertex is 60°). The 'corner-to-edge' distance along each triangular face is then:

Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle d_{face,c-e}=R\tan(60^{\circ })={\sqrt {3}}R}

This makes the height:

Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle {\begin{alignedat}{2}(d_{face,c-e})^{2}&=(H)^{2}+(R)^{2}\\H^{2}&=(d_{face,c-e})^{2}-(R)^{2}\\H&={\sqrt {({\sqrt {3}}R)^{2}-(R)^{2}}}\\&={\sqrt {3R^{2}-R^{2}}}\\&={\sqrt {2}}R\\\end{alignedat}}}

So that the pyramid face angle, is:

Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle {\begin{alignedat}{2}\tan(\alpha )&={\frac {H}{R}}\\\alpha &=\arctan \left({\frac {{\sqrt {2}}R}{R}}\right)\\&=\arctan({\sqrt {2}})\\&\approx 0.9553\\&\approx 54.75^{\circ }\end{alignedat}}}

The square base of the pyramid has edges of length 2R. The distance from the center of the square to any corner is H, such that:

Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle {\begin{alignedat}{2}\cos(45^{\circ })&={\frac {R}{H}}\\H&={\frac {R}{1/{\sqrt {2}}}}\\&={\sqrt {2}}R\end{alignedat}}}


Surface Area

For a non-truncated, regular pyramid, each face is an equilateral triangle (each vertex is 60°). So each face:

Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle {\begin{alignedat}{2}S_{face}&=2\times \left({\frac {RR\tan(60^{\circ })}{2}}\right)\\&=R^{2}{\sqrt {3}}\end{alignedat}}}

The base is simply:

Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle {\begin{alignedat}{2}S_{base}&=2R\times 2R\\&=4R^{2}\end{alignedat}}}

Total:

Volume

For a regular pyramid, the height and Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle \tan(\alpha )=H/R={\sqrt {2}}} :

Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle {\begin{alignedat}{2}V_{pyr}&={\frac {4}{3}}\tan(\alpha )\left[R^{3}-\left(R-{\frac {H}{\tan(\alpha )}}\right)^{3}\right]\\&={\frac {4}{3}}{\sqrt {2}}\left[R^{3}-\left(R-{\frac {{\sqrt {2}}R}{\sqrt {2}}}\right)^{3}\right]\\&={\frac {4{\sqrt {2}}}{3}}R^{3}\\\end{alignedat}}}