I have a very indirect proof of the following property involving a parametrized integral. If $a,a\_1,\ldots,a\_n\in\mathbb R^n$ (here $n\ge2$), let me denote $V(a,a\_1,\ldots,a\_n)$ the volume of the simplex spanned by these $n+1$ points. Then
$$\sup\_{a\_1,\ldots,a\_n}\int\_{\mathbb R^n}\left[\frac{V(x,a\_1,\ldots,a\_n)^2}{\left(\prod\_1^n|x-a\_i|\right)^{n+1}}\right]^{\frac1{n-1}}dx<+\infty.$$
Remark that the integral, $I(a\_1,\ldots,a\_n)$, is homogeneous of degree zero in the argument $(a\_1,\ldots,a\_n)$.


Mind also that the naive bound
$$\int\_{\mathbb R^n}\frac{dx}{\prod\_1^n|x-a\_i|}$$
diverges ; and the restriction of the latter integral to a ball $B(0;R)$, while finite, tends to infinity as all $a\_j$'s converge to $0$.



> 
> Is there an elementary proof of the boundedness of $(a\_1,\ldots,a\_n)\mapsto I$, or is it notoriously difficult ? If elementary, do you have an explicit reference ?
> 
> 
> 


Remark that if $n=2$, then $I(a\_1,a\_2)$ is a constant, because it is invariant under isometries and homogeneous of degree $0$, thus $\equiv I(0,\vec e\_1)$.



Integral inequality: an elementary proof?
By dyadic decomposition, it suffices to obtain bounds on the quantities
$$
\int\_{|x-a\_i| \sim R\_i \forall i} \left[ \frac{V(x,a\_1,\dots,a\_n)^2}{(R\_1 \cdots R\_n)^{n+1}}\right]^{\frac{1}{n-1}}\ dx$$
which, when summed over dyadic powers of two $R\_1,\dots,R\_n > 0$, becomes bounded by $O(1)$, uniformly over $a\_1,\dots,a\_n$.


By permuting we may assume that $R\_1 \leq \dots \leq R\_n$. The most important geometric feature of the set $\{a\_1,\dots,a\_n\}$ turns out to be its diameter; by rescaling we can normalize this diameter to be $1$. The triangle inequality then lets us assume $R\_n \gtrsim 1$. In fact the triangle inequality leaves us with just two scenarios: the "[far field](https://en.wikipedia.org/wiki/Near_and_far_field)" scenario
$$ R\_1 \sim \dots \sim R\_n \gtrsim 1 \quad (1)$$
and the "[near field](https://en.wikipedia.org/wiki/Near_and_far_field)" scenario
$$ R\_1 \leq \dots \leq R\_n \sim 1 \quad (2).$$


One could refine the range of possibilities of the near field further by inspecting the finer metric geometry of the $a\_1,\dots,a\_n$ beyond just using the diameter normalization, but fortunately we do not have to do so for this problem.


As implicitly observed by the OP, the volume $V(x,a\_1,\dots,a\_n)$ is of size at most $O(R\_1 \cdots R\_n)$ by multiplying all the lengths emenating from $x$. On the other hand, since $a\_1,\dots,a\_n$ lie in a diameter one set and $x$ is at distance $O(R\_1)$ from $a\_1$, we also have the bound $O(R\_1)$ by the base-times-height formula. [Note that this is a significantly superior bound in the far field case, which was identified by the OP as the case where the previous estimates were poor.] Finally, the condition $|x-a\_1| \sim R\_1$ restricts $x$ to a set of measure $O(R\_1^n)$. Thus the above integral can be bounded by
$$ \lesssim \left[\frac{\min(R\_1 \cdots R\_n, R\_1)^2}{(R\_1 \cdots R\_n)^{n+1}}\right]^{\frac{1}{n-1}} R\_1^n.$$


In the far field case (1), this simplifies to $O(R\_n^{-2})$ which is summable to $O(1)$. [NOTE: in a previous version I instead obtained the estimate $O(\frac{1}{R\_1^{\frac{n-3}{n-1}}R\_n^{\frac{n+1}{n-1}}})$, which was easily controlled for $n \geq 3$ but (as pointed out in comments) required a different argument to sum for $n=2$ if one did not exploit the additional property $R\_1 \sim R\_2$ of the far field.] In the near field case (2), we bound the above by
$$ \left[\frac{(R\_1 \cdots R\_n)^2}{(R\_1 \cdots R\_n)^{n+1}}\right]^{\frac{1}{n-1}} R\_1^n = \prod\_{j=1}^n \frac{R\_1}{R\_j} \leq \frac{R\_1}{R\_n} \sim R\_1$$
which is summable to $O(1)$ since $R\_1 \lesssim 1$ (and there are only $O(\log^{n-2} (1/R\_1))$ choices of the remaining scales $R\_2,\dots,R\_n$ for any fixed $R\_1$).


EDIT: one can refactor the above proof to avoid dyadic decomposition as follows.


We again normalize $\mathrm{diam}(a\_1,\dots,a\_n)=1$, and by translation we also normalize $a\_1=0$. We now split into the far field $|x| \geq 2$ and the near field $|x| < 2$. In the far field we have $V(x,a\_1,\dots,a\_n) \lesssim |x|$ by the base-times-height formula and $|x-a\_i| \sim |x|$ by the triangle inequality, so the contribution of the far field is
$$ \lesssim \int\_{|x| \geq 2} \left[ \frac{|x|^2}{|x|^{n(n+1)}}\right]^{\frac{1}{n-1}}\ dx = \int\_{|x| \geq 2} \frac{1}{|x|^{n+2}}\ dx = O(1).$$
In the near field we use the [Hadamard bound](https://en.wikipedia.org/wiki/Hadamard%27s_inequality) $V(x,a\_1,\dots,a\_n) \leq \frac{1}{n!} \prod\_{i=1}^n |x-a\_i|$ as in the OP to bound things by
$$ \lesssim \int\_{|x| \leq 2} \frac{dx}{\prod\_{i=1}^n |x-a\_i|}.$$
By the triangle inequality, at least one of the $|x-a\_i|$ is comparable to one. Deleting this factor and then applying the [AM-GM inequality](https://en.wikipedia.org/wiki/Inequality_of_arithmetic_and_geometric_means), one can bound this by
$$ \lesssim \int\_{|x| \leq 2} \sum\_{i=1}^n \frac{dx}{|x-a\_i|^{n-1}} = O(1)$$
and the claim follows. (Note how it was important to delete one of the factors in the denominator before applying AM-GM in order to obtain an integrable exponent of $n-1$ instead of the divergent exponent $n$. One can also use Holder's inequality in place of AM-GM if desired.)


This calculation and its refactoring illustrates how dyadic decomposition can be a systematic tool to estimate (up to constants) various unsigned geometric integrals (with every geometric fact about the integrand and the domain of integration reducible to some inequality regarding the various scales and amplitudes involved, which can then be combined and optimized by elementary algebra), but once one gains enough understanding of the dyadic decomposition argument, one can often streamline the proof to erase traces of the dyadic decomposition used, though sometimes at the cost of providing motivation as to where the proof came from.


