Use the curvature sign convention implicit in part (c): the Jacobi field equation is , so the operator in a positively curved normal direction has positive eigenvalue. With the alternative convention used for the catalog's Jacobi curvature operator, the paper's tensor is , and its operator is . This fixes the sign throughout Q1.
Let . Pair interchange for the Riemann curvature tensor gives
Thus is self-adjoint with respect to the Riemannian metric. The curvature sign does not affect this conclusion. The spectral theorem for real symmetric matrices therefore supplies the real orthonormal eigenbasis used in part (b). Also , so the tangent direction is a zero-eigenvalue direction.
Write and let denote parallel transport. Since is a geodesic, , while the transported vectors satisfy . The locally symmetric Riemannian manifold condition is . Hence the Levi-Civita connection product rule yields
Both this vector field and are parallel and agree at zero. Uniqueness of parallel transport therefore proves
The eigenvalues are independent of . The original PDF has here; the converted TeX omits the dot, which would incorrectly put a point rather than a tangent vector in the operator's subscript.
Expand a Jacobi field in the parallel orthonormal eigenbasis: . Part (b) reduces its equation to the scalar equations . Define
For arbitrary initial values and , the complete solution is
A point at parameter is a conjugate point of precisely when a nonzero Jacobi field vanishes at both and . The first condition makes every . The second requires for every , with at least one nonzero . For , the factors for zero and negative eigenvalues never vanish; positive eigenvalues give exactly the sine zeros. Therefore
The corresponding points are . If several positive eigenspaces have a zero at the same parameter, the conjugate multiplicity is the sum of their dimensions. This also shows that a locally symmetric Riemannian manifold has no conjugate points along a given geodesic when that geodesic's Jacobi curvature operator has no positive eigenvalues.
Parametrize the nonconstant closed geodesic at unit speed on . Its parallel transport around the loop fixes the tangent . Since the manifold is orientable, it preserves orientation on the full tangent space, and hence has determinant one on the normal space. This normal space has odd dimension because the manifold has even dimension.
An odd-dimensional special orthogonal transformation has a fixed vector. Indeed, nonreal eigenvalues of a real orthogonal transformation occur in conjugate pairs with product one; the real eigenvalues are , and determinant one makes the number of eigenvalues even. Odd dimension therefore requires a eigenvalue. Transport a unit fixed normal vector around the geodesic to obtain a periodic parallel field , with and .
The variation is defined uniformly for small , because the image of the original loop is compact. Periodicity gives closed curves and a homotopy through loops. With the standard positive-sectional-curvature Riemannian index form, the second variation of geodesic energy is
The first variation vanishes for a closed geodesic, since the endpoint terms cancel. Therefore for sufficiently small nonzero . The Cauchy-Schwarz inequality gives
Thus the closed geodesic is homotopic to a strictly shorter closed curve. This proves the instability of a closed geodesic in positive even-dimensional curvature; completeness of the manifold is unnecessary for this compact-loop variation. Constant loops are excluded by the usual nonconstant meaning of a closed geodesic in this assertion.
Interpret as a continuous finite-dimensional complex group representation, so its values lie in , and interpret approximation in the uniform norm. Assume the usual Hausdorff convention for the compact topological group. The following argument proves the required form of the Peter-Weyl theorem without first assuming that is a Lie group.
Start with normalized Haar measure. It is both left- and right-invariant on a compact group. Average any Hermitian inner product to make a finite-dimensional representation unitary. The left regular action is already unitary on .
Choose a continuous nonnegative function supported in a small symmetric identity neighbourhood, with and . Such functions are obtained from a nonzero continuous bump near the identity, symmetrizing and normalizing. Define
The convolution operator commutes with left translations. Its continuous kernel is Hermitian, hence it is self-adjoint. It is compact on : approximate the continuous kernel uniformly by finite sums of separated functions of and , using the Stone-Weierstrass theorem on the compact product. The corresponding operators have finite rank and converge in operator norm. It also maps into and satisfies
by the Cauchy-Schwarz inequality; continuity follows from the uniform continuity of the translated kernel.
The spectral theorem for compact self-adjoint operators decomposes the orthogonal complement of its kernel into finite-dimensional eigenspaces with nonzero eigenvalues. Every such eigenspace is left-translation invariant, and all its elements are continuous, since there. If is one of these invariant spaces, write . Evaluation at the identity gives
Thus is a matrix coefficient of the dual representation. This identifies each finite spectral sum with finite-dimensional representation data.
To get uniform rather than just approximation, let project onto increasing finite sums of the nonzero eigenspaces. For , is orthogonal to the kernel, so in . Applying the displayed smoothing bound gives
Meanwhile, choosing the support of sufficiently small makes uniformly close to , by uniform continuity under small translations. Since , it also makes uniformly close to . Every belongs to a finite sum of invariant eigenspaces, hence is a finite sum of matrix coefficients. This completes the convolution proof of uniform Peter-Weyl approximation.
A coefficient equals for the rank-one endomorphism . Finite sums of coefficients can be combined using a direct-sum representation and a block-diagonal . Therefore the conclusion has exactly the requested form:
For a class function , apply conjugation averaging on a compact group,
It fixes and is a contraction in the uniform norm, so averaging an approximation preserves its error bound. A trace coefficient averages to , where
This endomorphism commutes with the representation. A finite-dimensional unitary representation splits into irreducibles by repeatedly taking invariant orthogonal complements. On each irreducible diagonal block, Schur lemma makes its average a scalar multiple of the identity, namely its trace divided by the block dimension. Off-diagonal blocks do not contribute to the total trace. Thus the averaged coefficient is a finite linear combination of irreducible characters. Consequently irreducible complex characters span a uniformly dense subspace of the continuous class functions.
The within-group and between-group scatter decomposition separates noise from group separation: measures variation about each group sample mean, and measures separation of those means about the overall sample mean. In linear discriminant analysis, seek a projection maximizing between-group relative to within-group variation,
Normalize by ; a Lagrange multiplier calculation gives the generalized eigenvalue problem . The leading canonical discriminant directions have the largest eigenvalues, and subsequent directions are orthogonal in the inner product. There are at most nonzero separating directions. Projecting observations onto these directions gives a low-dimensional display of group separation.
For classification, estimate the common covariance matrix by and use the linear discriminant analysis scores
where are class prior probabilities. Assign to the largest score. The same within-group scatter sets the noise scale, while between-group scatter selects the directions in which the group means separate most strongly.
princomp(): principal component analysis. Supply observations in rows and numeric variables in columns. After centering, let denote the resulting matrix and its sample covariance matrix, with an explicitly stated normalization. Write
The first direction maximizes the Rayleigh quotient over unit vectors; successive directions maximize it subject to orthogonality to the earlier directions. The principal component scores are , with variances and zero pairwise sample covariance. Zero covariance alone does not imply independence. Each principal component loading describes a direction in variable space; its sign is arbitrary, and repeated eigenvalues leave the basis inside the corresponding eigenspace arbitrary.
The function returns principal component loadings, component standard deviations sdev and, when requested, principal component scores. By default it uses a covariance matrix; cor=T instead performs principal component analysis on a correlation matrix, standardizing nonconstant variables to variance one. This choice matters when variables have different measurement units. These interface details are documented in the S-Plus help for princomp.
A scree plot shows the ordered eigenvalues; a plot of the first two principal component scores can reveal groups, unusual observations or a dominant gradient. The first components explain of the total variance. Reconstructing each centered observation by its orthogonal projection onto these directions minimizes average squared reconstruction error among -dimensional linear subspaces. Thus PCA gives an optimal linear compression for squared reconstruction error, with scaling and the retained dimension chosen explicitly. It is an unsupervised learning method, rather than a classification rule fitted to supplied labels.
tree(): classification trees and regression trees. A call such as tree(response ~ predictors, data=dat) fits a recursive partitioning model. At a current node, test numerical predictors against candidate thresholds, or divide the levels of a categorical predictor into two sets. Select a split that most improves the fitting criterion, then repeat within its children. Numerical splits give axis-aligned regions. A categorical response produces a classification tree; a numeric response produces a regression tree, as specified in the S-Plus tree documentation.
For a regression tree, the constant prediction in a terminal region is its sample mean, minimizing the sum of squared residuals there. For a classification tree, estimate class probabilities by the observed proportions in a leaf and predict a most probable class under equal misclassification costs. The classification-tree deviance of a node is
with ; splitting maximizes the decrease from parent to children. This fitting criterion differs from the number of wrong classifications. The middle sketches show the same three-region partition as a binary tree: first split on variable one, then split the right-hand branch on variable two.
Growing to small leaves can overfit. Use cost-complexity tree pruning, balancing fitted loss against leaf count,
and select complexity by cross-validation. Plotting the fitted tree and annotating nodes displays the rules; prediction routes each new observation through the fitted splits. A tree gives an interpretable, piecewise-constant prediction rule; training fit alone does not measure prediction accuracy.
The rpart() alternative. This also fits classification and regression trees by recursive partitioning, with squared-error splitting for ordinary numeric responses and, by default for classification, Gini impurity . Classification settings can instead use information-based splitting and specify prior probabilities and a loss matrix. It supports surrogate splits to approximate a node's primary left/right assignment when the primary predictor is missing. Its complexity-parameter path and cross-validated errors support cost-complexity tree pruning; for example, select a complexity parameter and apply prune(fit, cp=...). The maintainer's rpart documentation describes the criteria and the distinction in handling surrogate predictors. The partition and tree sketches apply to either implementation; the fitted tree need not be the same because splitting criteria and controls can differ.
cmdscale(): classical multidimensional scaling. Its input is a dissimilarity matrix or stored distance object. It outputs coordinates for the observations in a chosen dimension, usually two for visualization. Let contain the squared entries of , and put
For centered Euclidean coordinates , the distance identity gives : multiplying by removes the row and column norm terms. Diagonalize this Gram matrix and retain its largest positive eigenvalues. The coordinates are
Thus classical scaling exactly recovers Euclidean distances when all positive eigenvalues are retained. Reducing dimension approximates the centered Gram matrix; it does not in general minimize squared discrepancies in the original distances. Negative eigenvalues identify a failure of the Euclidean distance matrix criterion, so their size matters when interpreting a display. The historical S-Plus cmdscale help describes the dimension argument, eigenvalue output and optional additive adjustment to dissimilarities.
Plot the two coordinate columns on equally scaled axes and label observations. Nearby plotted points suggest small dissimilarities only insofar as the dimension reduction fits the input well. The bottom sketches give an exact Euclidean example: a configuration and its independently recovered two-dimensional coordinates have identical pairwise distances, despite different coordinate axes. Classical scaling converts distances into geometry; translations, rotations and reflections do not change the represented distances. For Euclidean distances between observations, its coordinates agree with the corresponding principal component scores up to these harmless choices of axes.
Figure 1.
Illustrative PCA directions and scree plot, equivalent recursive partition and tree, and exact classical-scaling recovery of Euclidean distances
. Details.
Use a mostly-plus target metric tensor and the standard open-string mode expansion convention . The integer string oscillator level is , whose eigenvalues are nonnegative integers. Then , so the string mass-shell condition gives
The ground state is a tachyon, and the level-one state is massless. The physical string state conditions in old covariant string quantization are
Only the positive-mode Virasoro constraints annihilate a ket; their negative-mode partners annihilate the corresponding bra. Imposing both signs on every ket would contradict the nonzero Virasoro central extension.
The central charge comes from normal ordering the quadratic oscillator expressions. Classically their Poisson brackets have no central term. Quantum mechanically, moving annihilation operators through creation operators produces double contractions. From the oscillator commutator one obtains , which gives the noncentral part .
The coefficient is fixed by the free-boson Virasoro central term. On a formal zero-momentum vacuum, for ,
The two contraction pairings give
Here ; the time coordinate contributes one, just as each spatial coordinate does. Comparing with the Virasoro algebra proves
This counts covariant target coordinates. It is distinct from the physical transverse polarizations. Including the central charge of reparameterization ghosts gives ; cancellation of the worldsheet Weyl anomaly gives the usual critical dimension of the bosonic string .
A general level-one state is
The zero-mode constraint requires . Since , the remaining nontrivial positive-mode constraint is ; all with act trivially at this level. Thus the polarization vector is transverse to a nonzero null momentum.
There is also a null string state . It is physical when and orthogonal to every physical polarization. The null-state quotient of a string therefore identifies
Choose . Transversality sets , and the null-state identification removes this common component. The remaining spatial components have positive norm. Hence the physical level-one states form a massless vector with polarizations, or 24 at the critical dimension. The removed component is gauge redundancy. States violating the constraints, including independent timelike negative-norm polarizations, are unphysical; physical null states are removed by the quotient rather than counted as additional particles.
Start with the Clifford algebra . For the dimension count, complexify and rescale the timelike matrix so that the generators obey the Euclidean complex relation . This changes neither the complex representation dimension nor its decomposition under the complexified even algebra. In , pair the generators into
Direct substitution gives and . The commuting operators satisfy , so their occupation eigenvalues are zero or one. An irreducible module contains an all-zero occupation state , obtained by lowering the occupied entries of a simultaneous eigenvector. The states
are nonzero and independent by their distinct joint occupation eigenvalues. They span an invariant module; irreducibility makes this the full module. Conversely this Fock construction supplies a representation on exactly these states. Products of occupation projectors and raising/lowering operators give its matrix units, so the generated algebra is the full matrix algebra. Therefore the minimal complex Dirac spinor has
This module is irreducible under the full Clifford algebra, but the Spinor representation of the Lorentz group is generated by the even bilinears . The chirality matrix in the Euclideanized basis is
Its two eigenspaces therefore remain invariant under Lorentz transformations. Each has half the Fock states, since any gamma matrix bijects them. The even algebra connects every pair of occupation states of the same parity and has no further invariant subspaces within either one. Hence
These are the irreducible Weyl spinors. Globally they are ordinary representations of , the double cover; calling them representations of uses the conventional projective spinor terminology.
For , adjoin . It anticommutes with the old matrices and gives a Clifford module of the same dimension . The new even generators interchange the two old chirality eigenspaces, so the odd-dimensional Lorentz spinor is irreducible of complex dimension , with no Weyl splitting. The product of all odd-dimensional gamma matrices is central and acts as a scalar; it cannot be a new anticommuting chirality operator. Its two signs distinguish the two complex odd-Clifford modules, whose restrictions to Spin group are equivalent. Majorana spinor reality conditions, when the signature permits them, are additional conditions and should not be confused with complex chirality.
For the reduction, take the massless zero modes on a flat seven-torus with no flux or projections. The eleven-dimensional fields are the metric , a three-form and a Majorana gravitino . Split with and . The metric supplies , seven vectors and 28 scalar fields . The three-form supplies a four-dimensional three-form , seven two-forms , vectors and scalars .
A massless two-form in four dimensions has one physical particle polarization and can be dualized to a scalar by ; the seven two-forms thus supply seven more scalars. The four-dimensional three-form has no local particle polarization, since its four-form field strength is locally a constant times the volume form. The total bosonic content is consequently
The eleven-dimensional Majorana spinor has 32 real components, decomposing into eight four-dimensional Majorana spinors. Accordingly supplies eight gravitini, while supplies spin- fields after diagonalizing the gravitino/matter kinetic terms. The resulting theory has 32 real supercharges, hence is four-dimensional N=8 supergravity, with
The on-shell counts check the reduction: the metric has polarizations and the three-form has , totaling 128. In four dimensions bosonic polarizations match fermionic ones. Higher Kaluza–Klein levels, fluxes or projections would be additional choices beyond this massless toroidal reduction.
Interpret triviality as holomorphic triviality. Choose a global holomorphic bundle frame for . Let have as its columns the coordinates of these sections in the prescribed local frame over , with the convention that a fiber coordinate vector transforms as . Then
Each is holomorphic and invertible because its columns form a basis in every fiber. This proves the requested holomorphic splitting of transition functions on a trivial bundle. Smooth or topological triviality alone would not justify holomorphic splitting.
For the line bundle on projective twistor space, pull its patching function back to the twistor line using the incidence convention of this problem, , without an extra factor of . Triviality on each line gives nonvanishing linewise holomorphic functions with . The factors can be chosen as holomorphic functions of locally: a degree-zero scalar transition function has zero winding on the overlap annulus, hence admits a logarithm there; splitting its positive and negative Laurent series parts gives factors whose coefficients are Cauchy integrals that are holomorphic functions of . The explicit splitting is worked out below. Their only common linewise ambiguity is multiplication by a nonzero function : the ratio is a global holomorphic function on the compact complex projective line, hence constant on that line.
Take , and raise a two-component spinor by left multiplication with . The derivative along an alpha-plane is . It annihilates the unsplit twistor function, because
Thus
is holomorphic on the whole line and has projective weight one. A section of is linear in the homogeneous coordinates: in an affine coordinate, holomorphic regularity at infinity bounds its growth to at most first order, so the Cauchy integral formula excludes all higher powers. Therefore
The inverse splitting functions obey . Commuting these two equations proves
The spinor decomposition of gauge curvature is
where both two-component spinors are symmetric. Vanishing of the contraction for every forces . This is exactly the anti-self-dual Maxwell equations. Since the Abelian gauge curvature is , it is closed; its duality condition also makes it co-closed, giving the source-free Maxwell equations. Changing both splitting factors by changes by , the expected Abelian gauge transformation. This derives the rank-one Penrose-Ward correspondence in the needed local form.
For the exponential patching function, choose a representative of the Čech cohomology class, of homogeneous degree zero. On the overlap of the two line patches write . To compute the gauge potential explicitly, use a constant two-component spinor basis in which and the finite patch . Then , and
For a positively oriented contour in the overlap annulus, split its Laurent series as
The first is holomorphic inside the contour, the second outside including infinity; split . Since , coefficient comparison gives . Applying that operator separately to cancels all nonconstant terms and leaves
Thus yields and . Moreover
Substitution gives , which in homogeneous form is precisely
This is the Abelian Ward splitting by Laurent series. The factors have weights , respectively, so the integrand is invariant under a change of homogeneous representative. The contour avoids ; the basis chosen above places that point at infinity. One can keep the contour fixed for in a sufficiently small domain, allowing differentiation under the integral.
Finally differentiate the contour potential for an anti-self-dual Maxwell field. By the twistor incidence relation, . Consequently
The second equality uses the elementary two-spinor identity in the declared convention. The remaining two-component spinor is symmetric in , and no primed symmetric gauge curvature part survives. This directly verifies is anti-self-dual. In complexified Lorentzian signature duality has Hodge star eigenvalues ; here the vanishing primed two-component spinor fixes the stated ASD convention rather than incorrectly imposing the Euclidean equation on that signature.