(adsbygoogle = window.adsbygoogle || []).push({}); Cosmology problem

Ok i've been trying to work on this assignment and it's driving me crazy agh, i'm pretty sure i should be able to do this with my hands and not need the use of maple or anything, although i do think the others used maple in the end and got a result of something like 7000K, but ok here's the question from the book:

"Imagine that at the time of recombination, the baryonic portion of the unvierse consisted entirely of [tex]^4 He[/tex] (that is, helium with two protons and two neutrons in its nucleus). The ionization energy of helium (that is the energy required to convert neutral [tex]He[/tex] to [tex]He^+[/tex] is [tex]Q_{He}=24,6eV[/tex]. At what temperature would the fractional ionization of helium be [tex]X=1/2[/tex]? Assume that [tex]\eta=5.5 \cdot 10^{-10}[/tex] and that the number density of [tex]He^{++}[/tex] is negligably small. [The relevatn statistical weight factor for the ionization of helium is [tex]g_{He}/(g_{e}g_{He^+}=1/4[/tex]]."

Ok so here's what i've been doing so far:

The book already did the calculations on this assuming a Hydrogen only universe, so i redid the derivation with He and simply came to the same equation times 1/4, so the equation i'm assuming i should be working with is:

[tex]\frac{1-X}{X^2}=\frac{1}{4}3.84\eta \big( \frac{kT}{m_{e} c^2} \big)^{3/2} exp{\frac{Q}{kT}}[/tex]

Where k is boltzmans constant. Then putting X=1/2, and putting it all to the power of 2/3 and throwing the constants on the left side i get:

[tex]\big( \frac{2}{0.96\eta} \big)^(2/3) = \frac{kT}{m_{e}c^2} exp{\frac{2Q}{3kT}}[/tex]

And this is where i run into problems, i can't seem to isolate T properly. I at least can't do it without the help of maple or such as far as i can see, at least not in its unchanged form. So i tried to do some taylor expansions to see if i could solve it, but i've had no luck, i get a result but it's usually quite ludicrous numbers, first i tried expanding [tex]exp{\frac{2Q}{3kT}}[/tex] just to the second expansion, but i'm guessing the number becomes too big since it's divided by k, so then i tried to take the natural logarithm on both sides, and then trying to do a taylor expansion on [tex]ln(kT)[/tex] around the point kT=1 since that should be a relatively small number. I got a better result from that, i'll show my calculations so you can see if i go wrong somewhere:

ln on both sides

[tex]ln( \big( \frac{2}{0.96\eta} \big)^{2/3} )= ln( \frac{kT}{m_{e}c^2} exp{\frac{2Q}{3kT}})[/tex]

Using a few logarithm rulse

[tex]\frac{2}{3} ln(\frac{2}{0.96\eta})= ln(kT)-ln(m_{e}c^2) +\frac{2Q}{3kT}[/tex]

Putting in ln(x)=x-1 for a taylor expansion

[tex]\frac{2}{3} ln(\frac{2}{0.96\eta})= (kT-1)-ln(m_{e}c^2) +\frac{2Q}{3kT}[/tex]

isolating the T's (and collecting the two ln's on the left side into one).

[tex]\frac{2}{3} ln(\frac{2m_{e}c^2}{0.96\eta})+1= kT+\frac{2Q}{3kT}[/tex]

[tex]\frac{2}{3} ln(\frac{2m_{e}c^2}{0.96\eta})+1= \frac{3kT^2+2Q}{3kT}[/tex]

Define my constants:

[tex]\frac{2}{3} ln(\frac{2m_{e}c^2}{0.96\eta})+1=24,46=K_1[/tex]

[tex]2Q=49,6=K_2[/tex]

So i have:

[tex]K_1=\frac{3x^2+K_2}{x} -> 0=3x^2 - x \cdot K_1 + K_2 [/tex]

And solving for x:

[tex]x=\frac{K_1 \pm \sqrt{K_1^2 - 4K_2}}{2 \cdot 3} = \frac{24,5 \pm \sqrt{598,6 - 198,4}}{6} =\frac{24,5 \pm 20}{6} [/tex]

From this i get two solutions [tex]x_{1}=0.75[/tex] and [tex]x_{2}=7,4[/tex]. Now the first one gives a somewhat good answer, around 8700K, which is decent but it still seems to be quite a big % off from what people got solving it with the computer, the other one is ridiculously high, but i'm not sure i would argue that i use the other one, other then intuitively i wouldn't expect it to be correct. Also in this all i did take the logarithm of a number with units (i used [tex]m_{e}c^2 = 511000ev[/tex]), so i'm not sure how i could get around taht either. Any suggestions?

**Physics Forums | Science Articles, Homework Help, Discussion**

Join Physics Forums Today!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Kosmologi problem

**Physics Forums | Science Articles, Homework Help, Discussion**