BAAO 2020, Теория — Задача 3. Inspirals of Exoplanets Near Tidal Destruction

British Astronomy and Astrophysics Olympiad 2019-2020 — Astronomy & Astrophysics Competition Paper — Monday 20th January 2020 · 20 януари 2020 г. · 30 т.

Указания

This question paper must not be photographed or taken out of the exam room

Time: 3 hours (approximately 60 minutes per question).

Questions: All three questions should be attempted. The last question is worth slightly fewer marks than the rest.

Solutions: Answers and calculations are to be written on loose paper. Students should ensure their name and school is clearly written on all answer sheets and pages are numbered. EACH QUESTION ANSWERED must be started on a new page. A standard formula booklet with standard physical constants should be supplied.

Instructions: To accommodate students sitting the paper at different times, please do not discuss any aspect of the paper on the internet until 8 am Saturday 25th January.

Clarity: Solutions must be written legibly, in black pen (the papers are photocopied), and working down the page. Scribble will not be marked and overall clarity is an important aspect of this exam.

Eligibility: The International Olympiad will be held during September 2020; all sixth form students are eligible to participate, even if they will be attending university from October.

Calculators: Any standard calculator may be used, but calculators cannot be programmable and must not have symbolic algebra capability.

Training Dates and the IOAA (Bogota, Columbia, 13th to 22nd September 2020)

The best students taking this paper that are eligible to represent the UK at the IOAA will be invited to attend the Training Camp to be held in the Physics Department at the University of Oxford, (Friday 27th March to Tuesday 31st March 2020). Astronomy material will be covered; problem solving skills and observational skills (telescope and naked eye observations) will be developed. At the Training Camp a data analysis exam and a short theory paper will be sat. A team of five students (plus one reserve) will be selected for further training. From April there will be mentoring by email to cover some topics and problems, followed by additional training camps in the summer.

Important Constants
ConstantSymbolValue
Speed of lightcc3.00×108 m s−13.00 \times 10^{8}\ \mathrm{m\,s^{-1}}
Earth’s rotation period1 day24 hours
Earth’s orbital period1 year365.25 days
parsecpc3.09×1016 m3.09 \times 10^{16}\ \mathrm{m}
Astronomical Unitau1.50×1011 m1.50 \times 10^{11}\ \mathrm{m}
Semi-major axis of the Earth’s orbit1 au
Radius of the SunR⊙R_{\odot}6.96×108 m6.96 \times 10^{8}\ \mathrm{m}
Radius of the EarthR⊕R_{\oplus}6.37×106 m6.37 \times 10^{6}\ \mathrm{m}
Mass of the SunM⊙M_{\odot}1.99×1030 kg1.99 \times 10^{30}\ \mathrm{kg}
Mass of the EarthM⊕M_{\oplus}5.97×1024 kg5.97 \times 10^{24}\ \mathrm{kg}
Luminosity of the SunL⊙L_{\odot}3.85×1026 W3.85 \times 10^{26}\ \mathrm{W}
Stephan-Boltzmann constantσ\sigma5.67×10−8 J m−2 K−45.67 \times 10^{-8}\ \mathrm{J\,m^{-2}\,K^{-4}}
Gravitational constantGG6.67×10−11 m3 kg−1 s−26.67 \times 10^{-11}\ \mathrm{m^{3}\,kg^{-1}\,s^{-2}}
Boltzmann constantkBk_{\mathrm{B}}1.38×10−23 J K−11.38 \times 10^{-23}\ \mathrm{J\,K^{-1}}
Permittivity of free spaceε0\varepsilon_08.85×10−12 F m−18.85 \times 10^{-12}\ \mathrm{F\,m^{-1}}
Permeability of free spaceμ0\mu_04π×10−7 H m−14\pi \times 10^{-7}\ \mathrm{H\,m^{-1}}
Planck’s constanthh6.63×10−34 J s6.63 \times 10^{-34}\ \mathrm{J\,s}
Elementary chargeee1.60×10−19 C1.60 \times 10^{-19}\ \mathrm{C}
Proton rest massmpm_{\mathrm{p}}1.67×10−27 kg1.67 \times 10^{-27}\ \mathrm{kg}
Electron rest massmem_{\mathrm{e}}9.11×10−31 kg9.11 \times 10^{-31}\ \mathrm{kg}
Wien’s displacement lawλmax⁡T\lambda_{\max} T2.90×10−3 m K2.90 \times 10^{-3}\ \mathrm{m\,K}

