A sum involving derivatives of Vandermonde

Your expression is a polynomial $V^k(x_1,\ldots,x_n)$ that is still skew-symmetric. Therefore it is divisible by $V$. In addition, $V^k$ has the same degree as $V$. Thus the quotient $V^k/V$ is a constant. Hence the answer to your question is Yes.

By looking at the coefficient of the monomial $x_1^{n-1}x_2^{n-2}\cdots x_{n-1}$, one finds the constant $$c=\sum_{r=k}^{n-1}\frac{r!}{(r-k)!}.$$


Here is a direct way to obtain Denis Serre's formula: Just note that $x_i^k\frac{\partial^k}{\partial x_i^k}$ multiplies a monomial in the determinant by $\frac{r!}{(r-k)!}$ where $r$ is the power of $x_i$ in that monomial, provided that $r\geq k$. Otherwise, it acts by 0. That's all you need, since in each monomial, all powers between 1 and $n-1$ occur exactly once.


Perhaps a little Mathematica program will help us form a conjecture. For $k \geq n$ the answer is $0$, so we list your $c$ for $k < n$, as follows:

In[1]:= V[x_] := Product[x[[i]] - x[[j]], {i,1,Length[x]},{j,1,i-1}]
In[2]:= V[x/@Range[3]]
Out[2]= (-x[1]+x[2]) (-x[1]+x[3]) (-x[2]+x[3])
In[3]:= s[k_,x_] := Sum[x[[i]]^k*D[V[x],{x[[i]],k}],{i,1,Length[x]}]
In[4]:= Table[s[k,x/@Range[n]]/V[x/@Range[n]], {n,1,6},{k,0,n-1}]//Simplify//TableForm
Out[4]//TableForm=
1
2       1                               
3       3       2                       
4       6       8       6               
5       10      20      30      24      
6       15      40      90      144     120

Each row is showing a fixed $n$, and each column a fixed $k$. The first column is easy. The second one is triangular number $n(n+1)/2$, OEIS tells me the third one is $2 {n \choose 3}$ and the fourth one $n(n+1)(n+2)(n+3)/4$. After a bit of experimentation:

In[19]:= c[n_,k_] := Product[(n-k+i),{i,0,k}]/(k+1)
In[20]:= Table[c[n,k],{n,1,6},{k,0,n-1}]//TableForm
Out[20]//TableForm=
1
2       1                               
3       3       2                       
4       6       8       6               
5       10      20      30      24      
6       15      40      90      144     120

There we have it, your constant seems to be $n (n-1) \cdots (n-k)/(1+k)$. A true combinatorist should now be able to prove that this is really so. I already believe it.