Call your game Game $1$. Game $2$ is almost the same, except that we continue for the full $13$ rounds, and we win Game $2$ if we get one or more match. The probability of winning Game $2$ is the same as the probability of winning Game $1$. So to answer your question, it is enough to find the probability of winning Game $2$.
The probability of a match at position $i$ is $\frac{4}{52}$. Now sum over all $i$ from $1$ to $13$. We get $\binom{13}{1}\frac{4}{52}$. This is $1$, and is the expected (mean) number of matches in Game $2$. But of course it is not the probability of winning Game $2$.
The problem is that we have counted twice every situation in which we have a match at $i$ and also a match at $j$. For clarity of thought, though in fact it doesn't matter, assume that $i\lt j$. The probability of a match at $i$ is $\frac{4}{52}$. Given that we have a match at $i$, the probability of a match at $j$ is $\frac{4}{51}$. So the probability of a match at $i$ and $j$ is $\frac{4}{52}\cdot \frac{4}{51}$. There are $\binom{13}{2}$ ways of picking $i$ and $j$. So from $\binom{13}{1} \frac{4}{52}$ we subtract $\binom{13}{2}\frac{4}{52}\cdot\frac{4}{51}$.
But we have subtracted too much! For we have subtracted one too many times all situations in which we had three matches. Let $i\lt j\lt k$. The probability of matches at $i$, $j$, and $k$ is $\frac{4}{52}\cdot\frac{4}{51}\cdot \frac{4}{50}$. There are $\binom{13}{3}$ ways of choosing the triple $(i,j,k)$. So we will add back $\binom{13}{3} \frac{4}{52}\cdot\frac{4}{51}\cdot \frac{4}{50}$.
But we have added back too much, for we have added back too many times the situations in which we have $4$ matches. So we must subtract $\binom{13}{4}\frac{4}{52}\cdot\frac{4}{51}\cdot \frac{4}{50}\cdot\frac{4}{49}$.
Continue. We are using the Method of Inclusion/Exclusion.
For practical work, we could probably stop. Already the term $\binom{13}{4}\frac{4}{52}\cdot\frac{4}{51}\cdot \frac{4}{50}\cdot\frac{4}{49}$ is fairly small. But it is not hard to continue to the end and get an exact answer.
Remark: Your simulation gave an answer quite close to the truth. The terms I mentioned in the post give roughly $0.639$, and the next term in the Inclusion/Exclusion is an "add" term.