Magnitude Systems

Magnitude Systems

Magnitude is a measure of the brightness of a star, astronomical object, or other celestial objects like artificial satellites. We would like to be able to communicate and express how bright an object is. For this, one could use the luminosity of the object or the flux received from that object. However, these numbers can get very large or very small. As a result, a magnitude scale is used.

Magnitudes are denoted by the unit ‘mag’, or with a superscript m^m, although they have no dimensions. The scale is reverse logarithmic: the brighter an object is, the lower its magnitude number. A difference of 1.0 mag in magnitude corresponds to a brightness ratio of 1005\sqrt[5]{100}, or about 2.512. For example, a magnitude 2.0 star is 2.512 times brighter than a magnitude 3.0 star, 6.31 times brighter than a magnitude 4.0 star, and 100 times brighter than a magnitude 7.0 star.

Apparent magnitude is a measure of the brightness of an object as seen from Earth. It is related to how much flux is received from that object by the observer. It is usually denoted by mm.

Absolute magnitude, on the other hand, is a measure of the intrinsic brightness of an object, defined as the apparent magnitude it would have if it were located at a distance of 10 parsecs from the observer, without any extinction. The absolute magnitude of a star is directly related to its luminosity. It is usually denoted by MM.

👀

The Greek astronomer Hipparchus produced a catalogue which noted the apparent brightness of stars in the second century BCE. In the second century CE, the Alexandrian astronomer Ptolemy classified stars on a six-point scale, and originated the term magnitude. To the unaided eye, a more prominent star such as Sirius or Arcturus appears larger than a less prominent star such as Mizar, which in turn appears larger than a truly faint star such as Alcor.

After several centuries of progress, in 1856 Norman Pogson, an English astronomer, proposed the modern logarithmic scale for magnitudes. This is the scale we use today.

Pogson’s Equation

The dimmer an object appears, the higher the numerical value given to its magnitude, with a difference of 5 magnitudes corresponding to a brightness factor of exactly 100. Let a flux value of F0F_0 correspond to a magnitude of m0m_0. Then, the magnitude mm corresponding to an object having a flux of FF is given by Pogson’s equation, which reads

mm0=5log100(FF0)m - m_0 = -5 \log_{100} \left( \frac{F}{F_0} \right)

or equivalently,

mm0=2.5logFF0(2.2.1)\tag{2.2.1} \boxed{m - m_0 = -2.5 \log \frac{F}{F_0}}

where the base of the logarithm is 10.

Consider an object with luminosity LL, distance rr, apparent magnitude mm, and absolute magnitude MM. The flux FF received from the object at Earth is F=L/4πr2F = L/4\pi r^2. The flux F0F_0 at 10pc10 \mathrm{pc} from the object is F0=L/4π(10pc)2F_0 = L/4\pi (10 \mathrm{pc})^2. Ignoring extinction, we get

mM=2.5log(FF0)=2.5log(L/4πr2L/4π(10pc)2)=2.5log((10pc)2r2)=5logr10pc \begin{aligned} m - M &= -2.5 \log \left( \frac{F}{F_0} \right) \\ &= -2.5 \log \left( \frac{L/4\pi r^2}{L/4\pi (10 \, \mathrm{pc})^2} \right) \\ &= -2.5 \log \left( \frac{(10 \, \mathrm{pc})^2}{r^2} \right) \\ &= 5 \log \frac{r}{10 \, \mathrm{pc}} \\ \end{aligned}

Hence we get that

mM=5logr10pc(2.2.2)\tag{2.2.2} \boxed{m - M = 5 \log \frac{r}{10 \, \mathrm{pc}}}

If rr is measured in parsecs,

mM=5logr5m - M = 5 \log r - 5

The distance modulus is defined as μ=mM\mu = m - M.

📍

The difference in absolute magnitudes of two objects gives the ratio of their luminosities. If M1M_1 and M2M_2 are the absolute magnitudes of two objects, then

M1M2=2.5logL1L2M_1 - M_2 = -2.5 \log \frac{L_1}{L_2}

where L1L_1 and L2L_2 are the luminosities of the two objects.

If the two objects have radii R1R_1 and R2R_2, and temperatures T1T_1 and T2T_2, then we can find luminosity using the Stefan-Boltzmann law

L=4πR2σT4L = 4 \pi R^2 \sigma T^4

From this we get

M1M2=5logR1R210logT1T2M_1 - M_2 = -5 \log \frac{R_1}{R_2} - 10 \log \frac{T_1}{T_2}
What is the absolute magnitude of the Sun? Apparent magnitude of the sun as seen from Earth is m=26.8mm_\odot = -26.8^m.

Using Pogson’s equation,

M=m+5log1AU10pc=4.8mM_\odot = m_\odot + 5 \log \frac{1 \, \mathrm{AU}}{10 \, \mathrm{pc}} = \boxed{4.8^m}
Albireo is a star in the constellation of Cygnus. It is a double and the system has a total magnitude of 3.0m3.0^m. The fainter companion (Albireo B) has a magnitude of 5.1m5.1^m. Find out the apparent magnitude of the brighter component of the double, Albireo A.

The light which we see coming from a double is the sum of the light from both components. Let F0F_0 be the flux which corresponds to a magnitude of 0. Let LAL_A and LBL_B be the luminosities of Albireo A and B, respectively. The total luminosity of the system is thus LA+LBL_A + L_B.

The luminosity of Albireo B can be found using

