Let $\mathbb {X}$
bea Riemannian symmetric space of non-compact type. For a real-valued $f \in L^1(\mathbb {X})$
with $\int _{\mathbb {X}} f \ge 0$
, we prove that if f satisfies the iterated convolution inequality: $$ \begin{align*} f \ge \sum_{n=2}^N a_n \left(*^n f\right),\:\:\text{ a.e. on } \mathbb{X}, \end{align*} $$
where $N \ge 2$
is an integer and for $2 \le n \le N$
, $a_n$
are nonnegative integers with at least one of them positive, then f must be nonnegative a.e. and satisfy the nontrivial bound $\int _{\mathbb {X}} f \le t_{\mathcal {Q}}\:$
, where $\mathcal {Q}(t):=t-\sum _{n=2}^N a_n\:t^n$
and $t_{\mathcal {Q}}$
is the unique zero of $\mathcal {Q}'$
in $(0,\infty )$
. This result is new even for $\mathbb {R}^d$
and generalizes recent results of Carlen et al. (2021, Int. Math. Res. Not., 24, 18604–18612) and Nakamura and Sawano (2025, J. Geom. Anal., 35, 68) for two-fold and m-fold convolutions, respectively, to genuine polynomials. We then apply our result to obtain an a priori estimate for solutions of an integro-differential equation on symmetric spaces, which generalizes the equation on $\mathbb R^d$
, related to the physical problem of the ground state energy of the Bose gas, studied by Carlen et al. (2020, Pure Appl. Anal., 2, 659–684). We also obtain a surprising result on the nature of the extremizers of the convolution inequalities.