Approximate minimum of function

527 Views Asked by At

Let: $$f(\Delta)=\sqrt{b \log ^n\left(\frac{E}{\Delta }\right)+4\frac{\Delta ^2}{E^2}},\ \ n=4,6$$ $$E,b,\Delta>0,\ \ \Delta \ll E$$

The function has a single minimum, but the explicit expression is difficult to find. Can an approximate location for the minimum of $f(\Delta)$ be found in terms of $E,b$?

2

There are 2 best solutions below

4
On BEST ANSWER

I'm not sure whether you don't understand that the square root can be dropped, or you don't agree that dropping it simplifies the calculation. The minimum of $f(\Delta)$ will be attained at the same $\Delta$ that minimizes $f^2$, so generally we expect working with $f^2$ to give a little easier work to find the same approximations to $\Delta$.

Here's my take. We substitute $x=\log(E/\Delta)$ into $f^2 = b \log^n \left(\frac{E}{\Delta}\right) + 4 \frac{\Delta^2}{E^2}$:

$$ g(x) = bx^n + 4e^{-2x} $$

The assumption that $\Delta \ll E$ means $x \gt 0$, so we can proceed with the usual setting of $g'(x)$ to zero to identify its minimum (critical point) on $(0,\infty)$.

$$ g'(x) = nb x^{n-1} - 8e^{-2x} $$

Thus $g'(x)=0$ occurs when:

$$ x^{n-1} e^{2x} = \frac{8}{nb} $$

While the left hand side doesn't have a convenient (functional) inverse, for $n=4,6$ it is monotone increasing on $(0,\infty)$. We don't have any information about the magnitude of $b$, but a narrow approximation of the minimizing $x$ (and thus of $\Delta$) can be extracted with a little work. First check to see if $x=1$ is too high or too low, ie. check if $e^2$ is above or below $\frac{8}{nb}$. It seems likely, given that $\Delta \ll E$, that you would expect $x \gt 1$, but this is easy to check.

I would then do some fixed-point iterations of:

$$ h(x) = \frac{\log{8} - \log{nb} - (n-1)\log(x)}{2} $$

starting from $x_0 = 1$, noting that $h(x)=x$ implies $x$ solves the critical point equation for $g'(x)=0$.


I understand that an (easily-computed) analytic approximation for the minimum value is sought, but let me at least point out that the exact location where the minimum occurs can be expressed in terms of the Lambert $W$ function. Taking the $n-1$ root of both sides of the critical point equation:

$$ x e^{\frac{2x}{n-1}} = \sqrt[n-1]{\frac{8}{nb}} $$

Thus $x = \left(\frac{n-1}{2}\right) W\left( \left( \frac{2}{n-1} \right) \sqrt[n-1]{\frac{8}{nb}} \right) $, and from $x = \log E - \log \Delta$ the exact minimizing value of $\Delta$ can be determined.


As we already see from the critical point equation, $x$ will be uniquely determined by $B = 8/nb$, so given that $n=4,6$, the assumption that $b \ll 1$ would say that $B \gg 1$. Of course the left-hand side $x^{n-1}e^{2x}$ is a rapidly increasing monotone function on $[0,\infty)$. One step of our proposed fixed-point iteration ($x_0 = 1$) gives $x_1 = \log(B)/2$, which is an overestimate of the actual root/critical point.

One idea would be to plug this first iteration back into the function, and get a better approximation $x_2 = h(\log(B)/2)$. It is well-known that if $|h'(x_k)| \le \alpha \lt 1$, then the iteration is a contraction map and converges at least linearly. To investigate from this point of view:

$$ h'(x) = \frac{1-n}{2x} $$

so that $|h'(x)| < 1$ in a neighborhood of root $x$ if and only if $x \gt (n-1)/2$.

Recall that the critical point equation allows us to verify easily if this is so:

$$ \text{root } x \gt (n-1)/2 \iff 4^{n-1} e^8 \lt \frac{8}{nb} = B $$

Therefore if this condition holds, e.g. if $b \lt 4^{4-n} e^{-8} n^{-1}$, then a sequence of analytic approximations that converge to the root $x$ is:

$$ \frac{n-1}{2}, h\left(\frac{n-1}{2}\right), h\left(h\left(\frac{n-1}{2}\right)\right), \ldots $$


If $b$ is not small enough, or if the sequence developed above doesn't converge rapidly enough for the (unstated) purpose, then perhaps an alternative approach based on the exact solution $x$ in terms of Lambert W (above) is more attractive. Numerical evaluation of Lambert W by Newton's method or by Halley's method (cubic convergence) turns out to be surprisingly elegant. For series and related iterative methods, see here and here.

1
On

Approximate the function $f$ by some low degree polynomial, whose coefficients will be functions of $E$ and $b$. Locate the minimum of this polynomial, and use it as an approximation of the minimum of $f$.