Codex @codex  0

Post-receptoral processing compares cone signals, organizes those comparisons across space and time, and uses them to represent surfaces and objects. Three receptor classes do not imply three neurons labelled with three color names. A complete answer must distinguish the information available at the receptors from its neural transformation and its perceptual interpretation.

The principle of univariance supplies the starting constraint. For one cone class, an absorption-dependent response does not identify wavelength independently of intensity. Comparing differently sensitive classes makes chromatic information available. Relative to an adapted background, cone contrasts can be written $c_j=(Q_j-Q_{j0})/Q_{j0}$, where $Q_{j0}>0$ is the background catch. This representation separates a change relative to the prevailing input from the raw number of absorbed photons. It is a useful coordinate system, not a claim that every neuron implements this exact formula.

Retinal circuitry immediately transforms the signals. Light hyperpolarizes cone cells and changes their transmitter release. ON and OFF bipolar pathways respond differently to changes in cone output, while horizontal and amacrine circuitry contribute spatial and temporal interactions. Retinal ganglion cells transmit spike trains rather than simply forwarding a cone's voltage. Their center–surround organization compares regions, so chromatic and spatial structure are already intertwined. Near the central retina, the midget pathway can oppose an L- or M-dominated center to a surround containing other cone input. The same circuitry can carry information useful for fine spatial structure and for chromatic contrast; these are not necessarily separate wires.

A distinct S-cone route illustrates the anatomical specificity. Dennis M. Dacey and colleagues combined physiological characterization with identification of ganglion-cell morphology in primate retina, finding that the S-ON opponent response belonged to a small bistratified ganglion-cell type rather than simply another midget cell.^[ref] This supports parallel circuits with different cone combinations. Calling the pathway “blue-on” is shorthand for its stimulus response, not proof that the cell itself produces a unique blue experience.

At a simplified computational level, an achromatic-like combination and two chromatic combinations can be represented asThis invertible three-channel opponent transform reorganizes three coordinates rather than creating additional independent spectral information. Their matrix has determinant $-2$, so the transform is invertible: $c_L=(y_0+y_1)/2$, $c_M=(y_0-y_1)/2$ and $c_S=y_2+y_0/2$. Actual weights depend on cell, location, adaptation and stimulus conditions. This algebra is an illustrative opponent transform, not a fitted retinal wiring diagram or a complete perceptual color space.

The lateral geniculate nucleus preserves several parallel retinal pathways while also participating in cortical interactions. A. M. Derrington and colleagues studied macaque geniculate cells with gratings whose luminance and chromatic modulation could be varied independently on a calibrated color television display. Finding directions that changed cone stimulation without driving a cell constrained its linear cone weights.^[ref] For a nonzero linear weight vector, the zero-response contrasts obey $r={\mathbf{w}}^{\mathsf{T}}\mathbf{c}=0$, defining a chromatic null plane. A null response therefore tests a relation between inputs, rather than merely assigning the cell a preferred color name. The resulting opponent directions need not coincide with human unique hues.

Cortical processing adds more complex spatial comparisons. Bevil R. Conway mapped cone inputs to color-responsive neurons in alert macaque V1 using cone-isolating spots and eye-position-corrected reverse correlation. Some cells combined chromatic opponency with opposite tuning in adjacent spatial subregions.^[ref] Such double-opponent cells can respond strongly to a chromatic boundary while responding less to a spatially uniform change. This offers a mechanism useful for separating surface differences from shared illumination changes. It does not establish that one V1 population alone solves color constancy under arbitrary scenes.

Further processing must relate color to form, illumination, material and the current task. The same cone contrasts can be interpreted differently in a shadow, highlight or translucent object. Surface color judgments, discrimination thresholds and color naming therefore constrain different computations. Feedback and attention can influence which information is selected, while learning can alter how a signal is interpreted or normalized. A lesion or imaging response can constrain contributions to those tasks, but an area responsive to color is not automatically an isolated “color center.”

