Given a vector function $$f=(f_1,\ldots,f_n)\in L^2(\mathbb R,\mathbb R^n)$$ (for some $n\in\mathbb N$), let us define $$\Delta f:=(\Delta f_1,\ldots,\Delta f_n),$$ where $\Delta$ is the Laplacian operator, and let $Q:\mathbb R\to\mathbb R^{n\times n}$ be a potential taking values in symmetric $n\times n$ matrices. I'm interested in vector Schrödinger operators of the form $$Hf=-\Delta f+Qf,\qquad f\in L^2(\mathbb R,\mathbb R^n).$$
Question. Is there a Feynman-Kac type formula known for $H$'s semigroup in the case where $Q(x)$ is not necessarily diagonal?
(Note: I specify abote that I'm interested in the case where $Q$ is not diagonal; if $Q=\mathrm{diag}(Q_1,\ldots,Q_n)$, then $$Hf=(-\Delta f_1+Q_1 f_1,\ldots,-\Delta f_1+Q_1 f_1),$$ in which case we can simply apply the one-dimensional Feynman-Kac formula to each component.)