mB=2.5logLBF0    LB=F0100.4mBm_B = -2.5 \log \frac{L_B}{F_0} \implies L_B = F_0 \cdot 10^{-0.4 m_B}

The total luminosity of the system can be found using

mtot=2.5logLA+LBF0    LA+LB=F0100.4mtotm_\text{tot} = -2.5 \log \frac{L_A + L_B}{F_0} \implies L_A + L_B = F_0 \cdot 10^{-0.4 m_\text{tot}}

Hence the luminosity of Albireo A is

LA=F0100.4mtotF0100.4mBL_A = F_0 \cdot 10^{-0.4 m_\text{tot}} - F_0 \cdot 10^{-0.4 m_B}

The apparent magnitude of Albireo A is given by Pogson’s equation as

mA=2.5logLAF0=2.5log(100.4mtot100.4mB)m_A = -2.5 \log \frac{L_A}{F_0} = -2.5 \log \left(10^{-0.4 m_\text{tot}} - 10^{-0.4 m_B} \right)

Therefore we get that mA=3.2m\boxed{m_A = 3.2^m}

Color Indices

Usually, filters are used to allow only certain wavelengths to enter the detector. This gives more information about the object than measuring the flux through all wavelengths at once. One such system of filters used commonly is the UBV system. It consists of three filters:

  • The ultraviolet filter (U) is centered at 365 nm with an effective bandwidth of 66 nm
  • The blue filter (B) is centered at 445 nm with an effective bandwidth of 94 nm
  • The visual filter (V) is centered at 551 nm with an effective bandwidth of 88 nm

The apparent magnitudes measured through these filters are denoted by UU, BB and VV, or by mUm_U, mBm_B, and mVm_V. The absolute magnitudes are denoted by MUM_U, MBM_B, and MVM_V.

📍
Aside from the UBV system, there are other systems of filters used in astronomy, such as the Johnson-Cousins system, which includes additional filters like R (red), I (infrared) and J, H, K, L and M bands (which also lie in the infrared), and the SDSS system, which includes filters u, g, r, i, and z.

The color index is defined as the difference between the magnitudes measured through two different filters. For example, the UVU-V color index is the difference between its UU and VV magnitudes

(UV)=mUmV=MUMV=UV(U-V) = m_U - m_V = M_U - M_V = U - V

This gives the color of the star, which can be used to determine its temperature. The smaller the BVB-V index, the bluer the star is, hence the hotter it is. Color indices of a star are independent of its distance.

The sensitivity function S(λ)S(\lambda) is defined as the fraction of a star’s flux detected by a filter at a given wavelength. It depends on the optics of the system, filter bandwidth and response of the photometer. For example, the apparent UU magnitude is given by

U=2.5log0FλSU(λ)dλ+CU U = -2.5 \log \int_0^\infty F_\lambda S_U(\lambda)\, d\lambda + C_U

where CUC_U is a constant. The apparent BB and VV magnitudes are defined similarly.

The UBU-B color index is thus

(UB)=2.5log0FλSU(λ)dλ0FλSB(λ)dλ+CUB (U-B) = -2.5 \log \frac{\int_0^\infty F_\lambda S_U(\lambda)\, d\lambda}{\int_0^\infty F_\lambda S_B(\lambda)\, d\lambda} + C_{U-B}

where CUB=CUCBC_{U-B} = C_U - C_B.

Bolometric Magnitude

Apparent magnitude is defined in a specific bandpass, such as the visual band, and depends on the sensitivity of the detectors used. The apparent magnitude of an object in the visual band or ‘V’ band is called the apparent visual magnitude mVm_V of the object. The apparent bolometric magnitude mbolm_\text{bol} is defined as the apparent magnitude of an object if all its radiation were collected and measured over all wavelengths. This would correspond to an ideal detector which is able to measure at all wavelengths with complete efficiency.

The absolute bolometric magnitude MbolM_\text{bol} is defined as the absolute magnitude of an object if all its radiation were collected and measured over all wavelengths.

The bolometric correction BCBC is defined as the difference between the bolometric magnitude and the visual magnitude, given by

BC=mbolmV(2.2.3)\tag{2.2.3} \boxed{BC = m_\text{bol} - m_V}

We also have that BC=MbolMVBC = M_\text{bol} - M_V where MbolM_\text{bol} and MVM_V are the absolute bolometric and absolute visual magnitudes, respectively.

For the bolometric magnitude, the sensitivity function is unity everywhere: Sbol(λ)=1S_\text{bol} (\lambda) = 1.

mbol=2.5log0Fλdλ+Cbol\therefore m_\text{bol} = -2.5 \log \int_0^\infty F_\lambda\, d\lambda + C_\text{bol}

where CbolC_\text{bol} is a constant.

By definition, bolometric correction is zero for F5 spectral type stars.

Find an expression for the bolometric correction, assuming the visual sensitivity function is given by