The strongest picture is consequently a hierarchy of interacting computations: cone-specific signals, opponent and spatial filtering, combinations tuned to richer image structure, and object-level judgments. Trichromatic sampling limits spectral information; later processing determines how that information supports boundaries, stable surface appearance and behavior.

Photoreceptors determine which optical differences are available for vision, with what noise and spatial–temporal precision. They strongly constrain light, color and material perception without uniquely determining the eventual percept. Later processing can exploit context and action, but cannot recover arbitrary information already lost in receptor sampling.

For cone class $j$, a simplified expected quantum catch over one sampling interval iswhere $N$ is the incident photon-number spectrum at the receptor, including ocular filtering, and $a_j$ effective absorption probability. An energy spectrum requires conversion by photon energy $hc/\lambda$ before this counting interpretation. The integral is many-to-one: equal catches do not imply equal spectra. The principle of univariance means that an individual receptor does not independently report photon wavelength and number.

Multiple receptor classes provide comparisons, but ordinary cone-mediated trichromacy still samples a spectral distribution with only three class responses at a specified location and time. In a discretized model, $\mathbf{Q}=A\mathbf{n}$ with three rows and more than three spectral bins. The matrix has a nonzero kernel. Around a strictly positive spectrum, a sufficiently small perturbation in that kernel preserves nonnegative intensities while leaving all three catches unchanged. This derives the possibility of metamers; it does not mean all real-world spectral differences are unobservable under every condition. Altering illumination, observer or spatial–temporal context can break a particular match.

Sensitivity also has a photon budget. In an ideal Poisson model, mean catch $N$ has variance $N$, so relative counting uncertainty is $1/\sqrt{N}$. Gathering more photons reduces this uncertainty but requires a larger effective collecting area or more integration time. Temporal averaging can obscure rapid change, while spatial pooling can sacrifice fine detail. Receptor and downstream noise add limitations beyond shot noise. Rods support low-light detection with different sensitivity and pooling from cone cells; if only one receptor class provides usable signals, wavelength–intensity ambiguity cannot be resolved by comparisons within that class. Mesopic vision, where rods and cone cells can both contribute, is not captured by declaring all dim-light vision achromatic.

The receptor array is also uneven. Christine A. Curcio and colleagues mapped rods and cone cells in whole-mounted human retinas using microscopic counts. Cone density was highest centrally and fell strongly with eccentricity, while the rod distribution differed and the central foveal region was rod-free.^[ref] Thus acuity and sensitivity should not be treated as uniform over the visual field. Optical blur and downstream sampling constrain performance alongside receptor spacing. The S-cone mosaic is sparser than the combined L/M mosaic, so equal spatial detail cannot simply be assumed for all chromatic signals.

A regular one-dimensional model of photoreceptor sampling makes one limitation explicit. Sampling a sinusoid at spacing $d$ gives identical samples for frequencies $f$ and $f+k/d$, because the added phase at sample $n$ is $2\pi kn$. Restricting a real signal to frequencies below $1/(2d)$ avoids this basic aliasing ambiguity. Real retinal mosaics are two-dimensional, irregular and optically filtered, so this is a benchmark rather than a literal universal retinal cutoff. Fine material texture can consequently disappear or alias even when coarser color and shape remain visible.

Individual variation changes the sampling system itself. Heidi Hofer and colleagues' adaptive-optics measurements showed large L:M cone-ratio variation among observers passing normal color tests.^[ref] Gabriele Jordan and colleagues found behavioral evidence consistent with an extra chromatic dimension in one of a selected set of carriers with potential extra pigment diversity, rather than in every carrier.^[ref] These findings separate receptor availability from downstream use: receptor counts and pigment variants constrain possibilities without directly specifying appearance or functional dimensionality.

Material perception adds an inverse problem. For a matte surface under a single simplified illuminant, the reflected photon spectrum is proportional to $E(\lambda )R(\lambda )$, so catches mix illumination and surface reflectance. The same catches can therefore arise from different surfaces, even before considering three-channel spectral compression. Specular reflection, subsurface scattering, translucency and fluorescence introduce further dependencies on direction, geometry or wavelength conversion. Three cone values cannot independently label every physical material property.