Basic calculus formulae:

Chain ruleddxf(g(x))=f′(g(x))g′(x)\frac{\mathrm{d}}{\mathrm{d}x} f(g(x)) = f'(g(x))g'(x)
Product ruleddx(uv)=dudxv+udvdx\frac{\mathrm{d}}{\mathrm{d}x}(uv) = \frac{\mathrm{d}u}{\mathrm{d}x}v + u\frac{\mathrm{d}v}{\mathrm{d}x}
Quotient ruleddx(uv)=dudxv−udvdxv2\frac{\mathrm{d}}{\mathrm{d}x}\left(\frac{u}{v}\right) = \frac{\frac{\mathrm{d}u}{\mathrm{d}x}v - u\frac{\mathrm{d}v}{\mathrm{d}x}}{v^2}
Integration by parts∫udvdx dx=uv−∫vdudx dx\int u\frac{\mathrm{d}v}{\mathrm{d}x}\,\mathrm{d}x = uv - \int v\frac{\mathrm{d}u}{\mathrm{d}x}\,\mathrm{d}x
Standard integral∫1x dx=ln⁡∣x∣+C\int \frac{1}{x}\,\mathrm{d}x = \ln\lvert x\rvert + C
Important Formulae

You might find the diagram of an elliptical orbit below useful in solving some of the questions:

Elements of an elliptic orbit:

a=OA (=OP)a = \mathbf{OA}\ (= \mathbf{OP})semi-major axis
b=OB (=OC)b = \mathbf{OB}\ (= \mathbf{OC})semi-minor axis
e=1−b2a2e = \sqrt{1 - \frac{b^2}{a^2}}eccentricity
Ffocus
Pperiapsis (point nearest to F)
Aapoapsis (point furthest from F)

Kepler’s Third Law: For an elliptical orbit, the square of the period, TT, of an object about the focus is proportional to the cube of the semi-major axis, aa (as defined above), such that

T2=4π2GMa3 ,T^2 = \frac{4\pi^2}{GM}a^3\,,

where MM is the total mass of the system (typically dominated by the central object) and GG is the universal gravitational constant.

Vis-Viva Equation: For an elliptical orbit, the speed vv of an object at a distance rr from the focus is related to the semi-major axis, aa, total mass of the system, MM, and universal gravitational constant, GG, (as defined above), such that

v2=GM(2r−1a) .v^2 = GM\left(\frac{2}{r} - \frac{1}{a}\right)\,.

Magnitudes: The apparent magnitudes of two objects, m1m_1 and m0m_0, are related to their apparent brightnesses, b1b_1 and b0b_0, via the formula

b1b0=10−0.4(m1−m0) .\frac{b_1}{b_0} = 10^{-0.4(m_1 - m_0)}\,.

The absolute magnitude of an object, M\mathcal{M}, is the same as its apparent magnitude when viewed from 10 pc, hence the relationship between apparent and absolute magnitude and distance is

m−M=5log⁡(d10),m - \mathcal{M} = 5\log\left(\frac{d}{10}\right),

where dd is measured in parsecs.

Note for markers

BAAO 2019/20 Solutions and Marking Guidelines

Note for markers:

  • Answers to two or three significant figures are generally acceptable. The solution may give more in order to make the calculation clear. Units should be present on final answers when appropriate.
  • There are multiple ways to solve some of the questions; please accept all good solutions that arrive at the correct answer. If a candidate gets the final (numerical) answer then allow them all the marks for that part of the question (as indicated in red), so long as there are no unphysical / nonsensical steps or assumptions made.

Условие

Some of the very first exoplanets to be discovered in large surveys were dubbed ‘hot Jupiters’ as they were similar in mass to Jupiter (i.e. a gas giant) but were much closer to their star than Mercury is to the Sun (and hence are in a very hot environment). Planetary formation models suggest that they were unlikely to have formed there, but instead formed much further out from the star and migrated inwards, due to gravitational interactions with other planets in the system. Studies of ‘hot Jupiters’ show that there is an overabundance of them with periods of ∼3−4\sim 3 - 4 days, and very few with periods shorter than that. Since large, close-in planets should be the easiest to detect in all of the main methods of finding exoplanets, this scarcity is likely to be a real effect and suggests that exoplanets which are that close to their star are in a relatively rapid (by astronomical standards) inspiral towards destruction by their star.

