IOAA 2018, Анализ на данни — D1. Dust and Young Stars in Star-forming Galaxies

12th International Olympiad on Astronomy and Astrophysics — Data Analysis Competition · 75 т.

Указания
  1. The data analysis competition lasts for 5 hours and is worth a total of 150 points.

  2. Dedicated IOAA Summary Answer Sheets are provided for writing your answers. Enter the final

answers into the appropriate boxes in the corresponding Summary Answer Sheet. On each Answer
Sheet, please fill in

  • Student’s Code
  1. Graph Paper is required for your solutions. On each Graph Paper sheet, please fill in

    • Student’s Code
    • Question no.
    • Graph no. and total number of graph paper sheets used.
  2. There are Answer Sheets for carrying out detailed work/rough work. On each Answer Sheet, please

fill in

  • Student’s Code
  • Question no.
  • Page no. and total number of pages.
  1. Start each problem on a separate Answer Sheet. 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 marked, cross it out.

  1. Use as many mathematical expressions as you think may help the graders to better understand your

solutions. The graders may not understand your language. If it is necessary to explain something in words, please use short phrases (if possible in English).

  1. You are not allowed to leave your working desk without permission. If you need any assistance

(malfunctioning calculator, need to visit a restroom, need more Answer Sheets, Graph Paper etc.), please put up your hand to signal the invigilator.

  1. The beginning and end of the competition will be indicated by a long sound signal. Additionally,

there will be a short sound signal fifteen minutes before the end of the competition (before the final long sound signal).

  1. At the end of the competition you must stop writing immediately. Sort and put your Summary

Answer Sheets, Graph Papers, and Answer Sheets in one stack. Put all other papers in another stack.
You are not allowed to take any sheet of paper out of the examination area.

  1. Wait at your table until your envelope is collected. Once all envelopes are collected, your student

guide will escort you out of the competition room.

  1. A list of constants is given on the next page.
Table of constants
GroupQuantityValue
EarthMass5.98×1024 kg5.98\times10^{24}\ \mathrm{kg}
EarthRadius6.38×106 m6.38\times10^6\ \mathrm{m}
EarthAcceleration of gravity9.81 m s−29.81\ \mathrm{m\,s^{-2}}
EarthObliquity of Ecliptic23°27′
EarthLength of Tropical Year365.2422 mean solar days
EarthLength of Sidereal Year365.2564 mean solar days
EarthAlbedo0.39
MoonMass7.35×1022 kg7.35\times10^{22}\ \mathrm{kg}
MoonRadius1.74×106 m1.74\times10^6\ \mathrm{m}
MoonMean Earth-Moon distance3.84×108 m3.84\times10^8\ \mathrm{m}
MoonOrbital inclination with the Ecliptic5.14°
MoonAlbedo0.14
MoonApparent magnitude (mean full moon)−12.74
SunMass1.99×1030 kg1.99\times10^{30}\ \mathrm{kg}
SunRadius6.96×108 m6.96\times10^8\ \mathrm{m}
SunLuminosity3.83×1026 W3.83\times10^{26}\ \mathrm{W}
SunAbsolute Magnitude4.80 mag
SunSurface Temperature5772 K
SunAngular diameter at Earth30′
SunOrbital velocity in Galaxy220 km s−1220\ \mathrm{km\,s^{-1}}
SunDistance from Galactic center8.5 kpc
Physical constants1 au1.50×1011 m1.50\times10^{11}\ \mathrm{m}
Physical constants1 pc206265 au
Physical constantsGravitational constant (G)6.67×10−11 N m2 kg−26.67\times10^{-11}\ \mathrm{N\,m^2\,kg^{-2}}
Physical constantsPlanck constant (h)6.62×10−34 J s6.62\times10^{-34}\ \mathrm{J\,s}
Physical constantsBoltzmann constant1.38×10−23 J K−11.38\times10^{-23}\ \mathrm{J\,K^{-1}}
Physical constantsStefan-Boltzmann constant5.67×10−8 W m−2 K−45.67\times10^{-8}\ \mathrm{W\,m^{-2}\,K^{-4}}
Physical constantsHubble constant67.8 km s−1 Mpc−167.8\ \mathrm{km\,s^{-1}\,Mpc^{-1}}
Physical constantsSpeed of light in vacuum299792458 m/s
Physical constantsMagnetic Permeability of free space4π×10−7 H m−14\pi\times10^{-7}\ \mathrm{H\,m^{-1}}
Physical constants1 Jansky (Jy)10−26 W m−2 Hz−110^{-26}\ \mathrm{W\,m^{-2}\,Hz^{-1}}
Δlog⁡10(x)=Δxxln⁡10.\Delta\log_{10}(x)=\frac{\Delta x}{x\ln10}.
Бележки за изданието

