IOAA 2023, Анализ на данни — Задача 2. Isolated black hole

Data Analysis · 75 т.

Table of Constants
Fundamental constants
Speed of light in vacuumcc==2.998×108 m s−12.998 \times 10^{8}\ \mathrm{m\,s^{-1}}
Planck constanthh==6.626×10−34 J s6.626 \times 10^{-34}\ \mathrm{J\,s}
Boltzmann constantkBk_B==1.381×10−23 J K−11.381 \times 10^{-23}\ \mathrm{J\,K^{-1}}
Stefan-Boltzmann constantσ\sigma==5.670×10−8 W m−2 K−45.670 \times 10^{-8}\ \mathrm{W\,m^{-2}\,K^{-4}}
Elementary chargeee==1.602×10−19 C1.602 \times 10^{-19}\ \mathrm{C}
Universal gravitational constantGG==6.674×10−11 N m2 kg−26.674 \times 10^{-11}\ \mathrm{N\,m^{2}\,kg^{-2}}
Universal electric constantϵ0\epsilon_0==8.854×10−12 m−3kg−1s4A28.854 \times 10^{-12}\ \mathrm{m^{-3}kg^{-1}s^{4}A^{2}}
Universal gas constantRR==8.315 J mol−1 K−18.315\ \mathrm{J\,mol^{-1}\,K^{-1}}
Avogadro constantNAN_\mathrm{A}==6.022×1023 mol−16.022 \times 10^{23}\ \mathrm{mol^{-1}}
Wien’s displacement constantb=λmTb = \lambda_m T==2.898×10−3 m K2.898 \times 10^{-3}\ \mathrm{m\,K}
Mass of electronmem_\mathrm{e}==9.109×10−31 kg9.109 \times 10^{-31}\ \mathrm{kg}
Mass of protonmpm_\mathrm{p}==1.673×10−27 kg1.673 \times 10^{-27}\ \mathrm{kg}
Mass of neutronmnm_\mathrm{n}==1.675×10−27 kg1.675 \times 10^{-27}\ \mathrm{kg}
Mass of Helium nucleusmHem_\mathrm{He}==6.645×10−27 kg6.645 \times 10^{-27}\ \mathrm{kg}
Atomic mass unit (a.m.u., Dalton)==1.661×10−27 kg1.661 \times 10^{-27}\ \mathrm{kg}
Astronomical data
Hubble constantH0H_0==70 km s−1 Mpc−170\ \mathrm{km\,s^{-1}\,Mpc^{-1}}
North Ecliptic Pole (J2000.0J2000.0)(αE, δE)(\alpha_\mathrm{E},\ \delta_\mathrm{E})(18h00m00s,+66∘33′39′′)(18^\mathrm{h}00^\mathrm{m}00^\mathrm{s}, +66^{\circ}33'39'')
North Galactic Pole (J2000.0J2000.0)(αG, δG)(\alpha_\mathrm{G},\ \delta_\mathrm{G})(12h51m26s,+27∘07′42′′)(12^\mathrm{h}51^\mathrm{m}26^\mathrm{s}, +27^{\circ}07'42'')
1 jansky1 Jy==10−26 W m−2 Hz−110^{-26}\ \mathrm{W\,m^{-2}\,Hz^{-1}}
1 parsec1 pc==3.086×1016 m3.086 \times 10^{16}\ \mathrm{m}
206 265206\,265 au
3.2623.262 ly
1 astronomical unit (au)1 au==1.496×1011 m1.496 \times 10^{11}\ \mathrm{m}
1 sidereal dayTSDT_\mathrm{SD}==23.9344423.93444 h
23h56m04s23^\mathrm{h}56^\mathrm{m}04^\mathrm{s}
1 tropical year==365.2422365.2422 solar days
1 sidereal year==365.2564365.2564 solar days

Approximations

(1+x)n≈1+nx(1+x)^n \approx 1 + nx

(1+x)(1+y)≈1+x+y(1+x)(1+y) \approx 1 + x + y if x≪1x \ll 1 and y≪1y \ll 1

Gauss’s formulae

Spherical law of cosines: cos⁡a=cos⁡bcos⁡c+sin⁡bsin⁡ccos⁡A\cos a = \cos b \cos c + \sin b \sin c \cos A

Spherical law of sines: sin⁡Asin⁡a=sin⁡Bsin⁡b=sin⁡Csin⁡c\dfrac{\sin A}{\sin a} = \dfrac{\sin B}{\sin b} = \dfrac{\sin C}{\sin c}

The Sun
Solar luminosityL⊙L_\odot==3.826×1026 W3.826 \times 10^{26}\ \mathrm{W}
Apparent angular diameter of Sunθ⊙\theta_\odot==32’
Effective temperature of SunTeff,⊙T_{\mathrm{eff},\odot}==5778 K5778\ \mathrm{K}
Apparent visual magnitude==−26.75-26.75
Absolute visual magnitude==+4.82+4.82
Apparent bolometric magnitude==−26.83-26.83
Absolute bolometric magnitude==+4.74+4.74
Distance of the Sun from the Galactic centre≈\approx8 kpc8\ kpc
The Earth and Moon
Obliquity of the ecliptic (Earth)ϵ\epsilon==23.5∘23.5^{\circ}
Platonic year (period of precession of Earth’s axis)==25 76525\,765 years
Apparent visual magnitude of full Moon==−12.74-12.74
Apparent angular diameter of MoonθL\theta_\mathrm{L}==31′31'
Inclination of the lunar orbit to the ecliptic==05∘08′43′′05^{\circ}08'43''
Inclination of the lunar equator to its orbital plane==6.687∘6.687^{\circ}
Lunar sidereal monthTSLT_\mathrm{SL}==27.32166127.321661 d
655.71986655.71986 h
Synodic month==29.53058929.530589 d
Tropical month==27.32158227.321582 d
Anomalistic month (the average time between two successive Moon’s perigees)==27.55455027.554550 d
Draconic month (time between two successive transits of the Moon through the same node)==27.21222127.212221 d

The Solar System

ObjectMean radius [km]Mass [kg]Semimajor axis [au]Eccentricity
Sun695 700695\,7001.988×10301.988 \times 10^{30}——
Mercury2 4402\,4403.301×10233.301 \times 10^{23}0.3870.3870.2060.206
Venus6 0526\,0524.867×10244.867 \times 10^{24}0.7230.7230.0070.007
Earth6 3786\,3785.972×10245.972 \times 10^{24}1.0001.0000.016 7100.016\,710
Moon1 7371\,7377.346×10227.346 \times 10^{22}3.844×1053.844 \times 10^{5} km0.054 9000.054\,900 (range 0.026−0.0770.026 - 0.077)
Mars3 3903\,3906.417×10236.417 \times 10^{23}1.5241.5240.0930.093
Jupiter69 91169\,9111.898×10271.898 \times 10^{27}5.2035.2030.0480.048
Saturn58 23258\,2325.683×10265.683 \times 10^{26}9.5379.5370.0540.054
Uranus25 36225\,3628.681×10258.681 \times 10^{25}19.18919.1890.0470.047
Neptune24 62224\,6221.024×10261.024 \times 10^{26}30.07030.0700.0090.009
Data analysis - instructions
  • Do not touch envelopes until the start of the examination.
  • The data analysis examination lasts for 3 hours and is worth a total of 125 marks.
  • There are Answer Sheets for carrying out detailed work and Working Sheets for rough work, which are already marked with your student code and question number.
  • Use only the answer sheets for a particular question for your answer. Please write only on the printed side of the sheet. Do not use the reverse side. If you have written something on any sheet which you do not want to be evaluated, cross it out.
  • Use as many mathematical expressions as you think may help the evaluator to better understand your solutions. The evaluator may not understand your language. If it is necessary to explain something in words, please use short phrases (if possible in English).
  • You are not allowed to leave your work desk without permission. If you need any assistance (malfunctioning calculator, need to visit a restroom, etc.), please draw the attention of the supervisor.
  • The beginning and end of the examination will be indicated by the supervisor. The remaining time will be displayed on a clock.
  • At the end of the examination you must stop writing immediately. Put everything back in the envelope and leave it on the table.
  • Once all envelopes are collected, your student guide will escort you out of the examination room.
  • A list of constants and useful relations are included in the envelope.

Условие

In 2022, two independent groups reported the discovery of an isolated black hole based on observations of the gravitational microlensing event OGLE-2011-BLG-0462. In this problem, we will analyze data from the Hubble Space Telescope to reproduce their findings.

Gravitational microlensing occurs when the light of a distant star (the ‘source’) is bent and magnified by the gravitational field of an intervening object (the ‘lens’). The characteristic angular scale of gravitational microlensing events, called the angular Einstein radius θE\theta_E, depends on the mass MM and distance DℓD_\ell from the Earth to the lens:

θE=4GMc2Ds−DℓDsDℓ,\theta_\mathrm{E} = \sqrt{\frac{4GM}{c^2}\frac{D_s - D_\ell}{D_s D_\ell}},

where DsD_s is the distance to the source star. For typical microlensing events observed in the Milky Way, the source stars are in the Galactic bulge, near the Galactic center, so Ds≈8D_s \approx 8 kpc.

(a) Calculate the angular Einstein radius in milliarcseconds (mas) for an example lens of 1 M⊙1\ M_\odot located at a distance of 1 kpc. [2 т.]

Suppose that at time tt the lens and the source are separated by an angle θ≡u(t)θE\theta \equiv u(t)\theta_E on the sky. Two images of the source are created on a line through the positions of the source and the lens, at angular distances θ+\theta_+ and θ−\theta_- from the lens given by:

θ±=12(u±u2+4)θE.\theta_\pm = \frac{1}{2}\left(u \pm \sqrt{u^2+4}\right)\theta_\mathrm{E}.

These two images are magnified, relative to the unlensed brightness of the source. The absolute amplification of the images is:

A±=12(u2+2uu2+4±1).A_\pm = \frac{1}{2}\left(\frac{u^2+2}{u\sqrt{u^2+4}} \pm 1\right).

The image below shows the geometry of the event. The position of the lens is marked as LL, the unlensed position of the source is marked as SS, while A+A_+ and A−A_- mark the positions of the two images of the source. The dashed circle has a radius of one Einstein radius.

Geometry of the microlensing event: a dashed circle of one Einstein radius centred on the lens L; a straight line through L passes the unlensed source position S (small grey circle) inside the circle and reaches the image A+ (elongated ellipse) outside the circle on the left and the image A− (small ellipse) on the circle on the right; the segment A+–L is marked θ+ (blue), S–L is marked θ (orange), L–A− is marked θ− (green).
(A color version of this figure is available at a separate sheet.)

(b) Current telescopes cannot normally resolve this pair of images, but only measure the position of the image centroid, i.e. the brightness-weighted mean of the positions of the two images. Derive an expression for the angular separation θc\theta_c of the image centroid relative to the lens as a function of uu and θE\theta_E. [8 т.]

(c) Derive an expression for the source deflection Δθ\Delta\theta, i.e. the difference between the location of the centroid and the unlensed position of the source, as a function of uu and θE\theta_\mathrm{E}. What is the source deflection when the lens and the source are nearly perfectly aligned (u≈0u \approx 0)? [4 т.]

The source and lens are moving relative to each other in the sky. Thus, both the total magnification of the images and the position of the centroid changes with time, resulting in observable photometric and astrometric microlensing effects. For now, we assume that the source-lens relative motion is rectilinear.

The plot below shows the light curve of the gravitational microlensing event OGLE-2011-BLG-0462, discovered by the OGLE sky survey led by astronomers from the University of Warsaw. The solid line shows the best-fitting light curve model. The Einstein timescale of the event, i.e. the time needed for the source to move by one angular Einstein radius relative to the lens, was tE=247t_\mathrm{E} = 247 days. The event peaked on 21 July 2011 (HJD=2455763\mathrm{HJD} = 2455763). The minimal separation between the lens and the source was u0≈0u_0 \approx 0.

Light curve of OGLE-2011-BLG-0462: upper panel shows I magnitude (13 to 16, brighter upward) against HJD − 2450000 (days) from 5300 to 6300, with data points near 16.5 mag and a sharp peak reaching about 13.3 mag near 5763, and a blue best-fitting model line; lower panel shows residuals from −0.05 to 0.05 scattered around 0.00.

The table below shows the measured positions of the source star against the background objects in the East and North directions based on images from the Hubble Space Telescope.

HJDE position (mas)N position (mas)
2455765.22455765.22.58±0.132.58 \pm 0.137.29±0.167.29 \pm 0.16
2455865.72455865.72.32±0.122.32 \pm 0.125.44±0.245.44 \pm 0.24
2456179.72456179.70.46±0.140.46 \pm 0.141.62±0.081.62 \pm 0.08
2456195.82456195.80.88±0.360.88 \pm 0.361.56±0.771.56 \pm 0.77
2456426.22456426.2−1.02±0.21-1.02 \pm 0.21−0.94±0.12-0.94 \pm 0.12
2456587.72456587.7−2.04±0.07-2.04 \pm 0.07−1.88±0.40-1.88 \pm 0.40
2456956.62456956.6−4.54±0.25-4.54 \pm 0.25−5.16±0.29-5.16 \pm 0.29
2457995.22457995.2−11.14±0.12-11.14 \pm 0.12−15.14±0.17-15.14 \pm 0.17

(d) Plot the measured positions of the source star in the East and North directions as a function of time. [10 т.]

(e) The observed motion of the source star is the sum of two effects: rectilinear proper motion of the source and astrometric microlensing effects. Calculate the proper motion (in mas/year) of the source in the East and North directions and its uncertainty. [8 т.]

(f) After subtracting the effects of proper motion from the data, calculate and plot the total resultant astrometric deflection as a function of uu. Neglect uncertainty of the proper motion determination. [20 т.]

(g) Analyse the data to determine the angular Einstein radius θE\theta_\mathrm{E} of the event and its uncertainty. (Hint: it may be helpful to linearise the expression for Δθ\Delta\theta). [16 т.]

(h) For long-timescale events such as OGLE-2011-BLG-0462, the rectilinear approximation of the relative lens-source proper motion is not strictly true and the orbital motion of the Earth has to be taken into account. This allows measurement of a dimensionless quantity called the microlensing parallax, defined as πE=(πl−πs)/θE\pi_\mathrm{E} = (\pi_l - \pi_s)/\theta_\mathrm{E}, where πl\pi_l and πs\pi_s are parallaxes of the lens and the source, respectively.

For this event πE=0.095±0.009\pi_E = 0.095 \pm 0.009. Rearrange the expression for θE\theta_\mathrm{E} given earlier to calculate the mass of the lens in solar masses and its uncertainty. [7 т.]

Отговори

Решение


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