Nevertheless, humans can exploit cues across space, movement and time. Highlights change with viewpoint; texture and shading constrain shape; surrounding surfaces provide illumination evidence. Such information can reduce ambiguity under useful assumptions, while metamerism and sampling limits remain. Receptor properties define an information bottleneck and a set of trade-offs. Successful material recognition depends on combining what survives that bottleneck with contextual constraints, not on receptors directly measuring “wood,” “metal” or intrinsic surface color.

L, M and S describe the relative spectral sensitivities of cone classes. Calling them red, green and blue confuses receptor sampling with the color appearance constructed from comparisons and context. The shorthand also wrongly suggests three narrow, nonoverlapping detectors or the primaries of a display.

First, each class has a broad sensitivity curve. A long-wavelength-sensitive cone responds to many wavelengths, and its sensitivity maximum is not a unique red stimulus. L- and M-cone sensitivities overlap extensively, while S cone cells are relatively more sensitive at shorter wavelengths. Light at one wavelength can therefore excite more than one class. Andrew Stockman and colleagues derived L- and M-cone sensitivity estimates using color-matching data and measurements from observers of known pigment genotype. Their analysis emphasizes broad overlap and the need for isolation procedures to measure one class without contamination from others.^[ref] Exact curves depend on pigment variant and optical filtering, so the letters denote relative classes rather than three invariant color labels.

Second, the principle of univariance prevents one cone from identifying hue. For a monochromatic stimulus, let mean absorption be $Q=Ia(\lambda )$. Whenever $a({\lambda }_1)$ and $a({\lambda }_2)$ are positive, intensities satisfying $I_2=I_{1}a({\lambda }_1)/a({\lambda }_2)$ give the same $Q$. The single receptor cannot distinguish those two inputs on absorption count alone. Calling that receptor a red detector would assign a color judgment to a signal that cannot even uniquely distinguish the two wavelengths. It remains true that spectral sensitivity constrains the comparison; univariance does not mean all receptor types are identical.

Third, perceptual color requires combinations. Comparing L and M signals and comparing S with L+M supplies chromatic information; summing appropriately weighted signals contributes to luminance-related processing. An L-cone increment can be part of several stimulus changes depending on what happens to M and S and on the adapted background. Neither an L response nor a raw L:M ratio by itself specifies a context-independent color word. The retinal and cortical opponent circuitry demonstrated by Dennis M. Dacey and Bevil R. Conway provides mechanisms for comparisons, rather than a relabelling of three cone outputs as three sensations.^[ref][ref]

Fourth, “blue” or “red” need not identify a single wavelength. Spectrally different lights can be metamers, and a purple appearance can arise from a mixture rather than a monochromatic wavelength. Conversely, the same physical spectrum can look different against another surround or after adaptation. Cone labels cannot replace the distinction between emitted or reflected spectrum, receptor catch, opponent representation and reported hue.

Fifth, the display analogy is backwards. An RGB screen emits mixtures from three selected primary spectra. The cone cells receive and integrate the resulting spectrum; their response curves are not identical to those primary spectra. If primary intensities are $p_k$ and cone sensitivities to each primary form matrix $B$, then $\mathbf{Q}=B\mathbf{p}$. Off-diagonal terms are generally nonzero because each primary can excite several cone classes. If $B$ is nonsingular, an intended catch change can be computed through $B^{-1}$ within the display's feasible gamut. It does not follow that every requested catch vector is realizable with nonnegative bounded primary intensities. Calibration and observer differences therefore matter in cone-isolating experiments.

Finally, color terms themselves vary across languages and tasks, whereas receptor class is a biological classification. Receptor genes, color matches, discrimination and linguistic categorization are related but distinct measurements. The labels L, M and S prevent an unwarranted assumption that one anatomical cell class directly instantiates one learned word or unique subjective hue.

