# American Institute of Mathematical Sciences

## How to detect Wada basins

 1 Nonlinear Dynamics, Chaos and Complex Systems Group, Departamento de Física, Universidad Rey Juan Carlos, Tulipán s/n, 28933 Móstoles, Madrid, Spain 2 Department of Physics, Harvard University, Cambridge, Massachusetts 02138, USA 3 Department of Applied Informatics, Kaunas University of Technology, Studentu 50-415, Kaunas LT-51368, Lithuania

* Corresponding author: Alexandre Wagemakers

Received  April 2020 Revised  October 2020 Published  November 2020

We present a review of the different techniques available to study a special kind of fractal basins of attraction known as Wada basins, which have the intriguing property of having a single boundary separating three or more basins. We expose several approaches to identify this topological property that rely on different, but not exclusive, definitions of the Wada property.

Citation: Alexandre Wagemakers, Alvar Daza, Miguel A. F. Sanjuán. How to detect Wada basins. Discrete & Continuous Dynamical Systems - B, doi: 10.3934/dcdsb.2020330
##### References:

show all references

##### References:
Sketch of the condition 1 of the Nusse-Yorke method. The small disks represent areas of the basins $B_1$, $B_2$ and $B_3$. The unstable manifold of the unstable periodic orbit $P$ intersects the three basins. The preimages of the disks are stretched exponentially and asymptotically approach the stable manifold, which ultimately is the Wada boundary
Wada detection with the Nusse-Yorke method (a) Basins of attraction of the forced damped pendulum $\ddot{x}+0.2\dot{x}+\sin x = 1.66\cos t$, with the unstable manifold of a period three orbit (crosses on the basins). The unstable manifold intersects the three basins. (b) Basins of attraction of the forced damped pendulum $\ddot{x}+0.2\dot{x}+\sin x = 1.71\cos t$. There are four basins and we have found an accessible periodic orbit whose unstable manifold crosses only two basins. There is also a period-three periodic orbit similar to the case in (a). This basin is partially Wada
Sketch of the grid method. We set up a grid of boxes $box_j$ covering the whole disk. The center point of each box defines its color. In the first step, we see that $box_1$ belongs to the interior because its surrounding 8 boxes have the same color. On the other hand, $box_2$ and $box_3$ are in the boundary of two attractors, i.e., they are adjacent to boxes whose color is different. In the next step the algorithm classifies $box_2$ still in $G_2$ (boundary of two), while $box_3$ is now classified in $G_3$ (boundary of three). Ideally the process would keep on forever redefining the sets $G_1, G_2$ and $G_3$ at each step, though in practice we can impose some stopping condition. This plot constitutes an example of partially Wada basins
Wada detection with the grid method. (a) Basin of attraction of the forced damped pendulum $\ddot{x}+0.2\dot{x}+\sin x = 1.66\cos t$, (b) All $1000 \times 1000$ boxes are labeled either in the interior (white) or in the boundary of the three basins (black). (c) Histogram showing the number of points $N$ that take $q$ steps to be classified as boundary of three basins. (d) After reaching a maximum, there is an exponential decay of the computational effort related to the fractal structure of the basins. The log-plot reflects this tendency
Forced damped pendulum with eight basins. (a) The following damped forced pendulum $\ddot{x}+0.2\dot{x}+\sin x = 1.73\cos t$ shows eight basins of attraction mixed intricately. (b) Some boxes are classified to be in the boundary of eight basins (black dots), but not all of them (red dots), which is a clear example of a partially Wada basin. (c) The computational effort presents the usual shape for the Wada boundary, but the points which are not Wada keep refining until the algorithm meet the stop criterion at $q = 15$ (the red bar at rightmost represent the number of boxes not classified as Wada at this stage.). The grid method works best in systems with the Wada property. (d) Evolution of the proportion of boxes in the Wada boundary ($W_8$ in black) and proportion of boxes in a boundary which is not Wada ($W_{2 - 7}$) as a function of the $q$-step. The convergence of $W_8$ is used to determine the stopping rule
Graphical description of the merging of basins. In (a) upper left corner we have the original basin of the forced damped pendulum described by $\ddot{x}+0.2\dot{x}+\sin x = 1.66\cos t$. The other three panels are the modified basins with two merged basins. The colors above indicate which of the original basins have been merged. In (b) The case of the damped pendulum defined by $\ddot{x}+0.2\dot{x}+\sin x = 1.71\cos t$ is shown, which possesses four attractors. We have displayed only three of the four possible combinations of merging. However, these examples are enough to show that the boundaries are not identical
Interpretation of the Hausdorff distance. The figure represents two superimposed slim boundaries computed from two different merged basins. One of the boundaries is plotted with red pixels and the other one with green pixels. While it appears that most of the boundaries overlap, some parts of the red boundary do not coincide with the green boundary. The largest distance between the two boundaries is represented by a red circle of radius $max_d$ that corresponds to the Hausdorff distance between the two sets of points
Sketch of the saddle-straddle algorithm. Initially, two points are selected in such a way that each one lies on a different basin. Then, a bisection method is applied to reduce the distance between the two points to a desired accuracy. After that, the resulting points are iterated and the segment expands, so that the process must start over again. As a result, we obtain a set of arbitrarily small segments straddling the saddle
Computations of saddles with the saddle-straddle algorithm. (a) The picture represents the chaotic saddle embedded in the only boundary of the forced damped pendulum with equation $\ddot{x}+0.2\dot{x}+\sin x = 1.66\cos t$. (b) We have represented the computation of the saddle associated to the boundary between basins $B_2$ and $M_2$ of the forced damped pendulum with equation: $\ddot{x}+0.2\dot{x}+\sin x = 1.71\cos t$. In (c) we have the saddle corresponding to the boundary between basins $B_1$ and $M_1$. (d) shows the chaotic saddle of the boundary in the Hénon-Heiles Hamiltonian for the energy $E = 0.25$
Results of the computation of the Wada merging method for different systems with fractal basin boudaries. Some of these examples show a fractal basin according to the merging method. All the basins have been computed with a finite resolution of $1000 \times 1000$
 Dynamical system $max_d$ $min_d$ $(max_d - min_d)/ min_d$ Wada? Forced pendulum $N_A =3$ 0.0365 0.0219 0.667 YES Forced pendulum $N_A =4$ 0.368 0.0439 7.3826 NO Forced pendulum $N_A =8$ 0.3976 0.0655 5.0702 NO Hénon-Heiles Hamiltonian $E_0 = 0.2$ 0.0206 0.0168 0.2262 YES Hénon-Heiles Hamiltonian $E_0 = 0.3$ 0.0240 0.0236 0.0169 YES Newton method $N_A=3$ 0.0300 0.0240 0.2499 YES Newton method $N_A=4$ 0.0402 0.0350 0.1485 YES Newton method $N_A=5$ 0.0902 0.0420 1.1476 YES Newton method $N_A=6$ 0.0780 0.0566 0.3780 YES
 Dynamical system $max_d$ $min_d$ $(max_d - min_d)/ min_d$ Wada? Forced pendulum $N_A =3$ 0.0365 0.0219 0.667 YES Forced pendulum $N_A =4$ 0.368 0.0439 7.3826 NO Forced pendulum $N_A =8$ 0.3976 0.0655 5.0702 NO Hénon-Heiles Hamiltonian $E_0 = 0.2$ 0.0206 0.0168 0.2262 YES Hénon-Heiles Hamiltonian $E_0 = 0.3$ 0.0240 0.0236 0.0169 YES Newton method $N_A=3$ 0.0300 0.0240 0.2499 YES Newton method $N_A=4$ 0.0402 0.0350 0.1485 YES Newton method $N_A=5$ 0.0902 0.0420 1.1476 YES Newton method $N_A=6$ 0.0780 0.0566 0.3780 YES
