BFSS model is a theory of super-symmetric matrix quantum mechanics describing $N$ coincident $D0$-branes, defined by the action
$$S=\frac{1}{g^2}\int dt\ \text{Tr}\left\{ \frac{1}{2}(D_t X^I)^2 + \frac{1}{2}\psi_\alpha D_t \psi_\alpha + \frac{1}{4}[X^I,X^J]^2 + \frac{1}{2}i\psi_\alpha \gamma_{\alpha\beta}^I[\psi_\beta,X^I]\right\}.$$
You can read about the background and all the details regarding the indices, the gauge field, the spinor representations in a nice paper by Maldacena, To gauge or not to gauge? See section 2 specifically.
Now this is not a quantum field theory in $d+1$ dimensions. This is simply quantum mechanics in $0+1$ dimensions. In quantum mechanics we can similarly define Feynman diagrams for perturbative calculations, but this is not something known as commonly as Feynman diagram techniques in quantum field theory. See for example this nice review by Abbott, Feynman diagrams in quantum mechanics.
When I did a literature search, I couldn't find any results regarding Feynman rules for the BFSS model. Is this something that has never been considered? Or is it just too hard to find the propagators and vertex factors? How can we go about calculating the free propagator? Can we look at an expansion for the full propagator in the large-$N$ limit?
It feels as if we need to introduce the 1PI diagrams, and construct something like this Equation 7.22 in Peskin & Schroeder: writing the Fourier transform of a two-point function as a series of 1PI diagrams.
Is there an effective action for the BFSS? Maybe that can simplify the computations?