Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 115 1 d Solution Created 2026-10-03 Updated 2026-10-05
Under the standard smooth manifold hypotheses, the first answer is yes, although the isomorphism requires a choice. A partition of unity combines Euclidean inner products in manifold charts to produce a Riemannian metric . Nonnegative local weights summing to one preserve positivity, and local finiteness preserves smoothness. The musical isomorphism isIt is smooth and fiberwise invertible, with inverse having local matrix . It is therefore a vector bundle isomorphism. There is no preferred such map before a metric or equivalent additional structure is chosen.
The second answer is no. On , compare the trivial vector bundle with the Möbius line bundleBoth have rank one. A continuous section of a vector bundle of is represented on by a continuous function satisfying . If it were nowhere zero, the intermediate value theorem would prevent its endpoint signs from being opposite. Thus has no nowhere-zero section. The trivial bundle has the section , and a bundle isomorphism would carry that section to a nowhere-zero section of , a contradiction. Hence
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 131 2 Solution Created 2026-10-03 Updated 2026-10-05
For in the domain of the exponential map at , let be the unique geodesic with initial data and define . The domain consists of initial velocities whose geodesics exist on ; it is open, contains zero, and need not be the whole tangent space. The ordinary differential equation theorem used here says that a smooth vector field has unique maximal integral curves, with an open flow domain and smooth dependence on initial point and time. Applied to the geodesic equation on , it gives smoothness of the exponential map. Rescaling the parameter gives and hence . The inverse function theorem makes a diffeomorphism from a neighbourhood of zero to one of . Coordinates in an orthonormal tangent basis transported by this map are geodesic normal coordinates.
The geodesic sphere of radius is . For small , it is , a smooth compact hypersurface; larger distance spheres need not be smooth. The Gauss lemma statesIn particular radial and angular directions are orthogonal and the radial coordinate measures arc length.
For a piecewise smooth path , put . The Riemannian distance is over such paths from to . Connectedness makes it finite. Reversal and concatenation give symmetry and the triangle inequality; positivity and compatibility with the manifold topology follow locally from geodesic normal coordinates and the Gauss lemma.
Choose so that is a diffeomorphism on a neighbourhood of the closed tangent ball of radius . For , radial geodesics realize the distance from on that ball. Indeed, Gauss lemma bounds any path remaining in the normal ball below by its radial change, while a path leaving the radius- ball already has length at least . Thus the radius- sphere is the compact exponential image described above.
Given , take . The continuous function attains its minimum at some . Every path from to crosses that sphere, since the distance from is continuous. Splitting at a crossing gives length at least . Taking the infimum over paths, and using the triangle inequality in the other direction, proves the distance splitting through a small geodesic sphere equalityThis compact local argument does not assume completeness or the existence of a globally minimizing geodesic to .
For the Lie-group clause, a bi-invariant Riemannian metric has adjoint-invariant identity inner product. Differentiating this invariance gives . The Koszul formula for left-invariant vector fields has no metric-derivative terms, and this skew-adjointness reduces it toThis is the Levi-Civita connection of a bi-invariant metric.
Let be the integral curve through the identity of the left-invariant field with identity value . Uniqueness and left translation show wherever the local curves are defined. The field is complete: a fixed local existence interval about the identity translates to an equally long interval at every point of the group. Restarting the curve before any proposed finite endpoint extends it beyond that endpoint. Thus exists on , and uniqueness now gives the group law for every .
Since , this integral curve is a geodesic. Conversely any geodesic starting at the identity has the same initial data as one of these curves and agrees with it by uniqueness. Therefore the geodesics of a bi-invariant metric are one-parameter subgroups conclusion isThe zero velocity gives the constant, trivial subgroup. Here is the Exponential map of a Lie group; for this metric it agrees at the identity with the Riemannian exponential.
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 131 4 Solution Created 2026-10-03 Updated 2026-10-05
At each point, use the inverse metric on the cotangent space. It induces an inner product on -forms byEquivalently, wedges of distinct members of an orthonormal coframe form an orthonormal basis. Choose an orientation to define the ordinary global Hodge star operator by . It maps -forms to -forms and obeys . We use real forms; complex forms are obtained by complex-linear extension. Without an orientation an ordinary global star is not available without twisting the target, so orientation is implicit in this first part.
With the codifferential for and , define the Hodge Laplacian . On functions it is the nonnegative-sign Laplace-Beltrami operator, . A harmonic differential form satisfies .
For of degree , the displayed sign formula and star square giveApplying the identities again, at degrees and , proves the Hodge star commutes with the Hodge Laplacian relationSince star is invertible, is harmonic if and only if is harmonic. This part does not need compactness; zero-degree and top-degree terms are interpreted as zero when their degrees are outside the range.
For the remaining cohomological statements assume a compact oriented manifold without boundary. The Hodge decomposition theorem gives the -orthogonal decompositionwhere is finite-dimensional and consists of smooth forms. Formal adjointness gives , so harmonic differential forms are closed differential forms and coclosed differential forms.
If , formal adjointness and give . A continuous nonnegative function with zero integral vanishes, so . Then . Thus the a function with nonnegative Laplacian on a closed manifold is locally constant assertion isIt is globally constant when is connected. The PDF does not explicitly include connectedness here; on two disjoint circles one may take different constants, so literal global constancy needs that qualification. Boundaryless is the standard manifold convention in this argument; otherwise boundary conditions are necessary.
For a closed differential form , write its Hodge decomposition as . Since , we have . Formal adjointness yields , so . Thus represents its de Rham cohomology class. If two harmonic representatives differ by , their difference has . Therefore
Finally let . A connected Lie group admits a smooth path from the identity to . The maps give a smooth homotopy between the identity and . The homotopy invariance of de Rham cohomology says their pullbacks agree on cohomology. The maps preserve orientation, since their Jacobian signs cannot change from the identity along the path. As an orientation-preserving isometry, commutes on forms with star, exterior differentiation, codifferentiation and consequently the Hodge Laplacian. Hence is harmonic and represents the same cohomology class as a harmonic . Uniqueness now proves the harmonic forms are fixed by a connected isometric group action conclusionThe argument also applies componentwise when is disconnected.
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 215 2 b Solution Created 2026-10-03 Updated 2026-10-06
Take the real inner productThe stationary distribution is strictly positive for an irreducible Markov chain. Detailed balance gives , so is a self-adjoint operator. Choose a real orthonormal eigenbasis , where , , , , and for .
For fixed , apply the spectral decomposition to . Its coefficient against is , whereas . Thus the spectral kernel expansion of a reversible chain isThis holds for all integers ; at each multiplier is one, including when . Negative eigenvalues are allowed, and retain their signs for odd .
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 217 1 Solution Created 2026-10-03 Updated 2026-10-06
Put , the uniform bound on the supremum norm of the gradient. The fundamental theorem of calculus on a line segment, together with the inner product inequality , givesThus is Lipschitz continuous with respect to the L1 norm. The Gaussian random variables comprising and have finite exponential moments of their absolute values. More explicitly, for , Hölder's inequality gives . Consequently and every exponential expression below are finite, including when the covariance matrix is singular.
Use the convex function . Since is an independent copy of , Jensen inequality for the conditional expectation givesFor , define the Gaussian rotation of independent copiesWriting , both rotated covariance matrices are and their cross-covariance is zero. The pair has a centered multivariate normal distribution, so its components are independent and has the same probability law as . This calculation uses no inverse of and therefore also proves the assertion for degenerate multivariate normal distributions.
The path runs from to . The chain rule and the fundamental theorem of calculus give the pathwise identityApply Jensen inequality to the uniform probability measure on . Then use Tonelli theorem and the Gaussian rotation of independent copies to obtainThis proves the Gaussian rotation interpolation inequality in the required form:The stated almost-sure domination now gives an exponential moment of an absolute standard normal variable bound. If , completing the square in the standard normal density yieldswhere is the standard normal distribution function. Hence Markov inequality gives, for every ,The exponent is minimized at . Therefore the resulting Gaussian tail bound isNo independence of from or is needed: only the given almost-sure domination is used.
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 217 3 Solution Created 2026-10-03 Updated 2026-10-06
The Cameron-Martin space of a Gaussian random variable in a Banach space is most conveniently defined using the first Gaussian chaos. Work over the real separable Banach space , with continuous dual space . For a centered Gaussian random variable in a Banach space, defineThe set inside the closure is already a linear subspace. Every is a centered Gaussian random variable. The Bochner integral defining exists: Fernique's theorem states that for some , , and then Cauchy-Schwarz inequality gives .
The map is injective. Indeed, implies for all . Density of these observations in then implies . Define the Reproducing-kernel Hilbert space byThe injectivity makes the inner product unambiguous, and the completeness of makes a Hilbert space. The embedding into is continuous sinceFor put . The reproducing property isThis defines the Banach-space Gaussian RKHS even for a degenerate Gaussian measure.
For , write . There is a measurable coordinate on , obtained as a limit in L2 space of continuous linear observations, and has normal distribution . The Cameron-Martin theorem for a Gaussian measure, in its real separable Banach space form, states that if is the probability law of , then the probability law of is an equivalent probability measure to precisely when . For such a shift,For the two laws are mutually singular measures. The coordinate need not be a continuous functional on , and the sample need not belong to .
Apply this theorem to the shift . With , it givesChoose the measurable coordinate so that on a symmetric set of full -measure. This is possible by taking an almost-surely convergent subsequence of the approximating linear observations and intersecting its convergence set with its negative. Alternatively the joint law of is invariant under simultaneous negation, directly from those observations and their mean-square convergence. Since and the centered Gaussian measure is symmetric, the two exponential integrals with signs and agree. Consequentlybecause . The symmetric Gaussian translation lower bound is thereforeOnly symmetry and Borel measurability of are required; convexity is not needed.
There is one degenerate flaw in the last printed claim. The separable Banach space with has dense in , and all positive-radius balls have probability one. Nevertheless for every . Thus the statement needs the additional hypothesis . If the course convention excludes the zero space, this hypothesis is already implicit.
Under this necessary hypothesis, density of implies that contains a nonzero vector. Given , scale that vector to obtain with , and choose . The triangle inequality givesThe centered closed ball is a symmetric Borel set. Apply the symmetric Gaussian translation lower bound to getHence the Gaussian norm distribution is strictly increasing when is nonzero. This proves the strict increase of a Gaussian norm distribution without assuming that its distribution has a density or that spheres have zero probability.
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 301 1 d Solution Created 2026-10-03 Updated 2026-10-05
The oscillator-generated covariant photon Fock space has an indefinite Hermitian form, not a positive Hilbert space inner product. For a normalizable one-photon wave packet of polarization , its squared norm is proportional to . Thus gives a negative-norm photon state. The divergent of an unsmeared momentum eigenstate is a separate normalization issue, avoided by the wave packet.
Choose contravariant polarization vectors , , and two with for . The third spatial polarization is longitudinal polarization; the first two are transverse polarization. The timelike photon polarization is distinct from the longitudinal one.
The Gupta-Bleuler quantization condition sets the divergence of the positive-frequency part of a quantum field to zero on physical states:Restore the factors to the annihilation terms. Since and , Fourier transform gives the equivalent conditionIts sign depends on the chosen sign of the longitudinal polarization vector; the covariant condition does not.
To see its content for a general Fock state, temporarily discretize momentum and decompose one unphysical oscillator sector as , where . The temporal oscillator obeys , while . The condition therefore becomesEquivalently, the allowed finite-particle states use the transverse creation operators and only in the unphysical sector. Indeed the constraint acts on a polynomial of the two unphysical creation operators as , whose kernel consists of polynomials in their difference. Also , and a state containing is orthogonal to every constrained state because annihilates every such state. This argument applies mode by mode and extends by smearing to continuum momentum.
For one photon, is constrained only when , producing a null state. The Gupta-Bleuler null-state quotient removes these null directions. The condition excludes negative-norm physical states; quotienting its null states leaves the two positive-norm transverse photon polarizations. The condition alone gives a positive semidefinite Hermitian form, not yet a positive definite Hilbert space.
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 302 1 Solution Created 2026-10-03 Updated 2026-10-05
The two Lie groups have the same local infinitesimal structure, but different global topology. This difference determines which Lie algebra representations integrate to representations of each group.
An element of the SU(2) group is a unitary matrix of determinant one. Orthogonality of its columns and its determinant give the unique formThus its group manifold is the unit three-sphere in : this is SU(2) as the three-sphere. In particular it is compact, connected and simply connected. In terms of the Pauli matrices one may also write with real coefficients satisfying .
Differentiate and at . The tangent space consists of traceless skew-Hermitian matrices:The Lie bracket of a Matrix Lie group is the matrix commutator. With , the Pauli matrix multiplication law givesThis derives the SU(2) Lie algebra as a three-dimensional real Lie algebra. The Hermitian physics generators instead obey ; they are times the skew-Hermitian tangent generators, so these are consistent conventions.
The SO(3) group consists of real orthogonal matrices with determinant one. Differentiating at the identity givesThe determinant condition gives no additional infinitesimal constraint because a skew-symmetric matrix already has trace zero. Define . The cross product identity impliesTherefore the SO(3) Lie algebra has the same structure constants and is a Lie algebra isomorphism.
The global relation is the Adjoint double cover from SU(2) to SO(3). For , define byConjugation preserves the real space of traceless Hermitian matrices and its inner product . Hence is orthogonal. Continuity and connectedness, together with , put it in . Composition of conjugations makes a group homomorphism. If , then commutes with every Pauli matrix, hence is scalar; unitarity and determinant one leave precisely . The differential sends to , so it is an isomorphism. More concretely,induces rotation through angle about , by the Rodrigues rotation formula. Every three-dimensional rotation has such an axis and angle, proving surjectivity. ConsequentlyThe matrices and are antipodal points on the three-sphere, so the SO(3) group manifold is Real projective space . Equivalently, the closed axis-angle ball has opposite boundary points identified. The fundamental group is , whereas . A rotation lifts from to ; a rotation returns to . Thus the covering is the universal cover and the groups are not globally isomorphic.
For representation theory, specify finite-dimensional complex continuous representations. Compactness permits an invariant Hermitian inner product, obtained by averaging against Haar measure, and therefore complete reducibility. Complexifying either real Lie algebra gives the sl2 Lie algebra. Its finite-dimensional irreducibles are indexed by , have highest weight , and have dimension . By integration of a Lie-algebra representation, since is simply connected, every such Lie algebra representation integrates uniquely. The resulting homogeneous polynomial representation of SU2 isIn the spin angular momentum notation , its Hermitian eigenvalues are . The central matrix acts on the symmetric power by , so descent of an SU(2) representation to SO(3) occurs exactly when is even. HenceFor a reducible representation, every summand must satisfy the descent condition. The spin-one-half doublet is a genuine representation of the covering group but does not define a single-valued representation of ; the spin-one triplet does and is its vector representation. The distinction is topological, rather than a difference in their isomorphic Lie algebras.
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 302 3 Solution Created 2026-10-03 Updated 2026-10-05
The triple bond in the original Dynkin diagram gives the Cartan integerswhere is the long root and the short root. Their product is , and the inner product of distinct simple roots is nonpositive. The ratio of the two integers gives the squared-length ratio. ThereforeIt is convenient to normalize , and . Scaling the inner product does not affect the root system or the fundamental-weight relations.
For nonproportional roots , the root-string theorem states that the root string is consecutive:The endpoints are maximal, and reflection in reverses the string. The Cartan integer determines , not in general the total by itself. This result follows by restricting the Adjoint representation to the sl2 subalgebra associated with a root. If the roots are distinct simple roots, cannot be a root: its simple-root coefficients have opposite signs. Thus andHere length means number of roots; the number of intervals between them is one less. Distinctness matters. If , a reduced root system gives the set , with a missing zero between them, so the consecutive-string theorem and the displayed simple-root formula do not apply.
In the G2 root system, the initial strings areFor , . Since is not a root in a reduced root system, and : this generates . The remaining strings explain why the construction stops. The -strings through and are the same two-element string; the one through is a singleton since subtracting gives , and its Cartan integer is zero; the one through is the string . The -strings through , , and are the initial four-element string. Finally, is orthogonal to ; its -string is a singleton because is not a root. Apply the same reasoning to negatives. Using the permitted completeness of this procedure givesThe short positive roots are , of squared length two; the other three are long, of squared length six. Each root space is one-dimensional and the Cartan subalgebra has dimension two, soHere the dimension refers to the Lie algebra, with one Cartan generator per rank, not just the number of roots.
Write a prospective weight as . The pairings with simple coroots areThe fundamental weights are dual to those coroots. Solving the two linear systems gives, in the long-root-first numbering of this paper,The representation with Dynkin labels has highest weight , a short root. Numbering the short root first, as some references do, would call this the representation instead; the representation itself is unchanged.
The weight set of a finite-dimensional irreducible highest-weight representation is invariant under the Weyl group and lies in the convex hull of the orbit of its highest weight. All weights also differ from the highest weight by an element of the root lattice. The orbit of comprises the six short roots, so all six are weights. The lowering operators give the chainIn particular zero occurs: at weight , its pairing with is two, so the lowering operator is nonzero by the finite-dimensional sl2 Lie algebra representation theory. The -string through it is the usual three-weight string ; higher weight would lie outside the highest-weight convex hull.
There can be no further weights. Every point of that convex hull has squared norm at most , while a root-lattice point hasThe integer solutions of the bound are precisely zero and the six short roots: gives ; gives ; and gives . ThusThe final dimension uses the stipulated nondegeneracy of the weights. It also agrees with the Weyl dimension formula without that stipulation. The zero weight is not in the Weyl orbit of the nonzero weights, so this is not a minuscule representation.
In the diagram above, a coordinate label means . The left panel contains all twelve roots of a root system; the right panel contains the six short-root weights and the zero weight. The long-root-first convention is the same as in the calculations.
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 312 3 i Solution Created 2026-10-03 Updated 2026-10-06
Evaluate the metric inner products of the proposed orthonormal tetrad to first order in the vector cosmological perturbation:The orthonormal tetrad therefore has metric signature to the required order. Put and . The photon momentum components areThey satisfy the null vector condition to first order, with .
Use as parameter in the time component of the geodesic equation:The derivative of the direction is first order, since the background photon moves on a straight comoving line; is consequently second order. With ,The supplied Levi-Civita connection coefficients giveThe spatial gradients and all background expansion terms cancel. Since vanishes in the background, its product with is also second order. The photon energy redshift from a vector metric perturbation is thereforeAn independent check uses the covariant component . The covariant geodesic equation gives ; the conformal expansion term vanishes by the null vector condition and reproduces the same result.
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 314 3 Solution Created 2026-10-03 Updated 2026-10-06
Write for the constant isothermal sound speed squared and for the constant Alfvén speed squared. The horizontal magnetic field has no vertical magnetic tension; magnetostatic equilibrium isConsequently the magnetized isothermal atmosphere hasThe sign of specifies the orientation of the magnetic field. The adiabatic sound speed is , also constant.
For the proposed normal mode, let , , and temporarily write for its constant complex amplitude. The divergence amplitude is . Using and , the Eulerian and Lagrangian fluid perturbations giveSince , the magnetohydrodynamic total pressure amplitude isThe vertical derivative of is . For the magnetic tension perturbation, the terms proportional to from and cancel, leavingThus the linearization really reduces to the constant-coefficient eigenvalue problem
For a static ideal magnetohydrodynamics equilibrium with conservative gravitational acceleration, the magnetohydrodynamic energy principle gives a Hermitian operator for the restoring force in the mass density weighted inner product, provided the surface terms vanish. The factor removes the exponential mass density weight. This can be checked directly, without assuming stability: set , , and . Then , whereThis is a Hermitian matrix for real and real positive . Therefore every squared frequency is real; negative values give exponential growth, rather than an oscillation with a complex squared frequency. In an infinite atmosphere these are generalized plane waves, not finite-total-energy global modes.
At and , the equation, after cancelling its magnetic terms, gives . For nonzero magnetic field, the equation also gives . Substituting and then yieldsFor a nontrivial mode ; otherwise these equations and force every component to vanish. Substitution in the vertical equation givesFor example, the determinant provides an independent check:
To prove the claimed “if and only if”, rather than merely locating a zero eigenvalue, use the magnetic buoyancy energy criterionIf , every term is nonnegative, so all eigenvalues are nonnegative for every wavevector, including . If , choose , any fixed , and , , . The first two squares vanish; the remaining value iswhich is negative for sufficiently small nonzero . The Rayleigh quotient of this displacement is therefore negative, proving a negative eigenvalue and an unstable normal mode. This uses precisely the assumption that boundary conditions do not exclude the chosen wavelengths. Equality is a neutral threshold, not exponential growth.
Finally the plasma beta is . HenceThe magnetic field is assumed nonzero in the marginal-mode elimination and in the finite expression. The positive-sound-speed, nonmagnetic limit remains covered by the energy identity; for the usual it is stable. Parker instability occurs because magnetic pressure supports part of the atmosphere's weight while the gas can drain along gently bent magnetic field lines. The resulting buoyant displacement can release Newtonian gravitational potential energy; choosing long wavelengths along the magnetic field makes the stabilizing magnetic tension small. The completed-square calculation quantifies this competition, including the adiabatic process restoring force.
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 337 1 iii Solution Created 2026-10-03 Updated 2026-10-06
Eliminate the instantaneous Stokes flow velocity in favour of temperature. On a horizontal Fourier mode , the Stokes temperature-slaving operator maps to , where and . The temperature evolution has linear operator and bilinear map . Under the homogeneous thermal Dirichlet boundary conditions, is self-adjoint. Normalize its critical eigenfunction as and set ; the critical vertical velocity is .
At order , the critical eigenfunction equation gives . At order , the weakly nonlinear expansion contains the imposed second harmonic and the quadratic products of the critical mode: a horizontally uniform temperature correction proportional to and, in a general vertical-mode calculation, a second harmonic proportional to . These corrections are found by solving the noncritical boundary value problems, with homogeneous thermal data except for the imposed forcing.
At order , the method of multiple scales produces the slow derivative , the detuning term , and the two cross-advection terms involving first- and second-order fields. Project the component onto the adjoint eigenfunction using the vertical inner product. This is the solvability condition in the method of multiple scales: divide each resonant projection by . The detuning supplies with ; interactions of horizontal wavenumbers and permit with ; self-interaction through the slaved mean and second harmonic supplies . Other products have the wrong horizontal wavenumber. Reflection permits real coefficients with this cosine forcing. Thus the symmetry-allowed spatially forced convection amplitude equation isThere is a useful specialization that should not be silently missed. For the literal one-vertical-mode Stokes flow problem, the vanishing two-to-one forcing coefficient for Stokes convection makes at this order. To see this, write a positive second-harmonic forcing component as , incorporating the cosine's factor . Its coupling to the negative critical harmonic has projected integrand, apart from sign and its factor ,The integral vanishes because at both plates, even though is nonzero. This proves the cancellation without solving the forced profiles. The permitted coefficient is therefore zero times ; symmetry alone does not establish nonzero phase pinning for the equations actually supplied.
The same normalization makes the remaining coefficients explicit. Since , . The quadratic second harmonic cancels for , while the uniform correction is . Projecting gives . Thus for the literal model and this temperature normalization,A generic nonzero would require a nonvanishing projection in an amended physical model or a different forcing structure. It is still meaningful to classify the real-coefficient amplitude equation requested independently.
Write . Then and . These are a gradient flow for , so local minima give stable equilibrium points. At the origin the two eigenvalues are and . The origin has exponential asymptotic stability if , retains asymptotic stability with algebraic decay at , and is unstable if . At equality, obeys , since both linear coefficients are nonpositive. Integrating this inequality proves attraction even in the zero-eigenvalue direction.
For the stable nonzero equilibrium points are real; for they are imaginary:The real branch has Jacobian matrix eigenvalues ; the imaginary branch has . The oppositely aligned branch, when it exists, is a saddle equilibrium. No mixed real-imaginary nonzero equilibrium is possible when .
For , the origin is stable for , with algebraic decay at zero. If , the circle is radially attracting. Each point has Lyapunov stability but has a neutral phase direction, so it does not have individual asymptotic stability; the circle has orbital stability. This is the literal model's unpinned family. The general nonzero- branches instead exhibit phase locking to one of two phases separated by .
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 339 1 e Solution Created 2026-10-03 Updated 2026-10-06
For a convex cone in an inner product space, use the nonnegative-pairing conventionHere the pairing is the Frobenius inner product on real symmetric matrices. Let be the conic hull of the nonnegative rank-one matrices , and let . For , . This extends to conic combinations, and by continuity to their limits. Hence .
For the converse, if , separation from a closed convex cone supplies a symmetric withIn particular for every , so . The negative pairing then excludes from . ThereforeThe same generator test gives . This is the duality of copositive and completely positive cones; the next argument removes the closure.
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 340 2 c Solution Created 2026-10-03 Updated 2026-10-06
Use localization and cancellation of the assumed interval-adapted wavelet basis: at scale , every detail wavelet has support diameter at most , L1 norm at most , and annihilates polynomials of degree less than . These bounds also apply to the modified boundary wavelets; a restriction of an arbitrary whole-line wavelet without the cancellation-preserving boundary construction would not suffice. The finitely many coarse scaling functions can be treated separately.
Write with as above. Because , we have . Choose in the support of a function of , and let be the degree- Taylor polynomial of . The Hölder-Taylor remainder bound givesFor this is the Hölder seminorm bound. For , the integral Taylor remainder is bounded using , whose magnitude is at most .
The vanishing moments remove from the inner product. Applying the remainder estimate and the L1 norm bound on the localized support yields the wavelet coefficient decay for Hölder functionsThusHere depends on the fixed wavelet family, the exponent and boundary construction, but not on . The proof uses the assumed regularity and moments; it does not identify the minimal Daubechies wavelet order with its differentiability order, which need not coincide.
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 340 5 a Solution Created 2026-10-03 Updated 2026-10-06
Write and use the real pixelwise inner products. Let be the adjoint operator, defined by . For each nontrivial grid direction its one-dimensional formula iswith the analogous formula in for . Their sum is ; unused components and do not contribute. For , . If discrete divergence is defined, its sign convention is , as in the adjoint of a discrete forward gradient.
Define the compact convex set and its image . The image is compact and convex, hence closed, and contains zero. The finite-dimensional Euclidean duality formula , applied independently at every pixel, givesThus the penalty is the support function . Its convex conjugate is the indicator functional of a constraint set , equal to zero on and infinity elsewhere. For , the conjugate supremum is zero; for , the separation theorem supplies a direction with , and scaling that direction makes the supremum infinite.
The proximal operator is . It exists uniquely because the objective is continuous, coercive and strictly convex. The Moreau decomposition for a proper closed convex function gives . The second term is the Euclidean projection onto a convex set, so the projection residual for discrete total variation isAlternatively, the variational characterization of convex projection says exactly when for every . Hence attains the support function at , giving and the same optimality condition. This supplies a direct verification of the Moreau decomposition step in this instance.
The printed wording needs “the residual after a projection”, rather than “the projection itself”. In general this denoising map is not a projection onto any fixed closed convex set. A metric projection onto a closed convex set is idempotent, while denoising a sufficiently large jump twice shrinks it twice. Concretely, for and data with rows and , the minimizer has rows and . Reapplying the map gives rows and , so it is not idempotent. These formulas have a global certificate: take , all other dual entries zero, giving with rows and ; it saturates every positive row difference and yields the primal optimality condition. The qualification is therefore not merely a consequence of a restricted two-level ansatz.
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 341 4 ii Solution Created 2026-10-03 Updated 2026-10-06
The spatial shift in the PDF also requires a boundary value beyond the supplied right endpoint: if is the last interior index, the newest value is outside the interval. Thus the printed Dirichlet method has no specified boundary closure, and no unique stability range can be assigned to it as stated.
For completeness, the literal interior formula has an exact Von Neumann stability analysis on an infinite or periodic grid. Substitution of a Fourier mode givesFor the root of a polynomial tending to one as , expansion givesConsequently is unstable: sufficiently small nonzero frequencies have a root of a polynomial of modulus greater than one, and arbitrarily fine periodic grids contain such frequencies.
Conversely, if and , thenFor with ,Since , the complex quadratic Schur criterion yields both roots of a polynomial in the unit disk, strictly inside except at , where they are and . Uniform stability follows directly from the BDF2 discrete energy identity: the modal recurrence is , whose energy decreases because . Summing the modal energies proves a bound independent of the grid and even of ; no eigenvector-separation assumption is needed. This is the stability of a spatially shifted BDF2 stencil. Thus the literal periodic-grid result isThis cannot establish stability of an unspecified Dirichlet boundary condition closure. It also does not restore consistency of a numerical method.
For the intended unshifted backward differentiation formula, a spatial eigenvalue givesThe Schur stability criterion gives the root condition for a multistep method for every , with strict inequalities for . A useful uniform energy method proof, including the repeated interior roots that a bare separation argument would miss, is the BDF2 discrete energy identity. With andthe corrected Dirichlet discrete Laplacian scheme satisfiesThis follows by taking the inner product of with and using summation by parts. It controls both time levels uniformly, regardless of the mesh ratio. Hence the corrected method is stable for
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 341 6 Solution Created 2026-10-03 Updated 2026-10-06
A Runge-Kutta method advances through stagesThe Butcher tableau records , normally with . Linear absolute stability tests the amplification of a linear decay mode; nonlinear B-stability compares distances between two solutions of a dissipative vector field. Both describe propagation of perturbations, but they impose different conditions.
For the Dahlquist test equation, elimination of stages givesThe linear stability domain consists of for which the step is defined and . Poles or singular stage systems must be excluded, even if an output formula has a formal cancellation. For a scalar test equation, the proof is : the powers are bounded exactly when . For , the update is . If is a normal matrix, unitary diagonalization of a normal matrix makes the operator norm equal to the largest scalar amplification modulus. For a diagonalizable matrix, powers are bounded by the condition number of the eigenvector matrix times the maximal scalar power. These constants must remain uniform for mesh-dependent systems. A defective unit-modulus eigenvalue can instead produce a growing Jordan block; the scalar spectral test alone does not control a general matrix.
An A-stable method accepts the entire closed left half-plane. The Explicit Euler method has and a disk of stability , so it is not A-stable. No nonconstant polynomial stability function can be A-stable, since it is unbounded on the negative real axis; this excludes all consistent explicit Runge-Kutta methods. The Backward Euler method has , and when , proving A-stability. It is also L-stable, because for large left-half-plane .
The implicit midpoint rule and the trapezoidal rule have the same scalar stability function,The identity proves A-stability. On the imaginary axis their amplification has unit modulus, and along the negative real axis it tends to . Thus neither is L-stable; very stiff decay can persist as alternating numerical values. In contrast, the two-stage Radau IIA method has , with the explicit nonnegative modulus-gap proof in question 1, and is L-stable. This illustrates why L-stability matters beyond A-stability for a stiff differential equation.
For nonlinear problems, assume a real or complex Euclidean inner product and the dissipativity hypothesisfor all relevant states at each common time. Exact solutions satisfyA B-stable Runge-Kutta method reproduces this nonexpansion of distances whenever its stages exist. Put , , and . Then . Expansion of the output squared norm, followed by substituting in its linear terms, proves the Runge-Kutta contractivity identityIf and is a positive semidefinite matrix, the first sum is nonpositive by dissipativity, and the second is nonnegative: expand each coordinate of the stage vectors to express it as a sum of nonnegative quadratic forms. These are exactly the conditions for algebraic stability of a Runge-Kutta method. Hence algebraic stability implies B-stability. Applying B-stability to the linear dissipative problem , or its two-dimensional real form, also proves A-stability whenever stages are well defined.
The Backward Euler method has , hence ; the implicit midpoint rule has , , hence . Both are algebraically stable and B-stable. For the two-stage Radau IIA method,a positive semidefinite matrix, while both weights are positive. It therefore combines third-order accuracy, B-stability, and L-stability. The distinction between stability and solvability is essential: for Lipschitz continuous , the stage map is a contraction when , so the contraction mapping theorem proves small-step unique solvability. The preceding algebraic stability proof is an estimate for existing stages; it is not a theorem asserting stage existence for arbitrary step sizes and domains.
Linear A-stability is insufficient for nonlinear B-stability. The trapezoidal rule fails B-stability, even though its scalar stability function is identical to that of the B-stable implicit midpoint rule. Its Runge-Kutta method coefficients areFailure of this sufficient algebraic stability test alone would not prove failure of B-stability. An actual counterexample does: take and . Dissipativity follows fromThe unique step map obeys , because the left-hand side is strictly increasing and onto. At , , and implicit differentiation givesNearby initial states expand in distance. This proves a true nonlinear failure rather than merely failure of a coefficient criterion.
Finally, stability of a numerical method on finite intervals connects these contractivity properties to accuracy. If a one-step map is nonexpansive and its exact one-step defect has norm at most , then the error satisfies , hence for . A Lipschitz step bound gives the same order using the discrete Gronwall inequality. For a stiff differential equation, the value of B-stability is that the propagation estimate need not grow with a large negative dissipative rate. The accuracy constant still needs appropriate smoothness and uniform derivative bounds; stability alone does not establish a mesh-uniform error order.
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 341 7 Solution Created 2026-10-03 Updated 2026-10-06
For a fully discrete linear partial differential equation evolution, write in a specified discrete norm. Finite-time stability of a numerical method meanswith independent of the allowed meshes. For a multistep method, augment the state with its time-history values and include a stable starter. Bounds may grow as when the continuous problem has growth; demanding decay for all time would be a stronger assertion. The treatment of initial data, forcing, Dirichlet boundary conditions, and the norm is part of the hypothesis, not a detail supplied by an interior calculation.
For constant coefficients on an infinite or periodic uniform grid, Von Neumann stability analysis substitutes the Fourier mode . The discrete Fourier transform and the Parseval identity turn a uniform multiplier estimate into a discrete L2 norm estimate. In a one-step scalar scheme, proves contractivity, while gives finite-time stability. In a multilevel scheme, every root of a polynomial in its amplification polynomial of a multilevel finite difference scheme matters, as does uniform control of the associated companion matrix.
For the Forward Euler diffusion scheme with ,Thus for all frequencies exactly whenSufficiency follows since ; necessity follows from the highest-frequency mode on even periodic grids. This gives the usual parabolic mesh restriction . For the Backward Euler diffusion scheme,so every is stable. The Crank--Nicolson method hasagain of modulus at most one for all , but stiff modes approach rather than zero. These are PDE manifestations of the A-stability and L-stability distinctions for time integration.
For advection , assume and set . The upwind finite difference scheme is , withIt is contractive when . This also has a direct maximum-norm proof: each newest value is a convex combination of two old values. The resulting Courant–Friedrichs–Lewy condition expresses that the numerical domain of dependence covers the physical one. For negative , the upwind direction must be reversed. By contrast, Forward Euler method time stepping with the centered first finite difference has and for nonzero . Under fixed nonzero refinement it is unstable, since a fixed nontrivial mode grows geometrically over steps. This is not a claim of instability under every imaginable mesh coupling: instead bounds its spurious finite-time growth, since .
The leapfrog advection scheme illustrates a multilevel subtlety. Its amplification polynomial of a multilevel finite difference scheme isFor , both roots of a polynomial have unit modulus and separation at least . A uniformly conditioned eigenbasis of the companion matrix then proves power boundedness of a two-level Fourier scheme for arbitrary history data. At and a grid admitting , the roots coincide on the unit circle, and a Jordan block produces growth. Thus the endpoint fails the ordinary arbitrary-history root condition for a multistep method; checking only that both moduli equal one misses the instability. Restricting the starter or filtering the parasitic mode is an additional hypothesis.
The energy method handles variable coefficients and finite boundaries for which Fourier stability analysis may be unavailable. For the Backward Euler diffusion scheme with homogeneous Dirichlet boundary conditions, define . The discrete summation by parts identity isTake the inner product of with . The elementary identity givesThis proves unconditional contractivity including the actual boundary treatment. More generally, a dissipative operator gives of operator norm at most one: with , dissipativity implies , and the Cauchy-Schwarz inequality gives . In finite dimensions this also proves invertibility. For the second-order backward differentiation formula, the BDF2 discrete energy identity used in question 4 controls both time levels and proves unconditional diffusion stability. A repeated root strictly inside the disk is harmless here; the energy argument supplies the uniform bound without a singular eigenvector formula.
A second technique is eigenvalue stability analysis of a finite difference method. Under the method of lines, a time integrator advances by . If is a normal matrix, the scalar linear stability domain criterion on all controls its powers exactly in the corresponding L2 norm. For the centered Dirichlet discrete Laplacian, its eigenvalues lie between and zero; intersecting this interval with a time integrator's stability interval determines its mesh restriction. A uniform bound on the diagonalization of a matrix is needed for nonnormal diagonalizable systems. Eigenvalues alone can be misleading. For example,has both eigenvalues inside the disk for , but the upper-right entry of is . Taking , , this grows like . Hence stable scalar eigenvalues do not give mesh-uniform stability. Boundary closures can introduce precisely the extra growth that an interior Fourier symbol overlooks, so boundary stability of a finite-difference method must be checked separately.
For nonlinear spatial discretizations, a useful replacement for Fourier stability analysis is a convexity argument. If the Forward Euler method map is nonexpansive in a chosen norm for , the second-order strong stability preserving Runge-Kutta methodis also nonexpansive under the same restriction. Indeed, for two inputs, the triangle inequality gives . Its Taylor expansion is , proving order two. Such strong stability preserving Runge-Kutta methods transfer suitable forward-step bounds without relying on a linear spectral calculation.
Finally, consistency of a numerical method explains what stability buys. For a linear Hadamard well-posed problem, the Lax equivalence theorem equates convergence of a consistent discretization with its stability, under the stated approximation-space and norm hypotheses. Directly, if the error satisfies , iteration gives the discrete Duhamel principleThus vanishing normalized local truncation error and initial error yield convergence. With a nonlinear Lipschitz step estimate, the discrete Gronwall inequality plays the same role. A stable but inconsistent stencil, such as the literal defective formulas in questions 3 and 4 under their usual refinement, is not rescued by any amplification bound. The decisive checks are a mesh-uniform evolution bound, a valid boundary closure, and consistency with the actual PDE.
Past exam of the mathematics course of the University of Cambridge 2017 ii Paper 4 18G b Solution Created 2026-09-24 Updated 2026-10-05
The Clebsch-Gordan decomposition for SU2 statesRealize the tensor product as polynomials of bidegree in and . The determinant is invariant under the simultaneous SU2 action. For each , the vectorhas torus weight and is killed by . Applying produces nonzero vectors of distinct weights, with . The relations and give , so their span is an irreducible copy of .
These irreducible submodules have distinct highest weights and form a direct sum: a sum of earlier ones cannot contain an irreducible module inequivalent to all its summands. Equivalently use complete reducibility, obtained by averaging an inner product over the compact group. Assuming , their total dimension iswhich equals the dimension of the tensor product. Hence they exhaust it, proving the theorem.