Left: scatter plot with error bars of a/a_Roche (1.0 to 4.5) against M_star/M⊙ (0.6 to 1.6), a dotted horizontal line at 2.0; labelled points below it include WASP-43b, WASP-52b, WASP-4b, CoRoT-1b, OGLE-TR-56b, WASP-10b, HAT-P-32b, WASP-19b, WASP-12b and WASP-18b, with WTS-2b marked by a red square near 1.4. Right: logarithmic plot of T_shift / s against M_star / M⊙ (0.4 to 1.6) with points labelled WASP-18b, WASP-103b, WASP-19b, CoRoT-14b, OGLE-TR-56b, WASP-43b, WASP-78b, HAT-P-7b and WTS-2b, and a dotted horizontal detection line near 5 s.
Figure 6: Left: The orbital radius of several ‘hot Jupiters’ scaled by the Roche radius of the system (where tidal forces would destroy the planet). There is an expected pile up close to radii double the Roche radius (dotted line), and very few with radii smaller than that - those that are will inevitably spiral into the star and be destroyed by the tidal forces when they get too close. Credit: Birkby et al. (2014). Right: As the planets inspiral we should see a shift in when their transits occur. This figure shows the predicted size of the shift after a period of 10 years if the tidal dissipation quality factor Q⋆′=106Q'_{\star} = 10^{6}, as well as the current detection limit of 5 seconds (dotted line). Therefore measuring if there is any shift in the transit times over the course of a decade of observations can put stringent limits on the value of Q⋆′Q'_{\star}. Credit: Birkby et al. (2014).

The Roche radius, where a planet will be torn apart due to the tidal forces acting on it, is defined as

aRoche≈2.16RP(M⋆MP)1/3 ,a_{\mathrm{Roche}} \approx 2.16R_P\left(\frac{M_{\star}}{M_P}\right)^{1/3}\,,

where RPR_P is the radius of the planet, MPM_P is the mass of the planet and M⋆M_{\star} is the mass of the star. If a gas giant is knocked into a highly elliptical orbit (i.e. e≈1e \approx 1) that has a periapsis rperi<aRocher_{\mathrm{peri}} < a_{\mathrm{Roche}} then it will not survive. However, if the periapsis just grazes the Roche radius (rperi≈aRocher_{\mathrm{peri}} \approx a_{\mathrm{Roche}}) then the orbit will rapidly circularise. By conserving angular momentum, it can be shown that the circular orbit will have a radius a=2aRochea = 2a_{\mathrm{Roche}} (see the left panel of Fig 6). Exoplanets observed to be in an orbit with a radius less than that will be unstable and angular momentum will be transferred from the planet to the star, causing the star to spin more rapidly and the planet’s orbital radius to decrease. Eventually this will result in the planet’s orbit crossing the Roche radius and being destroyed by the tidal forces.

The duration of this inspiral will be dependent on how well the star can dissipate the orbital energy through frictional processes within the star, and can be parameterised by the tidal dissipation quality factor, Q⋆′Q'_{\star}. By looking for changes in the orbital period of the planet, detectable by shifts in the timing of transits by the planet in front of the star, we can determine an estimate of Q⋆′Q'_{\star}, which hence tells us about the internal structure of stars. These ‘hot Jupiters’ are the best laboratory we have for this, as they are the most likely to produce a measurable shift (i.e. ∼5\sim 5 s) in transit times within only ∼10\sim 10 years (see the right panel of Fig 6). We will try and reproduce these results in this question.

The WTS-2 system consists of a star of mass M⋆=0.820 M⊙M_{\star} = 0.820\,M_{\odot}, peak in its black-body spectrum at λmax⁡=580\lambda_{\max} = 580 nm, and distance from us of 1.03 kpc, with an orbiting planet (called WTS-2b) with a period P=1.0187P = 1.0187 days, mass of 1.12 MJM_J and radius 1.36 RJR_J. The mass and radius of Jupiter are MJ=1.90×1027M_J = 1.90 \times 10^{27} kg and RJ=7.15×107R_J = 7.15 \times 10^{7} m respectively.

