BPS polynomials and Welschinger invariants

We generalize Block-G\"ottsche polynomials, originally defined for toric del Pezzo surfaces, to arbitrary surfaces. To do this, we show that these polynomials arise as special cases of BPS polynomials, defined for any surface $S$ as Laurent polynomials in a formal variable $q$ encoding the BPS invariants of the $3$-fold $S \times \mathbb{P}^1$. We conjecture that for surfaces $S_n$ obtained by blowing up $\mathbb{P}^2$ at $n$ general points, the evaluation of BPS polynomials at $q=-1$ yields Welschinger invariants, given by signed counts of real rational curves. We prove this conjecture for all surfaces $S_n$ with $n \leq 6$.