Unveiling the fractal structure of Julia sets with Lagrangian descriptors

In this paper we explore by means of the method of Lagrangian descriptors the Julia sets arising from complex maps, and we analyze their underlying dynamics. In particular, we take a look at two classical examples: the quadratic mapping z(n+1) = z(n)(2) + c, and the maps generated by applying Newton's method to find the roots of complex polynomials. To achieve this goal, we provide an extension of this scalar diagnostic tool that is capable of revealing the phase space of open maps in the complex plane, allowing us to avoid potential issues of orbits escaping to infinity at an increasing rate. The simple idea is to compute the p-norm version of Lagrangian descriptors, not for the points on the complex plane, but for their projections onto the Riemann sphere in the extended complex plane. We demonstrate with several examples that this technique successfully reveals the rich and intricate dynamical features of Julia sets and their fractal structure. (C) 2020 Elsevier B.V. All rights reserved.