a. Show that WTS-2b has an orbital radius of ∼1.4 aRoche\sim 1.4\ a_{\mathrm{Roche}}. [4 т.]

b. (i) Given the apparent magnitude of WTS-2 in the visible is m=16.14m = 16.14 and the absolute magnitude of the Sun in the same part of the EM spectrum is M⊙=4.83\mathcal{M}_{\odot} = 4.83:

(i) Calculate the luminosity of the star, L⋆L_{\star}. Give your answer in units of L⊙L_{\odot}. [3 т.]

b. (ii) Hence work out the radius of the star, R⋆R_{\star}. Give your answer in units of R⊙R_{\odot}. [3 т.]

b. (iii) Show that the incident flux (in W m−2\mathrm{W\,m^{-2}}) on the planet is ∼103\sim 10^{3} times larger than what we receive on Earth from the Sun. [3 т.]

The change in the semi-major axis of the planet, aa, due to tidal forces is given by

∣a˙a∣=6k2ΔtMPM⋆(R⋆a)5n2\left|\frac{\dot{a}}{a}\right| = 6k_2\Delta t\frac{M_P}{M_{\star}}\left(\frac{R_{\star}}{a}\right)^5 n^2

where the dot notation is used to indicate the differential with respect to time (i.e. a˙≡da/dt\dot{a} \equiv \mathrm{d}a/\mathrm{d}t), k2k_2 is a constant related to the density structure of the star, Δt\Delta t is the (assumed constant) time lag between where the planet is in its orbit and the location of the tidal bulge on the star, and n=2π/Pn = 2\pi/P. By separating variables and integrating this equation, an expression can be derived for the time it takes for aa to decrease to zero. This is known as the inspiral time, τ\tau, and even though the planet will be destroyed when a=aRochea = a_{\mathrm{Roche}} the time to go from a=aRochea = a_{\mathrm{Roche}} to a=0a = 0 is negligible in comparison to the time to get to a=aRochea = a_{\mathrm{Roche}}, so τ\tau is a good estimate of the remaining lifetime of the planet.

c. Derive an equation for τ\tau, showing τ∝a8\tau \propto a^8. Hint: nn is a function of aa due to Kepler’s third law. [4 т.]

The tidal dissipation quality factor is defined as

Q⋆′=34k2ΔtnQ'_{\star} = \frac{3}{4k_2\Delta t n}

and so absorbs some of our unknown parameters about the stellar structure. Based upon studies of binary stars it seems a value of Q⋆′∼106−107Q'_{\star} \sim 10^{6} - 10^{7} is reasonable, although it is poorly constrained. The size of the shift in transit timings, TshiftT_{\mathrm{shift}}, over an observation period of TT seconds is

Tshift=12nT2(dndt) .T_{\mathrm{shift}} = \frac{1}{2n}T^2\left(\frac{\mathrm{d}n}{\mathrm{d}t}\right)\,.

d. (i) For the planet WTS-2b, assuming Q⋆′=106Q'_{\star} = 10^{6}:

(i) Show that τ∝Q⋆′nM⋆MP(aR⋆)5\tau \propto \frac{Q'_{\star}}{n}\frac{M_{\star}}{M_P}\left(\frac{a}{R_{\star}}\right)^5, finding the constant of proportionality (as a fraction), and hence the remaining lifetime of the planet. Give your answer in Myr. [1 Myr = 10610^{6} years] [5 т.]

d. (ii) Show that Tshift∝T2τT_{\mathrm{shift}} \propto \frac{T^2}{\tau}, finding the constant of proportionality, and hence verify that TshiftT_{\mathrm{shift}} is measurable (i.e. >5> 5 s) if T=10T = 10 years. Hint: Using the chain rule, dndt=dndaa˙\frac{\mathrm{d}n}{\mathrm{d}t} = \frac{\mathrm{d}n}{\mathrm{d}a}\dot{a}. [8 т.]

Отговори

Решение

Questions proposed by

END OF PAPER

Questions proposed by:

Dr Alex Calverley (Royal Grammar School, Guildford)

Alan Sun (University of Cambridge)

Navonil Neogi (University of Cambridge)


Оригинал в Архива: BAAO_2020.pdf · официални решения: BAAO_2020_solutions.pdf