Suppose we have a continuous strictly increasing function $f:\mathbb R^+\to\mathbb R^+$ such that $f(f(x))=\mathcal O(\exp(x))$ for $x\to\infty$.${^{[1]}}$ In other words, the function $f(x)$ exhibits (at most) half-exponential${^{[2]}}$$\!{^{[3]}}$$\!{^{[4]}}$ growth rate. It is known${^{[5]}}$ that there is no function that is built from arithmetical operators $(+,-,\times,/)$, $\exp$, $\log$, and real numbers that has exactly half-exponential growth rate. But, of course, there are elementary functions that are upper bounds for such functions, e.g. $\exp(x)$ is a very simple example of a (very loose) upper bound.
I'm interested in finding an elementary function $g:\mathbb R^+\to\mathbb R^+$ such that $f(x)=\mathcal O(g(x))$ that provides a much tighter upper bound (I mean an asymptotic upper bound, up to a coefficient). Ideally, I would like to find an infinite sequence of progressively tighter upper bounds that can approximate the half-exponential growth rate as close as we wish. If necessary, we can relax the requirement that the function must be elementary, and allow use of a well-known analytic special functions, such as Lambert $W$-function.
How can I approach this problem?