IOAA 2025, Анализ на данни — Задача 1. D01 · 30 Years of Exoplanets

18th IOAA, Mumbai 2025 — Data Analysis Examination · 90 т.

Внимание: Бележка към темата

The obtained Summary Answersheet contains D02 only. Referenced D01 response forms were unavailable; all D01 data graphs printed in the question and solution booklets are retained.

Условие

This problem explores some aspects of the two main methods of exoplanet detection: radial velocity and transit. Throughout this problem we shall consider a particular system of a single planet (P) in a circular orbit with radius aa around a solar-type star (S). We shall refer to this system as the “SP system”.

Stellar parameters

(D01.1) [8 marks] The V-band apparent magnitude of the star S is (7.65±0.03) mag(7.65\pm0.03)\ \mathrm{mag}, the parallax is (20.67±0.05)(20.67\pm0.05) milliarcsecond and the bolometric correction (BC) is −0.0650 mag-0.0650\ \mathrm{mag}. Thus the star has a higher bolometric luminosity than its luminosity in the V-band.

Estimate the mass of the star, MsM_s (in units of M⊙M_\odot), assuming a mass-luminosity (MM–LL) relation of the form L∝M4L\propto M^4. Also estimate the uncertainty in MsM_s. You may need ddxln⁡x=1x\frac d{dx}\ln x=\frac1x.

Radial Velocity method

The radial velocity method uses the Doppler shift δλ≡λobs−λ0\delta\lambda\equiv\lambda_{\mathrm{obs}}-\lambda_0 between the observed wavelength λobs\lambda_{\mathrm{obs}} and the rest wavelength λ0\lambda_0 of a known spectral line to detect an exoplanet and determine its characteristics.

The figure below shows the δλ\delta\lambda for the Fe I line (λ0=543.45×10−9 m)(\lambda_0=543.45\times10^{-9}\ \mathrm{m}) as a function of time as observed for the SP system.

Original observed Fe I wavelength-shift data, every point, full grid, wavelength axis and Days axis.

The radial velocity semi-amplitude KK is defined as K≡(vr,max⁡−vr,min⁡)/2K\equiv(v_{r,\max}-v_{r,\min})/2 where vr,max⁡v_{r,\max} and vr,min⁡v_{r,\min} are the minimum and maximum radial velocities, respectively. For a circular planetary orbit the semi-amplitude KK can be written as:

K=(2πGT)1/3Mpsin⁡i(Mp+Ms)2/3K=\left(\frac{2\pi G}T\right)^{1/3}\frac{M_p\sin i}{(M_p+M_s)^{2/3}}

where TT is the period, ii is the inclination of the planetary orbit (angle between the normal to the orbital plane of the planet and the line of sight of the observer), MpM_p and MsM_s are the masses of the planet and the star, respectively.

(D01.2) Use the graph given in the Summary Answersheet to answer the following.

(D01.2a) [2 marks] Draw a smooth curve associated with the observed data shown in the graph.

(D01.2b) [11 marks] Select appropriate points on your drawn curve and use suitable methods to determine TT and KK along with respective uncertainties. All data points used for the calculation of TT and KK must be shown in the table in the Summary Answersheet. Use the rest of the table to show your intermediate calculations, as needed, with appropriate headers.

(D01.2c) [5 marks] Find the minimum mass of the planet Mp,min⁡M_{p,\min} (in M⊙M_\odot), and its corresponding uncertainty, assuming Mp≪MsM_p\ll M_s.

(D01.2d) [4 marks] Using the value of Mp,min⁡M_{p,\min} estimated in part (D01.2c), calculate the minimum value of the semi-major axis of the planet’s orbit, amin⁡a_{\min}, in au and its uncertainty.

Transit method (without limb darkening)

The schematic diagram of a planet transit (not drawn to scale) is shown below. Initially, we shall assume the stellar disk to have a uniform average intensity with some intrinsic noise due to the star itself.