SV(λ)={1if λ1<λ<λ20otherwiseS_V(\lambda) = \begin{cases} 1 & \text{if } \lambda_1 < \lambda < \lambda_2 \\ 0 & \text{otherwise} \end{cases}

Where λ1=55144\lambda_1 = 551 - 44 nm and λ2=551+44\lambda_2 = 551 + 44 nm.

Consider the zero point of both the visual and bolometric magnitudes to be the same, i.e. CV=CbolC_V = C_\text{bol}. Assume the emitting body emits like a perfect blackbody of temperature 5800K5800 \, \mathrm{K}.

The visual magnitude is given by

mV=2.5log0FλSV(λ)dλ+CV m_V = -2.5 \log \int_0^\infty F_\lambda S_V(\lambda)\, d\lambda + C_V

The bolometric magnitude is given by

mbol=2.5log0Fλdλ+Cbol m_\text{bol} = -2.5 \log \int_0^\infty F_\lambda\, d\lambda + C_\text{bol}

Thus the bolometric correction is given by

BC=mbolmV=2.5log(0Fλdλ0FλSV(λ)dλ) BC = m_\text{bol} - m_V = -2.5 \log \left( \frac{\int_0^\infty F_\lambda\, d\lambda}{\int_0^\infty F_\lambda S_V(\lambda)\, d\lambda} \right)

Since FλF_\lambda follows a blackbody spectrum, we can write

Fλ=πR2D22πhc2λ51ehcλkT1 F_\lambda = \pi \frac{R^2}{D^2} \cdot \frac{2 \pi h c^2}{\lambda^5} \frac{1}{e^{\frac{hc}{\lambda k T}} - 1}

Where TT is the temperature of the object, RR is the radius of the object and DD is the distance to the object. Defining x=hcλkTx = \frac{hc}{\lambda kT}, and doing some algebra, we get that the expression for the bolometric correction simplifies to

BC=2.5log(15π4x2x1x3ex1dx) BC = 2.5 \log \left( \frac{15}{\pi^4} \int_{x_2}^{x_1} \frac{x^3}{e^x - 1} \,dx \right)

where x1=hcλ1kT4.901x_1 = \frac{hc}{\lambda_1 k T} \approx 4.901 and x2=hcλ2kT4.176x_2 = \frac{hc}{\lambda_2 k T} \approx 4.176. Numerical integration gives

BC2.37m \boxed{BC \approx -2.37^m}

The bolometric correction of the sun, which too has a temperature of about 5800K5800 \, \mathrm{K}, is about 0.7m-0.7^m. This large difference is due to the assumptions we made: Sun does not emit like a perfect blackbody, the zero points of the visual and bolometric magnitudes are not the same, and the sensitivity function is not a perfect rectangle, but rather shaped like a bell curve.

Lightcurve

A lightcurve is a plot of brightness versus time for a celestial object, such as a star.

Lightcurve
A lightcurve is a plot of brightness versus time for a celestial object, such as a star.

Tracking the lightcurve of a star can provide information about its properties. For example, periodic dips in brightness can indicate the presence of an exoplanet orbiting the star, as the planet passes in front of the star and blocks some of its light. Variations in brightness can also indicate changes in the star’s activity, such as flares or pulsations.

Extended Objects

For extended objects (which are not point-like), the total magnitude can be obtained by summing up the luminosity over the area of the object. The apparent magnitude of an extended object is defined as the magnitude of a point source that would produce the same total flux as the extended object.

Thus the total amount of light from a galaxy of 12m12^m would be the same as that from a star of 12m12^m. However, since the light from the galaxy is spread over a larger area, it will appear dimmer than the star.

For a comet, the total magnitude is the combined magnitude of the coma and the nucleus.

The surface brightness SS quantifies the apparent brightness of flux density per unit angular area of the spatially extended obejct. For an object with total magnitude mm extended over a visual area Aarcsec2A \mathrm{\: arcsec^2}, its surface brightness is defined as

S=m+2.5logA(2.2.4)\tag{2.2.4} \boxed{S = m + 2.5 \log A}

Surface brightness is measured in units of mag/arcsec2\mathrm{mag \,/\, arcsec^{2}}. This is a much more useful quantity when comparing brightness of extended objects.

Extinction

ISM

If the space between the observer and the radiating source is not empty and is filled with interstellar medium (ISM), which consists of gas and dust. As light passes through the ISM, some of the it will be absorbed or scattered. This reduction of intensity of light in the direction of propagation is called extinction.

Suppose a star is radiating a flux L0L_0 into a solid angle ω\omega. The flux will decrease with distance because of extinction

dL=αLdr(2.2.5)\tag{2.2.5} dL = - \alpha L \, dr

The opacity α\alpha quantifies how effectively the medium can obsecure radiation. For vacuum, α=0\alpha = 0. The dimensionless quanitity τ\tau, called the optical thickness or optical depth, is defined as

dτ=αdr(2.2.6)\tag{2.2.6} d\tau = \alpha \, {dr} dL=Ldτ    L=L0eτ(2.2.7)\tag{2.2.7} \therefore dL = -L \, d\tau \implies \boxed{L = L_0 e^{-\tau}}

Here τ=0rαdr=αr\tau = \int_0^r \alpha \, dr = \langle \alpha \rangle r is the optical thickness between the source and the observer. This exponential decay is known as the Beer-Lambert law.

If F0F_0 is the flux density at the star’s surface and F(r)F(r) is the flux density at a distance rr,

F(r)=F0(Rr)2eτ F(r) = F_0 \left( \frac{R}{r} \right)^2 e^{-\tau}

From this we get

mM=2.5log(F(r)F(10pc))=5logr10pc+A=5logr10pc+ar(2.2.8) \tag{2.2.8} \begin{aligned} m - M &= -2.5 \log \left( \frac{F(r)}{F(10\,\mathrm{pc})} \right) \\ &= 5 \log \frac{r}{10 \, \mathrm{pc}} + A \\ &= 5 \log \frac{r}{10 \, \mathrm{pc}} + ar \\ \end{aligned}

where A=ar0A = ar \ge 0 is the extinction in magnitudes due to the interstellar medium, as light goes from the star to the observer. a=2.5αloge1.1αa = 2.5 \langle \alpha \rangle \log e \approx 1.1 \langle \alpha \rangle gives the extinction in magnitudes per unit distance.

Extinction also causes reddening of light; blue light is scattered and absorbed more than red light, hence the color index BVB-V increases.

V=MV+5logr10pc+AVB=MB+5logr10pc+AB \begin{aligned} V &= M_V + 5 \log \frac{r}{10 \mathrm{\, pc}} + A_V \\ B &= M_B + 5 \log \frac{r}{10 \mathrm{\, pc}} + A_B \\ \end{aligned} (BV)=(BV)0+EBV(2.2.9)\tag{2.2.9} \therefore \boxed{(B-V) = (B-V)_0 + E_{B-V}}

where (BV)0=MBMV(B-V)_0 = M_B - M_V is the intrinsic color index of the star and the color excess is defined as EBV=ABAVE_{B-V} = A_B - A_V.

Studies show that the ratio of the visual extinction AVA_V to the (BV)(B-V) color excess EBVE_{B-V} is almost constant for all stars

R=AVEBV3.1(2.2.10)\tag{2.2.10} R = \frac{A_V}{E_{B-V}} \approx 3.1

A star’s apparent magnitude in the V band is 15.1m15.1^m and absolute magnitude is equal to 1.3m1.3^m. Extinction of the interstellar medium per kpc is 1m1^m.

a) Find the distance to the star.
b) If the star’s observed color index is (BV)=1.6m(B-V) = 1.6^m, find its intrinsic color.
c) If BC=0.6mBC = -0.6^m for the star, find its luminosity.

