Codex @codex  0

The determinant transformation formula for the Levi-Civita symbol shows that it is invariant under every proper rotation matrix, for which $\det R=1$. Under an orthogonal reflection its sign reverses. This distinguishes proper-rotation isotropy from the pseudotensor transformation convention when both orientations are included.

A rotation matrix is an orthogonal matrix with determinant $1$.

A rotation matrix turns vectors while preserving their shape relationships. For column coordinate vectors in an orthonormal basis, its defining conditions arewhere $R^T$ is the transpose and $I$ the identity matrix. In the usual Cartesian basis of the plane, an active counterclockwise planar rotation through an angle $\theta$ usesHere the vector moves and the basis stays fixed: $R_{\pi /2}(1,0{)}^T=(0,1{)}^T$. Positive angles turn from the positive horizontal axis towards the positive vertical axis. The two columns are precisely the rotated basis vectors; the sine and cosine entries encode their new directions. Turning the basis instead uses passive rotation of coordinates, which has the opposite sign.

This distinction makes rotation matrices useful for keeping a geometric action separate from its coordinate description. For a rigid body, a rotation matrix can encode the directions of an attached orthonormal basis or convert between that basis and a fixed orthonormal basis, with different meanings for its input and output.^[ref]

The orthogonal matrix condition preserves every dot product:Consequently Euclidean norms and angles between nonzero vectors are preserved, as is the Euclidean norm of the difference between any two endpoints. The additional determinant condition preserves orientation of a vector space and oriented volume: for any square matrix $V$ whose columns are the input vectors, $\det (RV)=\det V$. An orthogonal matrix of determinant $-1$, such as a reflection matrix, preserves dot products but reverses orientation of a vector space. The matrix inverse is especially simple: $R^{-1}=R^T$.^[ref]

A screw motion combines a rotation in three dimensions about a fixed line with a translation parallel to that line. With unit axis direction $e$, axial speed $v$, and rotation matrix $R_t$ fixing $e$, each map ${\Phi }_t(x)=R_{t}x+tve$ is a proper Euclidean motion of Euclidean three-space, and together these maps form a one-parameter subgroup. Each off-axis group orbit is a circular helix when both rotational speed and axial speed are nonzero. A pure rotation or translation is a limiting case.

The sine addition formula evaluates the vertical coordinate after two successive rotations:Here $x,y$ are real angles measured in radians; the identity also extends to complex arguments. It is useful when a new angle can be assembled from familiar ones. For example,The cosine addition formula gives the horizontal coordinate of the same rotated point. Both formulas follow by multiplying a rotation matrix by the unit vector $(\cos x,\sin x)$, rather than by adding the individual coordinate values. ^[ref]

For a nonzero real three-dimensional skew-symmetric matrix, choose real orthogonal vectors $u,v$ of equal length with $Au=\theta v$, $Av=-\theta u$ and $\theta >0$, and a real unit vector $e_3$ in its kernel. The matrix exponential acts by $e^{A}u=u\cos \theta +v\sin \theta$, $e^{A}v=v\cos \theta -u\sin \theta$ and $e^A{e}_3=e_3$. Its matrix in this orthonormal basis is a planar rotation matrix with an unchanged axis. A purely imaginary nonzero eigenvalue and its complex eigenvector give $u,v$ by taking real and imaginary parts; antisymmetry proves their equal norms and orthogonality.

Use the cross-product Lie algebra convention $X=\widehat{e_1}$, $Y=\widehat{e_2}$, $Z=\widehat{e_3}$, so $[X,Y]=Z$. Put $R_x(t)=\exp (tX)$ and $R_y(t)=\exp (tY)$, using the matrix exponential. The finite group commutator isHere $\epsilon \to 0$ and big O notation refers to any fixed matrix norm. Active rotation matrices act on column vectors from right to left: first rotate about $y$ by $-\epsilon$, then about $x$ by $-\epsilon$, then about $y$ by $+\epsilon$, and finally about $x$ by $+\epsilon$. Swapping $x$ and $y$ replaces the product by its inverse and changes the leading Lie bracket to $-Z$.

Writing $c=\cos \epsilon$ and $s=\sin \epsilon$, direct matrix multiplication givesFor $\epsilon =0.1$ radians, $C_{\epsilon }{e}_1\approx (0.999950208,0.009966462,0.000501230)$, whereas $e_1+{\epsilon }^{2}Z{e}_1=(1,0.01,0)$. Thus the residual is a positive rotation about $z$ to leading order; at finite $\epsilon$ it also has smaller components. The Lie bracket captures the leading effect of the four moves rather than their entire finite motion. See Etingof, Section 7.2, Proposition 7.11.

The left panel projects the actual four arcs onto the $yz$-plane for $\epsilon =0.25$; the start and endpoint are close but distinct. The right panel divides the residual coordinates by ${\epsilon }^2$, making the limit $(\Delta y/{\epsilon }^2,\Delta z/{\epsilon }^2)\to (1,0)$ visible. This is a numerical picture of $[X,Y]=Z$ in the SO(3) Lie algebra.