Original face-on and side-on transit geometry and full schematic lightcurve: all planet positions, radii, b, i, a, observer, Delta, tT and tF.

The lightcurve of the normalized intensity, II, as a function of time tt is shown in the schematic diagram of the transit above. The average stellar intensity outside the transit is taken as unity. The maximum decrease in the intensity is given by Δ\Delta in the normalized light curve. For a uniformly bright stellar disk, the radius of the planet, RpR_p, is related to Δ\Delta as

(RpRs)2=Δ,\left(\frac{R_p}{R_s}\right)^2=\Delta,

where RsR_s is the radius of the star.

The total duration of transit (when part or all of the planet covers the stellar disk) is given by tTt_T, while tFt_F gives the duration when the planet is fully in front of the stellar disk. The “impact parameter” bb is the projected distance between the planet and centre of the stellar disk at the mid-point of the transit, in units of the stellar radius, RsR_s.

For a nearly edge-on star-planet orbit, the impact parameter is given by the formula

b=[(1−Δ)2−(tF/tT)2(1+Δ)21−(tF/tT)2]1/2b=\left[\frac{(1-\sqrt\Delta)^2-(t_F/t_T)^2(1+\sqrt\Delta)^2}{1-(t_F/t_T)^2}\right]^{1/2}

(D01.3) [3 marks] For the SP system, the stellar radius is known to be Rs=1.20 R⊙R_s=1.20\,R_\odot, and the transit of the planet is indeed visible. Using the minimum orbital radius, amin⁡a_{\min}, estimated in part (D01.2d), find the minimum value, imin⁡i_{\min}, of the inclination angle.

Assuming a stellar disk of uniform brightness, the transit lightcurve would look like as shown below.

Original measured transit lightcurve, normalized-intensity and time axes, complete grid and data points.

(D01.4) Using the given lightcurve answer the following questions. For your reference the above lightcurve is also given in the Summary Answersheet.

(D01.4a) [3 marks] Estimate the values of tTt_T and tFt_F in days by marking appropriate readings on the graph.

(D01.4b) [2 marks] Estimate the mean value of Δ\Delta, by marking appropriate readings on the graph and hence find RpR_p in units of R⊙R_\odot.

(D01.4c) [2 marks] Determine the value of ii in degrees assuming the orbital radius to be amin⁡a_{\min}.

Introducing limb darkening

So far we have assumed the stellar disk to be uniformly bright. In reality, the observed brightness of the stellar disk is not uniform due to “limb darkening” — an optical effect where the central part of the stellar disk appears brighter than the edge, or the “limb”.

The limb darkening effect can be measured by the relative intensity J(θ)≡I(θ)I(0)J(\theta)\equiv\frac{I(\theta)}{I(0)}, where θ\theta is the angle between the normal to the stellar surface at a point and the line joining the observer to that point, I(θ)I(\theta) is the observed intensity of the stellar disk at that point (I(0)I(0) being the intensity at the centre of the stellar disk). For a distant observer, θ\theta varies from θ=0\theta=0 (centre of the disk) to θ≈90∘\theta\approx90^\circ (edge of the disk).

(D01.5) The table below gives measured J(θ)J(\theta) at a certain wavelength for the Sun. We shall assume that the same limb darkening profile holds for the star S.

θ\thetaJ(θ)J(\theta)θ\thetaJ(θ)J(\theta)θ\thetaJ(θ)J(\theta)θ\thetaJ(θ)J(\theta)
0∘0^\circ1.00020∘20^\circ0.97140∘40^\circ0.88370∘70^\circ0.595
10∘10^\circ0.99425∘25^\circ0.95050∘50^\circ0.79480∘80^\circ0.475
15∘15^\circ0.98430∘30^\circ0.94360∘60^\circ0.72490∘90^\circ0.312

The limb darkening profile can be modelled by a quadratic formula:

J(θ)=1−a1(1−cos⁡θ)−a2(1−cos⁡θ)2,J(\theta)=1-a_1(1-\cos\theta)-a_2(1-\cos\theta)^2,

