Past exam of the mathematics course of the University of Cambridge 2002 iii Paper 69 1 a Solution Created 2026-10-03 Updated 2026-10-14
Put for the membrane's internal dimension and assume . The small-slope expansion isThe zeroth-order term is the projected area and is independent of the height field. Discarding it changes only an overall normalization of the partition function. Thus the leading fluctuation energy and its statistical weight areOne obtains this Gaussian field theory by discretizing the height field on a microscopic mesh, integrating every height with its Boltzmann weight, and then writing the continuum measure. A microscopic cutoff is retained for quantities dominated by short wavelengths. Quadratic order here really means small slopes: a large rigid height shift is still free because the energy has no dependence on itself.
The uniform height mode has zero energy. Fix its mean, impose a boundary condition, or divide out the integral over global translations to normalize the partition function. Height differences are independent of that choice. The approximation is the tension-only part of small-slope elastic membrane energy; there is no bending rigidity or confining potential in this part.
Past exam of the mathematics course of the University of Cambridge 2002 iii Paper 69 1 Solution 2026-10-14
Spontaneous symmetry breaking occurs when the statistical Hamiltonian is invariant under a symmetry but a chosen equilibrium pure state is not. A small selecting field or a boundary condition chooses one of several states related by the symmetry. For a finite unconstrained system the symmetric ensemble generally has zero order parameter; spontaneous order is defined by taking the thermodynamic limit before removing the selecting field. A ferromagnet with a selected magnetization direction is the standard example.
If the broken symmetry is continuous, nearby symmetry-related states cost no energy for a uniform displacement along the manifold of states. Locality and a nonzero stiffness then make a slowly varying displacement cost gradient energy rather than a mass term. These gapless long-wavelength fluctuations are Goldstone modes. In the conventional relativistic internal-symmetry setting their count is the number of broken generators. Nonrelativistic dynamics can pair broken generators, as in a ferromagnet; gauge symmetries and broken space-time symmetries require additional qualifications. A discrete broken symmetry alone does not produce a Goldstone mode.
For an interface with a selected mean height, the continuous displacement is a rigid translation normal to the mean plane. The tension gives a stiffness for spatially varying translations. The following calculation concerns this static Goldstone sector. It determines the excitation energy and equilibrium fluctuations, but a frequency dispersion requires a separate inertial or dissipative dynamics. In low dimensions thermal roughness also limits the existence of a strictly localized flat interface in an infinite system.
Past exam of the mathematics course of the University of Cambridge 2002 iii Paper 69 4 b Solution Created 2026-10-03 Updated 2026-10-14
Work first on the infinite chain, or with a periodized exponential kernel. Put and use lattice Fourier momentum . Summing the geometric series gives the exponentially coupled Ising inverse kernel:ThusAlthough the original kernel is long-ranged, its inverse contains only a constant and a cosine and therefore is nearest-neighbor. Fourier orthogonality givesSubstitution into part (a) yieldsThis exact bulk identity does not require replacing by ; small is used later to motivate a smooth continuum field.
For a finite open chain the inverse has off-diagonal entries , interior diagonal entries and endpoint diagonal entries . Accordingly the displayed homogeneous potential acquires the endpoint correction . These boundary terms do not change the bulk thermodynamic conclusion. Stating infinite-chain or periodized boundary conditions is necessary if the translation-invariant expression is asserted exactly.
Past exam of the mathematics course of the University of Cambridge 2002 iii Paper 69 4 c Solution Created 2026-10-03 Updated 2026-10-14
With lattice spacing set to one and a smoothly varying auxiliary field, the discrete action becomesThis is precisely the Euclidean path integral of a quantum particle with , mass and potential . With fixed endpoints it gives the transition amplitude , while periodic classical boundary conditions give its trace. The corresponding operator isThis is classical Ising chain as a quantum particle. It is the continuum approximation to the exact positive transfer kernel, so its numerical spectrum is accurate when the relevant auxiliary configurations vary slowly. The transfer-kernel reasoning and the domain-wall argument below also support the qualitative conclusion without treating a continuum approximation as an exact finite-lattice identity.
At zero field write . Expansion givesA double well exists when ; the stationary points satisfy . For there is a single minimum at zero, with the quadratic curvature marginal at equality. Thus a double well is the low-effective-temperature regime of the model, not an unconditional property for every . The exact potential remains confining because at large .
Even when two wells exist, a quantum particle with finite mass and a finite barrier tunnels between them. Its one-dimensional confining Schrödinger operator has a unique, strictly positive ground-state wavefunction. At reflection symmetry makes that state even, so every odd observable has zero expectation. An antisymmetric first excited state is split from it by a positive gap. A spectral expansion of an odd-observable correlation has factors ; the smallest contributing gap therefore gives a finite classical correlation length. The original spin observable, after decoupling, is , which is odd at zero field. HenceA deep well can give an exponentially small quantum tunnelling gap and an extremely long correlation length, which can mimic order in a finite sample, but it does not produce true long-range order.
The classical counterpart is the finite cost of a domain wall. Because the pair sum is doubled, a sharp wall has dimensionless energyAt any finite temperature such walls have nonzero fugacity, divide a long chain into domains and destroy spontaneous magnetization. The range is large for but finite. The zero-temperature limit and an infinite-range limit are singular and must not be interchanged casually with the thermodynamic limit; an infinite-range interaction also requires an extensive normalization before comparing models.
A nonzero tilts the potential: to first order . A positive field favors the positive well and produces positive magnetization. For fixed finite parameters the ground state and classical free energy vary smoothly, so there is no zero-field discontinuity due merely to a double-well bare potential. A discrete broken-symmetry state appears only in an appropriate singular limit, not by ignoring quantum tunnelling.
More generally a local -dimensional classical statistical model can be sliced along one coordinate and represented by a transfer operator. That coordinate becomes imaginary time; the other coordinates become the spatial coordinates of a quantum system. Infinite extent in the transfer direction corresponds to its ground-state, or zero-quantum-temperature, limit. The classical Ising model is thus related to a quantum Ising system in spatial dimensions, with transverse dynamics generated by the transfer operator. For this is quantum mechanics of finitely many local variables, where quantum tunnelling restores symmetry. For the quantum system has infinitely many spatial degrees of freedom and can have distinct symmetry-broken ground states, allowing a classical finite-temperature ordering transition. Locality and anisotropy of the mapping determine the quantum couplings; the dimension correspondence alone does not fix them.
Past exam of the mathematics course of the University of Cambridge 2003 ib Paper 4 16C Solution Created 2026-09-24 Updated 2026-10-14
Irrotational flow means ; incompressible flow means . The given streamfunction automatically satisfies the second condition, since . For an oriented path with unit tangent , choose its right-hand unit normal . ThenThe opposite normal reverses the sign; the source's unspecified normal orientation must be chosen consistently. The volume flux between two streamlines is thus their streamfunction difference. A closed contour has zero net flux, expressing mass conservation for constant-density flow without a source inside.
For the channel, a rigid boundary has zero normal velocity and hence constant streamfunction. Choose zero on the lower wall. Uniform inflow at infinity must have horizontal velocity , so there . The upper wall then has constant . The unbroken left wall connects to it and also has constant . The jump by between the walls at the corner represents the withdrawn flux. Finally the vorticity is , so irrotationality gives .
Write . Then vanishes at and at infinity, while for . Separation of variables gives the harmonic sine mode in a half-strip . The Fourier sine series coefficients of the left boundary areTherefore the corner sink in a semi-infinite channel isFor the series and its differentiated series converge on compact subsets, making the sum harmonic. At , for , so the left-wall value is exactly . The other wall values and far-field limit follow directly.
The plus sign printed in the final series is inconsistent with the printed boundary conditions. It would give on the left wall, not . The negative sign above is required. As a local flux check, near the corner the wall values give ; the radial velocity is and its outward flux across a quarter-circle is , precisely withdrawal.
Past exam of the mathematics course of the University of Cambridge 2003 iii Paper 50 4 Solution Created 2026-10-03 Updated 2026-10-14
Write the integrand, without its common overall factor, as . Keep each constant Grassmann variable on the left. For the first supersymmetry component, varying the bosonic term givesVarying its fermion term gives . The plus sign in the second term comes from moving the odd parameter past . The middle terms cancel and the rest form . The second supersymmetry component similarly gives . Hence, off shell,This proves the local total-derivative identity for rigid worldsheet supersymmetry; the action is invariant when the corresponding boundary flux vanishes and the transformation preserves the allowed fields. A constant parameter preserves a periodic boson and a periodic fermion. In a Neveu–Schwarz sector, however, is antiperiodic for constant and does not preserve the periodic boson. Thus the bulk identity does not supply a global constant supercharge in every spin sector. This is the spin-structure obstruction to constant worldsheet supersymmetry; the local superconformal formulation handles the different sectors consistently.
For independent field variations, integration by parts gives bulk termsThe equations of motion are thereforeFor variations vanishing at the temporal ends, the spatial seam term is proportional to the difference between its values at and ofThe bosonic part cancels for closed string periodicity, also in a fixed winding sector. For a real fermion impose and the same relation on its variation. Seam cancellation requires , so the two untwisted possibilities are periodic, (Ramond sector), or antiperiodic, (Neveu–Schwarz sector). The two chiralities of a closed string can choose these signs independently, giving R–R, NS–NS, R–NS and NS–R sectors. This discussion assumes the ordinary Lorentz-invariant untwisted boundary conditions of the flat theory.
The chiral equations and boundary signs give the following Fourier serieswhere in R and in NS. Thus integer or half-integer modes follow directly from the Fourier series monodromy. With the conventional string oscillator normalization their anticommutation relations and reality conditions areFor an open string the endpoint gluing identifies the two sets of coefficients, leaving one R or NS tower.
In ten-dimensional light-cone gauge in string theory, the eight transverse bosonic modes have positive integer level and each NS fermionic creator has a positive half-integer level. For the unprojected low levels of the open RNS string, all states through in the NS sector areThe antisymmetric index restrictions are consequences of fermionic anticommutation; the boson–fermion pair has all 64 choices. The displayed partitions exhaust levels at most . Here the NS normal-ordering constant of a string is .
In the R sector , and the zero modes satisfy . Thus generate the Ramond zero-mode Clifford algebra. The unprojected ground states span a sixteen-dimensional transverse spinor space, . All states in the requested range areThere is no intervening half-level, since all nonzero R modes have integer level. Acting with the zero modes only moves within the ground spinor space and does not generate further independent copies.
The GSO projection keeps odd NS fermion-excitation parity, removing the tachyonic vacuum and the entire NS level. In R it keeps one eigenvalue of total fermion parity, whose zero-mode part is spinor chirality. It leaves eight massless spinor states. At R level one, the retained families are and , of opposite ground chirality, so there are states. The retained NS states are spacetime bosons and the R states spacetime fermions. Consequently the retained levels have at mass zero and at , with no tachyon.
Equality at every higher mass level follows from the standard Jacobi abstruse identity, or equivalently the equality of the transverse generating functionsThe left side counts the projected NS bosons, the right side the projected R fermions. Their common coefficients start . This consistent truncation yields spacetime supersymmetry; the bulk worldsheet symmetry alone would not remove the tachyon or impose the projection.
Past exam of the mathematics course of the University of Cambridge 2003 iii Paper 51 2 Solution Created 2026-10-03 Updated 2026-10-14
The stress-energy tensor generates translations. One may improve it by adding a divergence without changing the translation charges under suitable boundary conditions, and choose it symmetric in a relativistic theory. If is symmetric and conserved, a spacetime vector field defines the current . For a Conformal Killing vector field,Thus a traceless conserved stress-energy tensor produces all the conformal charges. The dilation current is , and the special conformal currents are . A trace that is merely the divergence of a virial current gives scale invariance; conformal invariance requires the additional improvement that removes the trace. These statements distinguish an actual traceless tensor from a trace discarded without justification.
For a free massless scalar field, the canonical stress-energy tensor has trace . For , the improved stress-energy tensor of a free massless scalar isIts improvement is identically conserved. On , the identity shows that its trace vanishes. In , no improvement is needed: the canonical stress-energy tensor is already traceless. In Lorentzian null coordinates, conservation and tracelessness imply that its two surviving components depend separately on and . After Euclidean continuation, these become the holomorphic stress-energy tensor and its antiholomorphic counterpart , withLocal conformal changes are arbitrary holomorphic and antiholomorphic coordinate changes, not just the three global Möbius transformation generators in either sector.
To make the free-field quantum normalization explicit, rescale the field to , so that its plane contraction is . Define , so . Normal ordering givesA single Wick contraction shows that is a weight-one primary operator. The two double contractions in produce ; single contractions, followed by Taylor expansion at , produce . Thus the operator product expansion isThe same holds independently for , with . The undifferentiated noncompact scalar has a logarithmic contraction and a zero-mode subtlety; its derivative current has the ordinary primary transformation law.
The holomorphic stress-energy tensor implements infinitesimal conformal transformations through contour charges . A primary field of weight obeysFor example the free-field vertex operator has holomorphic weight , as follows by contracting with the exponential twice. In correlators the contour may encircle each insertion, giving the conformal Ward identity . This also yields the global conformal Ward identities for chiral correlators.
Writing and integrating the stress-tensor operator product expansion gives the Virasoro algebraIts central charge is the quantum correction to the classical conformal algebra; the correction vanishes on . Equivalently, in the pullback convention a finite holomorphic change sends to , where is the Schwarzian derivative. This anomalous term vanishes for Möbius transformations. On a curved background the related trace anomaly is proportional to its curvature; it does not destroy the flat-space conformal symmetry of this massless theory.
Past exam of the mathematics course of the University of Cambridge 2003 iii Paper 51 4 Solution Created 2026-10-03 Updated 2026-10-14
For , let be a Conformal Killing vector field, put , and assign the massless scalar field scaling dimension . In a finite conformal transformation with metric scale , the transformation law is . Its infinitesimal fixed-coordinate form isFor , contract the Conformal Killing equation with . Direct differentiation then givesHere . Flat conformal Killing integrability gives , with the Lorentzian version obtained by replacing the flat metric consistently. Thus the action variation is a boundary integral, vanishing for suitable boundary conditions. This proves conformal invariance with field dimension , without using the field equation. The connected finite symmetry follows by integrating the infinitesimal generators. Equivalently its conformal currents follow from the improved stress-energy tensor of a free massless scalar.
In one dimension the metric conformality condition is vacuous. The same action-variation calculation contains , so the free action is invariant under the projective subgroup , rather than arbitrary reparametrizations. The usual field-theory conformal-invariance claim above concerns .
In , , so the field transforms simply as a scalar. There is no improvement term and no derivative of the conformal scale in the variation. All local independent reparametrizations of the two null coordinates preserve the action; after Euclidean continuation they are holomorphic and antiholomorphic maps. The conformal algebra is therefore infinite-dimensional, in contrast to .
Now use Lorentzian coordinates with spatial period . The canonical momentum is , and canonical quantization imposes , with the periodic delta distribution . Solving the wave equation and imposing reality gives a convenient oscillator normalization:Strict spatial periodicity excludes a winding term. The canonical commutation relation is equivalent towith and the corresponding tilded relation. Both oscillator families commute with . To verify the normalization, the zero-mode commutator supplies , while each oscillator family supplies . Their sum is exactly . The positive-index modes are annihilation operators and the negative-index modes are creation operators, giving two independent bosonic Fock spaces together with the free-particle zero mode.
Define and the normal ordering quadratic generatorsThese are the Fourier modes of the two chiral components of the stress-energy tensor. The oscillator operator commutators immediately give . Applying this to both factors of gives , except for a scalar term when , produced by reordering creators and annihilators. Its value is independent of the momentum sector. Evaluate it in the formal zero-momentum oscillator vacuum. For ,The two contractions of each pair explain the factor . Hence the free-boson Virasoro central term yieldsThe tilded modes give a second commuting Virasoro algebra with . The adjoint condition and the positive oscillator norms make this a unitary highest-weight Virasoro representation in each momentum sector. In a fixed momentum sector, , so the oscillator ground state is a highest-weight state of weight in either sector. For a noncompact field the zero mode is a free particle: definite momentum ground states are generalized states, rather than a normalizable zero-mode vacuum. The oscillator representation and central term remain well defined. With the standard cylinder vacuum-energy subtraction, the Hamiltonian isThe constant is the Casimir energy of a two-dimensional conformal field theory; omitting it amounts to a different additive Hamiltonian convention, not a change of the plane Virasoro algebra.
Past exam of the mathematics course of the University of Cambridge 2003 iii Paper 59 2 a Solution Created 2026-10-03 Updated 2026-10-14
Use the Galerkin method: substitute the retained trigonometric modes and project onto them with the inner product over the rectangle. Each mode satisfies the printed boundary conditions. Write , , and abbreviate the four mode prefactors by , , and . The Jacobian determinant in the advection terms is .
The stream function is a single Laplacian eigenfunction, so and exactly. For the temperature modes, direct multiplication givesFor the transverse-velocity modes, the retained products areFor example, and supply the signs of the two retained and couplings. The third harmonics are orthogonal to the retained modes and are discarded by this Galerkin method; they are not claimed to vanish in the full partial differential equations.
Before changing time, the five coefficient equations areDividing by gives the stated five-mode rotating-convection truncation, withThe squared rotation parameter is . The printed is proportional to rotation rate, so it is a square-root-type rotation parameter rather than the more usual Taylor number proportional to squared rotation rate; the algebra uses the coefficient as printed. At zero rotation the physical transverse velocity modes vanish; the normalized equations can still be understood through their continuous zero-rotation limit.
Past exam of the mathematics course of the University of Cambridge 2003 iii Paper 60 4 Solution Created 2026-10-03 Updated 2026-10-14
For a stationary finite-wavenumber instability, the relevant Landau–Ginzburg model is the real Ginzburg–Landau equation: its coefficients are real although its pattern amplitude is complex. The complex amplitude records both the strength and the position of the underlying pattern. Assume a simple critical pair of Fourier modes , a quadratic nondegenerate maximum of the linear growth rate, and no additional conserved or marginal mode. The last assumption is essential: the conserved field in the previous problem requires an additional evolution equation.
To derive the model, write a real evolution equation schematically as , with symmetric multilinear nonlinearities. Put , , and use the method of multiple scales expansionHere . At the next order the noncritical mean and second harmonic are slaved to the critical mode:where and . These inverses exist under the no-other-critical-mode hypothesis. Let be an adjoint null vector normalized by . At order , the Fredholm solvability condition requires the resonant forcing to be orthogonal to . The cubic resonant contribution isThis shows why the cubic coefficient includes feedback from slaved harmonics as well as any original cubic nonlinearity. The resulting amplitude equation, after naming its linear control coefficient , isHere is the critical linear dispersion relation for growth. The condition gives a supercritical bifurcation; calls for higher-order saturation, for example the quintic real Ginzburg-Landau equation.
A concrete check is the cubic Swift–Hohenberg equation . With , its linear operator contributes at resonant order, while the cubic resonance is . Hence . This directly realizes the perturbation-theory coefficients above.
Symmetry independently constrains the form. A fast spatial translation rotates by a constant phase, so every nondifferential monomial must transform like . The leading allowed nonlinear term is , and a quadratic term is forbidden. A spatial reflection maps to ; it makes the leading coefficients real and rules out a drift term in a reflection-symmetric rest frame. The vanishing first derivative of the linear growth rate at is consistent with the absence of a leading first spatial derivative. These arguments constrain terms but do not calculate or replace the solvability calculation. Breaking reflection or studying an oscillatory instability generally leads to the complex Ginzburg–Landau equation with genuinely complex coefficients.
The principal steady solutions are the zero state and the Ginzburg-Landau plane wavesTheir physical carrier wavenumber is . For the constant phase is arbitrary, giving a circle of uniform-amplitude states. The zero state has Fourier growth rate and is stable for . For , each existing plane wave is stable to spatially homogeneous amplitude perturbations, with radial eigenvalue , and neutral to a constant phase shift.
Steady modulated states are also possible. Writing gives the conserved spatial phase current andIts first integral is . Boundary conditions select nonuniform periodic amplitudes, fronts or defects from this spatial dynamics. In the invariant real subspace, for example,is a steady kink joining opposite real states. It is not stable in the unrestricted complex equation: the imaginary perturbation has eigenvalue under . This illustrates why a real-amplitude calculation can miss a phase instability. With real coefficients the Lyapunov functionalsatisfies and under compatible boundary conditions. It supplies a variational structure absent from a generic complex-coefficient model.
To test full linear stability of a plane wave, set with real perturbations . Let . The linear equations areFor perturbation wavenumber , the sideband spectrum of a real Ginzburg-Landau plane wave isThe phase branch at small isTherefore the robust stable band on the infinite line is , narrower than the existence band . Crossing its boundary gives the longitudinal Eckhaus instability. At the boundary the quadratic phase-diffusion coefficient vanishes and the leading nonzero term is negative of order . More exactly, a nonzero sideband grows only when ; a finite periodic domain for the wave envelope may exclude that interval. Constant-amplitude stability is understood modulo its neutral global phase, and the allowed perturbation wavenumbers must always be specified.
In a second spatial dimension, an anisotropic medium that selects one roll orientation permits an equation with . For positive coefficients its plane-wave analysis is similar, with the second-dimensional contribution adding damping. But simply replacing by an isotropic Laplacian is not a faithful extension for rotationally invariant stationary rolls. There the critical wavevectors form a circle. Near ,This forces the anisotropic transverse scaling of an isotropic roll envelope: , , . For the positive carrier convention used here the Newell–Whitehead–Segel equation isThe opposite fast-carrier convention reverses the accompanying imaginary sign and detuning convention. For a transverse phase perturbation of our convention gives . Thus permits the transverse zigzag instability, in addition to the longitudinal Eckhaus instability. The fourth transverse derivative and mixed derivative in the square encode the critical-circle geometry.
One complex amplitude still describes only a narrow cone of nearby roll orientations. Competing well-separated wavevectors require coupled amplitude equations; resonant triples can admit quadratic interactions and produce the hexagonal convection amplitude equations when the underlying sign symmetry allows them. Strong disorder, defect cores with large gradients, far-from-onset patterns, boundaries, mean flows and conserved fields may all lie outside the single-amplitude approximation. The method of multiple scales gives controlled leading behavior for small amplitudes and slow modulations, rather than a universal model of all two-dimensional pattern selection.
Past exam of the mathematics course of the University of Cambridge 2003 iii Paper 61 1 a Solution Created 2026-10-03 Updated 2026-10-14
Let . Expanding the proposed local relation gives . Choose the derivative terms in to cancel the mixed spectral terms:Indeed,so the resulting divergence is exactly . This establishes the identity for every complex spectral parameter , with the given differentiability. The cubic dispersion relation has the reverse sign from the usual convention; that sign affects both the spectral sectors and the number of required boundary conditions.
Past exam of the mathematics course of the University of Cambridge 2003 iii Paper 74 5 a i Solution Created 2026-10-03 Updated 2026-10-14
The negative drift moves the outflow boundary layer to . The first-order outer equation is . It must satisfy the left boundary condition, giving . This tends to two at the right endpoint.
Use . The leading inner equation is , and impose and on matching. Thus , and a leading uniform asymptotic approximation isFor higher accuracy, expand the drift as in the inner equation, solve the successive inhomogeneous equations, and match to outer corrections satisfying zero left endpoint data. There is no need for a left endpoint layer. Choosing the growing exponential at the wrong endpoint would prevent matching and is the reason the layer location matters.
Past exam of the mathematics course of the University of Cambridge 2003 iii Paper 74 5 a Solution 2026-10-14
The positive drift places the boundary layer at the left endpoint. Set and impose the right boundary condition order by order. The leading ordinary differential equation is , so . At the next order,which integrates toThe outer limit at the left endpoint is two, so it cannot satisfy the left boundary value without a rapidly varying correction.
For the inner expansion use and . Its equations and boundary values areMatching the leading constant gives . The first-order forcing is then . A particular solution is , and matching determines the constant to be . Enforcing givesIn the overlap region the common expansion is . Adding the inner and outer approximations and subtracting that common part gives the uniform asymptotic approximationIt has error uniformly for , apart from an exponentially small endpoint adjustment of the same asymptotic irrelevance. At it equals one exactly, and at it differs from one only exponentially. The polynomial multiplying the exponential remains bounded on the inner scale.
There is an independent check from a first integral. The equation is , so its exact solution has the formwhere is selected by the right boundary value. Endpoint expansion of this integral gives , the same outer coefficients and the same inner expansion. This check also establishes the uniform order of the composite approximation without treating an inner residual as an outer estimate.
Past exam of the mathematics course of the University of Cambridge 2004 iii Paper 60 3 a Solution Created 2026-10-03 Updated 2026-10-14
Write and use modes satisfying all boundary conditions,The quadratic temperature distortion is retained for the saturation calculation; it does not enter the linear onset. Substitution of the primary modes givesThe determinant vanishes whenThe trace is . Thus no conjugate eigenvalues can cross the imaginary axis away from zero, and the initial instability is stationary, not oscillatory. Indeed for the discriminant is . This is quasistatic vertical-field magnetoconvection: there is no independent magnetic evolution mode that could produce a magnetic-thermal Hopf bifurcation instability.
For a fixed box the other allowed horizontal wavenumbers are , and the actual first threshold is the minimum of their corresponding values. Vertical index replaces by and strictly increases the threshold at fixed horizontal wavenumber, so is selected. The displayed formula describes the prescribed fundamental roll; part (b) optimizes its width.
Past exam of the mathematics course of the University of Cambridge 2004 iii Paper 66 1 Solution Created 2026-10-03 Updated 2026-10-14
The separated disk modes. Away from the plane of a razor-thin disk, there is no volume source, so the Newtonian gravitational potential obeys the Laplace equation in cylindrical coordinates. Separation of variables givesHere is a radial wavenumber, with units of inverse length. Regularity on the axis selects the Bessel function of the first kind rather than the singular second solution. Decay away from the disk, continuity across it and reflection symmetry select . These conditions apply to each nonzero-wavenumber mode; a uniform sheet is a separate zero-mode problem.
Integrate the Poisson equation for Newtonian gravity across the disk to obtain the jump conditionFor the unit-amplitude separated mode, the jump is . ThusAn individual basis density is signed; positive physical surface densities are built by superposition.
An arbitrary axisymmetric surface density. Define its order-zero Hankel transform byMatching the coefficients of the mode densities gives the Hankel representation of thin-disk gravity,The boundary conditions exclude added external fields. For sufficiently decaying surface density, this is the isolated-disk potential; for infinite scale-free disks, only potential differences or regularized forces need be finite.
The Mestel disk. Put . The Bessel function integral , understood with a decaying regulator, gives . Differentiate the potential before removing the regulator and use :The last integral is , since . Therefore the circular speed and enclosed mass areConsequently . This equality follows from the specific Mestel disk calculation; the spherical shell theorem cannot be applied to general disks. Its total mass and absolute potential referenced to infinity diverge. A convenient finite midplane reference is .
The rotating distribution function. Use specific energy and . For the printed positive exponential, take and ; has units of squared velocity. If instead means ordinary total energy and , the velocity integral diverges, so that literal convention cannot define the requested galactic distribution function. Equivalently one can retain ordinary energy and choose the printed negative.
With the binding-energy convention, write the amplitude as to distinguish it from the function:Integration over both radial velocities and positive azimuthal velocities givesThe half-line integral follows by setting in the defining Gamma function integral and requires . Matching its radial power and normalization to the Mestel disk yields the rotating Mestel disk distribution function parametersSince , the required integrability condition is automatic. The amplitude uses the stated potential reference; changing the additive potential constant changes without changing the physical distribution.
Past exam of the mathematics course of the University of Cambridge 2004 iii Paper 69 2 a Solution Created 2026-10-03 Updated 2026-10-14
Write , , and assume are sufficiently smooth for the following Taylor expansions. The symmetric half-grid diffusion consistency calculation giveswhereThe normalized local truncation error of the full update is consequentlyThe PDF prints an advection equation, , whereas this stencil approximates conservative diffusion, . For the literal printed equation the leading defect is , which generally does not vanish: the method is inconsistent, with an normalized defect.
For the intended variable-coefficient conservative diffusion equation, the first two terms cancel and . Under the parabolic step scaling , this is , while the unnormalized one-step local truncation error is . These two orders use different normalizations. With numerical stability and compatible smooth initial and boundary conditions, the intended diffusion problem has global error on a fixed time interval. Positivity and boundedness of alone do not justify the displayed Taylor expansion; coefficient regularity is also needed for that error formula.
Past exam of the mathematics course of the University of Cambridge 2004 iii Paper 73 4 Circulation application Solution Created 2026-10-03 Updated 2026-10-14
Use the meridional overturning streamfunction conventionDifferentiating geostrophic and hydrostatic balance gives the thermal wind relation . Differentiate it in time and use the mean momentum and density equations. The resulting diagnostic elliptic boundary value problem isFor the transformed circulation it becomes the Eliassen equation for residual circulation, with signs appropriate to the convention above:At the rigid meridional walls, impermeability requires each streamfunction to be constant along each wall. Set those constants to zero, choosing the solution with no imposed net vertical throughflow. Since the prescribed density flux vanishes at both walls, the same conditions apply to .
There is a genuine data limitation: prescribing alone does not prescribe . The term remains in both circulation equations. To give the usual explicit step-flux circulation, assume vanishing meridional momentum-flux divergence, , with no independently imposed circulation. The formulas below are conditional on that assumption, not a unique consequence of the density flux alone.
Let , and , using the flux unit in which its lower plateau is . The bounded solution has . Its vertical equation isAt vertical infinity require bounded Eulerian flow and vanishing residual flow. Requiring the Eulerian vertical velocity to vanish also at lower infinity would be inconsistent with the continuing lower-region flux divergence. At , continuity of and follows because the forcing has a step but no delta function. The step-flux residual circulation in a stratified channel is thereforeHere , because below zero and vanishes above it. The residual streamfunction decays at both infinities, whereas the Eulerian streamfunction tends to below and to zero above.
For , the Eulerian velocity isThe flow descends near the wall and rises near , with a negative meridional return flow concentrated near . Below the dissipation level its vertical flow balances the continuing divergence of eddy density transport. There is no independently added throughflow, but opposing vertical transports at lower infinity are required by this idealized forcing.
Away from the discontinuity the residual meridional velocity equals the Eulerian value, whileThus the residual mean circulation forms two oppositely rotating cells, confined within vertical distance of the termination of the wave flux. Its streamfunction jumps by at zero. Distributionally this supplies the horizontal return transportThe Dirac delta sheet is a consequence of the discontinuous imposed flux, not an extra boundary condition. Smoothing the flux termination replaces it by a thin finite return current. Streamlines must not be drawn as continuously crossing the jump without that transport.
The physical Eulerian mean flow is continuous. In the transformed momentum equation, the sheet contribution to cancels the sheet in the flux divergence, leaving the smooth mean accelerationThus termination of the vertical Eliassen–Palm flux produces a finite column-integrated wave drag, redistributed over depth by the residual mean circulation.
Eulerian and residual circulations for zero momentum-flux divergence; the residual flux-termination sheet closes two opposite cells
. Details. For completeness, the general bounded response can be calculated without setting to zero. Expand in wall-compatible sine modes, with coefficients , and . Add to either displayed streamfunction the same correctionfor decaying momentum-flux forcing. The Green function sign follows from .
An explicit counterexample to uniqueness is the additional momentum flux , with . It vanishes at both walls and leaves the prescribed density flux unchanged. For it addsto both circulations. This nonzero bounded, wall-impermeable solution verifies that the original mean equations need an additional momentum-flux specification to select a unique sketch.
Past exam of the mathematics course of the University of Cambridge 2004 iii Paper 8 1 iii Solution Created 2026-10-03 Updated 2026-10-14
A field has infinitely many degrees of freedom: its initial state is a function rather than a finite list of coordinates. For a Hamiltonian evolution, the analogue of the finite-dimensional construction seeks an infinite sequence of independent conserved functionals whose Poisson brackets vanish. Their Hamiltonian vector fields then define commuting time evolutions. For example, a Poisson operator produces the field equation , where is the variational derivative.
A Lax equation packages many conservation laws into one operator identity. If , cyclicity of an appropriate trace givesFor a formal pseudodifferential operator the relevant trace is the Adler trace; suitable fractional powers of a monic differential operator produce further conserved functionals. A Lenard-Magri recursion can prove that such functionals are Poisson-commuting functions, rather than merely conserved quantities. A Lax equation by itself does not establish all the independence and completeness properties needed for integrability.
The Korteweg-De Vries equation illustrates the mechanism. In one time normalization take and . Direct composition givesIts spectral problem, together with the inverse scattering transform for decaying data or suitable spectral coordinates for periodic data, turns nonlinear evolution into simple evolution of spectral data. Solitary waves arise from discrete spectral data, while continuous data describe dispersive radiation.
Thus infinite-dimensional Hamiltonian integrability is more than an infinite list of formulas: one wants sufficiently complete commuting invariants and a reconstruction mechanism that solves the evolution. Boundary conditions and the function space are essential. A formal hierarchy need not converge, and an infinite family of conserved quantities does not automatically give a compact infinite-dimensional torus or an unrestricted analogue of the Arnold-Liouville theorem.
Past exam of the mathematics course of the University of Cambridge 2005 ia Paper 2 6B Solution Created 2026-09-24 Updated 2026-10-14
Substitution of the power series, with termwise differentiation, gives the coefficient identityThus and are arbitrary, and all power series solutions of the Hermite differential equation areAn empty product is one. For a nonterminating branch the ratio of consecutive terms in its own parity subsequence is asymptotic to , hence tends to zero for every fixed . The ratio test proves convergence everywhere; the resulting power series can be differentiated termwise and indeed solve the differential equation. A terminating branch is already a polynomial.
The boundary condition forces ; the solution is odd. Apart from the zero solution, which works for every , a polynomial occurs exactly when an odd index satisfies . Before that index the odd coefficients are nonzero, and after it they vanish. Thus the odd polynomial solutions of the Hermite differential equation haveEquivalently, comparison of leading polynomial coefficients gives , proving that no other nonzero polynomial is possible. Normalizing by , the polynomials of degree below six areEach can be multiplied by an arbitrary nonzero constant; these are normalized odd Hermite polynomials.
Let . The forcing in the final inhomogeneous linear differential equation is precisely the displayed degree-five eigenfunction , with . HenceDirect differentiation gives , and . Among polynomial solutions this choice is unique: a polynomial in the kernel of must be constant, and its value at zero eliminates that constant.
Past exam of the mathematics course of the University of Cambridge 2005 iii Paper 20 1 Solution Created 2026-10-03 Updated 2026-10-14
The central problem of spectral geometry is how much of a Riemannian metric is determined by the eigenvalues, with multiplicities, of its Laplace-Beltrami operator. For a membrane with fixed boundary, the squared vibration frequencies are proportional to the Dirichlet eigenvalues. Kac's 1966 paper brought this inverse problem into focus: the full sequence of frequencies might encode a domain's shape, but recovering a few geometric quantities is much weaker than recovering an isometry class.
The first major line of development extracts geometric information from spectral sums. On a smooth closed -dimensional manifold, the heat trace has the short-time expansionwhere are the nonnegative eigenvalues of and is scalar curvature. The coefficients are heat invariants, formed by integrating local curvature expressions. They recover dimension, Riemannian volume and total scalar curvature. For smooth planar domains with the Dirichlet boundary condition,Thus area and boundary length are audible. The constant term also contains topological information through Gauss-Bonnet theorem, with corner corrections when the boundary is polygonal. The Weyl law is the corresponding leading eigenvalue-counting result. Neither the leading term nor the whole list of local heat coefficients automatically reconstructs global geometry.
In the 1970s, the wave trace introduced a complementary link with global closed geodesics. The distributionis spectral, while its singularities are governed by periodic geodesic motion. The 1975 Duistermaat-Guillemin work made this relationship precise. Its trace formula places singularities at geodesic lengths and computes their coefficients under suitable clean or nondegenerate hypotheses. When a coefficient is nonzero, the corresponding length is audible. Possible cancellation must be addressed before identifying the entire length spectrum with the singular support. This approach led to spectral rigidity questions, rather than just the recovery of integral curvature quantities.
The second main line constructs counterexamples to metric determination. A precedent was the 1964 flat-torus example: different sixteen-dimensional Euclidean lattices could give the same Laplacian eigenvalues. Subsequent lattice constructions lowered the dimension; Conway and Sloane's 1992 work exhibited four-dimensional Euclidean lattices with matching lattice theta series. For a flat torus , the spectrum is the multiset for , so equality of lattice theta series of the dual lattices gives isospectrality even when the Euclidean lattices are not related by an orthogonal transformation.
A systematic geometric construction arrived with the Sunada theorem in 1985. For a common Riemannian cover with finite isometry group , almost conjugate subgroups have the same number of invariant vectors in each Laplacian eigenspace. Their freely acting quotients are therefore isospectral manifolds. This reduces an analytic equality to a finite-group representation identity. It produces many examples on surfaces and in higher dimensions, but nonconjugacy of the subgroups must be supplemented by an argument excluding extra quotient isometries. Transplantation theorem provides an explicit version: piecewise eigenfunctions are recombined by a fixed invertible linear transformation that respects the gluing and boundary conditions.
The planar problem itself was answered negatively by Gordon, Webb and Wolpert in 1992: noncongruent simply connected polygonal plane domains have equal Dirichlet spectra. This is stronger than counterexamples among abstract manifolds or flat tori, since these are genuine two-dimensional planar membranes. It does not settle every restricted version, such as asking for two domains within a specified smooth or convex class.
By the 1990s, isospectrality was also known to coexist with continuous metric variation and changes in local geometry. For example, a 1997 construction gives continuous isospectral families on products of spheres and tori whose members need not be locally isometric; in some examples the maximum scalar curvature changes. Conversely, rigidity and compactness results show that an isospectral class is not unrestricted. The 1988 surface compactness theorem makes fixed-spectrum sets of closed orientable surface metrics compact modulo diffeomorphisms. Compactness does not mean there is only one metric, or even finitely many metrics. By 2005, the subject therefore combined local heat invariants, global geodesic information, group-theoretic constructions, explicit transplantation and geometric rigidity. The spectrum determines substantial geometry, while generally failing to determine the full metric.