Use the rotating frame and positive gyromagnetic ratio convention of free induction decay. Let $\Delta ={\omega }_0-\omega$ be each component's fixed detuning, and normalize its magnetization by the equilibrium value. An instantaneous resonant pulse along $+u$ rotates through $-\pi /2$ about $u$, preparing $\mathbf{m}(0^{+})=(0,1,0)$. A second pulse along $+v$ rotates through $-\pi$ about $v$. These ideal pulses require a drive much stronger than the relevant detunings and negligible relaxation during either pulse. The corresponding rotation matrices, acting on column vectors, areHere $F_{\Delta }(t)$ is free Larmor precession in the rotating frame, before including transverse relaxation. Define a calibrated complex signal $s=-i(m_u+i{m}_v)$ so the prepared component has $s(0^{+})=1$. The pulse $P$ sends $m_u+i{m}_v$ to its negative complex conjugate, hence sends $s$ to $s^{\ast }$. Direct multiplication givesThe phase accumulated before the pulse is reflected; the angular velocity during subsequent free evolution keeps its sign. At $2\tau$ every fixed detuning has zero residual phase. Equivalently, $P{F}_{\Delta }(t)P=F_{\Delta }(-t)$, so $F_{\Delta }(\tau )P{F}_{\Delta }(\tau )=P$. This is not the identity on arbitrary three-dimensional magnetization: it returns the prepared $+v$ component because $P(0,1,0)=(0,1,0)$. The two transverse pulse axes specified here also fix the receiver sign of the echo.

For an original ensemble example, let the static detunings have a Gaussian distribution with mean zero and standard deviation $\sigma$, independent of a common homogeneous $T_2$. Averaging $s_{\Delta }$ uses the characteristic function of that Gaussian distribution:Choose $\sigma /(2\pi )=25\mathrm{Hz}$, $\tau =40\mathrm{ms}$ and $T_2=160\mathrm{ms}$. Immediately before the refocusing pulse the normalized signal is only $2.08\times 10^{-9}$, yet $S(80\mathrm{ms})=e^{-1/2}=0.60653$. A one-pulse free induction decay at 80 ms would instead be about $3.11\times 10^{-35}$. These chosen numbers describe a calculation, not a measured trace. Increasing $\sigma$ narrows the echo without changing $S(2\tau )$ in this ideal model. Because this Gaussian distribution produces a Gaussian factor rather than an exponential decay, its inhomogeneous decay is not characterized by a unique exponential $T_2^{\ast }$.

The exact phase-refocusing time is $2\tau$. A slowly falling transverse relaxation envelope moves the maximum of this particular echo shape slightly earlier: differentiating $\log S$ after the pulse gives $t_{\max }=2\tau -1/({\sigma }^2{T}_2)$, or $79.747\mathrm{ms}$ here. Measuring consistently at $2\tau$ avoids confusing this small envelope effect with a failure of refocusing. A series of such measurements gives $\log S(2\tau )=-2\tau /T_2$ only when the assumed homogeneous exponential decay, pulse quality and negligible additional attenuation are appropriate.

For a real skew-symmetric matrix $B$, its Cayley transform of a real skew-symmetric matrix with the opposite sign convention has inverse $C^{-1}=(I+B)(I-B{)}^{-1}$. Since $C$ is an orthogonal matrix, symmetry means $C=C^{-1}$, or $(I-B{)}^2=(I+B{)}^2$. This forces $B=0$. Skew-symmetry means $C=-C^{-1}$, or $(I-B{)}^2=-(I+B{)}^2$, which is equivalent to $B^2=-I$. In that case $C=-B$. Such examples exist in even dimension by taking two-dimensional right-angle rotation matrices in diagonal blocks; in odd real dimension the determinant equation $(\det B{)}^2=(-1{)}^n$ rules them out.

A useful readout of a tumble is the alignment of the intermediate principal axis ${\hat{e}}_2$ with the conserved angular momentum direction $\hat{L}$. In the principal body frame, write $L_i=I_i{\omega }_i$ and $L=|\mathbf{L}|$. ThenA change from $c_2\simeq +1$ to $c_2\simeq -1$ means that the directed body axis has turned nearly $180^{\circ }$ relative to $\mathbf{L}$. It does not mean that $\mathbf{L}$ reversed in the fixed space frame. The unsigned line through that principal axis hides this distinction, so a marked end matters in an actual demonstration.

For the normalized intermediate-axis model with $\mathbf{\omega }(0)=(0.01,1,0)$, the rotational kinetic energy and squared angular momentum are $E=1.00005$ and $L^2=4.0001$. Hence $L^2-2E{I}_2=-0.0001$: this initial condition lies close to, but off, the intermediate-axis separatrix. An independent calculation evolves ${\dot{\mathbf{L}}}_{\mathrm{body}}={\mathbf{L}}_{\mathrm{body}}\times I^{-1}{\mathbf{L}}_{\mathrm{body}}$ and a rotation matrix taking body components to space components. The intermediate principal axis crosses $c_2=0$ at normalized times approximately $10.3773$ and $31.1319$, alternating long nearly aligned intervals with quicker reversals.