where a1a_1 and a2a_2 are two constants.

We shall estimate the unknown coefficients a1a_1 and a2a_2 from the given data by making a plot with suitable variables.

(D01.5a) [2 marks] Choose a pair of variables (x1,y1)(x_1,y_1) which are suitable functions of θ\theta and JJ, that you want to plot along xx and yy axes, respectively, to determine a1a_1 and a2a_2. Write the expressions for x1x_1 and y1y_1.

If you need to define additional variables for additional plots, define them as (x2,y2)(x_2,y_2), etc.

(D01.5b) [4 marks] Tabulate the values necessary for your plots.

(D01.5c) [7 marks] Plot the newly defined variables on the given graph paper (mark your graph as “D01.5c”).

(D01.5d) [7 marks] Obtain a1a_1 and a2a_2 from the plot. Uncertainties on the values are not needed.

Transit in the presence of limb darkening

Now, we consider planetary transits across a limb darkened stellar disk. In the presence of limb darkening, which we shall model by the quadratic formula of J(θ)J(\theta) given above, the average observed intensity of the entire stellar disk (without any transit), ⟨I⟩\langle I\rangle, is given by:

⟨I⟩=(1−a13−a26)I(0)\langle I\rangle=\left(1-\frac{a_1}3-\frac{a_2}6\right)I(0)

Further, the dip in the light caused by the transiting planet now depends not only on the relative size of the planet and the star, (RpRs)\left(\frac{R_p}{R_s}\right), but also on the intensity profile of the stellar disk along the transit chord, which in turn, depends on the angle of inclination, ii.

The schematic diagram below (not drawn to scale) shows the configuration. Note that the brighter part of the star is shown in a darker shade, while the planet is shown as a black dot.

Original limb-darkened stellar transit face-on and side-on views, transit chord, C, i, thetaC, planet and observer.

Here the relation between (RpRs)\left(\frac{R_p}{R_s}\right) and the measured Δ\Delta from the light curve is

Δ=I(θC)⟨I⟩(RpRs)2,\Delta=\frac{I(\theta_C)}{\langle I\rangle}\left(\frac{R_p}{R_s}\right)^2,

where I(θC)I(\theta_C) is the intensity of the stellar disk at the midpoint of the transit chord (point C in the figure above), θC\theta_C being the angle between the line of sight and the normal to the surface at that point. From the above it is obvious that for a given star, the same value of Δ\Delta can be produced by many combinations of the planet size, RpR_p, and the inclination angle ii.

(D01.6) It is possible to uniquely determine both RpR_p and ii by using data from transit lightcurves at two wavelengths, say, λB\lambda_B (blue) and λR\lambda_R (red). The limb darkening coefficients for these two wavelengths are given below:

Wavelengtha1a_1a2a_2
λB\lambda_B0.820.05
λR\lambda_R0.240.20

(D01.6a) [2 marks] Choose the correct statement among the following that describes the relation between the maximum depth of the transit (Δ\Delta) for λB\lambda_B and the inclination angle (ii) of the orbit and tick (✓) it in the Summary Answersheet.

(A) Δ\Delta increases with decreasing ii.

(B) Δ\Delta decreases with decreasing ii.

(C) Δ\Delta is independent of ii.

(D01.6b) [4 marks] The maximum depth of the transit (Δ\Delta) for the “SP system” was measured to be 0.0182 and 0.0159 for λB\lambda_B and λR\lambda_R, respectively.

Draw schematic transit light curves for both λB\lambda_B and λR\lambda_R on the given grid and label the curves by “B” and “R” respectively. Assume that the total transit duration is same for both wavelengths. The curves need not be to scale, but should represent the shapes of the light curves correctly.

(D01.7) We shall use a graphical method to find the values of RpR_p and ii for the SP system using measurements of Δ\Delta at λB\lambda_B and λR\lambda_R.

(D01.7a) [6 marks] Write an appropriate expression connecting the relevant variables that are to be plotted. (Hint: You may consider ii or bb, and RpR_p among the relevant variables.)

