La fórmula de Riemann-Liouville para la semiderivada es la siguiente:
$${}_0 D_x^{(1/2)} f(x)=\frac{f(0)}{\sqrt{\pi x}}+\frac1{\sqrt\pi}\int_0^x \frac{f^\prime(t)}{\sqrt{x-t}}\mathrm dt$$
Tomando $f(x)=\exp(-\alpha x^2+\beta x)$ entonces tenemos
$$\begin{align*} {}_0 D_x^{(1/2)} f(x)&=\frac1{\sqrt{\pi x}}+\frac1{\sqrt\pi}\int_0^x \frac{(\beta-2\alpha t)\exp(-\alpha t^2+\beta t)}{\sqrt{x-t}}\mathrm dt\\ &=\frac1{\sqrt{\pi x}}+\frac1{\sqrt\pi}\left(\beta\int_0^x \frac{\exp(-\alpha t^2+\beta t)}{\sqrt{x-t}}\mathrm dt-2\alpha\int_0^x \frac{t\exp(-\alpha t^2+\beta t)}{\sqrt{x-t}}\mathrm dt\right) \end{align*}$$
Por desgracia, ninguna de las dos integrales parece tener una evaluación fácil en términos de funciones conocidas, a menos que $\alpha=0$ o $\beta=0$ En ese caso,
$$\begin{align*} {}_0 D_x^{(1/2)} \exp(-\alpha x^2)&=\frac1{\sqrt{\pi x}}-\frac{8\alpha x^{3/2}}{3\sqrt\pi}{}_2 F_2\left({{1,\frac32}\atop{\frac54,\frac74}}\mid -\alpha x^2\right)\\ {}_0 D_x^{(1/2)} \exp(\beta x)&=\frac1{\sqrt{\pi x}}+\sqrt{\beta}\exp(\beta x)\mathrm{erf}(\sqrt{\beta x}) \end{align*}$$
De hecho, existe una versión semiderivada de la fórmula de Leibniz, pero, como ya he dicho, parece algo difícil de manejar:
$${}_0 D_x^{(1/2)}f(x)g(x)=\sum_{j=0}^\infty \binom{\frac12}{j}\left({}_0 D_x^{(1/2-j)}f(x)\right)g^{(j)}(x)$$
Si tomamos $f(x)=\exp(\beta x)$ y $g(x)=\exp(-\alpha x^2)$ obtenemos
$$\begin{split}&{}_0 D_x^{(1/2)} \exp(-\alpha x^2+\beta x)=\\&\quad\exp(-\alpha x^2+\beta x)\sum_{j=0}^\infty \binom{\frac12}{j}\frac{(-1)^j \alpha^{j/2} \beta^{1/2-j}}{\Gamma(j-1/2)}\gamma(j-1/2,\beta x) H_j(x\sqrt\alpha)\end{split}$$
donde $\gamma(a,x)$ es una función gamma incompleta (inferior), y $H_j(x)$ es un polinomio de Hermite. Sin embargo, no conozco ninguna forma más sencilla de expresar esta suma infinita.