Yes, there is!
The derivation, which involves a lot of computer-assisted manipulations, is a little bit long, so I won’t put it here just now, but the answer is
where C is the Catalan constant, γ is Euler’s gamma, and Li3 is a polylogarithm, and (edited) the last term ψ(5/4)=4−γ−π2−log8 is there because everything below computes the sum over n≥0, so I need to subtract n=0 from the answer below to get the sum over n≥1.
Edit to explain the calculation.
The way to derive this is to write the sum as
where I’ve introduced the notation
The intermediate sums can be found in closed form mainly using the integral representation of the digamma function, in the form
a sum representation in the form
and using also the definition of the Lerch transcendent and its relation to the hypergeometric function and the incomplete beta function:
The intermediate sum S(54;14,14) can be found using
Using the above relationship between incomplete beta function, the fact that this hypergeometric is a repeated integral of the incomplete beta function (which you know because of the ones on top and matching twos on the bottom) and the fact that there is a nice closed form for the incomplete beta function with a rational first parameter and zero second parameter, like so:
it can be simplified to
The second intermediate sum can be simplified using partial fractions on 1(n+14)(n+12), and the fact that
where I've expanded difference of digamma functions as an infinite sum over k and performed the sum over n. This expression can be handled using the integral representation for ψ mentioned above, giving
Putting everything together, and using Mathematica to evaluate the easier sums, gives the expression I gave above.
Basically, the most important steps in the derivation are the integral representation for ψ(x)−ψ(y) and the explicit formula for B(z,β,0) when β is a rational number. By luck all the integrals then simplify to something Mathematica can do in closed form.
I don't know any standard reference for this; the whole question is basically just computing
and in general derivatives of hypergeometric functions are given by Kampe de Feriet functions, which don't have a closed form in terms of hypergeometric functions. So in general there should be no closed form, but in this case the parameters are just right that all integrals can be done.
When Q(n) is some rational function having a partial fraction expansion ∑kqk(n+αk)−s, the sum
can be expressed in terms of sums
where Φ′ is the Lerch transcendent's derivative. If the sum on the left diverges, it is necessary to consider the asymptotic expansion
and only take the non-divergent term S(β,α) as the value of the sum at z→1 (the divergent terms should be independent of α and cancel out when the partial fractions of Q(n) are added together).
So the problem of evaluating the sum
reduces to calculating asymptotic expansions of integrals of the form