search

ladsantos / p-winds

Python implementation of Parker wind models for exoplanet atmospheres
https://p-winds.readthedocs.io/
MIT License
10 stars 5 forks source link

Exact calculation of radiative processes. #27

Closed ladsantos closed 3 years ago

ladsantos commented 3 years ago

As pointed out by @astroshrey:

The current version of the code uses the approximation in Equation (9) of Oklopcic & Hirata (2018), where the photoionization rate Phi_prime is approximated as the photoionization rate at null optical depth Phi multiplied by exp(-tau_0). While this approximation works well in the optically thin outer part of the outflow, I have found that it can overestimate the photoionization rate by a factor of several near the ionization front. Therefore I have written an exact version of the photoionization calculation which can be accessed in the top-level "ion_fraction" function by setting the argument "exact_phi = True". I have left it as False by default so it does not change code behavior, but the default could be updated in the future.

ladsantos commented 3 years ago

The PR #25 may have some conflicts with the newest commits of hydrogen.py. I will work on correcting them.