a) The distance modulus is given by

mM=5logr10×103+arm - M = 5 \log \frac{r}{10 \times 10^{-3}} + ar

Where rr is measured in kpc, and a=1mkpc1a = 1^m \, \mathrm{kpc^{-1}}. Substituting the values, we get

3.8=logr+r3.8 = \log r + r

Solving it numerically, we get r3.28kpc\boxed{r \approx 3.28 \, \mathrm{kpc}}.

b) The intrinsic color index is given by

(BV)0=(BV)EBV(B-V)_0 = (B-V) - E_{B-V}

We can use the relation between AVA_V and EBVE_{B-V} to find the color excess. The extinction in visual band is

AV=ar=3.28mA_V = ar = 3.28^m

Hence the color excess is

EBV=AVR1.06mE_{B-V} = \frac{A_V}{R} \approx 1.06^m

Thus the intrinsic color index of the star is

(BV)0=1.6m1.06m=0.54m\boxed{(B-V)_0 = 1.6^m - 1.06^m = 0.54^m}

c) The bolometric magnitude of the star is

Mbol=MV+BC=1.3m0.6m=0.7mM_\text{bol} = M_V + BC = 1.3^m - 0.6^m = 0.7^m

Comparing with the bolometric magnitude of the Sun gives the luminosity of the star

MbolM,bol=2.5logLLM_\text{bol} - M_{\odot, \text{bol}} = -2.5 \log \frac{L}{L_\odot}

    L=L100.4(MbolM,bol)\implies L = L_\odot \cdot 10^{-0.4(M_\text{bol} - M_{\odot, \text{bol}})}

The absolute bolometric magnitude of the Sun is 4.75m4.75^m. Thus we get that the luminosity of the star is

L41.7L\boxed{L \approx 41.7 L_\odot}

Atmospheric Extinction

Earth’s atmosphere also causes extinction. The observed magnitude mm depends on the location of the observer and the zenith distance of the object. If the zenith distance zz is not too large, the atmosphere can be modelled by a plane layer of constant thickness. The airmass is defined as

X=secz(2.2.11)\tag{2.2.11} \boxed{X = \sec z}
Atmospheric extinction due to airmass

This approximation holds good for z70z \lesssim 70^\circ.

The increase in magnitude is

m=m0+kX(2.2.12)\tag{2.2.12} m = m_0 + kX

where kk is the extinction coefficient, which depends on the wavelength of light and the atmospheric conditions. The extinction coefficient is defined as the increase in magnitude per unit airmass. m0m_0 is the true apparent magnitude, devoid of extinction effects.

Atmospheric extinction is usually negligible for objects at zenith, but can be significant for objects near the horizon.

Since kk depends on the wavelength, the color index is also affected by atmospheric extinction.

Problems

Vega is known to have an apparent magnitude of about 00 (in reality it is about 0.030.03 and fluctuates). Given that it has an absolute magnitude of about 0.580.58, find its distance.

Using Pogson’s equation, we have

mM=5logr10pcm - M = 5 \log \frac{r}{10 \, \mathrm{pc}}

Substituting the values, we get

00.58=5logr100 - 0.58 = 5 \log \frac{r}{10}

    r=10100.1167.7pc\implies r = 10 \cdot 10^{-0.116} \approx 7.7 \, \mathrm{pc}

Hence the distance to Vega is approximately 7.7pc\boxed{7.7 \, \mathrm{pc}}.

Note that since Vega is so close to us, the extinction effects are negligible, and we can ignore them in this case.