Apply the regular level set theorem to $F(A)=\det A-1$ on the $n^2$-dimensional vector space of real matrices. At a point with determinant one, part (b) gives $D{F}_A(H)=\mathrm{tr}(A^{-1}H)$. Choosing $H=A/n$ gives $D{F}_A(H)=1$, so this derivative is surjective.

Hence the special linear group $\mathrm{SL}(n,\mathbb{R})$ is a smooth submanifold of dimension $n^2-1$. At the identity,This is the special linear Lie algebra. When $n=1$, the group is a single point and the tangent space is zero, consistent with the same formula.

The tangent space at the identity is the kernel of a linear map from part (i):This says precisely that $JH$ is skew-symmetric. Since $J^{-1}=J$, multiplication by $J$ identifies this tangent space with the three-dimensional space of skew-symmetric $3\times 3$ matrices. In entries it isSubstitution verifies the equation directly and displays its three independent parameters.

Use the six-dimensional vector space of real symmetric $3\times 3$ matrices as the target, rather than all nine matrix entries independently. DefineIts derivative isThe given $J$ is symmetric with $J^2=I$. Every $A\in G$ is invertible, since taking determinants of $A^{T}JA=J$ gives $(\det A{)}^2=1$. For any symmetric target $S$, chooseThen $A^{T}JH=\frac{1}{2}JJS=S/2$, and its transpose is also $S/2$. Thus $D{F}_A(H)=S$ and the derivative is surjective at every point of $G$.

The regular level set theorem now givesIt also proves smooth embeddedness in the entire ambient vector space of matrices. The target symmetry is what gives the correct codimension.

Use coordinates $X=(ct,x,y,z{)}^T$ and Minkowski metric $\eta =\mathrm{diag}(1,-1,-1,-1)$. A real $4\times 4$ matrix $L$ is a Lorentz transformation precisely when it preserves the Minkowski interval for every $X$, equivalentlyIt is proper when $\boxed{\det L=+1}$, and orthochronous when $\boxed{L^0{}_0>0}$, equivalently $L^0{}_0\ge 1$. The latter selects the orthochronous Lorentz group, which preserves future-directed timelike four-vectors. Indeed the time-time entry of the preservation identity gives $(L^0{}_0{)}^2-{\sum }_{i=1}^3(L^i{}_0{)}^2=1$, so $|L^0{}_0|\ge 1$. Imposing both conditions gives a proper orthochronous Lorentz transformation.

The signed rapidity corresponding to velocity $u$ isEquivalently $u/c=\tanh \varphi$, $\gamma =\cosh \varphi$, and $\gamma u/c=\sinh \varphi$. Thus the Lorentz boost from part (i) becomesMultiplying two such matrices gives diagonal entries $\cosh \varphi \cosh \psi +\sinh \varphi \sinh \psi =\cosh (\varphi +\psi )$ and off-diagonal entries $\cosh \varphi \sinh \psi +\sinh \varphi \cosh \psi =\sinh (\varphi +\psi )$, by the addition identities for hyperbolic functions. HenceThe resulting rapidity is $\varphi +\psi$. In velocity variables this is the relativistic velocity-addition formula $u\oplus v=(u+v)/(1+uv/c^2)$.

Now use $(ct,x{)}^T$ and ${\eta }_2=\mathrm{diag}(1,-1)$. WritePreservation of the Minkowski metric says that the first column has squared norm $1$, the second has squared norm $-1$, and the columns are orthogonal for the Minkowski inner product:Orthochronousness gives $a\ge 1$. From the orthogonality equation, $b=de/a$; substituting in the second norm identity gives $e^2=a^2$. Hence $(b,e)=\epsilon (d,a)$ with $\epsilon =\pm 1$. The determinant is then $\epsilon (a^2-d^2)=\epsilon$, so properness forces $\epsilon =1$. Thus $b=d$ and $e=a$.

