We initiate a new study of differential operators with symmetries and combine this with the study of branching laws for Verma modules of reductive Lie algebras. By the criterion for discretely decomposable and multiplicity-free restrictions of generalized Verma modules (T.
Kobayashi (2012) [22]), we are brought to natural settings of parabolic geometries for which there exist unique equivariant differential operators to submanifolds. Then we apply a new method (F-method) relying on the Fourier transform to find singular vectors in generalized Verma modules, which significantly simplifies and generalizes many preceding works.
In certain cases, it also determines the Jordan-Hölder series of the restriction for singular parameters. The F-method yields an explicit formula of such unique operators, for example, giving an intrinsic and new proof of Juhl's conformally invariant differential operators (Juhl (2009) [16]) and its generalizations to spinor bundles.
This article is the first in the series, and the next ones include their extension to curved cases together with more 57applications of the F-method to various settings in parabolic geometries.