The dynamics of the annihilation of a vortex-antivortex pair is investigated. The pair is activated magnetically during the run of a simulated hysteresis loop on a square mesoscopic superconducting cylinder with an antidot inserted at its center. We study the nucleation of vortices and antivortices by first increasing the magnetic field, applied parallel to the axis of the sample, from zero until the first vortex is created. A further increase in the field pulls the vortex in, until it reaches the antidot. As the polarity of the field is reversed, an antivortex enters the scene, travels toward the center of the sample, and eventually the pair is annihilated. Depending on the sample size, its temperature, and Ginzburg-Landau parameter, the vortex-antivortex encounter takes place at the antidot or at the superconducting sea around it. The position and velocity of the vortex and antivortex singularities were evaluated as a function of time. The current density, magnetization, and orderparameter topology were also calculated. Achieving a deep understanding of the nucleation and propagation of vortices in real superconductors is a truly complex task, since these entities interact with almost everything: first, with the surface of the specimen, to surpass it; upon entrance, with other vortices that might have already penetrated, and also with defects, which might attract them and even act as pinning centers. Additional difficulties to emulate the problem arise from the fact that vortices generate heat while propagating, what can be harmful to the robustness of the superconducting properties, if not catastrophic, as is the case of vortex avalanches observed in some superconducting films. [1][2][3][4][5][6] It is quite common, however, that the existence of pinning potentials represent a beneficial feature, since vortices can thus be prevented from undergoing dissipative motion. An interesting approach to the problem, which enables one to address most specificities without excessive complexity, is to work in the small universe of mesoscopic samples. In such an ambient, one can accommodate the essential ingredients: relatively important surface-tovolume ratio, only a few vortices on scene, and a number of defects-the so-called antidots-usually arranged in a regular pattern. Furthermore, one can study the interaction of an individual vortex-antivortex ͑V-AV͒ pair and, eventually, witness their mutual annihilation.Recently, there have been many studies about V-AV configurations in mesoscopic superconductors ͑see for instance Refs. 7-11͒. The authors of these references have found that vortices and antivortices may coexist in equilibrium in configurations which look like a V-AV molecule. A somewhat common approach is to assume an a priori configuration and minimize the free energy in terms of some relevant parameter for which the V-AV molecule is a stable configuration. Here, we will focus in a rather different approach concerning more with the dynamics of a V-AV encounter. The aim of the present work is to elucidate the de...