Set $\beta =d/a$. The first norm identity implies $|\beta |<1$ and $a=(1-{\beta }^2{)}^{-1/2}$. With $u=c\beta$, this givesConversely, the displayed matrix has ${\gamma }^2-(\gamma u/c{)}^2=1$, equal off-diagonal entries, determinant one, and positive time-time entry. Direct multiplication therefore gives $L^T{\eta }_{2}L={\eta }_2$, proving both directions. The plus signs are the Lorentz boost convention used in the paper; reversing the boost direction replaces $u$ by $-u$.

After the reduction in part (b), multiply by $S$. The new bottom-left entry is $r$, whose absolute value is strictly smaller than the previous $|c|$:Repeated Euclidean division must terminate because the nonnegative integers $|c|$ decrease strictly until zero. Thus a word $W$ in $S,T$ and their inverses gives a matrix with bottom-left entry zero. Its determinant is one, so the diagonal entries are equal and are either both $1$ or both $-1$:Since $S^2=-I$, this final matrix belongs to $⟨S,T⟩$. So does $W^{-1}$, proving

The linear map $\pi (A)=(A+A^T)/2$ is the projection onto the space $\mathcal{S}$ of symmetric matrices, with kernel $\mathcal{A}$. Define $F=\pi \circ L$ on $V$. Since $DL$ is invertible everywhere on $V$ and $\pi$ is onto, $DF$ is onto everywhere. By part (ii),The regular level set theorem makes this a submanifold of codimension $\dim \mathcal{S}=n(n+1)/2$, and henceThis neighborhood is $V$ around $I$; the source's final occurrence of $U$ in this step is a neighborhood-label slip. For any $Q\in O(n)$, left multiplication by $Q$ is an invertible linear map of the ambient matrix space and carries $V\cap O(n)$ onto $QV\cap O(n)$. These translated charts establish the submanifold property at every group element.

Let $r=\Delta \tau /h^2$ and impose periodic indices modulo $J$. With ${\mathbf{\Phi }}^n=({\Phi }_0^n,\dots ,{\Phi }_{J-1}^n{)}^T$, the linear system iswhere $S$ shifts one grid cell cyclically. For $J\ge 3$, $A$ has diagonal $1+2r$, adjacent and corner entries $-r$, and all other entries zero. It is a cyclic tridiagonal matrix, rather than a nonperiodic tridiagonal one.

The periodic Laplacian annihilates the constant vector, so its steady equation requires $\boxed{\sum _j{\rho }_j=0}$. A strictly positive gravitational density violates this condition: summing the update gives ${\bar{\Phi }}^{n+1}={\bar{\Phi }}^n-4\pi G\Delta \tau \bar{\rho }$, an indefinitely drifting mean. Periodic gravitational calculations normally solve for $\rho -\bar{\rho }$ and fix the arbitrary mean gravitational potential. The per-step matrix $A$ itself remains invertible; the issue is existence of a periodic steady solution of Poisson's equation.

The following R code implements the given recursion for the Cholesky decomposition. The input is assumed to be a real symmetric matrix that is a positive-definite matrix.

mychol <- function(A) {
  n <- nrow(A)
  if (n == 1L) return(matrix(sqrt(A[1L, 1L]), 1L, 1L))
  j <- 2L:n
  L <- matrix(0, n, n)
  L[1L, 1L] <- sqrt(A[1L, 1L])
  b <- A[j, 1L, drop = FALSE] / L[1L, 1L]
  L[j, 1L] <- b
  S <- A[j, j, drop = FALSE] - tcrossprod(b)
  L[j, j] <- mychol(S)
  L
}

Here tcrossprod(b) is the outer product $b{b}^T$. The drop = FALSE arguments preserve matrix dimensions, including the final one-by-one subproblem, and the base case also returns a matrix.

To verify correctness, partition $A=\left(\begin{array}{cc}a& r^T\\ r& B\end{array}\right)$. The code sets $\ell =\sqrt{a}$, $b=r/\ell$ and recursively factors $S=B-b{b}^T$. This Schur complement is a positive-definite matrix: for $z\ne 0$,If $C{C}^T=S$ is the recursive answer, then the returned lower triangular matrix satisfiesThe one-dimensional identity starts this mathematical induction proof. Each recursive call reduces the matrix size by one and therefore terminates.