The independently archived question and solution PDFs are the source editions used here. D1.4 remains in the questions, answer form and worked key, but the key explicitly marks its 12 marks “discarded”; printed totals of 75 per problem and 150 overall are retained. Section 1.3 calls F1600F_{1600} “unattenuated” although section 1.2 defines it as observed; the source wording is preserved. The answer form labels the mass-function field D2.2.2 and omits a separate orbital-radius field. D2.3.3 calls the sin⁡i+Δsin⁡i\sin i+\Delta\sin i curve an upper limit and the sin⁡i−Δsin⁡i\sin i-\Delta\sin i curve a lower limit; the official figure places them in the opposite vertical order, as the mass function requires. In the key, equation (19) prints a strict inequality at the radius bound, equation (32) leaves the cgs-to-SI conversion implicit, and the D2.3.4 rubric assigns only four of its printed five marks. These source inconsistencies have not been silently corrected.

The solution PDF repeats an older wording of the questions: BB replaces QQ, f(β)f(\beta) replaces S(β)S(\beta), D1.1.3 explicitly requests least-squares fitting, Table 1 is called UV photometry, and the dex explanation is written dex(x)=10x\mathrm{dex}(x)=10^x. It additionally gives the 21 µJy ALMA limit and explicit redshift-factor hint reproduced with D1.4. Its APOGEE velocities are rounded to 34.85 and 42.57 km/s; its D2.3.2 options are (1) Blue Giant (2) Blue Dwarf (3) Red Giant (4) Red Dwarf (5) Main-Sequence Star (6) White Dwarf. The question edition’s data and five options are retained above; the official answer remains (3) Red Giant.

Условие

As a by-product of the star-forming process in a galaxy, interstellar dust can significantly absorb stellar light in ultraviolet (UV) and optical bands, and then re-emit in far-infrared (FIR), which corresponds to a wavelength range of 1010–300 μm300\,\mu\mathrm m.

1.1.

In the UV spectrum of a galaxy, the major contribution is from the light of the young stellar population generated in recent star-formation processes, thus the UV luminosity can act as a reliable tracer of the star-formation rate (SFR) of a galaxy. Since the observed UV luminosity is strongly affected by dust attenuation, extragalactic astronomers define an index called the UV continuum slope (β\beta) to quantify the shape of the UV continuum:

fλ=Q⋅λβf_\lambda=Q\cdot\lambda^\beta

where fλf_\lambda is the monochromatic flux of the galaxy at a given wavelength λ\lambda (in the unit of W m−3\mathrm{W\,m^{-3}}) and QQ is a scaling constant.

(D1.1.1) AB magnitude is a specific magnitude system. The AB magnitude is defined as:

mAB=−2.5log⁡fν3631 Jym_{\mathrm{AB}}=-2.5\log\frac{f_\nu}{3631\ \mathrm{Jy}}

The AB magnitude of a typical galaxy is roughly constant in the UV band. What is the UV continuum slope of this kind of galaxy? (Hint: fνΔν=fλΔλf_\nu\Delta\nu=f_\lambda\Delta\lambda) [6 т.]

(D1.1.2) Table 1 presents the observed IR photometry results for a z=6.60z=6.60 galaxy called CR7. Plot the AB magnitude of CR7 versus the logarithm of the rest-frame wavelength on graph paper and labelled as Figure 1. [12 т.]

