Subsections

Torsionless pseudo-ellipse equations

Figure 16.19: The torsionless pseudo-ellipse model simulated and calculated in-frame Daeg(1) frame order matrix elements. In these plots, θX corresponds to the cone opening half-angle θx and θY to the cone opening half-angle θy. When the half-angle is not varied, the angle is fixed to either θx = π/4 or θy = 3π/8. Frame order matrix values have been calculated every 10 degrees. The first angle for the calculated elements is set to 0.01 degrees as a pseudo-ellipse cone opening angle of 0.0 cannot be correctly handled by the numerical integration.
 \includegraphics[width=.5\textwidth]{images/frame_order_matrix/Sij_pseudo-ellipse_torsionless_in_frame_theta_x_ens1000000.eps} \includegraphics[width=.5\textwidth]{images/frame_order_matrix/Sij_pseudo-ellipse_torsionless_in_frame_theta_x_calc.eps}  
      
 \includegraphics[width=.5\textwidth]{images/frame_order_matrix/Sij_pseudo-ellipse_torsionless_in_frame_theta_y_ens1000000.eps} \includegraphics[width=.5\textwidth]{images/frame_order_matrix/Sij_pseudo-ellipse_torsionless_in_frame_theta_y_calc.eps}  

Figure 16.20: The torsionless pseudo-ellipse model simulated and calculated in-frame Daeg(2) frame order matrix elements. In these plots, θX corresponds to the cone opening half-angle θx and θY to the cone opening half-angle θy. When the half-angle is not varied, the angle is fixed to either θx = π/4 or θy = 3π/8. Frame order matrix values have been calculated every 10 degrees. The first angle for the calculated elements is set to 0.01 degrees as a pseudo-ellipse cone opening angle of 0.0 cannot be correctly handled by the numerical integration.
 \includegraphics[width=.5\textwidth]{images/frame_order_matrix/Sijkl_pseudo-ellipse_torsionless_in_frame_theta_x_ens1000000.eps} \includegraphics[width=.5\textwidth]{images/frame_order_matrix/Sijkl_pseudo-ellipse_torsionless_in_frame_theta_x_calc.eps}  
      
 \includegraphics[width=.5\textwidth]{images/frame_order_matrix/Sijkl_pseudo-ellipse_torsionless_in_frame_theta_y_ens1000000.eps} \includegraphics[width=.5\textwidth]{images/frame_order_matrix/Sijkl_pseudo-ellipse_torsionless_in_frame_theta_y_calc.eps}  

Figure 16.21: The torsionless pseudo-ellipse model simulated and calculated out-of-frame Daeg(1) frame order matrix elements. In these plots, θX corresponds to the cone opening half-angle θx and θY to the cone opening half-angle θy. When the half-angle is not varied, the angle is fixed to either θx = π/4 or θy = 3π/8. Frame order matrix values have been calculated every 10 degrees. The first angle for the calculated elements is set to 0.01 degrees as a pseudo-ellipse cone opening angle of 0.0 cannot be correctly handled by the numerical integration.
 \includegraphics[width=.5\textwidth]{images/frame_order_matrix/Sij_pseudo-ellipse_torsionless_out_of_frame_theta_x_ens1000000.eps} \includegraphics[width=.5\textwidth]{images/frame_order_matrix/Sij_pseudo-ellipse_torsionless_out_of_frame_theta_x_calc.eps}  
      
 \includegraphics[width=.5\textwidth]{images/frame_order_matrix/Sij_pseudo-ellipse_torsionless_out_of_frame_theta_y_ens1000000.eps} \includegraphics[width=.5\textwidth]{images/frame_order_matrix/Sij_pseudo-ellipse_torsionless_out_of_frame_theta_y_calc.eps}  

Figure 16.22: The torsionless pseudo-ellipse model simulated and calculated out-of-frame Daeg(2) frame order matrix elements. In these plots, θX corresponds to the cone opening half-angle θx and θY to the cone opening half-angle θy. When the half-angle is not varied, the angle is fixed to either θx = π/4 or θy = 3π/8. Frame order matrix values have been calculated every 10 degrees. The first angle for the calculated elements is set to 0.01 degrees as a pseudo-ellipse cone opening angle of 0.0 cannot be correctly handled by the numerical integration.
 \includegraphics[width=.5\textwidth]{images/frame_order_matrix/Sijkl_pseudo-ellipse_out_of_frame_theta_x_ens1000000.eps} \includegraphics[width=.5\textwidth]{images/frame_order_matrix/Sijkl_pseudo-ellipse_out_of_frame_theta_x_calc.eps}  
      
 \includegraphics[width=.5\textwidth]{images/frame_order_matrix/Sijkl_pseudo-ellipse_out_of_frame_theta_y_ens1000000.eps} \includegraphics[width=.5\textwidth]{images/frame_order_matrix/Sijkl_pseudo-ellipse_out_of_frame_theta_y_calc.eps}  

Torsionless pseudo-ellipse rotation matrices

Setting the torsion angle σ to zero in the full torsion-tilt rotation matrix of equation 12.74c, the matrix becomes

R(θ, φ) = $\displaystyle \begin{pmatrix}
\cos^2\phi \cos\theta + \sin^2\phi & \cos\phi \si...
...\
- \cos\phi \sin\theta & - \sin\phi \sin\theta & \cos\theta \\
\end{pmatrix}$. (16.55)


Torsionless pseudo-ellipse frame order matrix

The frame order matrix is

