Proving that $\sum_{n=1}^\infty \frac{\sin\left(n\frac{\pi}{3}\right)}{(2n+1)^2}=\frac{G}{\sqrt 3} -\frac{\pi^2}{24}$
\begin{align*} I&=\int_0^1 \frac{\sqrt{x}\ln x} {x^2-x+1}\,dx \end{align*}
Perform the change of variable $y=\sqrt{x}$,
\begin{align*} J&=4\int_0^1 \frac{x^2\ln x} {x^4-x^2+1}\,dx\\ &=\left[\left(\frac{1}{\sqrt{3}}\ln\left(\frac{x^2-\sqrt{3}x+1}{x^2+\sqrt{3}x+1}\right)+2\arctan\left(\frac{x}{1-x^2}\right)\right)\ln x\right]_0^1-\\ &\int_0^1 \frac{1}{x\sqrt{3}}\ln\left(\frac{x^2-\sqrt{3}x+1}{x^2+\sqrt{3}x+1}\right)\,dx-2\int_0^1\frac{1}{x}\arctan\left(\frac{x}{1-x^2}\right)\,dx\\ &=-\int_0^1 \frac{1}{x\sqrt{3}}\ln\left(\frac{x^2-\sqrt{3}x+1}{x^2+\sqrt{3}x+1}\right)\,dx-2\int_0^1\frac{1}{x}\arctan\left(\frac{x}{1-x^2}\right)\,dx\\ \end{align*}
\begin{align*} L&=\int_0^1 \frac{1}{x}\ln\left(\frac{x^2-\sqrt{3}x+1}{x^2+\sqrt{3}x+1}\right)\,dx\\ &=\int_0^1 \frac{1}{x}\ln\left(\frac{\left(x^2-\sqrt{3}x+1\right)\left(x^2+\sqrt{3}x+1\right)}{\left(x^2+\sqrt{3}x+1\right)^2}\right)\,dx\\ &=\int_0^1 \frac{1}{x}\ln\left(x^4-x^2+1\right)\,dx-2\int_0^1 \frac{1}{x}\ln\left(x^2+\sqrt{3}x+1)\right)\,dx\\ &=\int_0^1 \frac{x}{x^2}\ln\left(\frac{1+x^6}{1+x^2}\right)\,dx-2\int_0^1 \frac{1}{x}\ln\left(x^2+\sqrt{3}x+1)\right)\,dx\\ \end{align*}
In the first integral perform the change of variable $y=x^2$,
\begin{align*}L&=\frac{1}{2}\int_0^1 \frac{1}{x}\ln\left(\frac{1+x^3}{1+x}\right)\,dx-2\int_0^1 \frac{1}{x}\ln\left(x^2+\sqrt{3}x+1\right)\,dx\\ &=\frac{1}{2}\int_0^1 \frac{x^2}{x^3}\ln\left(1+x^3\right)\,dx-\frac{1}{2}\int_0^1 \frac{1}{x}\ln\left(1+x\right)\,dx-2\int_0^1 \frac{1}{x}\ln\left(x^2+\sqrt{3}x+1\right)\,dx\\ \end{align*}
In the first integral perform the change of variable $y=x^3$,
\begin{align*}L&=\frac{1}{6}\int_0^1 \frac{1}{x}\ln\left(1+x\right)\,dx-\frac{1}{2}\int_0^1 \frac{1}{x}\ln\left(1+x\right)\,dx-2\int_0^1 \frac{1}{x}\ln\left(x^2+\sqrt{3}x+1\right)\,dx\\ &=-\frac{1}{3}\int_0^1 \frac{1}{x}\ln\left(1+x\right)\,dx-2\int_0^1 \frac{1}{x}\ln\left(x^2+\sqrt{3}x+1\right)\,dx\\ &=\frac{1}{3}\int_0^1 \frac{1}{x}\ln\left(1-x\right)\,dx-\frac{1}{3}\int_0^1 \frac{1}{x}\ln\left(1-x^2\right)\,dx-2\int_0^1 \frac{1}{x}\ln\left(x^2+\sqrt{3}x+1\right)\,dx\\ \end{align*}
In the second integral perform the change of variable $y=x^2$,
\begin{align*}L&=\frac{1}{3}\int_0^1 \frac{1}{x}\ln\left(1-x\right)\,dx-\frac{1}{6}\int_0^1 \frac{1}{x}\ln\left(1-x\right)\,dx-2\int_0^1 \frac{1}{x}\ln\left(x^2+\sqrt{3}x+1\right)\,dx\\ &=\frac{1}{6}\int_0^1 \frac{1}{x}\ln\left(1-x\right)\,dx-2\int_0^1 \frac{1}{x}\ln\left(x^2+\sqrt{3}x+1\right)\,dx\\ \end{align*}
\begin{align*}M&=\int_0^1 \frac{1}{x}\ln\left(x^2+\sqrt{3}x+1\right)\,dx\\\end{align*}
Define the function $F$ on $\left[0;\sqrt{3}\right]$ by,
\begin{align*}F(a)=\int_0^1 \frac{1}{x}\ln\left(x^2+ax+1\right)\,dx\\\end{align*}
Observe that $F(\sqrt{3})=M$ and,
\begin{align*}F(0)&=\int_0^1 \frac{1}{x}\ln\left(x^2+1\right)\,dx\\ &=-\frac{1}{4}\int_0^1 \frac{1}{x}\ln\left(1-x\right)\,dx\\ \end{align*}
Since $0<a<2$,
\begin{align*}F^\prime(a)&=\int_0^1 \frac{1}{x^2+ax+1}\,dx\\ &=\left[\frac{2}{\sqrt{4-a^2}}\arctan\left(\frac{2x+a}{\sqrt{4-a^2}}\right)\right]_0^1\\ &=\frac{2}{\sqrt{4-a^2}}\arctan\left(\frac{2+a}{\sqrt{4-a^2}}\right)-\frac{2}{\sqrt{4-a^2}}\arctan\left(\frac{a}{\sqrt{4-a^2}}\right)\\ &=\frac{2}{\sqrt{4-a^2}}\arctan\left(\sqrt{\frac{2-a}{2+a}}\right)\\ \end{align*}
\begin{align*}F(\sqrt{3})-F(0)=\int_0^{\sqrt{3}}\frac{2}{\sqrt{4-a^2}}\arctan\left(\sqrt{\frac{2-a}{2+a}}\right)\,da\\ \end{align*}
Perform the change of variable $y=\sqrt{\frac{2-a}{2+a}}$,
\begin{align*}F(\sqrt{3})-F(0)&=4\int_{2-\sqrt{3}}^1 \frac{\arctan y}{1+y^2}\,dy\\ &=2\Big[\arctan^2(x)\Big]_{2-\sqrt{3}}^1\\ &=2\times \frac{\pi^2}{4^2} -2\times \frac{\pi^2}{12^2}\\ &=\frac{\pi^2}{9} \end{align*}
Since,
\begin{align*}\int_0^1 \frac{1}{x}\ln\left(1-x\right)\,dx=-\frac{\pi^2}{6}\end{align*}
Then, \begin{align*}M&=F(\sqrt{3})\\ &=F(0)+\frac{\pi^2}{9}\\ &=-\frac{1}{4}\times -\frac{\pi^2}{6}+\frac{\pi^2}{9}\\ &=\frac{11\pi^2}{72} \end{align*}
Therefore,
\begin{align*}L&=-\frac{1}{6}\times -\frac{\pi^2}{6}-2M\\ &=\frac{1}{6}\times -\frac{\pi^2}{6}-2\times \frac{11\pi^2}{72}\\ &=-\frac{\pi^2}{3} \end{align*}
\begin{align*} K&=\int_0^1 \frac{\arctan\left(\frac{x}{1-x^2}\right)}{x}\,dx\\ \end{align*}
Perform the change of variable $x=\tan\left(\frac{t}{2}\right) $,
\begin{align*} K&=\int_0^{\frac{\pi}{2}} \frac{\arctan\left(\frac{1}{2}\tan t\right)}{\sin t}\,dt \end{align*}
Define the function $H$ on $\left[\frac{1}{2};1\right]$ by,
\begin{align*}H(a)&=\int_0^{\frac{\pi}{2}} \frac{\arctan\left(a\tan t\right)}{\sin t}\,dt\end{align*}
Observe that $K=H\left(\dfrac{1}{2}\right)$ and,
\begin{align*}H(1)&=\int_0^{\frac{\pi}{2}} \frac{t}{\sin t}\,dt\\ &=\Big[t\ln\left(\tan\left(\frac{t}{2} \right)\right)\Big]_0^{\frac{\pi}{2}}-\int_0^{\frac{\pi}{2}}\ln\left(\tan\left(\frac{t}{2} \right)\right)\,dt\\ &=-\int_0^{\frac{\pi}{2}}\ln\left(\tan\left(\frac{t}{2} \right)\right)\,dt\\ \\\end{align*}
Perform the change of variable $x=\dfrac{t}{2}$,
\begin{align*}H(1)&=-2\int_0^{\frac{\pi}{4}}\ln\left(\tan\left(t \right)\right)\,dt\\ &=2\text{G} \\\end{align*}
\begin{align*}H^\prime (a)&=\int_0^{\frac{\pi}{2}} \frac{\cos x}{1-(1-a^2)\sin^2 x}\,dt\\ &=\left[\frac{1}{2\sqrt{1-a^2}}\ln\left(\frac{1+\sin(x)\sqrt{1-a^2}}{1-\sin(x)\sqrt{1-a^2}}\right)\right]_0^{\frac{\pi}{2}}\\ &=\frac{1}{2\sqrt{1-a^2}}\ln\left(\frac{1+\sqrt{1-a^2}}{1-\sqrt{1-a^2}}\right) \end{align*}
Therefore,
\begin{align*}H(1)-H\left(\frac{1}{2}\right)&=\int_{\frac{1}{2}}^1 \frac{1}{2\sqrt{1-a^2}}\ln\left(\frac{1+\sqrt{1-a^2}}{1-\sqrt{1-a^2}}\right)\,da\end{align*}
Perform the change of variable $y=\arctan\left(\sqrt{\dfrac{1+\sqrt{1-x^2}}{1-\sqrt{1-x^2}}}\right)$
\begin{align*}H(1)-H\left(\frac{1}{2}\right)&=-2\int_{\frac{\pi}{12}}^{\frac{\pi}{4}} \ln\left(\tan y\right)\,dy\\ &=-2\int_0^{\frac{\pi}{4}} \ln\left(\tan y\right)\,dy+\int_0^{\frac{\pi}{12}} \ln\left(\tan y\right)\,dy \end{align*}
But, it's well known that:
\begin{align*}\int_0^{\frac{\pi}{4}} \ln\left(\tan y\right)\,dy=-\text{G}\\ \int_0^{\frac{\pi}{12}} \ln\left(\tan y\right)\,dy=-\frac{2}{3}\text{G}\\ \end{align*}
(see Integral: $\int_0^{\pi/12} \ln(\tan x)\,dx$ )
Therefore,
\begin{align*}H(1)-H\left(\frac{1}{2}\right)&=2\text{G}+2\times -\frac{2}{3}\text{G}\\ &=\frac{2}{3}\text{G} \end{align*}
Thus,
\begin{align*} K&=H\left(\frac{1}{2}\right)\\ &=H(1)-\frac{2}{3}\text{G}\\ &=2\text{G}-\frac{2}{3}\text{G}\\ &=\frac{4}{3}\text{G} \end{align*}
\begin{align*}I&=-\frac{1}{\sqrt{3}}L-2K\\ &=-\frac{1}{\sqrt{3}}\times -\frac{\pi^2}{3}-2\times \frac{4}{3}\text{G}\\ &=\boxed{\frac{\sqrt{3}\pi^2}{9}-\frac{8}{3}\text{G}} \end{align*}
PS:
Actually,
\begin{align}\int_0^1 \frac{\arctan\left( \frac{x}{1-x^2}\right)}{x}\,dx-\int_0^1 \frac{\arctan x}{x}\,dx=\int_0^1 \frac{\arctan \left(x^3\right)}{x}\,dx\end{align}
In the latter integral perform the change of variable $y=x^3$,
\begin{align}\int_0^1 \frac{\arctan\left( \frac{x}{1-x^2}\right)}{x}\,dx-\int_0^1 \frac{\arctan x}{x}\,dx=\frac{1}{3}\int_0^1 \frac{\arctan x}{x}\,dx\end{align}
therefore,
\begin{align}\int_0^1 \frac{\arctan\left( \frac{x}{1-x^2}\right)}{x}\,dx&=\frac{1}{3}\int_0^1 \frac{\arctan x}{x}\,dx+\int_0^1 \frac{\arctan x}{x}\,dx\\ &=\frac{4}{3}\int_0^1 \frac{\arctan x}{x}\,dx\\ &=\frac{4}{3}\text{G} \end{align}
This is a major revision of the last few lines. The conclusion is that
$\sum_{n=1}^\infty \frac{\sin\left(n\frac{\pi}{3}\right)}{(2n+1)^2} =\frac{\sqrt{3}}{72}\psi^{(1)}(\frac56)-\frac{\sqrt{3}\pi^2}{144} $
where $\psi^{(1)}$ is a polygamma function (reference below).
$\begin{array}\\ \sum_{n=1}^\infty \frac{\sin\left(n\frac{\pi}{3}\right)}{(2n+1)^2} &=\sum_{n=1}^\infty \left(\frac{\sin\left((3n-2)\frac{\pi}{3}\right)}{(2(3n-2)+1)^2}+\frac{\sin\left((3n-1)\frac{\pi}{3}\right)}{(2(3n-1)+1)^2}+\frac{\sin\left(3n\frac{\pi}{3}\right)}{(2(3n)+1)^2}\right)\\ &=\sum_{n=1}^\infty \left(\frac{\sin\left(-2\frac{\pi}{3}\right)}{(6n-3)^2}+\frac{\sin\left(-\frac{\pi}{3}\right)}{(6n-1)^2}+\frac{\sin\left(n\pi\right)}{(6n+1)^2}\right)\\ &=\sum_{n=1}^\infty \left(-\frac{\sqrt{3}/2}{(6n-3)^2}+\frac{\sqrt{3}/2}{(6n-1)^2}\right)\\ &=\frac{\sqrt{3}}{2}\sum_{n=1}^\infty \left(-\frac1{(6n-3)^2}+\frac1{(6n-1)^2}\right)\\ &=\frac{\sqrt{3}}{2}\sum_{n=1}^\infty \left(\frac1{(6n-1)^2}-\frac1{(6n-3)^2}\right)\\ &=\frac{\sqrt{3}}{2}\left(\sum_{n=1}^\infty \frac1{(6n-1)^2}-\frac19\sum_{n=1}^\infty\frac1{(2n-1)^2}\right)\\ &=\frac{\sqrt{3}}{2}\left(\frac1{36}\sum_{n=1}^\infty \frac1{(n-1/6)^2}-\frac19(1-\frac14)\zeta(2)\right)\\ &=\frac{\sqrt{3}}{2}\left(\frac1{36}\sum_{n=0}^\infty \frac1{(n+5/6)^2}-\frac19\frac34\frac{\pi^2}{6}\right)\\ &=\frac{\sqrt{3}}{72}\psi^{(1)}(\frac56)-\frac{\sqrt{3}\pi^2}{144}\\ \end{array} $
$\psi^{(1)}(z)$ is a polygamma function: https://en.wikipedia.org/wiki/Polygamma_function
$\psi^{(m)}(z) =(-1)^{m+1}m!\sum_{n=0}^{\infty} \dfrac1{(z+n)^{m+1}} $ so that $\psi^{(1)}(z) =\sum_{n=0}^{\infty} \dfrac1{(z+n)^2} $
Also
$\begin{array}\\ \zeta(m) &=\sum_{n=1}^{\infty} \dfrac1{n^m}\\ &=\sum_{n=1}^{\infty} \dfrac1{(2n-1)^m}+\sum_{n=1}^{\infty} \dfrac1{(2n)^m}\\ &=\sum_{n=1}^{\infty} \dfrac1{(2n-1)^m}+\dfrac1{2^m}\sum_{n=1}^{\infty} \dfrac1{n^m}\\ &=\sum_{n=1}^{\infty} \dfrac1{(2n-1)^m}+\dfrac1{2^m}\zeta(m)\\ \end{array} $
so
$\zeta(m)(1-2^{-m}) =\sum_{n=1}^{\infty} \dfrac1{(2n-1)^m} $.
Integrate by parts \begin{align} \int_0^1 \frac{\sqrt x \ln x}{x^2-x+1}dx =& \int_0^1 \ln x \>d\left( \tan^{-1} \frac{\sqrt x} {1-x} -\frac1{\sqrt3 } \tanh^{-1} \frac{\sqrt {3x}} {1+x} \right)\\ =&\frac1{\sqrt3 }I_1 - I_2\tag1 \end{align} where \begin{align} I_1&=\int_0^1 \frac{\tanh^{-1} \frac{\sqrt {3x}} {1+x}}x dx\\ I_2 & = \int_0^1 \frac{\tan^{-1} \frac{\sqrt {x}} {1-x}}x dx = \underset{\sqrt x\to x}{\int_0^1 \frac{\tan^{-1} \sqrt x}x dx } + \underset{\sqrt {x^3}\to x}{\int_0^1 \frac{\tan^{-1} \sqrt {x^3}}x }dx \\ &= \left(2+ \frac23 \right) \int_0^1 \frac{\tan^{-1} x}x dx = \frac83G\tag2 \end{align} Evaluate $I_1$ with $J(a) = \int_0^1 \frac{\tanh^{-1} \frac{2a\sqrt {x}} {1+x}}x dx$ $$J’(a) =\int_0^1 \frac{2(\frac1{\sqrt x}+\sqrt x)dx}{(x+1)^2-(2a\sqrt x)^2} =\frac{2 \tan^{-1}\frac{\sqrt{(1-a^2)x}}{1-x}\bigg|_0^1} {\sqrt{1-a^2}}=\frac\pi{\sqrt{1-a^2}} $$ Then $$I_1 = J(\frac{\sqrt3}2)=\int_0^{ \frac{\sqrt3}2}J’(a)da =\int_0^{ \frac{\sqrt3}2} \frac\pi{\sqrt{1-a^2}} da =\frac{\pi^2}3\tag3 $$ Plug (2) and (3) into (1) to obtain $$\int_0^1 \frac{\sqrt x \ln x}{x^2-x+1}dx = \frac{\pi^2}{3\sqrt3}-\frac83G$$