Is Re-ionisation of Hydrogen due to Lyman-alpha emitters (LAEs)?
GitHub Repository: https://github.com/JP-Carr/ReHILAE
Group members and roles:
James Gold – Coordinator/PI
Connor Donovan – Write-up lead/blog
Jack Bowden – Theory lead
James Carr – Code lead
Joe Philips – Administrator
Our group, ReHILAE, is investigating whether or not LAEs are the dominating LyC (ionising radiation) sources during re-ionisation and therefore the main driving force behind the epoch of re-ionisation.
We will be building on the work of Faisst, A. (2016), initially we will be running and improving on existing python scripts which represent the equations on page 7 of this paper. Following on from this we focus on improving the availability of ionising photons, which is currently given by:
This formula can be improved by considering the work of Sobral, D. and Matthee, J. (2019). Using equation (9) in this paper (page 8) we can replace the in the equation above with:
Week 1 Research Update:
Today, we created our own Python scripts describing all the equations needed for our model of re-ionisation. These equations are shown below:
Where and . The final factor we need to consider before calculating the fraction of ionised hydrogen is the recombination time for ionised hydrogen.
Where and are the primordial mass fractions of helium (25%) and hydrogen (75%) respectively. The term is the recombination coefficient and the term is the hydrogen number density, both of which are described below.
Where is the temperature of the IGM. We then plotted each of these equations to illustrate the relationships between these variables.
We also plotted the UV density luminosity function as a function of redshift:
The main and final equation used in our model is a differential equation that describes the evolution of the fraction of ionised hydrogen as function of redshift or cosmic time, which is shown below.
The code below shows how this equation was implemented and solved in Python. The odeint function returns an array containing the fraction of ionised hydrogen for each inputted value of time (which is immediately converted to redshift).
def dQ_dt(Q,t): z= ((((28./(t))-1.)**(1./2.)-1.)) trec = 1./(C*alpha_B*(1.+Y_p/(4.*X_p))*n_H*(1.+z)**(3.)) nion = f_esc0*xi_ion*p(z)/ (2.938e+73) return ((nion)/n_H-Q/trec)*3.1536e+16 ts = np.linspace(0.051,14,10000000) # time in Gyr zs= ((((28./(ts))-1.)**(1./2.)-1.)) # conversion from Gyr to redshift Q = odeint(dQ_dt, 0., ts)
This array can be plotted, to visualise the epoch of re-ionisation:
In the coming weeks we will also explore how varying our initial assumptions effect the plot above.