f[n_] := Sum[j^n*3^j/Binomial[2*j, j], {j, 1, Infinity}];
a[n_] := FindIntegerNullVector[{Pi/Sqrt[3], 1, N[-f[n], 20]}][[2]];
Table[Expand[FunctionExpand[FullSimplify[Sum[j^n*3^j/Binomial[2*j, j], {j, 1, Infinity}]]]][[1]], {n, 0, 20}] (*
Vaclav Kotesovec, May 14 2020 *)
S[k_, z_] := Sum[n!*(z/(4 - z))^n* StirlingS2[k + 1, n]*(1/n + Sum[(-1)^p*Pochhammer[1/2, p]/(p + 1)!* Binomial[n - 1, p]*(4/z)^(p + 1)*(Sqrt[z/(4 - z)]*ArcSin[Sqrt[z]/2] - 1/2*Sum[Gamma[l]/Pochhammer[1/2, l]*(z/4)^l, {l, 1, p}]), {p, 0, n - 1}]), {n, 1, k + 2}]; Table[Expand[Simplify[S[j, 3]]][[1]], {j, 0, 20}] (*
Vaclav Kotesovec, May 15 2020 *)