\begin{subequations}\begin{align}
\textrm{Daeg}^{(n)}&= \left. \int_S R^{\otimes...
...a \,\mathrm{d}\phi \right / \int_S \,\mathrm{d}S .
\end{align}\end{subequations}

The surface normalisation factor is

\begin{subequations}\begin{align}
\int_S \,\mathrm{d}S &= \int_{-\pi}^{\pi} \int...
...thrm{d}\phi , \\
&= \pec (\theta_x , \theta_y ) .
\end{align}\end{subequations}

Torsionless pseudo-ellipse 1 degree frame order

The 1 degree frame order matrix with tensor rank-2 consists of the following elements

\begin{subequations}\begin{align}
\textrm{Daeg}_{00} &= \frac{1}{2\pec (\theta_x...
...pi}
\sin^2\theta_{\textrm{max}}
\,\mathrm{d}\phi .
\end{align}\end{subequations}

As the trigonometric functions of θmax cannot be symbolically integrated, these components must be numerically integrated.

Torsionless pseudo-ellipse 2 degree frame order

The 2 degree frame order matrix with tensor rank-4 consists of the following elements, using Kronecker product double indices from 0 to 8

\begin{subequations}\begin{flalign}
\begin{split}
\textrm{Daeg}_{00} = \frac{1}{...
...rm{max}}
\,\mathrm{d}\phi
\Bigg] ,
\end{split} &
\end{flalign}\end{subequations}\begin{subequations}\begin{flalign}
\begin{split}
\textrm{Daeg}_{11} = \frac{1}{...
...rm{max}}
\,\mathrm{d}\phi
\Bigg] ,
\end{split} &
\end{flalign}\end{subequations}\begin{subequations}\begin{flalign}
&\textrm{Daeg}_{22} = \frac{1}{6\pec (\theta...
...theta_{\textrm{max}}
\,\mathrm{d}\phi
\Bigg] , &
\end{flalign}\end{subequations}\begin{subequations}\begin{flalign}
&\textrm{Daeg}_{33} = \textrm{Daeg}_{11} , &
\end{flalign}\end{subequations}\begin{subequations}\begin{flalign}
\begin{split}
\textrm{Daeg}_{44} = \frac{1}{...
...rm{max}}
\,\mathrm{d}\phi
\Bigg] ,
\end{split} &
\end{flalign}\end{subequations}\begin{subequations}\begin{flalign}
&\textrm{Daeg}_{55} = \frac{1}{6\pec (\theta...
...theta_{\textrm{max}}
\,\mathrm{d}\phi
\Bigg] , &
\end{flalign}\end{subequations}\begin{subequations}\begin{flalign}
&\textrm{Daeg}_{66} = \textrm{Daeg}_{22} , &
\end{flalign}\end{subequations}\begin{subequations}\begin{flalign}
&\textrm{Daeg}_{77} = \textrm{Daeg}_{55} , &
\end{flalign}\end{subequations}\begin{subequations}\begin{flalign}
&\textrm{Daeg}_{88} = \frac{1}{3\pec (\theta...
...theta_{\textrm{max}}
\,\mathrm{d}\phi
\Bigg] , &
\end{flalign}\end{subequations}\begin{subequations}\begin{flalign}
\begin{split}
\textrm{Daeg}_{04} = \frac{1}{...
...rm{max}}
\,\mathrm{d}\phi
\Bigg] ,
\end{split} &
\end{flalign}\end{subequations}\begin{subequations}\begin{flalign}
&\textrm{Daeg}_{40} = \textrm{Daeg}_{04} , &
\end{flalign}\end{subequations}\begin{subequations}\begin{flalign}
&\textrm{Daeg}_{08} = \frac{1}{3\pec (\theta...
...theta_{\textrm{max}}
\,\mathrm{d}\phi
\Bigg] , &
\end{flalign}\end{subequations}\begin{subequations}\begin{flalign}
&\textrm{Daeg}_{80} = \textrm{Daeg}_{08} , &
\end{flalign}\end{subequations}\begin{subequations}\begin{flalign}
&\textrm{Daeg}_{48} = \frac{1}{3\pec (\theta...
...theta_{\textrm{max}}
\,\mathrm{d}\phi
\Bigg] , &
\end{flalign}\end{subequations}\begin{subequations}\begin{flalign}
&\textrm{Daeg}_{84} = \textrm{Daeg}_{48} , &
\end{flalign}\end{subequations}\begin{subequations}\begin{flalign}
&\textrm{Daeg}_{13} = \textrm{Daeg}_{04} , &
\end{flalign}\end{subequations}\begin{subequations}\begin{flalign}
&\textrm{Daeg}_{31} = \textrm{Daeg}_{04} , &
\end{flalign}\end{subequations}\begin{subequations}\begin{flalign}
&\textrm{Daeg}_{26} = -\textrm{Daeg}_{80} , &
\end{flalign}\end{subequations}\begin{subequations}\begin{flalign}
&\textrm{Daeg}_{62} = -\textrm{Daeg}_{80} , &
\end{flalign}\end{subequations}\begin{subequations}\begin{flalign}
&\textrm{Daeg}_{57} = -\textrm{Daeg}_{48} , &
\end{flalign}\end{subequations}\begin{subequations}\begin{flalign}
&\textrm{Daeg}_{75} = -\textrm{Daeg}_{48} . &
\end{flalign}\end{subequations}    

Torsionless pseudo-ellipse frame order matrix simulation and calculation

The frame order matrix element simulation script from Section 16.2, page [*] was used to compare the implementation of equations 16.58 and 16.59 above. Frame order matrix Daeg(1) and Daeg(2) values were both simulated and calculated, both within and out of the motional eigenframe. The in-frame Daeg(1) values are shown in figure 16.19 and Daeg(2) in figure 16.20. The out-of-frame Daeg(1) values are shown in figure 16.21 and Daeg(2) in figure 16.22.

The relax user manual (PDF), created 2020-08-26.