(D1.1.3) Calculate CR7’s UV slope, plot the best-fit UV continuum on Figure 1 and make a comparison with the results you obtained in (D1.1.1). Is it dustier than the typical galaxy in (D1.1.1)? Please answer with [YES] or [NO]. (Hint: Express mABm_{\mathrm{AB}} as a function of λ\lambda and m1600m_{1600}, where m1600m_{1600} is the AB magnitude at λ0=160 nm\lambda_0=160\ \mathrm{nm} (1600 Å))

Table 1. (Observed Frame) IR Photometry of CR7 at z = 6.60

BandYJHK
Central Wavelength (µm)1.051.251.652.15
AB Magnitude24.71 ± 0.1124.63 ± 0.1325.08 ± 0.1425.15 ± 0.15

[5 т.]

1.2.

Under the assumption that dust grains in the galaxy absorb the energy of UV photons and re-emit it by blackbody radiation, the relation between the UV continuum slope (β\beta), UV brightness (at 1600 Å) and FIR brightness could be established:

IRX≡log⁡(FFIRF1600)=S(β)\mathrm{IRX}\equiv\log\left(\frac{F_{\mathrm{FIR}}}{F_{1600}}\right)=S(\beta)

where FFIRF_{\mathrm{FIR}} is the observed far-infrared flux and F1600F_{1600} is the observed flux at rest-frame wavelength 160 nm (1600 Å) (The “flux” FλF_\lambda is defined as Fλ=λ⋅fλF_\lambda=\lambda\cdot f_\lambda). Table 2 presents 20 measurements of β\beta, FFIRF_{\mathrm{FIR}} and F1600F_{1600} in nearby galaxies (Meurer et al. 1999).

Table 2. UV slope, flux and FIR flux of 20 nearby galaxies

Galaxy NameUV Slope β\betalog⁡(F1600/(10−3 W m−2))\log(F_{1600}/(10^{-3}\ \mathrm{W\,m^{-2}}))log⁡(FFIR/(10−3 W m−2))\log(F_{\mathrm{FIR}}/(10^{-3}\ \mathrm{W\,m^{-2}}))
NGC4861-2.46-9.89-9.97
Mrk 153-2.41-10.37-10.92
Tol 1924-416-2.12-10.05-10.17
UGC 9560-2.02-10.38-10.41
NGC 3991-1.91-10.14-9.80
Mrk 357-1.80-10.58-10.37
Mrk 36-1.72-10.68-10.94
NGC 4670-1.65-10.02-9.85
NGC 3125-1.49-10.19-9.64
UGC 3838-1.41-10.81-10.55
NGC 7250-1.33-10.23-9.77
NGC 7714-1.23-10.16-9.32
NGC 3049-1.14-10.69-9.84
NGC 3310-1.05-9.84-8.83
NGC 2782-0.90-10.50-9.33
NGC 1614-0.76-10.91-8.84
NGC 6052-0.72-10.62-9.48
NGC 3504-0.56-10.41-8.96
NGC 4194-0.26-10.62-8.99
NGC 32560.16-10.32-8.44

(D1.2.1) Based on the data given in Table 2, plot the IRX−β\mathrm{IRX}-\beta diagram on graph paper and labelled as Figure 2 and find a linear fit to the data. Write down your best-fit equation (i.e. IRX=a⋅β+b\mathrm{IRX}=a\cdot\beta+b) by the side of your diagram. [14 т.]

(D1.2.2) Quantify the dispersion (in ‘units’ of dex, where for example, log⁡(109)−log⁡(104)=5 dex\log(10^9)-\log(10^4)=5\ \mathrm{dex}) between the observed IRXobs\mathrm{IRX}_{\mathrm{obs}} and predicted IRXpred\mathrm{IRX}_{\mathrm{pred}} using the following equation:

σ=∑(ΔIRXi)2N−1(unit: dex),ΔIRXi=IRXi,obs−IRXi,pred.\sigma=\sqrt{\frac{\sum(\Delta\mathrm{IRX}_i)^2}{N-1}}\quad\text{(unit: dex)},\qquad\Delta\mathrm{IRX}_i=\mathrm{IRX}_{i,\mathrm{obs}}-\mathrm{IRX}_{i,\mathrm{pred}}.

