  We study the interaction and dynamics of two half-quantized vortices in
two-component Bose- Einstein condensates. Using the Pade approximation for the
vortex core profile, we calculate the intervortex potential, whose asymptotic
form for a large distance has been derived by Eto et al. [Phys. Rev. A, 83,
063603 (2011)]. Through numerical simulations of the two-dimensional
Gross-Pitaevskii equations, we reveal different kinds of dynamical trajectories
of the vortices depending on the combinations of signs of circulations and the
intercomponent density coupling. Under the adiabatic limit, we derive the
equations of motion for the vortex coordinates, in which the motion is caused
by the balance between Magnus force and the intervortex forces. The initial
velocity of the vortex motion can be explained quantitatively by this point
vortex approximation, but under- standing the long-time behavior of the
dynamics needs more consideration beyond our model.
