Given a graph $G$ on $n$ vertices, its chromatic polynomial $P(G,x)$ is a function that gives the number of proper colorings of G using $x$ colors.
When $P(G,x)$ is written using the basis $\{x, \ldots , x^n\}$, it is known that the coefficients of $x, \ldots , x^n$ alternate in signs; when $P(G,x)$ is written using the basis $\{ (x)_0, (x)_1, \ldots, (x)_n \}$, where $(x)_k = x(x-1)(x-2) \cdot \ldots \cdot (x-(k-1))$ denotes the falling factorial, then its coefficients are nonnegative.
I verified at SageMath that, up to $9$ vertices, when $P(G,x)$ is written using the basis $\{ x^{n-1}(x)_1, x^{n-2}(x)_2, \ldots, (x)_n \}$, then its coefficients are also nonnegative, but I could not find an example of graph with at least one negative coefficient and I have not found any proof of this fact in the literature.
Is anyone aware of this result? This question arises from work in progress by Petter Brändén and Leonardo Saud.