There seems to exist an asymptotic line
$$\displaystyle a+bx\sim \frac{x e^x}{p_{n}-x e^x}\; ,\;n=\lfloor e^x\rfloor\tag1$$
Which suggests an estimation
$$\displaystyle g(n)=\Big(1+\frac{1}{a+b\ln n}\Big)\cdot n \ln n\tag2$$ of $p_n$.
With $a=4,8461520091$ and $b=0,2473078665$, from the linear regression in the diagram above, $|p_n-g(n)|<|p_n-n\cdot(\ln n+\ln\ln n -1)|$ for $30456<n\leq5761455$, that is up to all primes $<100,000,000$ (the present limit of my prime number table). For $1941$ numbers $n\leq 30456\;\; g_n$ is less accurate.
Is this estimation known?
Are there methods to calculate values of $a$ and $b$?
Can the asymptotic equivalence in $(1)$ be proved?
No, I guess that Vincenzo Oliva on MSE is right. The accuracy comes from a local approximation and b must be a function of n that very slowly approach to zero when $n\to \infty$.