What is the optical thickness of a layer of fog, if the Sun seen through the fog seems as bright as a full moon in a cloudless sky?

The apparent magnitudes of the Sun and the Moon are 26.8−26.8 and 12.5−12.5, respectively. Thus the total extinction in the cloud must be A=14.3A = 14.3. Since

A=(2.5loge)τ A = (2.5 \log e) \tau

where τ\tau is the optical thickness of the fog, we can rearrange to get

τ=A2.5loge=14.31.08613.2 \tau = \frac{A}{2.5 \log e} = \frac{14.3}{1.086} \approx \boxed{13.2}
Assume that all stars have the same absolute magnitude and stars are evenly distributed in space. Let N(m)N(m) be the number of stars brighter than m magnitudes. Find the ratio N(m+1)/N(m)N(m + 1)/N(m).

Let the number density of stars be nn, and their absolute magnitude be MM. Let the apparent magnitude of the star at a distance rr pc from the observer be mm. mm and rr are related by the equation

mM=5logr5 m - M = 5 \log r - 5

Rearranging gives

r=100.2(mM)+1 r = 10^{0.2(m - M) + 1}

All the stars brighter than magnitude mm will lies within rr pc distance to the observer. Thus the volume of the sphere upto radius rr is

V=43πr3=43π(100.2(mM)+1)3 V = \frac{4}{3} \pi r^3 = \frac{4}{3} \pi \left( 10^{0.2(m - M) + 1} \right)^3

The number of stars brighter than magnitude mm is given by

N(m)=nV=n43π(100.2(mM)+1)3 N(m) = n V = n \frac{4}{3} \pi \left( 10^{0.2(m - M) + 1} \right)^3

The number of stars brighter than magnitude m+1m + 1 is given by

N(m+1)=n43π(100.2(m+1M)+1)3 N(m + 1) = n \frac{4}{3} \pi \left( 10^{0.2(m + 1 - M) + 1} \right)^3

The ratio of the two is given by

N(m+1)N(m)=(100.2(m+1M)+1100.2(mM)+1)3=100.63.98 \frac{N(m + 1)}{N(m)} = \left( \frac{10^{0.2(m + 1 - M) + 1}}{10^{0.2(m - M) + 1}} \right)^3 = 10^{0.6} \approx \boxed{3.98}

An old planetary nebula, with a white dwarf (WD) in its center, is located 50 pc away from Earth. Exactly in the same direction, but behind the nebula, lies another WD, identical to the frist, but located at 150 pc from the Earth. Consider that the two WDs have absolute bolometric magnitude +14.2 and intrinsic color indexes BV=0.300B - V = 0.300 and UV=0.330U - V = 0.330. Extinction occurs in the interstellar medium and in the planetary nebula.

When we measure the color indices for the closer WD (the one who lies at the center of the nebula), we find the values BV=0.327B - V = 0.327 and UB=0.038U - B = 0.038. In this part of the Galaxy, the interstellar extinction rates are 1.50, 1.23 and 1.00 magnitudes per kiloparsec for the filters U, B and V, respectively. Calculate the color indices as they would be measured for the second star.

Let B0B_0, V0V_0, U0U_0 be the intrinsic absolute magnitudes of the white dwarf. Let aBa_B, aVa_V, aUa_U be the interstellar extinction coefficients per kpc. Let Aneb,BA_{neb,B}, Aneb,VA_{neb,V}, Aneb,UA_{neb,U} be the total extinction for the planetary nebula around the first WD, for light travelling from the centre of the nebula to the edge.

(UB)0=(UV)0(BV)0=0.3300.300=0.030(U-B)_0 = (U-V)_0 - (B-V)_0 = 0.330 - 0.300 = 0.030

For the closer WD,

(UV)1=(UB)1+(BV)1=0.038+0.327=0.365(U-V)_1 = (U-B)_1 + (B_V)_1 = 0.038 + 0.327 = 0.365

(BV)1=(BV)0+(aBaV)r+(Aneb,BAneb,V)(B-V)_1 = (B-V)_0 + (a_B - a_V) r + (A_{neb,B} - A_{neb,V})

(Aneb,BAneb,V)=(BV)1(BV)0(aBaV)r=0.0155\begin{align*} \therefore (A_{neb,B} - A_{neb,V}) &= (B-V)_1 - (B-V)_0 - (a_B - a_V) r \\ &= 0.0155 \end{align*}

Similarily,

(Aneb,UAneb,V)=(UV)1(UV)0(aUaV)r=0.0100\begin{align*} (A_{neb,U} - A_{neb,V}) &= (U-V)_1 - (U-V)_0 - (a_U - a_V) r \\ &= 0.0100 \end{align*}

(Aneb,UAneb,B)=(UB)1(UB)0(aUaB)r=0.0055\begin{align*} (A_{neb,U} - A_{neb,B}) &= (U-B)_1 - (U-B)_0 - (a_U - a_B) r \\ &= -0.0055 \end{align*}

For second white dwarf, the extinction in planetary nebula will be doubled as extinction the light from this pulsar has to cross the entire diameter of the nebula. Overall distance of this pulsar is 3 times larger.

(BV)2=(BV)0+3(aBaV)r+2(Aneb,BAneb,V)(B-V)_2 = (B-V)_0 + 3(a_B-a_V)r + 2(A_{neb,B} - A_{neb,V})

    (BV)20.366\implies \boxed{(B-V)_2 \approx 0.366}

