







Study with the several resources on Docsity
Earn points by helping other students or get them with a premium plan
Prepare for your exams
Study with the several resources on Docsity
Earn points to download
Earn points by helping other students or get them with a premium plan
Advanced problems and solutions from various fields including dynamics, functional analysis, electromagnetism, graph theory, stochastic financial models, computational statistics, quantum mechanics, general relativity, statistical physics, transport processes, theoretical geophysics, and numerical analysis. It covers a wide range of topics such as lorentz force law, poisson brackets, magnetic monopoles, wellfounded relations, extremal graph theory, mean-variance approach, inference in multi-parameter exponential families, asymptotic properties of maximum-likelihood estimators, bootstrap inference, time-dependent schrödinger equation, density of photons, einstein a and b coefficients, weak gravitational fields, bose-einstein distribution, hot big bang theory, fluid flows, advection-diffusion equation, phase speed and group velocity, internal gravity waves, perturbation theory, and computation of eigenvalues and eigenvectors.
Typology: Exams
1 / 13
This page cannot be seen from the preview
Don't miss anything!








Thursday 7 June 2001 1.30 to 4.
Before you begin read these instructions carefully.
Candidates must not attempt more than FOUR questions.
The number of marks for each question is the same. Additional credit will be given for a substantially complete answer.
Write legibly and on only one side of the paper.
Begin each answer on a separate sheet.
At the end of the examination:
Tie your answers in separate bundles, marked A, B, C,... , L according to the letter affixed to each question. (For example, 3A, 9A should be in one bundle and 4C, 10C in another bundle.)
Attach a completed cover sheet to each bundle.
Complete a master cover sheet listing all questions attempted.
It is essential that every cover sheet bear the candidate’s examination number and desk number.
1D Markov Chains
Write an essay on the convergence to equilibrium of a discrete-time Markov chain on a countable state-space. You should include a discussion of the existence of invariant distributions, and of the limit theorem in the non-null recurrent case.
2H Principles of Dynamics
(i) Consider a particle of charge q and mass m, moving in a stationary magnetic field B. Show that Lagrange’s equations applied to the Lagrangian
mr˙^2 + q r˙ · A(r),
where A is the vector potential such that B = curl A, lead to the correct Lorentz force law. Compute the canonical momentum p, and show that the Hamiltonian is H = 12 mr˙^2.
(ii) Expressing the velocity components r˙i in terms of the canonical momenta and co-ordinates for the above system, derive the following formulae for Poisson brackets:
(a) {F G, H} = F {G, H} + {F, H} G, for any functions F , G, H;
(b) {m r˙i, m r˙j } = qijkBk;
(c) {m r˙i, rj } = −δij ; (d) {m r˙i, f (rj )} = − (^) ∂r∂i f (rj ).
Now consider a particle moving in the field of a magnetic monopole,
Bi = g
ri r^3
Show that {H, J} = 0, where J = mr ∧ r˙ − gqˆr. Explain why this means that J is conserved.
Show that, if g = 0, conservation of J implies that the particle moves in a plane perpendicular to J. What type of surface does the particle move on if g 6 = 0?
Paper 4
5J Electromagnetism
Write down the form of Ohm’s Law that applies to a conductor if at a point r it is moving with velocity v(r).
Use two of Maxwell’s equations to prove that ∫
C
(E + v ∧ B) · dr = −
d dt
S
B · dS ,
where C(t) is a moving closed loop, v is the velocity at the point r on C, and S is a surface spanning C. The time derivative on the right hand side accounts for changes in both C and B. Explain briefly the physical importance of this result.
Find and sketch the magnetic field B described in the vector potential
A(r, θ, z) = (0, 12 brz, 0)
in cylindrical polar coordinates (r, θ, z) , where b > 0 is constant.
A conducting circular loop of radius a and resistance R lies in the plane z = h(t) with its centre on the z-axis.
Find the magnitude and direction of the current induced in the loop as h(t) changes with time, neglecting self-inductance.
At time t = 0 the loop is at rest at z = 0. For time t > 0 the loop moves with constant velocity dh/dt = v > 0. Ignoring the inertia of the loop, use energy considerations to find the force F (t) necessary to maintain this motion.
[In cylindrical polar coordinates
curl A =
r
∂Az ∂θ
∂Aθ ∂z
∂Ar ∂z
∂Az ∂r
r
∂r
(rAθ ) −
r
∂Ar ∂θ
6K Dynamics of Differential Equations Write a short essay about periodic orbits in flows in two dimensions. Your essay should include criteria for the existence and non-existence of periodic orbits, and should mention (with sketches) at least two bifurcations that create or destroy periodic orbits in flows as a parameter is altered (though a detailed analysis of any bifurcation is not required).
7C Geometry of Surfaces Write an essay on the Gauss-Bonnet theorem. Make sure that your essay contains a precise statement of the theorem, in its local form, and a discussion of some of its applications, including the global Gauss-Bonnet theorem.
Paper 4
8B Logic, Computation and Set Theory
What is a wellfounded relation, and how does wellfoundedness underpin wellfounded induction?
A formula φ(x, y) with two free variables defines an ∈-automorphism if for all x there is a unique y, the function f , defined by y = f (x) if and only if φ(x, y), is a permutation of the universe, and we have (∀xy)(x ∈ y ↔ f (x) ∈ f (y)).
Use wellfounded induction over ∈ to prove that all formulæ defining ∈-automorphisms are equivalent to x = y.
9A Graph Theory Write an essay on extremal graph theory. Your essay should include proofs of at least two major results and a discussion of variations on these results or their proofs.
10C Number Theory
Attempt one of the following:
(i) Discuss pseudoprimes and primality testing. Find all bases for which 57 is a Fermat pseudoprime. Determine whether 57 is also an Euler pseudoprime for these bases. (ii) Write a brief account of various methods for factoring large numbers. Use Fermat factorization to find the factors of 10033. Would Pollard’s p − 1 method also be practical in this instance?
(iii) Show that
1 /pn is divergent, where pn denotes the n-th prime.
Write a brief account of basic properties of the Riemann zeta-function. State the prime number theorem. Show that it implies that for all sufficiently large positive integers n there is a prime p satisfying n < p 6 2 n.
11E Algorithms and Networks
State the optimal distribution problem. Carefully describe the simplex-on-a-graph algorithm for solving optimal distribution problems when the flow in each arc in the network is constrained to lie in the interval [0, ∞). Explain how the algorithm can be initialised if there is no obvious feasible solution with which to begin. Describe the adjustments that are needed for the algorithm to cope with more general capacity constraints x(j) ∈ [c−(j), c+(j)] for each arc j (where c±(j) may be finite or infinite).
Paper 4 [TURN OVER
14E Computational Statistics and Statistical Modelling
(i) Assume that independent observations Y 1 ,... , Yn are such that
Yi ∼ Binomial(ti, πi), log
πi 1 − πi
= βT^ xi for 1 6 i 6 n ,
where x 1 ,... , xn are given covariates. Discuss carefully how to estimate β, and how to test that the model fits.
(ii) Carmichael et al. (1989) collected data on the numbers of 5-year old children with “dmft”, i.e. with 5 or more decayed, missing or filled teeth, classified by social class, and by whether or not their tap water was fluoridated or non-fluoridated. The numbers of such children with dmft, and the total numbers, are given in the table below:
dmft Social Class Fluoridated Non-fluoridated I 12/117 12/ II 26/170 48/ III 11/52 29/ Unclassified 24/118 49/
A (slightly edited) version of the R output is given below. Explain carefully what model is being fitted, whether it does actually fit, and what the parameter estimates and Std. Errors are telling you. (You may assume that the factors SClass (social class) and Fl (with/without) have been correctly set up.)
Call:
glm(formula = Yes/Total ∼ SClass + Fl, family = binomial, weights = Total)
Coefficients:
Estimate Std. Error z value (Intercept) -2.2716 0.2396 -9. SClassII 0.5099 0.2628 1. SClassIII 0.9857 0.3021 3. SClassUnc 1.0020 0.2684 3. Flwithout 1.0813 0.1694 6.
(Dispersion parameter for binomial family taken to be 1)
Null deviance: 68.53785 on 7 degrees of freedom
Residual deviance: 0.64225 on 3 degrees of freedom
Number of Fisher Scoring iterations: 3
Here ‘Yes’ is the vector of numbers with dmft, taking values 12, 12 ,... , 24 , 49, ‘Total’ is the vector of Total in each category, taking values 117, 56 ,... , 118 , 104, and SClass, Fl are the factors corresponding to Social class and Fluoride status, defined in the obvious way.
Paper 4 [TURN OVER
15F Foundations of Quantum Mechanics
(i) The two states of a spin- 12 particle corresponding to spin pointing along the z axis are denoted by | ↑〉 and | ↓〉. Explain why the states
| ↑, θ〉 = cos
θ 2
| ↑〉 + sin
θ 2
| ↓〉, | ↓, θ〉 = − sin
θ 2
| ↑〉 + cos
θ 2
correspond to the spins being aligned along a direction at an angle θ to the z direction.
The spin-0 state of two spin- 12 particles is
Show that this is independent of the direction chosen to define | ↑〉 1 , 2 , | ↓〉 1 , 2. If the spin of particle 1 along some direction is measured to be 12 ℏ show that the spin of particle 2 along the same direction is determined, giving its value.
[The Pauli matrices are given by
σ 1 =
, σ 2 =
0 −i i 0
, σ 3 =
(ii) Starting from the commutation relation for angular momentum in the form
[J 3 , J±] = ±ℏJ±, [J+, J−] = 2ℏJ 3 ,
obtain the possible values of j, m, where mℏ are the eigenvalues of J 3 and j(j + 1)ℏ^2 are the eigenvalues of J^2 = 12 (J+J− + J−J+) + J 32. Show that the corresponding normalized eigenvectors, |j, m〉, satisfy
J±|j, m〉 = ℏ ((j ∓ m)(j ± m + 1))^1 /^2 |j, m± 1 〉,
and that 1 n!
J (^) −n |j, j〉 = ℏn
(2j)! n!(2j − n)!
|j, j−n〉, n ≤ 2 j.
The state |w〉 is defined by
|w〉 = ewJ−/ℏ|j, j〉,
for any complex w. By expanding the exponential show that 〈w|w〉 = (1 + |w|^2 )^2 j^. Verify that e−wJ−/ℏJ 3 ewJ−/ℏ^ = J 3 − wJ−,
and hence show that
J 3 |w〉 = ℏ
j − w
∂w
|w〉.
If H = αJ 3 verify that |eiαt〉 e−ijαt^ is a solution of the time-dependent Schr¨odinger equation.
Paper 4
18J Statistical Physics and Cosmology
(i) Given that g(p)dp is the number of eigenstates of a gas particle with momentum between p and p + dp, write down the Bose-Einstein distribution ¯n(p) for the average number of particles with momentum between p and p + dp, as a function of temperature T and chemical potential μ.
Given that μ = 0 and g(p) = 8π V p
2 h^3 for a gas of photons, obtain a formula for the energy density ρT at temperature T in the form
ρT =
0
T (ν)dν,
where T (ν) is a function of the photon frequency ν that you should determine. Hence show that the value νpeak of ν at the maximum of T (ν) is proportional to T.
A thermally isolated photon gas undergoes a slow change of its volume V. Why is ¯n(p) unaffected by this change? Use this fact to show that V T 3 remains constant.
(ii) According to the “Hot Big Bang” theory, the Universe evolved by expansion from an earlier state in which it was filled with a gas of electrons, protons and photons (with ne = np) at thermal equilibrium at a temperature T such that
2 mec^2 kT B ,
where me is the electron mass and B is the binding energy of a hydrogen atom. Why must the composition have been different when kT 2 mec^2? Why must it change as the temperature falls to kT B? Why does this lead to a thermal decoupling of radiation from matter?
The baryon number of the Universe can be taken to be the number of protons, either as free particles or as hydrogen atom nuclei. Let nb be the baryon number density and nγ the photon number density. Why is the ratio η = nb/nγ unchanged by the expansion of the universe? Given that η 1, obtain an estimate of the temperature TD at which decoupling occurs, as a function of η and B. How does this decoupling lead to the concept of a “surface of last scattering” and a prediction of a Cosmic Microwave Background Radiation (CMBR)?
Paper 4
19H Transport Processes
Fluid flows in the x-direction past the infinite plane y = 0 with uniform but time- dependent velocity U (t) = U 0 G (t/t 0 ), where G is a positive function with timescale t 0. A long region of the plane, 0 < x < L, is heated and has temperature T 0 (1 + γ (x/L)n), where T 0 , γ, n are constants [γ = O(1)]; the remainder of the plane is insulating (Ty = 0). The fluid temperature far from the heated region is T 0. A thermal boundary layer is formed over the heated region. The full advection–diffusion equation for temperature T (x, y, t) is
Tt + U (t)Tx = D (Tyy + Txx) , (1)
where D is the thermal diffusivity. By considering the steady case (G ≡ 1), derive a scale for the thickness of the boundary layer, and explain why the term Txx in (1) can be neglected if U 0 L/D 1.
Neglecting Txx, use the change of variables
τ =
t t 0
, ξ =
x L
, η = y
U (t) Dx
= γ
( (^) x L
)n f (ξ, η, τ )
to transform the governing equation to
fηη +
ηfη − nf = ξfξ +
Lξ t 0 U 0
Gτ 2 G^2
ηfη +
fτ
Write down the boundary conditions to be satisfied by f in the region 0 < ξ < 1.
In the case in which U is slowly-varying, so = (^) t 0 LU 0 1, consider a solution for f in the form f = f 0 (η) + f 1 (ξ, η, τ ) + O
Explain why f 0 is independent of ξ and τ.
Henceforth take n = 12. Calculate f 0 (η) and show that
f 1 =
Gτ ξ G^2
g(η) ,
where g satisfies the ordinary differential equation
g′′^ +
ηg′^ −
g =
−η 4
η
e−u
(^2) / 4 du.
State the boundary conditions on g(η).
The heat transfer per unit length of the heated region is −DTy |y=0. Use the above results to show that the total rate of heat transfer is
T 0 [DLU (t)]^1 /^2
γ 2
π −
Gτ G^2
g′(0) + O
Paper 4 [TURN OVER
21H Mathematical Methods
The equation Ax = λx ,
where A is a real square matrix and x a column vector, has a simple eigenvalue λ = μ with corresponding right-eigenvector x = X. Show how to find expressions for the perturbed eigenvalue and right-eigenvector solutions of
Ax + b(x) = λx , || 1 ,
to first order in , where b is a vector function of x. State clearly any assumptions you make.
If A is (n×n) and has a complete set of right-eigenvectors X(j), j = 1, 2 , ...n, which span Rn^ and correspond to separate eigenvalues μ(j), j = 1, 2 , ...n, find an expression for the first-order perturbation to X(1)^ in terms of the
X(j)
and the corresponding left- eigenvectors of A.
Find the normalised eigenfunctions and eigenvalues of the equation
d^2 y dx^2
with y(0) = y(1) = 0. Let these be the zeroth order approximations to the eigenfunctions of d^2 y dx^2
with y(0) = y(1) = 0 and where b is a function of y. Show that the first-order perturbations of the eigenvalues are given by
λ(1) n = −
0
sin(nπx) b
2 sin nπx
dx.
22K Numerical Analysis
Write an essay on the computation of eigenvalues and eigenvectors of matrices.
Paper 4