(D01.7b) [5 marks] Tabulate the appropriate quantities that are to be plotted.

(D01.7c) [7 marks] Draw a suitable graph and mark it as “D01.7c”.

(D01.7d) [4 marks] Estimate the values of RpR_p (in R⊙R_\odot) and ii (in degrees) from the graph.

(D01.8) [2 marks] Based on the results obtained in this problem, indicate whether the planet P is “ROCKY” or “GASEOUS” by ticking (✓) the appropriate box in the Summary Answersheet.

Решение

D02 Summary Answersheet — SA-1

Summary Answersheet — FRONT

Contestant Code: ______ – ______

Q. No. D02 · Page No. SA-1

Final answers

(D02.1a) Marks (3.0)

CME Nameuu (km s−1\mathrm{km\,s^{-1}})vv (km s−1\mathrm{km\,s^{-1}})τ\tau (h)
CME-A80447074.5
CME-B247360127.5
CME-C523396103.5
CME-D83041571.0
CME-E665400104.5
CME-F347350101.5
CME-G44637599.5
CME-H15536097.0
CME-I101651567.0
CME-J68341054.0

(D02.1b) Marks (15.0)

Graph number: D02.1b

m=m= ______

α=\alpha= ______

D02 Summary Answersheet — SA-2

Summary Answersheet — BACK

Contestant Code: ______ – ______

Q. No. D02 · Page No. SA-2

Final answers

(D02.1c) Marks (4.0)

CME Nameamodela_{\mathrm{model}} (m s−2\mathrm{m\,s^{-2}})
CME-A
CME-B
CME-C
CME-D
CME-E
CME-F
CME-G
CME-H
CME-I
CME-J

δarms=\delta a_{\mathrm{rms}}= ______

(D02.2a) Marks (4.0)

τ1,m=\tau_{1,\mathrm m}= ______ h

τ2,m=\tau_{2,\mathrm m}= ______ h

(D02.2b) Marks (2.0)

VALIDNOT VALID
CME-1
CME-2
D02 Summary Answersheet — SA-3

Summary Answersheet — FRONT

Contestant Code: ______ – ______

Q. No. D02 · Page No. SA-3

Final answers

(D02.3a) Marks (6.0)

C1:

δRD=\delta R_D= ______ R⊙R_\odot

C2:

δRD=\delta R_D= ______ R⊙R_\odot

(D02.3b) Marks (4.0)

Data Pointtt (h)C3: RD(t)R_D(t) (R⊙R_\odot)C4: RD(t)R_D(t) (R⊙R_\odot)
P521.0
P650.0
P785.0
P8111

(D02.3c) Marks (10.0)

Graph number: D02.3c

D02 Summary Answersheet — SA-4

Summary Answersheet — BACK

Contestant Code: ______ – ______

Q. No. D02 · Page No. SA-4

Final answers

(D02.3d) Marks (4.0)

C3: ∣δτ∣=|\delta\tau|= ______

C4: ∣δτ∣=|\delta\tau|= ______

(D02.3e) Marks (1.0)

TRUE □  FALSE □

(D02.4) Marks (7.0)

CME Nameuu (km s−1\mathrm{km\,s^{-1}})vv (km s−1\mathrm{km\,s^{-1}})τ\tau (h)VsV_s (km s−1\mathrm{km\,s^{-1}})
CME-A80447074.5
CME-B247360127.5
CME-C523396103.5
CME-D83041571.0
CME-E665400104.5
CME-F347350101.5
CME-G44637599.5
CME-H15536097.0
CME-I101651567.0
CME-J68341054.0

Vs,avg=V_{s,\mathrm{avg}}= ______

Total: ______

End of Summary Answersheet


Оригинал в Архива: Data-Analysis-Questions.pdf · официални решения: Data-Analysis-Solutions.pdf · допълнителни данни: Data-Analysis-Summary-Answersheet.pdf