(UV)2=(UV)0+3(aUaV)r+2(Aneb,UAneb,V)(U-V)_2 = (U-V)_0 + 3(a_U-a_V)r + 2(A_{neb,U} - A_{neb,V})

    (UV)20.425\implies \boxed{(U-V)_2 \approx 0.425}

(UB)2=(UB)0+3(aUaB)r+2(Aneb,UAneb,B)(U-B)_2 = (U-B)_0 + 3(a_U-a_B)r + 2(A_{neb,U} - A_{neb,B})

    (UB)20.060\implies \boxed{(U-B)_2 \approx 0.060}

A UBV photometric observation of a star gives U=8.15U = 8.15, B=8.50B = 8.50, and V=8.14V = 8.14. Based on the spectral class, one gets the intrinsic colour (UB)0=0.45(U − B)_0 = −0.45. If the star is known to have radius of 2.3R2.3R_\odot, absolute bolometric magnitude of 0.25−0.25, and bolometric correction (BC) of 0.15−0.15, determine:

a) the intrinsic magnitudes U, B, and V of the star
b) the effective temperature of the star
c) the distance to the star

Take the ratio of total to selective extinction RV=3.2R_V = 3.2 and the colour excess in (BV)(B − V) to be about 72% of the colour excess in (UB)(U − B).

a) Let intrinsic magnitudes of the star be U0U_0, B0B_0, and V0V_0 respectively.

EUB=(UB)(UB)0=(3.158.50)(0.45)=0.1E_{U-B} = (U-B) - (U - B)_0 = (3.15 - 8.50) - (-0.45) = 0.1

Since EUB=0.72EBVE_{U-B} = 0.72 \, E_{B-V}, EBV=0.10.72=0.14E_{B-V} = \frac{0.1}{0.72} = 0.14.

Using the relation AV=RVEBVA_V = R_V \, E_{B-V}, we get

AV=3.20.14=0.45A_V = 3.2 \cdot 0.14 = 0.45

V0=V+AV=7.69\therefore \boxed{V_0 = V + A_V = 7.69}

The intrinsic B0B_0 and U0U_0 color can be found using

(BV)0=(BV)EBV=(8.508.14)0.14=0.22(B-V)_0 = (B-V) - E_{B-V} = (8.50 - 8.14) - 0.14 = 0.22

    B0=V0+(BV)0=7.91U0=B0+(UB)0=7.46\begin{align*} \implies &\boxed{B_0 = V_0 + (B-V)_0 = 7.91} \\ &\boxed{U_0 = B_0 + (U-B)_0 = 7.46} \end{align*}

b) Comparing the star’s bolometric magnitude to that of the Sun,

MbolM,bol=2.5logLL    M,bolMbol=5logRR+10logTT\begin{align*} &M_\text{bol} - M_{\odot, \text{bol}} = -2.5 \log \frac{L}{L_\odot} \\ \implies &M_{\odot, \text{bol}} - M_\text{bol} = 5 \log \frac{R}{R_\odot} + 10 \log \frac{T}{T_\odot} \end{align*}

Substituting the values, we get

T=2.07T=12000K\boxed{T = 2.07 T_\odot = 12000 \, \mathrm{K}}

c) The absolute visual magnitude can be found using the bolometric correction

MV=Mbol+BC=0.25+(0.15)=0.40M_V = M_\text{bol} + BC = -0.25 + (-0.15) = -0.40

Hence the distance can be found

mVMV=5logr10pcm_V - M_V = 5 \log \frac{r}{10 \, \mathrm{pc}}

    7.69(0.40)=5logr10\implies 7.69 - (-0.40) = 5 \log \frac{r}{10}

    r415pc\implies \boxed{r \approx 415 \, \mathrm{pc}}
A pulsar, located 1000 pc far from Earth and 10,000 times more luminous than our Sun, emits radiation only from its two opposite poles, creating a homogeneous emission beam shaped as double cone with opening angle α=4\alpha = 4^\circ. Assuming the angle between the rotation axis and the emission axis is θ=30\theta = 30^\circ, and assuming a random orientation of the pulsar beams in relation to an observer on Earth, what is the probability of detecting the pulses? In case we can see it, what is the apparent bolometric magnitude of the pulsar?

As the pulsar rotates, its beam sweeps out a conical surface. The solid angle into which radiation is emitted (by both poles) in one rotation is

Ω=2[2π(1cos(θ+α2))2π(1cos(θα2))]=8πsinθsinα2\begin{align*} \Omega &= 2 \cdot \left[ 2 \pi (1 - \cos (\theta + \frac{\alpha}{2})) - 2 \pi (1 - \cos (\theta - \frac{\alpha}{2})) \right] \\ &= 8 \pi \sin \theta \sin \frac{\alpha}{2} \end{align*}

To detect the pulsar, Earth must lie inside this solid angle. Hence the probability is

p=Ω4π=2sinθsinα40.035p = \frac{\Omega}{4 \pi} = 2 \sin \theta \sin \frac{\alpha}{4} \approx \boxed{0.035}

The pulsar emits radiation into a cone having solid angle

Ω=24πsin2α4\Omega = 2 \cdot 4 \pi \sin^2 \frac{\alpha}{4}

Thus the flux recieved on Earth is

F=LA=Lωr2=L8πr21sin2α4F = \frac{L}{A} = \frac{L}{\omega r^2} = \frac{L}{8 \pi r^2} \frac{1}{\sin^2 \frac{\alpha}{4}}

