Genus $g$ Virasoro Correlation Functions for Vertex Operator Algebras

For a simple, self-dual, strong CFT-type vertex operator algebra (VOA) of central charge $c$, we describe the Virasoro $n$-point correlation function on a genus $g$ marked Riemann surface in the Schottky uniformisation. We show that this $n$-point function determines the correlation functions for all Virasoro vacuum descendants. Using our recent work on genus $g$ Zhu recursion, we show that the Virasoro $n$-point function is determined by a differential operator $\mathcal{D}_{n}$ acting on the genus $g$ VOA partition function normalised by the Heisenberg partition function to the power of $c$. We express $\mathcal{D}_{n}$ as the sum of weights over certain Virasoro graphs where the weights explicitly depend on $c$, the classical bidifferential of the second kind, the projective connection, holomorphic 1-forms and derivatives with respect to any $3g-3$ locally independent period matrix elements. We also describe the modular properties of $\mathcal{D}_{n}$ under a homology base change.