--- title: Globally Coupled Logistic Maps url: https://www.emergentmind.com/topics/globally-coupled-logistic-maps-gcm type: topic --- # Globally Coupled Logistic Maps Globally coupled logistic maps (GCM), often also discussed under the broader label of globally coupled map lattices when the local map is logistic or quadratic-logistic, are discrete-time dynamical systems in which each degree of freedom is updated by a nonlinear logistic map together with a mean-field interaction involving all units. In canonical form one writes \[ x_{n+1}(i)=(1-\varepsilon)f\bigl(x_n(i)\bigr)+\frac{\varepsilon}{N}\sum_{j=1}^{N}f\bigl(x_n(j)\bigr), \] with either the quadratic logistic form \(f_a(x)=1-a x^2\) or the standard logistic form \(f(x)=\alpha x(1-x)\); these parameterizations are topologically conjugate in the unimodal quadratic class and underpin much of the modern theory of synchronization, clustering, collective chaos, and chaotic itinerancy in high-dimensional nonlinear dynamics [2507.22643] [1112.4744]. ## 1. Canonical formulation and state-space organization A standard Kaneko-type GCM uses the quadratic local map \[ f_a(x)=1-a x^2, \qquad x_{n+1}(i)=(1-\varepsilon)f_a(x_n(i))+\frac{\varepsilon}{N}\sum_{j=1}^{N}f_a(x_n(j)), \] with \(a\) controlling local nonlinearity and \(\varepsilon\) the global coupling strength [2507.22643] [1112.4744]. A widely used alternative writes the local dynamics as \[ f(x)=\alpha x(1-x), \] so that the globally coupled logistic system evolves on \([0,1]^N\) under \[ x_{n+1}(i)=(1-\varepsilon)f(x_n(i))+\frac{\varepsilon}{N}\sum_{j=1}^{N}f(x_n(j)), \] or, in a closely related variant, \[ x_{n+1}(i)=(1-\varepsilon)f(x_n(i))+\frac{\varepsilon}{N-1}\sum_{j\neq i}f(x_n(j)), \] which excludes the self-term from the mean field and replaces \(N\) by \(N-1\) in the normalization [2510.00538] [2508.21721]. The distinction is structurally minor in large systems but relevant when comparing precise definitions. The mean-field structure produces a large family of invariant subspaces corresponding to cluster partitions. In the \([0,1]^N\) formulation, a permutation \(\sigma\in S_N\) defines an invariant set \[ H_\sigma=\{\mathbf{x}\in[0,1]^N\mid x(i)=x(\sigma(i)),\ i=1,\dots,N\}, \] so cluster states can be treated as symmetry-restricted dynamics on lower-dimensional manifolds [2510.00538]. A local “effective dimension” \[ \mathrm{ED}(x,\delta) := \inf_{\sigma\in S_N}\left\{\dim H_\sigma \mid H_\sigma\cap B_\delta(x)\neq\emptyset\right\} \] then provides a coarse-grained count of the number of clusters near a given state [2510.00538]. In this sense, GCM phase space is naturally stratified by synchronization patterns rather than by geometry in the usual lattice sense. The two-map case already contains the essential mean-field symmetry. For \(N=2\), \[ x_{n+1}=(1-\epsilon)f(x_n)+\epsilon f(y_n),\qquad y_{n+1}=\epsilon f(x_n)+(1-\epsilon)f(y_n), \] which is exactly the globally coupled form specialized to two units and has been used as an analytically tractable prototype for symmetry breaking, invariant manifolds, and bifurcation analysis [2005.08698] [1905.12977]. ## 2. Collective regimes, clustering, and synchronization Classical GCM phenomenology is organized into a small number of collective regimes. In the Kaneko classification, one encounters a **coherent phase**, an **ordered phase**, a **partially ordered phase**, and a **desynchronized phase** [2008.08204]. In the quadratic-map formulation these appear as **coherent chaos**, **two-clustered periodic states**, **periodic-cluster regimes**, and **turbulent states** as the coupling is decreased at fixed local nonlinearity [1112.4744]. The large-\(\varepsilon\) end is characterized by full synchronization, so the entire ensemble collapses onto the dynamics of a single map. At weaker coupling, cluster states emerge, and at still smaller coupling the system enters turbulent many-cluster dynamics. A central structure is the two-clustered periodic regime. There the population splits into two synchronized groups, and the cluster-population ratio \[ \vartheta=\frac{N_+}{N_++N_-} \] acts as an internal bifurcation parameter. In GCML phase diagrams, the balanced case \(\vartheta=0.5\) yields a maximally symmetric two-cluster attractor, while departures from balance produce a successive period-doubling cascade of the cluster orbit, the so-called \(\vartheta\)-bifurcation [1112.4744]. This makes cluster population itself a dynamical control parameter, not merely a state descriptor. Closely related globally coupled logistic-type systems confirm that synchronization can also be induced by size. In a multiplicative mean-field model for symbiotic logistic-type maps, the case \(\lambda=1.07\), \(\epsilon=0.01\) changes qualitatively when \(N\) is increased from \(100\) to \(1000\): the mean field develops a clear period-2 oscillation and the single-unit dynamics becomes entrained to it, an explicit example of size-induced synchronization [1908.08000]. This suggests that finite-population effects in GCMs are not limited to fluctuation amplitudes; they can reorganize the attractor itself. ## 3. Chaotic itinerancy and attractor-ruins One of the most distinctive dynamical signatures of GCMs is chaotic itinerancy: trajectories wander among quasi-stable collective states, stay there for irregular residence times, and then depart through chaotic transition episodes. In an entropy-and-clustering analysis of the quadratic GCM \[ x_{n+1}(i)=(1-\varepsilon)f_a(x_n(i))+\frac{\varepsilon}{N}\sum_{j=1}^N f_a(x_n(j)), \qquad f_a(x)=1-a x^2, \] local Shannon entropy and local permutation entropy were used to identify parameter regions likely to support itinerancy, and density-based clustering then separated dense “attractor-ruin” regions from sparse transition states [2507.22643]. For \(N=3\), a broad band in the \((a,\varepsilon)\)-plane running approximately from \((1.6,0)\) to \((2,0.2)\) was identified as the region where chaotic itinerancy is most pronounced, matching Kaneko’s intermittent or partially ordered phase [2507.22643]. At a representative point \(N=3\), \(a=2\), \(\varepsilon=0.2574\), DBSCAN identifies 12 clusters that organize into three 4-cycles, and these are naturally grouped into three attractor ruins. Time series segments assigned to those ruins are ordered and low-dimensional, whereas unassigned points correspond to chaotic transition states [2507.22643]. Statistical tests on the sequence of visited ruins and on the residence times were then used to show that the switching process is irregular but stationary rather than simply periodic [2507.22643]. A complementary description uses cluster statistics and optimal transport. In the standard logistic-map formulation, \[ x_{n+1}(i)=(1-\varepsilon)f(x_n(i))+\frac{\varepsilon}{N}\sum_{j=1}^{N}f(x_n(j)), \qquad f(x)=\alpha x(1-x), \] an attractor-ruin is treated as the remnant of a destabilized cluster attractor and is associated with a Milnor attractor [2510.00538]. Cluster configurations are compressed into a probability vector over cluster sizes, and an optimal-transport distance between successive cluster distributions detects when the system is stationary in cluster space. The entropy of the effective-dimension distribution at zero-transport times then measures the “strength” of attractor-ruins. This strength is maximal in the partially ordered phase, especially in the regime labeled Partially Ordered Phase II, where chaotic itinerancy is strongest [2510.00538]. These analyses correct a common oversimplification of chaotic itinerancy as a merely visual phenomenon. In GCMs it can be formalized through state-space clustering, entropy variance, residence-time statistics, and transport-based cluster diagnostics [2507.22643] [2510.00538]. ## 4. Rigorous results, invariant structures, and finite-size scaling Fully rigorous results for standard logistic GCMs remain limited, but two complementary lines of work are established. First, the two-site symmetrically coupled logistic map admits a detailed invariant-set analysis. For \[ F_{\mu,\epsilon}(x,y)=\bigl((1-\epsilon)f_\mu(x)+\epsilon f_\mu(y),\ (1-\epsilon)f_\mu(y)+\epsilon f_\mu(x)\bigr), \qquad f_\mu(t)=\mu t(1-t), \] the diagonal \(\Delta=\{(x,x)\}\) is invariant and carries the usual logistic dynamics [1905.12977]. In explicit parameter ranges there exists a positively invariant Jordan curve \(\Gamma\) containing \(F^{-1}(O)\) that bounds the immediate basin of infinity and organizes the bounded set. For small coupling and suitable \((\mu,\epsilon)\), every bounded orbit synchronizes to the diagonal; in other regimes off-diagonal fixed points, Cantor sets, annular basin components, and numerically observed fat attractors appear [1905.12977]. This two-dimensional prototype already exhibits the synchronization manifold, symmetry breaking, and fractal basin geometry that recur in larger GCMs. Second, rigorous continuum-limit results are known for globally coupled expanding circle maps. In that setting the state is a probability density, the weak-coupling regime admits a unique absolutely continuous invariant distribution that attracts all admissible initial densities exponentially fast in the \(BV\) norm, and the invariant density depends Lipschitz-continuously on \(\varepsilon\); for sufficiently strong coupling, a broad class of initial measures converges in Wasserstein distance to a moving Dirac mass, interpreted as synchronization in a chaotic state [1711.05461]. These theorems apply to uniformly expanding circle maps rather than logistic maps, but they provide the most complete rigorous mean-field picture currently available within the coupled-map framework. At the finite-size level, the largest Lyapunov exponent in turbulent GCMs is not universal. For globally coupled maps with positive multipliers, the asymptotic law \[ \lambda_\infty-\lambda(N)\simeq \frac{c}{N^\gamma}, \qquad 0<\gamma\le 1, \] replaces the previously proposed universal logarithmic correction, and the exponent \(\gamma\) depends on the multiplier statistics [2110.01949]. For strongly dissimilar multipliers, the convergence becomes even slower and no universal law is supported numerically [2110.01949]. Since standard logistic maps have derivatives of both signs, the positive-multiplier analysis does not directly resolve the classical logistic case, but it establishes that universal finite-size scaling of \(\lambda(N)\) cannot be assumed even in paradigmatic GCMs [2110.01949]. ## 5. Delays, adaptive couplings, and low-dimensional bifurcation theory Several important generalizations preserve the GCM mean-field spirit while changing the effective collective dynamics. A delayed, distance-dependent logistic-map lattice introduces a coupling range parameter \(\alpha\) and a delay parameter \(\beta\), with the global-coupling limit at \(\alpha=0\). In the completely synchronized state, \[ x_{t+1}^\star=(1-\varepsilon)f(x_t^\star)+\varepsilon f\bigl(\beta x_{t-1}^\star+(1-\beta)x_t^\star\bigr), \] so the synchronized manifold reduces to a two-dimensional non-delayed map independent of \(\alpha\) [1802.10539]. Moderate delay regularizes chaotic orbits and makes synchronization possible even for short-range coupling, while stronger delay introduces genuinely two-dimensional phenomena such as Neimark–Sacker bifurcations and quasiperiodic synchronized states [1802.10539]. The globally coupled case is therefore not just a stronger version of the undelayed model; delay changes the internal dimensionality of collective motion. At the smallest nontrivial scale, strong mean-field coupling in two globally coupled logistic maps produces explicitly computable local bifurcations. The asymmetric fixed points \((x^*,y^*)\) of the symmetric two-map system undergo a reverse flip bifurcation and a Neimark–Sacker bifurcation as the coupling \(\epsilon\) is varied [2005.08698]. In the reverse flip, unstable period-2 orbits collide with the symmetry-broken fixed point and a stable 1-cycle is recovered as coupling increases; on the Neimark–Sacker curve, the asymmetric fixed point loses stability to an invariant closed curve, giving quasiperiodic motion and then more complicated attractors [2005.08698]. These results show that coupling strength can act as a primary bifurcation parameter, not merely as a synchronizing perturbation. Global coupling can also be generalized from a static complete graph to adaptive mean-field networks. In a delayed Hebbian extension of coupled logistic maps, \[ x_i(n+1)=f\!\left((1-c)x_i(n)+c\sum_j w_{ij}(n)x_j(n)\right), \] uniform weights \(w_{ij}=1/(N-1)\) recover a standard GCM, while adaptive weights generate ten distinct network states, including pacemaker networks, paired layers, loop networks, hidden modular states, and hidden randomly connected states [2008.08204]. The familiar coherent, ordered, partially ordered, and desynchronized dynamical phases remain visible, but the coupling architecture becomes part of the emergent state [2008.08204]. A different logistic-type extension uses a multiplicative global mean field that modulates the effective growth parameter rather than adding an additive mean-field term, again producing synchronization and size-induced collective order [1908.08000]. ## 6. Experimental realizations and empirical hysteresis Electronic realizations have made GCM-like dynamics experimentally accessible with high precision. A modular analog implementation of the logistic map based on an AD633 analog multiplier plus a sample-and-hold block realizes \[ x_{n+1}=r x_n(1-x_n) \] with parameter control at \(\sim 10^{-3}\) resolution and reproduces the single-map bifurcation diagram and Lyapunov exponent with strong agreement between experiment and simulation [1512.08947]. A coupling block then implements \[ x_{n+1}^{(i)}=(1-\epsilon)f(r_i,x_n^{(i)})+\epsilon\sum_j \frac{A_{ij}}{d_i}f(r_j,x_n^{(j)}), \] so the same architecture extends directly to arbitrary networks and, in the complete-graph case, to finite GCMs [1512.08947]. A six-map experimental system with configurable diffusive coupling realized 52 distinct network topologies, including the complete graph, and operated with a signal-to-noise ratio of approximately \(10^6\) [2006.11378]. The isolated circuits reproduced Feigenbaum scaling with less than 8% deviation in \(\delta\), while the coupled system showed robust multistable regions and hysteresis even in the presence of measured parameter mismatch and unavoidable electronic noise [2006.11378]. In complete and ring-like configurations, hysteresis was observed under slow sweeps of the coupling strength \(\varepsilon\), and across all 52 configurations a strong-coupling hysteresis around the period-2 to fixed-point transition occurred near \(\varepsilon_c\approx 0.8\) [2006.11378]. The experiment therefore confirms that multistability in coupled logistic maps is not a fine-tuned numerical artifact; it survives realistic imperfections. Experimental work also sharpens the distinction between synchronization and identicality. In the six-map setup, the effective parameters \(r_i\) differed slightly across circuits, so the exact synchronization manifold was no longer strictly invariant; nonetheless, the average pairwise variance \[ \overline{\sigma^2}=\frac{2}{N(N-1)}\sum_{i