Gravitational lensing of massive particles in the charged NUT spacetime

In astronomy, gravitational lensing of light leads to the formation of multiple images, arcs, Einstein rings, and, most important, the shadow of black holes. Analogously in the vicinity of a massive compact object massive particles, following timelike geodesics, are gravitationally lensed. So far gravitational lensing of massive particles was mainly investigated in the weak and strong field limits. In this paper we will, for the first time, investigate exact gravitational lensing of massive particles using the example of the charged Newman-Unti-Tamburino (NUT) metric (and its special cases) which contains three physical parameters, the mass parameter m, the electric charge e, and the gravitomagnetic charge n. We will first discuss and solve the equations of motion for unbound timelike geodesics using elementary and Jacobi's elliptic functions and Legendre's elliptic integrals. Then we will introduce an orthonormal tetrad to relate the z component of the angular momentum and the Carter constant to the energy E of the particles along the timelike geodesics and latitude-longitude coordinates on the celestial sphere of a stationary observer in the domain of outer communication. We will use these relations to derive the angular radius of the particle shadow of the black hole, to formulate an exact lens equation, and to derive the travel time of the particles in terms of the time coordinate and the proper time. Finally, we will discuss the impact of the physical parameters and the energy of the particles on observable lensing features. We will also comment on how we can use these features alone and in a multimessenger context together with the corresponding features for light rays to determine if an astrophysical black hole can be described by the charged NUT metric or one of its special cases.