Consider the elliptic second-order PDE on some bounded domain $D$ (in any dimension)
$-\Delta_D u + \alpha u = 0$
subject to the constraint
$\gamma^i \nabla_i u + \beta u = 0$,
where $\Delta_D$ is the Laplacian on $D$ (with respect to some metric), and $\alpha$, $\beta$, and $\gamma^i$ are real smooth functions on $D$.
A zeroth question, to address user254433's comment below, is what compatibility conditions on $\alpha$, $\beta$, $\gamma^i$, and presumably the metric on $D$ are required to ensure that the system above has solutions? user254433 derived such a condition in one dimension.
Now my main question. Assume the system satisfies the above compatibility conditions so that it admits solutions, and say if any (real) solution to the above system is positive on the boundary $\partial D$, it must be positive everywhere in $D$. What does this imply on the coefficients $\alpha$, $\beta$, and $\gamma^i$?
Note that here I'm asking for necessary conditions on the coefficients. Indeed, I do know of a sufficient condition: if $\alpha \geq 0$ everywhere, this ensures by standard minimum principles that $u$ cannot have any negative local minima (since otherwise at the minimum we'd have $-\Delta_D u + \alpha u < 0$), and so if it's positive on $\partial D$ it must be positive everywhere. However, if necessary conditions on the coefficients are unknown, I'm also interested in learning about sufficient conditions that are weaker than $\alpha \geq 0$ everywhere.
Note that these conditions (necessary or sufficient) need not be local; for instance I'm happy to have conditions on integrals of these coefficients over $D$, etc.