- #1
dthomas
- 4
- 0
I am trying to calculate the ionization curve for hydrogen at a given matter density (in this case it happens to be 55 atoms/cm^3)
From the Saha equation and putting the concentrations in terms of reaction coordinates. It looks something like this:
Log[x[t]^2/(1 - x[t])] = 1/(k t)E - 3/2 Log[-(1/(k t))E] + 3/2 Log[(-2pimE)/(n^(2/3) h^2)]
k is the Boltzmann constant, h is Planck's constant, and E is a given energy value
t is Temperature
Basically I need to get Mathematica to plot the reaction coordinate (x) vs. t
If anyone has some advice as to how to do this, it would be greatly appreciated
From the Saha equation and putting the concentrations in terms of reaction coordinates. It looks something like this:
Log[x[t]^2/(1 - x[t])] = 1/(k t)E - 3/2 Log[-(1/(k t))E] + 3/2 Log[(-2pimE)/(n^(2/3) h^2)]
k is the Boltzmann constant, h is Planck's constant, and E is a given energy value
t is Temperature
Basically I need to get Mathematica to plot the reaction coordinate (x) vs. t
If anyone has some advice as to how to do this, it would be greatly appreciated