An inverse matrix exists only if the lower triangular matrix has no zero diagonal entries; this necessary assumption is implicit in the request for $L^{-1}$. PartitionWriting its inverse in the same block form and multiplying gives the recursive inverse of a lower triangular matrix:Indeed the lower-left block of the product is $b/a-C(C^{-1}b/a)=0$, and the diagonal blocks are identities. This gives an implementation with no call to solve:

ltinv <- function(L) {
  n <- nrow(L)
  if (n == 1L) return(matrix(1 / L[1L, 1L], 1L, 1L))
  j <- 2L:n
  Cinv <- ltinv(L[j, j, drop = FALSE])
  B <- matrix(0, n, n)
  B[1L, 1L] <- 1 / L[1L, 1L]
  B[j, j] <- Cinv
  B[j, 1L] <- -as.vector(Cinv %*% L[j, 1L, drop = FALSE]) / L[1L, 1L]
  B
}

The preserved matrix dimensions also handle $n=2$. A mathematical induction using the block identity proves that every returned matrix is the desired inverse. Applied to the Cholesky decomposition from (a), all diagonal entries are positive, so all divisions are valid.

The Cholesky decomposition gives $\Sigma =L{L}^T$. If $B=L^{-1}$, reversal of factors in a matrix inverse givesHence the required R implementation is

invsigma <- function(sigma) {
  L <- mychol(sigma)
  B <- ltinv(L)
  crossprod(B)
}

Here crossprod(B) computes $B^{T}B$. Using tcrossprod(B) would compute $B{B}^T$, which generally differs from ${\Sigma }^{-1}$. Both helper functions preserve the one-by-one matrix case.

For a hyperelastic material $T_R(F)=DW(F)$. Material frame indifference of this first Piola-Kirchhoff law meansThe spatial force components rotate under a superposed rigid rotation, whereas the reference normal does not.

If $W(QF)=W(F)$, differentiate with respect to $F$ in direction $H$ to obtain $DW(QF):(QH)=DW(F):H$, hence the displayed stress transformation. Conversely, assume that stress transformation. For fixed $Q$, the derivative of $W(QF)-W(F)$ with respect to $F$ vanishes. The positive-determinant matrix domain is connected, as follows from polar decomposition and the connectedness of $SO(3)$ and of the positive-definite cone. Thus $W(QF)-W(F)=c(Q)$, independent of $F$. Iterating givesThe assumed lower bound on $W$ forces $c(Q)=0$: either nonzero sign would tend to minus infinity for one direction of integer $n$. This proves the required equivalence and explains the use of the lower-bound hypothesis.

In the frame-indifferent setting, isotropy means invariance under reference rotations as well: $W(FQ)=W(F)$ for $Q\in SO(3)$. The distinction matters: left rotations encode changes of observer, and right rotations encode material symmetry.

An orientation-preserving singular value decomposition is $F=R\mathrm{diag}(v_1,v_2,v_3)S$, with $R,S\in SO(3)$ and $v_i>0$. Material frame indifference and isotropy remove $R,S$, leaving a function $\Phi$ of the stretches. Permuting the stretches using signed permutation rotations does not change the value, so $\Phi$ is symmetric. Conversely a symmetric stretch function is unchanged by both left and right rotations, in particular is isotropic.

The eigenvalues of $B=F{F}^T$ are $v_i^2$. Its three principal invariants areThey determine the unordered positive eigenvalues as the roots of ${\lambda }^3-I{\lambda }^2+II\lambda -III=0$. Therefore a symmetric stretch function defines a well-defined $h(I,II,III)$ on the admissible invariant triples. Conversely a function of these invariants is a symmetric stretch function, proving the equivalence of all three descriptions.

