May I have some further discussion on
$$\sum_{n=0}^{\infty} {\frac{(-1)^{n+1}H_{2n}}{2n+1} \big(\psi(n+1)-\psi(n+\tfrac1{2})\big)} = \frac{\pi^3}{96} - \frac{\pi}{8}\ln^2\!2 \tag1$$
where $H_{n}$ is harmonic-number and $\psi(x)$ is digamma function.
Formerly, I deduce the result by this elementary approach, also the only solution I have due to the asset for me is quite insufficient.
Start with famous series
$$\frac1{2} \arctan x \ln (1+x^2) = \sum_{n=0}^{\infty} {\frac{(-1)^{n+1}H_{2n}}{2n+1} x^{2n+1}}$$
by which somehow I find (1) actually equal to
$$\sum_{n=0}^{\infty} {\frac{(-1)^{n+1}H_{2n}}{2n+1} \big(\psi(n+1)-\psi(n+\tfrac1{2})\big)} = \int_{0}^{1} {\frac{\arctan x \ln(1+x^2)}{x(1+x)} \mathrm{d}x}$$
and this integral is not impossible to solve.
Question 1: Any other possible access to summation as (1)?
For digamma function has strong relation with harmonic-number, here $\psi(n+1)-\psi(n+\tfrac1{2})$ can also be substitution for $H_{2n}-H_{n}$ with some other constants. So what if I want a part of (1) like
$$\sum_{n=0}^{\infty} {\frac{(-1)^{n+1}H_{2n}^{2}}{2n+1}} \quad \text{or} \quad \sum_{n=0}^{\infty} {\frac{(-1)^{n+1}H_{2n}H_{n}}{2n+1}} \tag2$$
but, of course, I have no idea how to deal with such summation involving higher order of harmonic-number.
Question 2: How to find the closed form for summation in (2)?