Let us call
$$L = -\beta + \sum_{i=1}^n b_i\frac{\partial}{\partial x_i};\qquad L^\ast = -\beta - \sum_{i=1}^n b_i\frac{\partial}{\partial x_i}.$$
That a distribution $T$ is a fundamental solution of $L$ means $LT = \delta$, where $\delta$ is the Dirac distribution. With the dual $L^\ast$ of $L$, since by definition
$$\langle LT,\varphi\rangle = \langle T, L^\ast \varphi\rangle,$$
that means $\langle T, L^\ast\varphi\rangle = \varphi(0).$
So you have to compute
$$-\int_0^\infty \left(\beta\varphi(tb) + \sum_{i=1}^n b_i\frac{\partial\varphi}{\partial x_i}(tb)\right)e^{\beta t}\,dt$$
and see that the result is $\varphi(0)$. Integration by parts after seeing that
$$\sum_{i=1}^n b_i\frac{\partial\varphi}{\partial x_i}(tb) = \frac{d}{dt} \varphi(tb)$$
looks like a very promising approach.
With regard to a point brought up by mkl314, the formula
$$\langle T,\varphi\rangle = \int_0^\infty \varphi(tb)e^{\beta t}\,dt$$
does not define a tempered distribution for $\beta > 0$, the integral generally does not converge for $\varphi \in \mathscr{S}(\mathbb{R}^n)$ then. When treating tempered distributions, the cases $\beta < 0$ and $\beta > 0$ must be distinguished ($\beta = 0$ fits both ways), and lead to different fundamental solutions of $L$. However, in $\mathscr{D}'(\mathbb{R}^n)$, no such distinction is necessary; since $\varphi \in \mathscr{D}(\mathbb{R}^n)$ has compact support, the integral converges for all $\beta \in \mathbb{R}$, and defines a distribution. And
$$\langle T_1, \varphi\rangle = -\int_{-\infty}^0 \varphi(tb) e^{\beta t}\,dt$$
gives another fundamental solution of $L$ in $\mathscr{D}'(\mathbb{R}^n)$.