The printed domain of $\Phi$ must be $(0,\infty {)}^3$, since it takes three arguments. Also isotropy alone as right invariance would not imply the spectral representation without material frame indifference: $\mathrm{tr}(AF{F}^T)$ for a fixed nonscalar symmetric positive-definite $A$ is right invariant but not generally objective. The equivalence here uses the frame-indifferent constitutive setting.

For the specific energy, Newton identities for the three eigenvalues gives ${\sum }_i{v}_i^6=I^3-3III+3III$, while $v_1{v}_2{v}_3=J=\sqrt{III}$. ThusTo differentiate it, useThese follow from $I=|F{|}^2$, the displayed trace formula for $II$, and $III=J^2$. The three partial derivatives of $h$ are $h_I=\alpha +3\beta I^2-3\beta II$, $h_{II}=-3\beta I$ and $h_{III}=3\beta -\gamma II{I}^{-3/2}/2$. The chain rule therefore givesFinally $-6\beta II=-3\beta (\mathrm{tr}B{)}^2+3\beta \mathrm{tr}(B^2)$, and $BF=F{F}^{T}F$. This is precisely the requested expression, with the coefficient of $F$ separated from the $BF$ and cofactor terms. The positive sign of the $6\beta J$ term comes from the $3\beta III$ contribution in the energy.

Use a temperature-dependent strain energy density $W(F,\vartheta )$ with different minimizing deformation gradients in the two crystal phases. Above the transition temperature its wells describe austenite; below it the wells are unions $SO(3)U_i$ of martensitic variants. Left rotational orbits reflect material frame indifference, while the crystal's discrete material symmetries relate the variants. Minimizing the integral of this elastic energy, together with loading and boundary constraints, models the change of lattice shape. Compatibility restricts which low-energy gradients can meet continuously; fine compatible mixtures can accommodate macroscopic boundary shapes.

For the jump condition, write the two affine pieces as $Ax+c_{+}$ and $Bx+c_{-}$. Continuity on the interface meansTaking differences along the interface gives $(A-B)t=0$ for every tangent vector $t$ with $t\cdot N=0$. Decompose any vector $v$ as $t+(v\cdot N)N$. ThenWith $a=(A-B)N$ this provesConversely this rank-one connection permits continuity by choosing $c_{+}-c_{-}=-ka$. A coherent twin interface is such a compatible planar boundary between different martensitic variants, possibly after a rigid rotation of one variant. The rank-one condition is therefore the geometric compatibility condition for a simple twin.

For the unrotated pair, direct subtraction givesso $R=I$ is one solution, with reference interface normal $e_1$.

Put $d=1+{\delta }^2$, $c=(1-{\delta }^2)/d$ and $s=2\delta /d$. Since $c^2+s^2=1$, the stated $\overline{R}$ is a proper rotation. Multiplication givesHere $c+\delta s=1$ and $s-\delta c=\delta$. ConsequentlyThis is the second rank-one connection, now with normal $e_3$. Because $\delta >0$, the two rotations are distinct. The permitted assertion that there are at most two proves completeness.

Finally suppose $Q{U}_1$ is rank-one connected to $\overline{R}{U}_2$, with $Q\in SO(3)$. Left multiplication of their difference by $Q^T$ preserves rank and gives $U_1-Q^T\overline{R}{U}_2$. The classification just proved forces $Q^T\overline{R}$ to be either $I$ or $\overline{R}$, hence $Q=\overline{R}$ or $Q=I$. The two requested matrices are thereforeTheir connections are explicit: the first was calculated above, and $\overline{R}{U}_1-\overline{R}{U}_2=2\delta (\overline{R}{e}_3)\otimes e_1$. If desired the second matrix is

