Let $f=\sum_{n\geq 1} a_n q^n$ be a normalized eigenform which is supersingular and crystalline at a prime $p$ and let $V_f$ be the associated crystalline representation, then it follows from the work of Scholl and Faltings that the characteristic polynomial of the $\Phi$ operator on $D_{cris}(V_f^*)$ is $X^2-a_p X+p^{k-1}$. This is alluded to in "Explicit reduction modulo p of certain $2$-dimensional crystalline representations"- Buzzard and Gee.
I've been trying to find the location of this result by looking through some of Scholl's papers, in which paper can it be found?