Fix a bounded tile in Euclidean space with labelled flat boundary faces. For each face label , encode an assembly of copies by an involutive real symmetric signed permutation matrix : a paired interface contributes a two-by-two swap, an unglued face with a Neumann boundary condition contributes on the diagonal, and an unglued face with a Dirichlet boundary condition contributes . Use the same face-identification convention on each copy. If a second assembly has matrices and there is an invertible constant matrix satisfying
then the map is an isomorphism between corresponding Laplace-Beltrami operator eigenspaces. This is the transplantation theorem.
To prove it, restrict an eigenfunction to the tile copies and identify each copy with . Along a face, let be its vector of boundary values and its vector of outward normal derivatives. At a glued interface the values agree and the outward derivatives are opposite. At an unglued face with a Neumann boundary condition the derivative vanishes, and at a face with a Dirichlet boundary condition the value vanishes. All three cases are encoded by
The intertwining identities give and , so the transplanted function satisfies the second assembly's interface and boundary conditions. On each tile its Laplace-Beltrami operator is the same constant linear combination of the original Laplacians, hence it has the same eigenvalue. Across an interface, matching values and opposite normal derivatives make it a weak solution; elliptic regularity, or the allowed reflection principle, makes it a genuine smooth solution there. At corners the operator is understood by its finite-energy, Friedrichs extension; the same matching argument applies to its quadratic form domain. Applying reverses the construction. Thus dimensions of all eigenspaces agree, including multiplicities, and the two specified boundary-value spectra coincide.
The requested pair exists with explicitly specified mixed boundary conditions and polygonal boundary. Take four copies of an octagonal regular polygon in Euclidean space. Label four alternate sides in cyclic order; the four intervening sides always have a Neumann boundary condition. The labelled sides are disjoint even at their endpoints. Parameterize every side in the positive boundary direction of the reference octagon. When pairing same-labelled sides, identify equal parameters. This produces a locally flat surface with boundary: a neighbourhood of an interior point of an interface consists of two half-disks glued along their diameters. All tile vertices remain boundary points, with boundary angles or ; there are no interior cone points.
Here are the complete gluing and boundary conditions. An entry pairs the corresponding sides of . The indicated unpaired sides have a Dirichlet boundary condition; every other unpaired side has a Neumann boundary condition.
Both assemblies are compact and connected. The Laplace-Beltrami operator is the nonnegative operator associated with the finite-energy quadratic form , with zero Sobolev trace on the specified Dirichlet arcs. Thus the mixed boundary condition is defined without requiring classical differentiability at a corner or at a Dirichlet–Neumann junction.
The four signed permutation matrices for are
For they are
Use the Hadamard matrix
The intertwining relations can be checked on all columns at once. Writing the columns of as , direct application of the displayed gives
For the four intervening faces both assemblies have the identity face matrix, so these intertwining relations hold as well. Also . This four-tile mixed-boundary transplantation between a disk and a nonorientable surface gives a bijection of every pair of eigenspaces, preserving eigenvalues and multiplicities.
There is also a direct argument on the quadratic form domains that handles all boundary corners. The orthogonal matrix acts on tile restrictions, preserves their L2 inner product and the sum of their Dirichlet energies, and commutes with taking a Sobolev trace. The condition on each face encodes equality of the two traces at a paired interface and vanishing trace at an unpaired Dirichlet face. Intertwining takes these conditions exactly to . Consequently maps the whole form domain onto the other form domain and preserves its quadratic form. The associated self-adjoint operators are unitarily equivalent, proving isospectrality for the stated mixed boundary conditions without any assumption of smoothness at a Dirichlet–Neumann junction.
We now verify the actual topology using the topology of disk tiles glued along disjoint boundary arcs. Because the gluing arcs in each tile are disjoint, the resulting surface with boundary is a thickening of its gluing graph, with one vertex for each closed disc tile and one edge for each paired arc. One can construct this homotopy equivalence by taking a small central closed disc in each tile and narrow strips to its paired arcs. What remains is a collar neighbourhood of the boundary, which retracts onto this connected system of closed discs and strips. Shrinking each central closed disc and each strip then gives the gluing graph. Equivalently, begin with a spanning tree of the gluing graph: its closed disc tiles joined along boundary arcs form one closed disc, and every remaining pair of boundary arcs adds one band. Thus for tiles and paired arcs the fundamental group is a free group of rank , and the Euler characteristic is .
For , the gluing graph is the path
Gluing a new closed disc to a closed disc along one boundary arc again gives a closed disc; applying this three times proves that is homeomorphic to a closed disc. In particular
For , the five paired arcs give the graph edges . It is connected, and hence
Here is the free group on two generators, so is not simply connected. To decide orientation of a surface, give every reference tile its counterclockwise orientation. Our same-parameter gluing preserves boundary directions, so a global orientation requires opposite orientation signs on paired tiles: the two induced boundary orientations at an interior interface must be opposite. Such signs exist exactly when the gluing graph is bipartite. In the cycle
has length three. Following it forces the sign to reverse three times and return to its starting tile with the opposite sign. This proves that is a nonorientable surface. It has nonempty boundary because every intervening octagon side remains unpaired. Thus is a simply-connected orientable surface with boundary, whereas its isospectral partner is neither simply connected nor orientable, exactly as required, for the specified mixed boundary problem.
The boundary qualification is essential. With smooth boundary and the same pure Dirichlet boundary condition or pure Neumann boundary condition on both surfaces, the following existing obstruction applies. Polygonal corners and the mixed boundary conditions in the construction above place it outside that obstruction.
With the interpretation of compact smooth Riemannian surfaces with the same pure Dirichlet boundary condition or pure Neumann boundary condition, the printed simply-connected/nonorientable pair cannot exist. The two-dimensional heat trace has the following expansion (the smooth-boundary coefficients are given in the heat-coefficient formulas, Theorems 5.1–5.2):
where the boundary sign is minus for Dirichlet and plus for Neumann. The Gauss-Bonnet theorem identifies the constant term with . Isospectrality therefore forces equal Euler characteristics. A compact connected simply-connected surface with boundary is a disk and has . A connected nonorientable surface with crosscaps and boundary components has . Thus no such smooth pure-Dirichlet or pure-Neumann pair exists. This is the Euler characteristic from the heat trace of a bordered surface.
Piecewise-smooth boundaries introduce corner terms, and mixed boundary conditions change the heat coefficients, as the explicit construction above demonstrates. Moreover, in the usual pure-Neumann polygonal reflection-transplantation setup, the pure Neumann reflection transplantation preserves Euler characteristic argument applies even with corners. The numbers of tiles agree. For each face label the number of edge classes is , so these agree too. For each tile-vertex label, vertex classes are orbits of the subgroup generated by incident-face involutions. Their number is the dimension of invariant vectors in that permutation representation, which is unchanged by the intertwiner. Hence the numbers of vertices agree, proving equal . Thus this standard tiled construction also cannot supply a disk on one side and a nonorientable surface on the other.
For comparison, the following existing construction gives an orientable and a nonorientable pair with the same spectrum for the pure Neumann boundary condition. Both have nontrivial fundamental groups; this is a separate example, not the construction used to meet the simply-connected requirement above. This eight-tile Neumann transplantation across orientability makes the comparison explicit. The finite seed is the affine modulo-eight example in the original orientability construction; the intertwiner and the topology checks are derived below. Take eight copies of a regular Euclidean hexagon and label three pairwise nonadjacent sides . Every remaining side is a side with a Neumann boundary condition. Index the copies by . For the first assembly pair same-labelled faces according to
for the second, use
with all arithmetic modulo eight. A fixed point means an unglued side with a Neumann boundary condition. Pair faces by the reflection convention, matching the same endpoint parameter in the reference tiles. The distinguished sides are disjoint, so no interior vertex singularities are created; the resulting surfaces have polygonal boundary corners.
Let on the eight-dimensional tile-index space and set
It commutes with the permutation matrices for and , because its coefficients are invariant under reversing cyclic order. For , write , where , , and all other vanish. The required relation is exactly , verified from these four nonzero entries. Thus all three intertwining relations hold.
The cyclic Fourier vectors diagonalize . Its eigenvalues are
namely . None is zero, so is invertible and the two assemblies are isospectral manifolds for the Neumann boundary condition.
We prove their topology as well. Their gluing graphs are connected: the first contains the spanning path , and the second contains the spanning cycle . Each has eight tile vertices and ten paired-face edges. Since tiles are disks and gluing arcs are disjoint, each surface retracts up to homotopy to this multigraph; thus its Euler characteristic is and its fundamental group is a free group of rank three.
An orientation assigns a sign to each reference tile, and a reflection gluing requires opposite signs at the ends of every graph edge. The second graph admits this assignment: take sign on and sign on . Every paired face joins opposite signs, proving orientability. In the first graph the cycle
has odd length. Transporting orientation around it reverses the initial sign, proving nonorientability. Hence the first surface is also not simply connected. The second surface is orientable but likewise not simply connected; the disk required in the question is supplied instead by the four-tile mixed-boundary construction above.
The four-tile construction resolves the printed simply-connected requirement under explicit mixed boundary conditions. The eight-tile comparison only concerns orientation of a surface under a pure Neumann boundary condition.
We use the Dirichlet density form of the Chebotarev density theorem, the version directly proved by the permitted information at . If is a finite Galois extension of number fields with group and is a conjugacy class, let be the unramified prime ideals of whose Frobenius conjugacy class is . Then
The same conclusion holds for a union of conjugacy classes, replacing by the size of that union. The denominator may equivalently be the sum over all primes of . Here and below is real. This specifies the density being proved; an unweighted prime-counting asymptotic is a stronger formulation and is not being inferred merely from a limit at one.
First establish the necessary pole for a general base field. The Dedekind zeta function as a permutation Artin L-function gives it using exactly the allowed assumptions. Take a finite normal closure of , and let be the permutation representation on . Its trivial constituent has multiplicity one: an invariant function on the transitive set of embeddings is constant. At a rational prime, a basis of consists of the sums over inertia orbits of embeddings. Frobenius automorphism permutes these orbit sums. Within the decomposition-group orbit belonging to a prime of , the resulting cycle has length , because the residue embeddings are permuted cyclically by finite-field Frobenius. A cycle of length has determinant . Consequently the local factor of an Artin L-function gives
including at ramified primes. Decomposing gives
All factors after are holomorphic and nonzero at one by the allowed hypothesis. Thus has a simple pole there and
This argument also explains why the trivial relative Artin L-function is , rather than always the Riemann zeta function.
Now let be an irreducible character of , of degree . At an unramified prime the Euler product has local factor , where is a Frobenius representative. The eigenvalues of are roots of unity, since has finite order. Expanding the logarithm for therefore gives
where the finitely many ramified local factors contribute a function bounded as . They are nonzero and finite near one because their eigenvalues have modulus one while .
The terms are uniformly bounded. Indeed , and for and ,
There are at most primes of above each rational prime and each has norm at least that prime. Hence . Isolating proves the character prime sums in the Chebotarev proof identity
For nontrivial , the allowed holomorphy and nonvanishing at one give a local holomorphic logarithm and a bounded right side. The logarithm specified by the Euler series differs from that local branch by a constant multiple of on the real interval near one, so it too is bounded. For the trivial character, the right side is by the pole of . Thus
Expand the conjugacy-class indicator by character orthogonality. For ,
To see the coefficient, take its inner product with : . The irreducible characters form an orthonormal basis of class functions, so this is the full expansion. Substituting the prime-sum estimates leaves only the trivial character's logarithmic term:
Division by the diverging logarithm proves the theorem. The same trivial-character calculation shows that all primes of have density one; adding or removing finitely many primes does not change a Dirichlet density.
There is a necessary qualification for the final polynomial application: the printed claim is false for degree one. For example, is monic and irreducible in , yet has the root zero modulo every prime. We prove the intended assertion for degree .
Let be the splitting field of the monic irreducible polynomial over , and let be its roots. The Galois group acts transitively on . There is a derangement in a transitive group action: count the pairs satisfying . The orbit-stabilizer theorem gives
If every fixed at least one root, the identity, which fixes , would make this sum at least , a contradiction. Therefore some fixes no root, and neither does any conjugate of .
For completeness, relate this group action to reduction without assuming the desired splitting result. All roots are algebraic integers, since the polynomial is monic. Exclude the finitely many primes ramifying in or dividing the nonzero polynomial discriminant. At a prime over any remaining , the roots reduce to distinct elements of and give all roots of the reduced polynomial. Arithmetic Frobenius permutes them by . A reduced root belongs to exactly when it is fixed by this map. Because reduction is injective on the original root set, this happens exactly when the Frobenius permutation fixes an original root. This proves the Frobenius permutation and roots modulo a prime criterion.
The Chebotarev density theorem assigns the conjugacy class of our fixed-point-free the positive Dirichlet density . All but finitely many of these primes have rootless reduction, so there are infinitely many. More precisely, applying the same argument to the union of all fixed-point-free classes gives
In a linear regression model, a large number of predictors or nearly dependent columns of the design matrix can make some eigenvalues of small. The ordinary least squares estimate has covariance matrix , so these directions have large variance, even though the estimator is unbiased. Ridge regression reduces this variance by shrinking the estimate towards zero, at the cost of introducing bias. This bias-variance tradeoff can improve the mean squared error when the signal is not too large in the shrunken directions.
For the usual preprocessing with a regression intercept, separate the constant column from the predictors and leave the intercept unpenalized. Subtract the sample mean from the response and from each nonconstant predictor, then divide each centered predictor column by its positive scale, for example its Euclidean norm. Thus all predictor columns have squared norm one; using squared norm instead simply changes the convention for the penalty. Centering is performed on the slope columns, not on the constant column. A full-rank design containing an intercept retains full rank on its centered slope columns. A model without an intercept should retain its specified mean structure unless an intercept is deliberately introduced.
Let be the resulting scaled, centered slope design matrix, let be the centered response, and let denote the coefficients in these standardized coordinates. The ridge regression problem is
Here is the number of slopes; it is one less than the original number of columns when a constant column was present. The gradient of this objective is . Its Hessian matrix is a positive-definite matrix, giving the unique closed-form ridge regression estimator
The regression intercept in the original coordinates is recovered as , and the original slopes are , where is the scale used for column . This is the unpenalized intercept in ridge regression convention. For a fixed full-rank coefficient design written simply as , with response and coefficient , the same formula reads
The risk comparison below first uses the fixed coefficient coordinates in which this quadratic penalty is imposed, and then addresses reversal of the scaling.
Assume at least one penalized coefficient, zero-mean errors with covariance matrix , and . No normal distribution assumption is needed for this mean squared error calculation. Put and . The covariance and bias of a ridge regression estimator follow directly from its linear expression:
For the centered slope model, the response errors have covariance matrix , where . Since , the same coefficient covariance matrix formula still holds with .
By the spectral theorem for real symmetric matrices, write , with orthogonal and every , and put . The bias-variance decomposition of mean squared error now gives
The expression at is the ordinary least squares mean squared error. Its right derivative is
Therefore, for every fixed coefficient vector, all sufficiently small positive penalties strictly improve mean squared error.
An explicit sufficient choice follows by subtracting the ordinary least squares risk term by term:
Since , every summand is negative whenever
When , every positive penalty works. This is a pointwise result: the permitted penalty may depend on the true coefficient vector. A fixed positive penalty cannot improve risk for all arbitrarily large signals, because its squared bias grows quadratically with the signal. The same distinction appears in uniform directional risk improvement by ridge regression and unbounded directional risk of fixed ridge shrinkage. In centered coordinates the unpenalized intercept contributes the same to both coefficient risks, so it does not affect the improvement. Reversing a fixed predictor scaling measures coefficient error with a fixed positive-definite matrix instead of the identity matrix. For this weighted loss the squared bias still has derivative zero at , while the variance term has derivative
The inequality follows because . Thus sufficiently small penalties also improve the original-coordinate mean squared error, though the explicit sufficient bound above was for the unweighted coordinates. Recovering an uncentered intercept adds a fixed nonnegative quadratic form to the slope loss, and its sample-mean error is uncorrelated with the centered slope errors, so the same weighted argument applies. A model with only an unpenalized intercept gives identical estimators. If the noise variance were zero, the strict improvement assertion would fail; positive noise variance is essential.
A Lie algebra representation is a linear map preserving brackets: . Differentiating a smooth unitary representation of a Lie group gives , an anti-Hermitian operator for each real algebra element . In the physics convention one uses Hermitian generators and writes . Conversely, integration of a Lie-algebra representation gives a representation of the connected simply connected group. A representation of another group with the same algebra additionally has to be trivial on the kernel of its covering map.
For the angular momentum commutation relations, take Hermitian and on a positive-definite inner-product space. A highest-weight vector has and . Set . Commuting through each lowering operator gives . Furthermore,
so acting on the highest-weight state and summing the coefficients gives
Taking the inner product with produces the unitary highest-weight termination recursion
The first factor requires . If is not an integer, the first integer would make the norm negative, while every preceding norm is positive. Thus
For these values, have strictly positive norms and distinct eigenvalues, while has zero norm and hence is zero. They therefore span a space of dimension , invariant under all three generators. Its weights are . This proves integer or half-integer spin and finite dimension . The positive-definite unitary assumption matters: an unrestricted algebraic highest-weight module need not terminate. Normalizing the ladder states yields
This is the normalized highest-weight lowering formula with the usual positive ladder phases and units .
In the ordered spin-one-half basis , the matrices are
Thus , where the are Pauli matrices.
Choose the active-rotation convention. The SU(2) lift of an axis-angle rotation is
The Pauli matrix multiplication law shows that conjugating by this matrix produces , with
This is the ordinary three-dimensional rotation. In any spin- representation the corresponding unitary operator and its matrix elements are
For half-integer spin, in this notation includes a choice of lift to SU(2); it does not define a single-valued representation of SO(3).
For the specified Euler-angle order, the rightmost rotation acts first:
Here the Wigner D-matrix middle factor is . It can be calculated explicitly by realizing spin as the symmetric product of spin-one-half factors. Write and . The fundamental rotation sends to and to . Expanding a normalized symmetric state with plus signs gives
where . The exponent counts initially positive spinors that become negative, while counts initially negative ones that become positive; their minus signs and binomial coefficients give the displayed expression. Every factorial argument is an integer even for half-integer spin.
For spin one-half this reduces to the explicit matrix
Replacing by reverses both sine and cosine, so every entry changes sign. In particular, a pure rotation gives , while a rotation gives .
Finally, the Adjoint double cover from SU(2) to SO(3) is surjective by the axis-angle construction. Its kernel consists of matrices commuting with all Pauli matrices, hence scalar matrices; within SU(2) these are . Thus
Their Lie algebras are isomorphic, but their global topology and representations differ. In spin , the central element acts as , so exactly the integer-spin representations descend to ordinary SO(3) representations. This is the descent of an SU(2) representation to SO(3) criterion.
Let be the defining SU(3) representation. A mixed SU(3) tensor representation belongs to , with transformation law
The unrestricted tensor has dimension . The upper and lower counts are independent; setting them equal would lose part of the printed generality.
For an irreducible tensor, impose the appropriate Young symmetrizer on each set of indices and subtract all upper-lower contractions. Antisymmetric pairs can be converted using the invariant alternating tensor, since and . A full alternating triple is a singlet because . Traces use the invariant Kronecker delta and yield lower-rank tensors. Repeating these operations gives invariant irreducible subspaces, with multiplicities retained.
An explicit algorithm for the complete decomposition starts at , tensors with exactly times, and then with exactly times. At each step apply
omitting terms with negative labels and adding multiplicities of coincident terms. The first is the sl3 highest-weight tensor rule: a Young diagram of row lengths gains a box in one of its three rows. The third-row term loses a determinant column of height three, leaving labels . The second rule follows by conjugation, which exchanges and . Thus the procedure determines all constituents of the general tensor, not just one selected symmetry type.
The separately symmetric, traceless component has labels . Before imposing traces its dimension is . The contraction map has target the separately symmetric tensors of degrees and is onto: in polynomial coordinates it is , the adjoint of multiplication by in the monomial-factorial inner product. Multiplication is injective, so the adjoint is surjective. Hence the symmetric traceless SU(3) tensor representation has dimension
For or the subtracted term is zero. The highest tensor is trace-free and has highest weight . The three positive-root factors in the Weyl dimension formula are , and , giving exactly the kernel dimension above. Thus the whole symmetric traceless kernel is that irreducible constituent.
There is a genuine distinction between this highest-weight constituent and the literal largest-dimensional constituent. For example, contains , obtained by antisymmetrizing one pair, with dimension , whereas the fully symmetric has dimension . Thus cannot be the answer to “largest dimension” for every unrestricted pair of counts.
The maximal dimension in mixed SU(3) tensor powers can nevertheless be determined exactly. Every constituent highest weight differs from by with nonnegative integers , where and . Its labels therefore satisfy
If both are positive, subtracting from both increases the labels to , whose dimension is strictly larger. The resulting boundary labels occur in the tensor: use disjoint antisymmetric upper pairs or disjoint antisymmetric lower pairs, and then take the separately symmetric traceless component. Consequently the exact answer is the finite formula
This supplies a practical closed finite computation for arbitrary . To reduce it to at most four evaluations, put , and . The first continuous maximum is at , clipped to ; test its two neighbouring allowed integers. The second is similarly at , clipped to . Indeed differentiating the first cubic gives a derivative with the sign of , which changes sign once at that root. For equal upper and lower counts, both maxima are at zero and the familiar result is .
For three quarks, the unconstrained flavour product is
The brackets indicate permutation symmetry: symmetric, mixed, and alternating. The three-quark colour singlet is proportional to and is alternating. The fermion wavefunction must be alternating overall. In the usual spatially symmetric ground-state S-wave, the spin-flavour part must therefore be symmetric.
The three-spin-one-half space has a symmetric spin- sector and mixed-symmetry spin- sectors; it has no fully alternating sector because . Symmetric flavour can pair with symmetric spin, giving the baryon decuplet of spin . The flavour-octet multiplicity space and the spin- multiplicity space both transform as the two-dimensional standard representation of the symmetric group. Their product contains exactly one symmetric combination: its characters at the identity, a transposition and a three-cycle are , so the symmetric multiplicity is . This gives one baryon octet of spin . Alternating singlet flavour would need alternating spin to make the spin-flavour part symmetric, and that sector does not exist. Therefore the Pauli constraint on three-quark flavour multiplets yields
The spin-flavour dimension check is , the symmetric spin-flavour SU6 representation.
The spatial-symmetry qualification matters if the wording is interpreted as total alone. Rotationally scalar excited orbital functions can have mixed permutation symmetry. For example, the two independent differences among , multiplied by a symmetric radial factor, span such a space while remaining rotational scalars. Combining this standard permutation representation with the spin- standard representation gives an alternating spin-orbital combination; multiplying by alternating singlet flavour makes it symmetric, and the colour factor restores overall fermionic antisymmetry. Thus total without a symmetric ground-state orbital assumption would not by itself exclude a singlet.
For a quark-antiquark pair, colour singletness selects the contraction , and the independent flavour product gives
Quark and antiquark are distinguishable, so there is no three-identical-quark exclusion of either flavour multiplet or either spin. Both the octet and singlet admit spin singlets and spin triplets. Their baryon number is zero. The flavour octet has rows , while the singlet has . The octet's zero-isospin state is , and the singlet is . Each row contains all , and electric charge is .
For ground-state mesons with , opposite intrinsic quark-antiquark parities give . The neutral self-conjugate states have charge conjugation , so both flavour multiplets occur as pseudoscalar mesons and vector mesons, forming a flavour nonet in each channel. Charged and open-flavour states have no individual charge-conjugation eigenvalue. The meson request does not restrict orbital excitation: in general , , and the same formula for applies where defined. Thus the flavour representations stay , with these additional orbital-spin quantum numbers.
Use logarithms to base two, so all Von Neumann entropies and Shannon entropies are measured in bits. Write the nonzero eigenvalues of as . Orthogonality of the supports makes the mixture a block diagonal matrix, with eigenvalues . Therefore
Each has trace one. Terms with zero probability or zero eigenvalue vanish under the convention . This proves the entropy of an orthogonal quantum mixture.
Introduce the orthonormal local vectors and . Reexpressing the second state's factors in this basis gives its Schmidt decomposition
The squared Schmidt coefficients are therefore , whereas those of the first state are .
To see directly why these distinguish the states, a local unitary operation changes the reduced density matrix to
This is local-unitary invariance of reduced-state spectra: the unitary on the traced subsystem cancels, and conjugation on leaves its eigenvalues unchanged. Here the reductions are
Their spectra differ. Hence no pair of local unitary operations carries the first state to the second. More generally, normalized bipartite pure states are locally unitarily equivalent exactly when their squared Schmidt coefficients agree, including multiplicities; sufficiency follows by mapping one pair of Schmidt bases to the other.
An upper bound follows directly from the Schmidt-tail entanglement monotone
where is the reduced density matrix of a bipartite pure state. On these two-qubit states, this is the smaller squared Schmidt coefficient, so and .
Here is the average-monotonicity argument needed for all LOCC strategies. The largest eigenvalue is convex: for any unit vector , , and maximizing the left side gives the claim. Hence is concave. If Alice makes a local measurement, Bob's conditional reductions average to his original reduction. Concavity therefore gives . If Bob measures, use Alice's reductions instead. The two reductions of a pure state have the same nonzero eigenvalues by the Schmidt decomposition. Refining all measurement outcomes into individual Kraus operators keeps conditional states pure; iterating the inequality covers adaptive rounds of classical communication. The same bound persists for limits of such protocols.
If success produces the target with probability , its contribution to the final average is . Failure branches have nonnegative . Thus
This is a bound on arbitrary local protocols, not merely on a particular filtering measurement.
Interpret as the minimally coupled scalar positive Laplace-Beltrami operator on the smooth closed manifold. Writing for its metric volume measure (the Riemannian volume form on an oriented chart), integration by parts gives . The compact elliptic spectral theorem supplies an orthonormal basis with , real nonnegative eigenvalues tending to infinity and finite multiplicities. Solve this eigenvalue problem to obtain the spectral zeta function
The zero eigenvalues are excluded, not assigned inverse powers. Their multiplicity is the number of connected components: implies , so is constant on each component.
An equivalent calculation uses the spectral expansion of the Riemannian heat kernel,
Since for , the Mellin transform gives
The reduced trace decays exponentially at large . At small , subtract a finite number of terms of the heat kernel expansion, integrate those terms explicitly, and then integrate the remainder. Increasing the number of subtractions extends the formula successively to further left half-planes, producing the meromorphic continuation of rather than substituting in its divergent defining series.
For completeness, the needed curvature coefficient comes from the heat-kernel transport equations. In normal coordinates centered at , the Riemannian volume form has density . The leading Gaussian amplitude is therefore
Substitution of the Gaussian amplitude series into the heat equation gives the first diagonal transport relation . The curvature coefficient of the scalar heat kernel is consequently
In dimension two this implies
Set and . Splitting the Mellin integral at and subtracting these two heat terms gives, near ,
where is holomorphic for . As , only the second term survives at zero. Thus the zero-mode correction to the spectral zeta value is
The Gauss-Bonnet theorem, with twice the Gaussian curvature, makes this . For a connected surface , so the answer is ; a sphere gives , and a flat torus gives .
The purely local heat coefficient is . Sometimes a formal formula for quotes this coefficient with zero modes left implicit. For the actual positive-spectrum spectral zeta function defined above, the compact scalar problem necessarily has the additional subtraction . A strictly positive shifted operator is a different problem and has its own heat coefficients.
A Killing spinor is a nonzero spinor parameter for which the fermionic supersymmetry variations vanish on a bosonic background. In minimal four-dimensional supergravity without a cosmological term this means
It is an overdetermined equation and expresses an unbroken supersymmetry of the background. In the cosmological deformation the corresponding equation is ; in theories with additional bosonic fields their contributions also enter the fermion variations.
For the undeformed theory the connection preserves both the spinor adjoint and gamma matrices. Consequently the Dirac current
obeys . In particular , so is a Killing vector, and in fact is parallel. For the real cosmological deformation, and instead give , whose symmetric part again vanishes.
The causal property uses an ordinary commuting spinor. Take an orthonormal frame with , , and . Then . Reverse its time orientation and write , so
The matrices are Hermitian and obey . For every spatial unit vector , , whence . Choose along the spatial part of to obtain
A nonzero Killing spinor cannot vanish at an isolated point: its defining connection equation transports its value along every curve. Thus on a connected background this is a nonzero, nowhere-spacelike causal Dirac spinor current. The estimate does not need the Majorana condition, but it does need commuting components; Grassmann bilinears have no such ordinary positivity order.
For the positive energy theorem, take a complete asymptotically flat spacelike hypersurface with a spin structure , with suitable decay, Einstein constraint equations, and matter satisfying the dominant energy condition. Assume there are no inner boundaries, or that any inner-boundary contribution has the required nonnegative sign. The auxiliary commuting spinor approaches an arbitrary constant at spatial infinity. Define the real Nester two-form
Antisymmetry removes second derivatives except through the Riemann curvature tensor. Differentiation and the three-gamma spinor curvature identity used above give
The first term comes from differentiating the two spinor factors. For the second term, replace by half their commutator; its spin-curvature contraction is from each conjugate term. This is the reason the Einstein tensor, and hence the matter energy density, enters the boundary-charge identity.
In an orthonormal frame adapted to the future normal of , put . Since ,
Here is the spacetime spin connection pulled back to , including its extrinsic curvature. Choose the auxiliary spinor to solve the elliptic Witten spinor equation
The standard asymptotic-flatness and completeness hypotheses permit this boundary-value problem; the positivity identity also excludes a decaying homogeneous zero mode. It is this hypersurface equation that removes the negative square. A spacetime Killing spinor is not assumed for a general initial-data set.
Stokes' theorem and now give
The dominant energy condition makes the second integrand nonnegative because and are future causal. Evaluating the two-form with the asymptotic spin connection identifies its surface integral with the ADM energy and momentum:
For example its energy term uses to first order in the asymptotic spatial metric tensor perturbation. With this is exactly the ADM normalization. Since has eigenvalues , positivity for every boundary spinor proves
Killing spinors describe saturation: their spatial derivative term vanishes, and the contracted matter term must also vanish for the corresponding charge to be zero. They characterize preserved supersymmetry rather than being prerequisites for positivity. Under the usual rigidity hypotheses, zero total four-momentum makes the nonnegative integral vanish for a full basis of asymptotic spinors; their parallel extensions force a vanishing Riemann curvature tensor and give the Minkowski vacuum. A single Killing spinor with null current alone does not justify asserting that every component of the total momentum vanishes.
Use and the exterior mean-motion resonance condition . Then the resonant semi-major axis is , neglecting the small planet-to-star mass correction. For the geometric construction, take , , , and start the particle at pericentre. Solve Kepler's equation with . Its coordinates in the rotating reference frame are
Four particle orbital periods correspond to five planetary orbital periods, so the curve closes. At a pericentre, the mean longitude and the true longitude both equal , and the resonant argument becomes . Thus the four pericentre directions are modulo . The figure chooses and marks the corresponding angle to the planet; other phases rotate this pattern. Its second panel shows an instantaneous conjunction at apocentre, where the angle from pericentre to the common longitude is directly : here both true and mean anomaly equal . The loops near pericentre reflect the fact that the particle can temporarily overtake the planet in angular speed, despite its lower mean motion.
Figure 1.
Exterior 5:4 resonance with eccentricity 0.2 and resonant phase pi
. Details.
For the evolution equations, must be the mean longitude. The PDF calls it true longitude, but differentiating the literal true longitude gives the rapidly varying , not . To first order , so the true longitude version of has an extra orbital ripple. The stated Lagrange planetary equations and slow complex eccentricity equation are interpreted consistently as orbit-averaged equations in mean longitude. The pericentre geometry of an exterior resonance above is unchanged by this distinction.
Let be perturbing forces, and let be the particle's mass. The angular momentum about the star satisfies
The central gravitational force and the radial perturbation have zero torque. If the displayed denotes acceleration instead, the identical formula applies to specific angular momentum.
For a weak kick, the specific orbital energy changes by . A planet leading the particle pulls it forwards, while a trailing planet pulls it backwards. For a nearly circular exterior Kepler orbit, a forward tangential kick raises angular momentum and semi-major axis, lowers mean motion, and makes the next conjunction occur sooner in the synodic cycle. Its longitude retreats relative to the unperturbed conjunction sequence. A backward kick produces the opposite shift. At exact alignment the leading force is radially inward: before pericentre, and an inward kick does positive work; after pericentre, and it does negative work. These likewise raise or lower semi-major axis. At an apsis, a purely radial infinitesimal kick has no first-order energy change but changes the eccentricity vector, so angular momentum alone does not determine every timing change.
The timing sign follows explicitly from the near-circular synodic period, . At the nominal resonance, successive conjunction longitudes agree modulo , since . A persistent small change gives
This is the accumulated mean timing shift; an actual impulsive encounter can also change the instantaneous orbital phase and longitude of periapsis.
Retaining the given first-order resonant orbital perturbations, eccentricity damping and leading Keplerian detuning,
Here Kepler's third law gives . The approximation discards higher-order resonant contributions to and fast orbital terms. Differentiate the rotating complex eccentricity :
Since ,
The apparent singularity in at disappears in this complex eccentricity equation.
Take , and the conventional positive resonant coefficient . At a fixed point, the two balance conditions and give
Substitution yields
Thus fixed points exist only if . With strict inequality, put and . Then
respectively. The locations must also satisfy for the detuning expansion to apply. At equality the phases merge at and . If is assigned the opposite sign, shift the phase convention by ; the invariant stability distinction is the sign of .
For locally constant , the stable equilibrium is the branch in the leading slow model. At fixed , increasing on that branch increases the positive resonant torque, because . The resulting increase in semi-major axis lowers and hence lowers : the subsequent conjunctions shift back, restoring the phase. On the branch the feedback has the opposite sign and drives the phase farther away. Eccentricity damping removes the oscillatory disturbance on the restoring branch.
One can verify that this feedback stabilizes the coupled variables rather than just the phase. Treat as locally constant, and write . Then , and . At a fixed point . The characteristic polynomial of the linear stability analysis is
For , all coefficients are positive and the cubic Routh-Hurwitz stability criterion holds, since
All three eigenvalues have negative real parts. For the constant term is negative, so the polynomial has a positive real root and that branch is unstable. The merged case has a zero eigenvalue. The unrestricted functions printed in the question do not imply this unconditional stability classification. If their derivatives are retained, has, at equilibrium, the derivative
The trace of the full three-variable linearization is . As a concrete counterexample, take and , choosing their values to satisfy the displayed fixed-point conditions on the negative-cosine branch. Its trace is , so at least one eigenvalue has positive real part. Both coefficients are smooth. Thus the restoring conjunction feedback identifies the stable branch of the usual locally constant model, but additional assumptions on coefficient variation are necessary to assert stability for every allowed function.
Use von Neumann stability analysis on the infinite grid. For the Fourier mode , the semidiscrete eigenvalue is
Its real part is nonpositive for every frequency. The mode evolves by , whose modulus does not grow. The discrete Fourier transform and Parseval identity therefore yield
with a bound independent of , where . Hence the semidiscrete method is stable. Frequencies with are nondissipative rather than unstable.
The same sign can be seen directly by summation by parts: the centred operator is skew-adjoint, and for the unitary shift ,
This is an energy version of the dissipative second-order forward advection stencil. A subsequent time integrator would need its own stability condition; semidiscrete stability does not assert that every fully discrete method is stable.
The Padé approximant of the exponential is . Replacing each exponential in its original order gives the split Crank-Nicolson method
with the analogous definition for . Operationally, first solve the substep and then the substep. The matrices , are invertible for every , since are symmetric negative definite.
An orthonormal eigenbasis of diagonalizes . Its eigenvalue for is
Therefore , and similarly . Submultiplicativity gives the contractivity of split Crank-Nicolson diffusion:
No commutativity or simultaneous diagonalization is needed for this product bound. This proves unconditional discrete stability, not positivity of every stencil weight or unconditional maximum-norm monotonicity. Nor do the individually second-order rational substeps make the ordered splitting second order: when , its leading splitting defect remains .
A-stability removes the scalar linear decay restriction on the step size; it does not replace accuracy, nonlinear solvability or a suitable system norm. Its importance is clearest for a stiff differential equation, where rapidly decaying modes coexist with slowly varying quantities of interest. On the Dahlquist test equation , the exact solution decays for . A stable time integrator should not turn that decay into numerical growth.
For a one-step method write , where and is its stability function. Its linear stability domain consists of parameters with and a well-defined update. The method is A-stable when this domain contains the closed left half-plane. A consistent order- one-step method has near zero, but that local agreement alone says nothing about stability at large negative .
For a Runge-Kutta method with stage matrix , weights and , eliminating the stages on the test equation gives
An explicit finite-stage method has a polynomial . A nonconstant polynomial is unbounded along the negative real axis, so no consistent explicit Runge-Kutta method is A-stable. The explicit Euler method, with , is stable in the disk and on the negative real axis only for . A stiff mode with a large negative eigenvalue therefore forces an unnecessarily small step even after its exact contribution has almost disappeared.
The Backward Euler method has . For , , so it is A-stable. The trapezoidal rule has
and is also A-stable. The implicit midpoint method has the same scalar stability function, although its nonlinear update differs from the trapezoidal rule. The first method has order one and the latter two order two. Their implicit solves trade greater work per step for freedom from the scalar stiff stability restriction.
Figure 1.
Absolute stability domains of explicit Euler, backward Euler and the trapezoidal rule; shaded regions satisfy modulus of the amplification factor at most one
. Details.
Damping of very stiff modes separates A-stability from L-stability. For a one-step method, the latter additionally requires as tends to infinity in the left half-plane. Backward Euler method satisfies this. The trapezoidal and midpoint rules instead have along the negative real axis: a very stiff mode can alternate in sign with little damping even though its amplitude remains bounded. L-stability is often preferable for diffusion and rapid relaxation, while retaining oscillatory modes can be desirable in conservative dynamics.
For a linear multistep method, write
Its test-equation amplification roots satisfy . Absolute stability requires every root to lie in the unit disk or on its boundary with all unit-modulus roots simple; all parasitic modes matter. At this is the zero-stability root condition. Consistency requires and , and the Dahlquist equivalence theorem combines consistency and zero-stability to characterize convergence. Thus a claim of stability obtained by canceling a common factor is not a proof about the original recurrence and its arbitrary consistent starting data.
Explicit linear multistep methods cannot be A-stable either. If and , some lower coefficient of becomes unbounded as , while the leading coefficient stays fixed. If all roots remained in the unit disk, every coefficient-to-leading-coefficient ratio, an elementary symmetric function of those roots, would remain bounded. This contradiction forces an unstable root.
An important implicit example is BDF2 method, whose characteristic equation is . A unit root gives
The leading coefficient vanishes only at , outside the left half-plane. Near a small negative , both roots are inside the disk; no unit-circle crossing can occur in the connected open left half-plane. On the imaginary axis the boundary equality can occur only at , where the roots are . Hence BDF2 is A-stable and its two roots tend to zero for very large negative , giving strong stiff-mode damping. This is a complete boundary-locus argument, not just a boundary plot.
The Second Dahlquist barrier says that a convergent A-stable ordinary linear multistep method has order at most two. This theorem does not limit implicit Runge-Kutta methods or multiderivative multistep methods. In particular the A-stable third-order two-step multiderivative method in question 1 uses and satisfies different order conditions. Among collocation Runge-Kutta methods, the standard family results give Gauss collocation methods order with diagonal Padé approximants to as stability functions, and Radau IIA methods order with subdiagonal approximants. Both families are A-stable; the Gauss functions have a nonzero limit at infinity, whereas Radau IIA is L-stable. The corrected two-stage right-Radau method, for example, has
which matches the exponential through cubic order and tends to zero at infinity. Its poles lie in the right half-plane. For , , direct subtraction gives
proving its A-stability directly.
Sectorial A-alpha stability relaxes the requirement to a sector around the negative real axis, useful when the spectrum of a stiff problem remains in that sector; higher-order backward differentiation formulas illustrate the usefulness of this weaker requirement. For a normal linear system , unitary diagonalization reduces the numerical update to its scalar factors . For a nonnormal system, a scalar eigenvalue test alone does not give a uniform bound in the chosen norm. Nonlinear stability is different again: B-stability asks for contractivity when , and algebraic stability of a Runge-Kutta method is a standard sufficient criterion. Thus A-stability addresses a precise and valuable linear test, while step selection still has to meet accuracy requirements and the actual implicit equations must be solved on the relevant branch.