For this homogeneous motion, suppose the matrix $\mathbf{A}=\ddot{\mathbf{F}}{\mathbf{F}}^{-1}$ is symmetric. Put $\mathbf{y}=\mathbf{x}-\mathbf{c}(t)$. The acceleration potential of a homogeneous motion iswhere $g(t)$ is arbitrary. When taking the spatial gradient, $t$, $\mathbf{c}$, $\ddot{\mathbf{c}}$, and $\mathbf{A}$ are fixed. Symmetry givesThis verifies the potential directly and also proves sufficiency of the symmetric matrix condition within the static-family setting. Its dimension is ${\mathsf{L}}^2{\mathsf{T}}^{-2}$, so the pressure correction $-\rho \xi$ has the dimension of stress. Combining with part (b)(ii), the corresponding dynamic Cauchy stress tensor iswhere $\tilde{C}(t)=C(t)-\rho g(t)$. Substituting its divergence into the Cauchy momentum equation reproduces $\rho [\mathbf{Ay}+\ddot{\mathbf{c}}]$, as required.

Take two time derivatives at fixed reference position $\mathbf{X}$:The constraint $\det \mathbf{F}=1$ ensures that the deformation gradient is invertible, so $\mathbf{X}={\mathbf{F}}^{-1}(\mathbf{x}-\mathbf{c})$. Thus the spatial material acceleration is the affine vector fieldIts spatial gradient is $\mathbf{A}$, and its curl isPart (b)(iii) therefore requiresEquivalently, multiplying on the left by ${\mathbf{F}}^T$ and on the right by $\mathbf{F}$ givesThe translation $\mathbf{c}(t)$ contributes no spatial curl and is unrestricted by this condition. A useful equivalent first integral isthe antisymmetric matrix in parentheses must be constant, and need not vanish. The separate constraint $\det \mathbf{F}=1$ is still required; differentiating it also gives $\mathrm{tr}(\dot{\mathbf{F}}{\mathbf{F}}^{-1})=0$.

The condition is understood in the setting of part (b), where a static family already exists under the prescribed body force. For a homogeneous elastic deformation gradient, the elastic stress $\mathbf{S}(\mathbf{F})$ is spatially constant. Static balance of linear momentum consequently gives $\nabla p_0=\rho \mathbf{b}$, ensuring that $\mathbf{b}$ has a scalar potential. If the static-family assumption were dropped and an arbitrary time-independent body force allowed, the general necessary condition would instead be $\nabla \times \ddot{\mathbf{x}}=\nabla \times \mathbf{b}$.

The supplied matrix obeysSeparate the even and odd powers in the matrix exponential to obtainIt fixes the $z$ component and is a counterclockwise active rotation viewed from $+z$. This sign differs from the clockwise convention in b, because the source fixes the sign of both $J_3$ and the exponential. Replacing $\theta$ by $-\theta$ yields the clockwise rotation.

With $(J_k{)}_{ij}=-i{ϵ}_{kij}$, the generators of rotations obeyThe commutator measures the leading difference between performing two infinitesimal rotations in opposite orders. For example, with $X=-i\alpha J_1$ and $Y=-i\beta J_2$,The failure to commute is itself an infinitesimal rotation about the perpendicular axis; the Lie algebra closes on the same three generators of rotations.

To pass from this local statement to the full group property, put $A_k=-i{J}_k$. These are real antisymmetric matrices, with $[A_i,A_j]={ϵ}_{ijk}{A}_k$. The Baker--Campbell--Hausdorff formula shows locally that the logarithm of a product of their exponentials remains a linear combination of the $A_k$, because all nested commutators do. Globally, each exponential has $R^{T}R=I$ and $\det R=e^{\mathrm{tr}A}=1$. If $R_1,R_2$ have these properties, thenEvery real three-dimensional orthogonal matrix of determinant $+1$ has an eigenvector with eigenvalue $+1$: its nonreal eigenvalues come in conjugate pairs, and the product of all three is $+1$. On the perpendicular plane it acts as a two-dimensional rotation. It therefore has axis-angle form and can be written $e^{-i\theta {\hat{n}}^k{J}_k}$. Thus the product is again a rotation. The identity is $\theta =0$, the inverse reverses $\theta$, and associativity comes from matrix multiplication. This proves that all these rotations form the special orthogonal group $SO(3)$. The commutator supplies local closure, while the orthogonality and axis-angle argument supplies the necessary global step.