Comparing the flux with the solar flux gives the apparent magnitude of the pulsar

m=M2.5logFF=M2.5logL8πr2sin2α44πr02L=M2.5logLL+5logr10pc+5log(2sinα4)=4.728.04=3.32\begin{align*} m &= M_\odot - 2.5 \log \frac{F}{F_\odot} \\ &= M_\odot - 2.5 \log \frac{L}{8 \pi r^2 \sin^2 \frac{\alpha}{4}} \frac{4 \pi r_0^2}{L_\odot} \\ &= M_\odot -2.5 \log \frac{L}{L_\odot} + 5 \log \frac{r}{10 \, \mathrm{pc}} + 5 \log (\sqrt{2} \sin \frac{\alpha}{4}) \\ &= 4.72 - 8.04 = \boxed{-3.32} \end{align*}

A star has an apparent magnitude mU=15.0m_U = 15.0 in the U-band. The U-band filter is ideal, i.e., it has perfect (100%) transmission within the band and is completely opaque (0% transmission) outside the band. The filter is centered at 360 nm, and has a width of 80 nm. It is assumed that the star also has a flat energy spectrum with respect to frequency. The conversion between magnitude, mm, in any band and flux density, ff, of a star in Jansky is given by

f=3631×100.4mJyf = 3631 \times 10^{-0.4m} \, \mathrm{Jy}

a) Approximately how many U-band photons, N0N_0, from this star will be incident normally on a 1m21 \, \mathrm{m^2} area at the top of the Earth’s atmosphere every second?

This star is being observed in the U-band using a ground based telescope, whose primary mirror has a diameter of 2.0m2.0 \, \mathrm{m}. Atmospheric extinction in U-band during the observation is 50%. You may assume that the seeing is diffraction limited. Average surface brightness of night sky in U-band was measured to be 22.0mag/arcsec222.0 \, \mathrm{mag/arcsec^2}.

b) What is the ratio, RR, of number of photons received per second from the star to that received from the sky, when measured over a circular aperture of diameter 22''?
c) In practice, only 20% of U-band photons falling on the primary mirror are detected. How many photons, NtN_t, from the star are detected per second?

a) The U-band is defined as (360±40360 \pm 40) nm. Thus, the minimum, maximum and average frequencies of the band are

νmin=cλmax=7.495×1014Hzνmax=cλmin=9.369×1014Hzνavg=νmin+νmax2=8.432×1014HzΔν=νmaxνmin=1.874×1014Hz\begin{align*} \nu_\text{min} &= \frac{c}{\lambda_\text{max}} = 7.495 \times 10^{14} \, \mathrm{Hz} \\ \nu_\text{max} &= \frac{c}{\lambda_\text{min}} = 9.369 \times 10^{14} \, \mathrm{Hz} \\ \nu_\text{avg} &= \frac{\nu_\text{min} + \nu_\text{max}}{2} = 8.432 \times 10^{14} \, \mathrm{Hz} \\ \Delta \nu &= \nu_\text{max} - \nu_\text{min} = 1.874 \times 10^{14} \, \mathrm{Hz} \end{align*}f=3631×100.415.0Jy=3.631×1029WHz1m2\therefore f = 3631 \times 10^{-0.4 \cdot 15.0} \, \mathrm{Jy} = 3.631 \times 10^{-29} \, \mathrm{W \, Hz^{-1} \, m^{-2}}

Now, N0×hνavg=ΔνfAΔtN_0 \times h \nu_\text{avg} = \Delta \nu \cdot f \cdot A \cdot \Delta t, where A=1m2A = 1 \, \mathrm{m^2} and Δt=1s\Delta t = 1 \, \mathrm{s}.

N0=fAΔνΔthνavg12180 \therefore N_0 = \frac{f A \Delta \nu \Delta t}{h \nu_\text{avg}} \approx \boxed{12180}

b) Let us call sky flux per square arcsec as Φ\Phi and total sky flux for the given aperture as ϕsky\phi_\text{sky}. Let total star flux be ϕst\phi_\text{st}.

ϕsky=ΦA=π(1)2Φ\phi_\text{sky} = \Phi A = \pi (1'')^2 \Phi

Using the formula for surface brightness, we get that the magnitude corresponding to this flux is

msky=22.0+2.5logΦϕsky=20.76mm_{sky} = 22.0 + 2.5 \log \frac{\Phi}{\phi_\text{sky}} = 20.76^m

As extinction is 50%,

R=0.5ϕstϕsky=0.5×100.4(1520.76)100R = \frac{0.5 \phi_\text{st}}{\phi_\text{sky}} = 0.5 \times 10^{-0.4(15 - 20.76)} \approx \boxed{100}

c) Equating the number of photons,

NT×1m=N0×0.5×0.2×AtN_T \times 1 \, \mathrm{m} = N_0 \times 0.5 \times 0.2 \times A_t

    NT=12180×0.5×0.2×(2.02)23813\implies N_T = 12180 \times 0.5 \times 0.2 \times \left( \frac{2.0}{2} \right)^2 \approx \boxed{3813}

