4

Let $p\in \mathbb{R}$, I would like to investigate the asymptotic behavior of the following integral:

$$\int_0^1 e^{x(tp-\frac{1}{2}t^2)}dt$$ as $x\to\infty$. In particular, I would like to know how it depends on $p$ and $x$?.

If $p\in (0,1)$, then the maximum of $t\mapsto at-\frac{1}{2}t^2$ over $[0,1]$ is achieved at $p\in (0,1)$. Hence applying Laplace's method (e.g., here), $$ \int_0^1 e^{x(tp-\frac{1}{2}t^2)}dt\sim \sqrt{\frac{2\pi}{x}}e^{x\frac{p^2}{2}} $$ as $x\to \infty$. In this case, is it possible to quantify how the error term $$\left|\int_0^1 e^{x(tp-\frac{1}{2}t^2)}dt-\sqrt{\frac{2\pi}{x}}e^{x\frac{p^2}{2}} \right|$$ depends on $p$?

Alternatively, if $p\not \in (0,1)$, then the maximum of $t\mapsto at-\frac{1}{2}t^2$ over $[0,1]$ is obtained at boundary, and the first order derivative at the maximizer is not zero. In this case, I am not sure how to apply the classical result.

Could someone help to explain how to study the asymptotic behavior of the integral for different $p$?

John
  • 13,738

2 Answers2

2

I shall mainly address your first question. Your integral can be re-written in the form $$ \exp \left( {\frac{1}{2}p^2 x} \right)\int_0^1 {\exp \left( { - \frac{x}{2}(t - p)^2 } \right)dt} = \sqrt {\frac{{2\pi }}{x}} \exp \left( {\frac{1}{2}p^2 x} \right)(1 + R(x,p)), $$ where $$ R(x,p) = - \frac{1}{2}\left( {\operatorname{erfc}\left( {p\sqrt {\frac{x}{2}} } \right) + \operatorname{erfc}\left( {(1 - p)\sqrt {\frac{x}{2}} } \right)} \right), $$ with $\operatorname{erfc}$ being the complementary error function. By $\S7.12(\text{i})$, we have that $$ \left| {R(x,p)} \right| \le \frac{1}{{\sqrt {2\pi x} }}\left( {\frac{1}{p}\exp \left( { - \frac{1}{2}p^2 x} \right) + \frac{1}{{1 - p}}\exp \left( { - \frac{1}{2}(1 - p)^2 x} \right)} \right) $$ provided $x>0$ and $0<p<1$. If $p\to 0$ or $p\to 1$, the saddle coalesces with an endpoint. In that case, the asymptotics is given by the exact form of $R(x,p)$. See $\S2.3(\text{v})$ for more details. Outside of $(0,1)$, one may use $\operatorname{erfc}(-w)=2-\operatorname{erfc}(w)$.

Gary
  • 36,640
  • Thanks for your answer. Are you suggesting that for $p\not \in (0,1)$, we should also factor out $\exp(\frac{1}{2}p^2x)$ in the integral? For $p\in (0,1)$ it makes sense, as Laplace's method predicts the leading term is $\exp(\frac{1}{2}p^2x)$. But if $p\not \in (0,1)$, say $p>1$, formally I think the dominating part of the integral is near the maximiser $1$? In this case, shall we factor out $\exp(\frac{1}{2}1^2x)$? – John Feb 08 '24 at 07:22
  • Here I focused mainly on the case $0<p<1$. If $p>1$, $$ R(x,p) = - 1 + \frac{1}{2}\left( \operatorname{erfc}\left( {(p - 1)\sqrt {\frac{x}{2}} } \right)-\operatorname{erfc}\left( p\sqrt {\frac{x}{2}} \right) \right). $$ You can use this together with the asymptotic behaviour of $\operatorname{erfc}$ for large positive argument. – Gary Feb 08 '24 at 07:32
  • It is not clear to me how to use the asymptotic behavior of erfc when $p>1$, as both terms go to zero. Could you elaborate it further in your answer? – John Feb 08 '24 at 08:09
  • 1
    The second erfc term is exponentially smaller than the first one. Just substitute the asymptotics $(7.12.1)$ given in the DLMF. – Gary Feb 08 '24 at 08:36
1

$$I=\int_0^1e^{x \left(p t-\frac{t^2}{2}\right)}\,dt=\sqrt{\frac{\pi }{2x}}\,\,e^{\frac{p^2 x}{2}} \Bigg(\text{erf}\left(\frac{(1-p) \sqrt{x}}{\sqrt{2}}\right) -\text{erf}\left(\frac{p \sqrt{x}}{\sqrt{2}}\right)\Bigg)$$ If $x\to \infty$ and $p>1$, a detailed asymptotic is $$e^{\frac{(2 p-1) x}{2} }\Bigg(\frac{1}{(p-1) x}-\frac{1}{(p-1)^3 x^2}+O\left(\frac{1}{x^3}\right)\Bigg)-\Bigg(\frac{1}{p x}-\frac{1}{p^3 x^2}+O\left(\frac{1}{x^3}\right) \Bigg)$$

Using $p=3$ and $x=10$ the above gives $$\frac{39 }{800}e^{25}-\frac{89}{2700}=3.51024\times 10^9$$ to be compared to the exact $$\frac{1}{2} \sqrt{\frac{\pi }{5}}\, e^{45}\left(\text{erf}\left(3 \sqrt{5}\right)-\text{erf}\left(2 \sqrt{5}\right)\right)=3.51627\times 10^9$$

I suppose that you notice how small is $\frac{89}{2700}$ compared to the term on his left.

So, I presume that $$I \sim \frac {e^{\frac{(2 p-1) x}{2} } }{(p-1)x }$$ would be more than sufficient.

  • Thank you for elaborating the details for $p>1$. I eventually accepted Gary's answer as his answer pointed out the connection to the complementary error function first. – John Feb 09 '24 at 14:29