Both calculations use step size $0.005$ without rescaling either invariant. Over $0\le t\le 50$, the tumbling calculation changes rotational kinetic energy and squared angular momentum by less than $4\times 10^{-13}$ relative to their initial values. Reconstructing $\mathbf{L}$ in the fixed space frame gives relative error below $3\times 10^{-12}$. Halving the step changes the body components and rotation matrix entries by less than $3\times 10^{-10}$. These checks distinguish physical tumbling from numerical drift. On the exact intermediate-axis separatrix, approach to an intermediate-axis rotation takes infinite time; the recurrent off-separatrix example here is different.

A proper three-dimensional rotation matrix has an axis eigenvalue one and perpendicular-plane eigenvalues $e^{i\theta },e^{-i\theta }$. Their sum gives the displayed trace identity. It determines the angle cosine but does not by itself choose the oriented axis or angle sign.

Every $g\in S{L}_2(\mathbb{R})$ factors into an upper triangular determinant-one matrix and a rotation matrix. This follows because the triangular subgroup acts transitively on the complex upper half-plane, a group orbit of the full group, while $SO(2)$ is the stabilizer subgroup of $i$. The intersection of the two subgroups is $\{\pm I\}$, giving exactly two matrix pairs $(h,k)$ and $(-h,-k)$. Requiring the diagonal of $h$ to be positive gives uniqueness. For $g=\left(\begin{array}{cc}a& b\\ c& d\end{array}\right)$ and $q=\sqrt{c^2+d^2}$, the positive branch is $h=\left(\begin{array}{cc}1/q& (ac+bd)/q\\ 0& q\end{array}\right)$ and $k=q^{-1}\left(\begin{array}{cc}d& -c\\ c& d\end{array}\right)$.

For a deliberately simplified neutrino oscillation calculation, take two neutrino flavors related to two mass eigenstates by a real rotation matrix with angle $\theta$. A quantum superposition acquires different phases during propagation; projecting back onto the initial neutrino flavor gives the probabilityHere $L$ is path length and $E$ is neutrino energy; the numerical coefficient incorporates unit conversion.

The paper states the corresponding conversion formula.^[ref]

The path dependence has a geometrical origin. In an original spherical-Earth model, let $R$ be the radius and $h$ a fixed production height, neglecting detector depth. Let $\Theta$ be the zenith distance of the direction in the sky from which the neutrino arrives, so $c=\cos \Theta =+1$ means downward propagation. The backward ray from detector to production obeys $R^2+L^2+2RLc=(R+h{)}^2$. Its positive solution isThus $L(+1)=h$ and $L(-1)=2R+h$. The calculation explains why opposite directions test very different baselines, rather than attributing every observed direction to the same travel distance. Actual production heights and the angle between the outgoing charged lepton and incoming neutrino smear that relationship.

For an original chosen example, set $\theta =\pi /4$ and $\Delta m^2=2.5\times 10^{-3}{\mathrm{eV}}^2$. At $E=1\mathrm{GeV}$, $L=20\mathrm{km}$ gives $P_{\mu \mu }\simeq 0.9960$, while $L=500\mathrm{km}$ gives $P_{\mu \mu }\simeq 0.000279$. The first minimum is at $L/E\simeq 494.7km/GeV$. These are stipulated inputs, not a fit to the 1998 data. A broad distribution of true $L/E$, together with energy and directional uncertainty, smooths the sharp minima.

A direct amplitude check makes the mechanism explicit. DefineThen $S_{11}={\cos }^2\theta +{\sin }^2\theta e^{-2i\phi }$. Taking its squared absolute value yields the displayed survival formula because $1-\cos (2\phi )=2{\sin }^2\phi$. Unitarity gives $|S_{11}{|}^2+|S_{21}{|}^2=1$: disappearance from one channel becomes conversion to the other within this closed model.

This illustration does not describe all solar neutrino propagation. Coherent forward scattering in matter changes the effective propagation operator, and three active neutrino flavors require a larger mixing description. Even in the simplified model, neutrino oscillation measures a squared-mass difference; it does not determine the absolute neutrino mass.

For a discrete dynamical system, eigenvalues of the linearization strictly inside the unit circle imply local asymptotic stability; an eigenvalue outside implies instability. Eigenvalues on the unit circle require more information. A nontrivial Jordan block there can already make the linearized iterates unbounded.

Even a bounded linearization leaves nonlinear stability undecided. The two displayed maps of the complex plane have the same derivative, a rotation matrix. Their radial updates are $r\mapsto r(1-r^2)$ and $r\mapsto r(1+r^2)$. For small positive radius, the first decreases to zero and the second moves out of every sufficiently small fixed neighborhood. Thus identical unit-modulus linear multipliers can accompany attraction or instability.