The star β\beta-Doradus is a Cepheid variable star with a pulsation period of 9.84 days. We make a simplifying assumption that the star is brightest when it is most contracted (radius being R1R_1) and it is faintest when it is most expanded (radius being R2R_2). For simplicity, assume that the star maintains its spherical shape and behaves as a perfect black body at every instant during the entire cycle. The bolometric magnitude of the star varies from 3.463.46 to 4.084.08. From Doppler measurements, we know that during pulsation the stellar surface expands or contracts at an average radial speed of 12.8kms112.8 \, \mathrm{km \, s^{-1}}. Over the period of pulsation, the peak of thermal radiation (intrinsic) of the star varies from 531.0 nm to 649.1 nm.

a) Find the ratio of radii of the star in its most contracted and most expanded states (R1/R2R_1/R_2).
b) Find the radii of the star (in metres) in its most contracted and most expanded states (R1R_1 and R2R_2).
c) Calculate the flux of the star, F2F_2, when it is in its most expanded state.
d) Find the distance to the star, DstarD_\text{star}, in parsecs.

a) The ratio of the fluxes is

m1m2=2.5logF1F2    F1F2=100.4(m1m2)=1.77m_1 - m_2 = -2.5 \log \frac{F_1}{F_2} \implies \frac{F_1}{F_2} = 10^{-0.4(m_1 - m_2)} = 1.77

The flux is given by F=4πR2σT44πD2F = \frac{4 \pi R^2 \sigma T^4}{4 \pi D^2}, if the star is at a distance DD from the observer.

Hence we get that

F1F2=R12T14R22T24=1.77\frac{F_1}{F_2} = \frac{R_1^2 T_1^4}{R_2^2 T_2^4} = 1.77

Using Wien’s displacement law, T1T2=λ2λ1\frac{T_1}{T_2} = \frac{\lambda_2}{\lambda_1}.

    R1R2=F1F2(λ1λ2)2=0.890\implies \frac{R_1}{R_2} = \sqrt{\frac{F_1}{F_2}} \left( \frac{\lambda_1}{\lambda_2} \right)^2 = \boxed{0.890}

b) The radial velocity causes the change in radius of the star

R2R1=vP2=5.441×109mR_2 - R_1 = v \cdot \frac{P}{2} = 5.441 \times 10^9 \, \mathrm{m}

(10.890)R2=5.441×109m(1-0.890)R_2 = 5.441 \times 10^9 \, \mathrm{m}

R2=4.95×1010mR1=4.40×1010m\begin{align*} \therefore \: &\boxed{R_2 = 4.95 \times 10^10 \, \mathrm{m}} \\ &\boxed{R_1 = 4.40 \times 10^{10} \, \mathrm{m}} \end{align*}

c) To get the absolute value of flux (F2F_2) we must compare it with observed flux of the Sun.

m2m=2.5logF2Fm_2 - m_\odot = -2.5 \log \frac{F_2}{F_\odot}

    F2=F100.4(m2m)\implies F_2 = F_\odot \cdot 10^{-0.4(m_2 - m_\odot)}F2=1361Wm2100.4(4.08+26.72)=6.51×1010Wm2F_2 = 1361 \, \mathrm{W \, m^{-2}} \cdot 10^{-0.4(4.08 + 26.72)} = \boxed{6.51 \times 10^{-10} \, \mathrm{W \, m^{-2}}}

d) The temperature of the star when it is fully expanded is

T2=bλ2=4465KT_2 = \frac{b}{\lambda_2} = 4465 \, \mathrm{K}

The distance can be found using the observed flux

F2=4πR22σT244πDstar2F_2 = \frac{4 \pi R_2^2 \sigma T_2^4}{4 \pi D_\text{star}^2}

    Dstar=R2T22σF2=298pc\implies D_\text{star} = R_2 T_2^2 \sqrt{\frac{\sigma}{F_2}} = \boxed{298 \, \mathrm{pc}}
Given the intrinsic (BV)(B-V) color index of a star, find its color temperature TcT_c.

Assuming the star to be a blackbody, the spectral flux density in the Wien approximation of the star is propotional to

Fλλ5ehc/λkTc F_\lambda \propto \lambda^{-5} e^{-hc/\lambda kT_c}

The color index is defined as

(BV)=2.5log(FBFV) (B-V) = -2.5 \log \left( \frac{F_B}{F_V} \right)

    (BV)=2.5log(λB5ehc/λBkTcλV5ehc/λVkTc)\implies (B-V) = -2.5 \log \left( \frac{\lambda_B^{-5} e^{-hc/\lambda_B kT_c}}{\lambda_V^{-5} e^{-hc/\lambda_V kT_c}} \right)

    (BV)=12.5log(λB5λV5)+2.5hckTc(1λB1λV)\implies (B-V) = 12.5 \log \left( \frac{\lambda_B^5}{\lambda_V^5} \right) + 2.5 \frac{hc}{kT_c} \left( \frac{1}{\lambda_B} - \frac{1}{\lambda_V} \right)

    Tc=hck2.5(1λB1λV)(BV)+2.5log(λB5λV5)\implies T_c = \frac{hc}{k} \cdot \frac{2.5 \left( \frac{1}{\lambda_B} - \frac{1}{\lambda_V} \right)}{(B-V) + 2.5 \log \left( \frac{\lambda_B^5}{\lambda_V^5} \right)}

Substituting the values λB=440nm\lambda_B = 440 \, \mathrm{nm} and λV=550nm\lambda_V = 550 \, \mathrm{nm}, we get

Tc7000K(BV)+0.47\boxed{T_c \approx \frac{7000 \mathrm{K}}{(B-V) + 0.47}}