We are trying to estimate the cardinality $K(n,p)$ of so-called Kuratowski monoid with $p$ positive and $n$ negative linearly ordered idempotent generators. In particular, we are interesting in the case $n=p$.
This notion arose in an intersection of general topology and abstract algebra. We calculated the exact value of
$$K(n,p)=\sum_{i=0}^n\sum_{j=0}^p \binom {i+j}i\binom{i+j}j.$$
We try to simplify this formula, but we stuck. The problem is that many different binomial sums are known, but I cannot remember among them sums of binomial “polynomials” of degree greater than 1, containing the summation index above.
If the summation diapason was a triangle $\{i,j\ge 0:i+j\le 2n\}$ instead of a rectangle $[0;n]\times [0;p]$, a sum should be simplified as follows:
$$\sum_{i,j=0}^{i+j\le 2n} \binom {i+j}i^2=\sum_{k=0}^{2n}\sum_{i,j=0}^{i+j=k}\binom ki^2=\sum_{k=0}^{2n}\binom {2k}k.$$
A rough asymptotic of $K(n,n)$ should be $\log K(n,n)\sim n\log n+O(\log n)$, which follows form the bounds:
$$\binom {2n}n^2\le K(n,n)= \sum_{i=0}^n\sum_{j=0}^p \binom {i+j}i\binom{i+j}j\le 2(n+1)^2\binom {2n}n^2.$$
Since we are topologists, it is complicated for us to obtain a highly estimated asymptotic, so we decided to ask specialists about that.
The argument you've written down shows that
$$K(n, n) \le \sum_{k=0}^{2n} {2k \choose k}.$$
Stirling's formula gives the asymptotic
$${2k \choose k} \sim \frac{4^k}{\sqrt{\pi k}}$$
so this sum behaves approximately like a geometric series, and we expect that
$$\sum_{k=0}^{2n} {2k \choose k} \sim \frac{4}{3} \left( \frac{4^{2n}}{ \sqrt{2 \pi n} } \right).$$
We also have the (easily improved) lower bound
$${2n \choose n}^2 \le K(n, n)$$
and Stirling's formula here gives
$${2n \choose n}^2 \sim \frac{4^{2n}}{\pi n}$$
so we get upper and lower bounds that are within a multiplicative factor of $O(\sqrt{n})$ of each other this way. You can tighten the lower bound using the refined entropy formula.