In my partial answer to this question: Integral involving polylogarithm and an exponential, I arrive at the integral
$$ \int_0^{a^2} \mathrm{Ei} (-s) \frac{1 - e^s}{s} ds , ~~~~ (\ast) $$
where $a \in \mathbb{R}$ and $\mathrm{Ei}$ is the standard exponential integral defined by $$ \mathrm{Ei} (z) = \int_{-\infty}^z \frac{e^t}{t} d t , $$ where the Cauchy principal value is taken for $z > 0$. The integral ($\ast$) is very similar to the integrals #28, #30, and #31 on page 198 of https://ia600303.us.archive.org/1/items/jresv73Bn3p191/jresv73Bn3p191_A1b.pdf. I haven't been able to complete the connection. Does anyone know how to express ($\ast$) in terms of known special functions?