Geometric singularities of Feynman integrals

We provide a new method to calculate the full microlocal description of the singularities of Feynman integrals. This is done by associating a unique constructible function to the system of partial differential equations (PDEs) annihilating the integral and from this function the singularities can directly be read-off. This function can be constructed explicitly even if the system of PDEs is unknown and describes both the location of the singularities and the number of (generalized) master integrals on them. Our framework is flexible enough to perform the calculation in any of the Lee-Pomeransky, Feynman, or momentum representations.