Comparison of the principal procedures to test if a basin of attraction has the Wada property. The time noted with $^*$ refers only to the computation time and does not take into account the previous study of the system
 Name Type of system Dim. Computation What we need time Nusse-Yorke method [21] ODEs Hamiltonians Maps 2D 1$^*$ It requires a detailed knowledge of the basin and the boundaries (accessible unstable periodic orbit embedded in the basin boundary). Grid method [10] Any dynamical system n-D 100 It requires the basins and the dynamical system to compute parts of the basin at a higher resolution. Merging method [8] Any dynamical system n-D 0.01 It needs to know the basins, but not the dynamical system. Saddle-straddle method [30] ODEs Hamiltonians Maps 2D 1 It needs to know the dynamical system, but not the basins.
 Name Type of system Dim. Computation What we need time Nusse-Yorke method [21] ODEs Hamiltonians Maps 2D 1$^*$ It requires a detailed knowledge of the basin and the boundaries (accessible unstable periodic orbit embedded in the basin boundary). Grid method [10] Any dynamical system n-D 100 It requires the basins and the dynamical system to compute parts of the basin at a higher resolution. Merging method [8] Any dynamical system n-D 0.01 It needs to know the basins, but not the dynamical system. Saddle-straddle method [30] ODEs Hamiltonians Maps 2D 1 It needs to know the dynamical system, but not the basins.
 [1] Predrag S. Stanimirović, Branislav Ivanov, Haifeng Ma, Dijana Mosić. A survey of gradient methods for solving nonlinear optimization. Electronic Research Archive, 2020, 28 (4) : 1573-1624. doi: 10.3934/era.2020115 [2] Xuefei He, Kun Wang, Liwei Xu. Efficient finite difference methods for the nonlinear Helmholtz equation in Kerr medium. Electronic Research Archive, 2020, 28 (4) : 1503-1528. doi: 10.3934/era.2020079 [3] Xin Guo, Lei Shi. Preface of the special issue on analysis in data science: Methods and applications. Mathematical Foundations of Computing, 2020, 3 (4) : i-ii. doi: 10.3934/mfc.2020026 [4] Zexuan Liu, Zhiyuan Sun, Jerry Zhijian Yang. A numerical study of superconvergence of the discontinuous Galerkin method by patch reconstruction. Electronic Research Archive, 2020, 28 (4) : 1487-1501. doi: 10.3934/era.2020078 [5] Vieri Benci, Marco Cococcioni. The algorithmic numbers in non-archimedean numerical computing environments. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020449 [6] Wenbin Li, Jianliang Qian. Simultaneously recovering both domain and varying density in inverse gravimetry by efficient level-set methods. Inverse Problems & Imaging, , () : -. doi: 10.3934/ipi.2020073 [7] Leilei Wei, Yinnian He. A fully discrete local discontinuous Galerkin method with the generalized numerical flux to solve the tempered fractional reaction-diffusion equation. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020319 [8] Marion Darbas, Jérémy Heleine, Stephanie Lohrengel. Numerical resolution by the quasi-reversibility method of a data completion problem for Maxwell's equations. Inverse Problems & Imaging, 2020, 14 (6) : 1107-1133. doi: 10.3934/ipi.2020056 [9] Christopher S. Goodrich, Benjamin Lyons, Mihaela T. Velcsov. Analytical and numerical monotonicity results for discrete fractional sequential differences with negative lower bound. Communications on Pure & Applied Analysis, , () : -. doi: 10.3934/cpaa.2020269 [10] Omid Nikan, Seyedeh Mahboubeh Molavi-Arabshai, Hossein Jafari. Numerical simulation of the nonlinear fractional regularized long-wave model arising in ion acoustic plasma waves. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020466

2019 Impact Factor: 1.27

## Tools

Article outline

Figures and Tables