[6 т.]

1.3.

Under the previous assumption of the energy transfer process, the ratio of FFIRF_{\mathrm{FIR}} to F1600F_{1600} can be expressed as:

FFIRF1600≈100.4A1600−1\frac{F_{\mathrm{FIR}}}{F_{1600}}\approx10^{0.4A_{1600}}-1

Where F1600F_{1600} is the unattenuated flux and AλA_\lambda is the dust absorption (in magnitudes) as a function of wavelength λ\lambda.

(D1.3.1) Express A1600A_{1600} as a function of IRX. [6 т.]

(D1.3.2) Based on Table 2 data and the function of A1600(IRX)A_{1600}(\mathrm{IRX}) you derived above, plot the A1600−βA_{1600}-\beta diagram on graph paper and label it as Figure 3 and find a linear fit to the data. Write down your best-fit equation (i.e. A1600=a′⋅β+b′A_{1600}=a^\prime\cdot\beta+b^\prime) by the side of your diagram. [12 т.]

(D1.3.3) If your linear model in (D1.3.2) is correct, what is the expected UV continuum slope β0\beta_0 of a dust-free galaxy? [2 т.]

1.4.

After establishing the local relation between UV continuum slope and IRX, we could probably test this empirical law in the high-redshift universe. In 2016, researchers obtained an Atacama Large Millimeter / submillimeter Array (ALMA) observation of CR7, and the FIR continuum corresponded to a 3σ3\sigma upper limit of an FIR flux of 1.5×10−19 W m−21.5\times10^{-19}\ \mathrm{W\,m^{-2}}.

(D1.4.1) Calculate the IRX of CR7. Is it an upper limit or lower limit?

Hint: here F1600F_{1600} should be written in the form of:

F1600=λ0⋅f1600F_{1600}=\lambda_0\cdot f_{1600}

where λ0=160 nm\lambda_0=160\ \mathrm{nm} (1600 Å) and f1600f_{1600} is the observed flux in the rest-frame [6 т.]

(D1.4.2) Is the current observation long enough to show any deviation of CR7 from the IRX−β\mathrm{IRX}-\beta relationship you just derived in the local universe? Please answer with [YES] or [NO] on the summary answer sheet, give the IRX difference and show the working used to calculate it on the answer sheet. [6 т.]

Summary Answer Sheet — (D1) Dust and Young Stars in Galaxies

Student’s Code: □□□–□□

(D1.1.1) The UV slope of such a typical galaxy is □;

(D1.1.2) See Figure 1;

(D1.1.3) The UV slope of CR7 is □; [YES] or [NO]? □;

(D1.2.1) See Figure 2; The best-fit Equation is: □;

(D1.2.2) The dispersion is □ dex;

(D1.3.1) A1600A_{1600} as a function of IRX: □;

(D1.3.2) See Figure 3; The Best-fit Equation is: □;

(D1.3.3) Expected UV slope of a dust-free galaxy: □;

(D1.4.1) The IRX of CR7 is □; Here the value is a □ limit (please fill with [upper] or [lower]);

(D1.4.2) Is the current observation deep enough ([YES] or [NO])? □.

Graph Paper

Student’s Code: □□□–□□

Blank linear graph paper sheet with Student’s Code boxes in the header.
Graph Paper
Blank logarithmic graph paper sheet; the left scale is labelled 10^3.00 to 10^3.50 in steps of 0.05, the right scale 1000 to 3200 in steps of 200.
Graph Paper(log)

Решение

Теми и трудност

Трудност: неоценена.

Светлина и астрономически измервания · Звездни величини · Междузвездно поглъщане

Предпоставки: Разширена олимпиадна подготовка


Оригинал в Архива: Data Analysis Problems.pdf · официални решения: Data Analysis Solutions.pdf · допълнителни данни: 18_da.pdf · допълнителни данни: 18_da_anssheet.pdf · допълнителни данни: 18_da_graph.pdf · допълнителни данни: 18_da_graph_log.pdf