I was playing with an approximation I found in this question: $$\ln(1+e^x) \approx \begin{cases} \frac{x}{1-e^{-\frac{x}{\ln(2)}}},\quad x\neq 0\\ \ln(2),\quad x=0\end{cases}$$
When comparing $f(x) = \ln(1+e^x)$ with $g(x) = x$ you could see that $f(x)$ behaves as a smooth transition from zero to the straight line $g(x)$ as depicted in the plot below by Wolfram-Alpha:
So then I tried to integrate them and see how close their integrals are, which intuition says should approximate a quadratic function $G_{q}(x) = \frac{x^2}{2}$. By taking the integral in Wolfram-Alpha shows it matches $$F(x) = \int\ln(1+e^x)dx = -\text{Li}_2(-e^x)$$ with $\text{Li}_n(x)$ the Polylogarithm function with parameter $n=2$ (later through the answers I realized the function is named Dilogarithm).
When comparing the plots of $F(x)$ with the quadratic $G_q(x)$ I noticed that for big values of $x$ there was always an offset difference, so I took the limit of the difference in Wolfram-Alpha and result it was coincidentally the result of the Basel problem: $$\lim\limits_{x\to\infty} \left[-\text{Li}_2(-e^x)-\frac{x^2}{2}\right] = \frac{\pi^2}{6}$$
and indeed by using $G(x)=\frac{x^2}{2} + \frac{\pi^2}{6}$ the function matches $F(x)$ for big values of $x$ as one could see in Wolfram-Alpha:
With this I found an smooth transition function between zero and a quadratic function through $F(x)$.
But now, my question is: Is this just a coincidence? Or there is a relation to the Basil problem I cannot see?
Added later
Furthermore, $$\int\limits_{-\infty}^0 2\ln(1+e^x)\ dx = \frac{\pi^2}{6}$$
As also by the approximation, $$\int\limits_{-\infty}^0 \dfrac{x}{\ln^2(2)(1-e^{-\frac{x}{\ln(2)}})}\ dx = \frac{\pi^2}{6}$$
And, $$\int\limits_{-\infty}^0 \dfrac{x}{(1-e^{-x})}\ dx = \frac{\pi^2}{6}$$
All of them associated to $\operatorname{Li}_2(e^x)$ or $\operatorname{Li}_2(-e^x)$.

