REVIEW 3 major objections 6 minor 14 references
Multifractal analysis of intermingled basins and blowout bifurcations in a parametetric family of skew product maps
T0 review · 3 major / 6 minor · reviewed 2026-08-08 · deepseek-v4-flash
Pith's one-line read For an open region of parameters, a two-parameter skew product family has two chaotic attractors whose basins intermingle, separated by a single invariant graph whose stability-index level sets have dimensions given by a Legendre transform.
desk verdict A solid explicit example of intermingled basins, but the multifractal spectrum is not proved as written; send to review with Section 5 flagged. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the skew product map $F_{a,b}(x,y)=(f(x),g_{a,b,x}(y))$ with base map $f(x)=2x\bmod 1$ and fiber maps $g_a(y)=y+ay(1-y)$ for $x\in[0,1/2)$ and $g_b(y)=y-by(1-y)$ for $x\in[1/2,1]$. The two constant graphs $\phi_0(x)=0$ and $\phi_1(x)=1$ are invariant, and the proof of Theorem 1 rests on the normal Lyapunov exponents of the SRB measures on them: for $\Gamma_0$, $\Lambda_{\mathrm{SRB}}(A_0)=\tfrac12\log(1+a)+\tfrac12\log(1-b)<0$, and for $\Gamma_1$, $\Lambda_{\mathrm{SRB}}(A_1)=\tfrac12\log(1-a)+\tfrac12\log(1+b)<0$. Riddling is obtained by applying a cited criterion that combines negative SRB normal exponent with a dense set of periodic points having positive normal exponent; the separating graph $\phi^*$ is identified with the common value of the lower boundary function $\ell_{a,b}(x)=\sup\{y:\lim_n g^n_{a,b,x}(y)=0\}$ and the upper boundary function $r_{a,b}(x)=\inf\{y:\lim_n g^n_{a,b,x}(y)=1\}$, whose equality is imported from a cited proposition. The stability-index and multifractal results are carried by the pressure functions $p_{i,\psi}(t)$ and by the implicitly defined curves $T(q)$, $S(q)$, whose Legendre transforms give the Hausdorff dimensions of the level sets.
What would settle it
Compute $\ell_{a,b}(x)$ and $r_{a,b}(x)$ numerically for a dense grid of $x$ and for fine parameter grids inside $\Gamma$; if for some set of positive $\nu_{\mathrm{ac}}$-measure the two functions differ, the separating graph $\phi^*$ of Theorem 1 fails and the dimension formulas collapse. A cheaper check is to measure the stability index at points just above and below the numerically located boundary for a fixed parameter in $\Gamma$ and compare the values with the constants predicted by Theorem 4.
Extended reading notes
Core claim
For $F_{a,b}$ defined by (2.2)--(2.3), if $(a,b)$ lies in $\Gamma=\{0<a<1/2,\,0<b<1/2,\,b>a/(1+a),\,b<a/(1-a)\}$, then $F_{a,b}$ has two chaotic essential Milnor attractors $A_0\subset I\times\{0\}$ and $A_1\subset I\times\{1\}$ whose basins are locally riddled and intermingled (Theorem 1). The basins are separated by a $\nu_{\mathrm{ac}}$-almost everywhere invariant graph $\phi^*$ with $0<\phi^*<1$. Crossing $b=a/(1+a)$ or $b=a/(1-a)$ turns one of the attractors into a chaotic saddle and produces a blowout bifurcation, one supercritical and one subcritical (Theorem 2 and Corollary 2.13). Under condition (H1), the stability index $\sigma_\nu(x,y)$ is constant $\nu_\psi$-almost everywhere on each side of $\phi^*$, with values written in terms of positive zeros $t^*_0,t^*_1$ of pressure functions $p_{i,\psi}(t)=P(\psi+t\log dg_{a,b}(\phi_i))$ (Theorems 3 and 4). The level sets $A_{i,\nu}(\sigma)$ have Hausdorff dimension $T(q)+q\sigma(q)$ or $S(q)+q\sigma(q)$, where $T$ and $S$ are defined implicitly by pressure equations, so the multifractal spectrum of the stability index is the Legendre transform of a thermodynamic function (Theorems 5 and 6).
Load-bearing premise
The load-bearing premise is that for parameters in $\Gamma$ the lower and upper boundary functions $\ell_{a,b}(x)$ and $r_{a,b}(x)$ coincide almost everywhere, so that one separating graph $\phi^*$ exists, and that condition (H1) holds: a Gibbs measure $\nu_\psi$ with negative normal Lyapunov exponents on both invariant lines; the paper verifies the equality only by importing a cited proposition, and verifies (H1) explicitly only for the SRB potential.
Editorial extensions
If this is right
- For every $(a,b)\in\Gamma$, almost every point of the square is eventually attracted to one of the two invariant lines $y=0$ or $y=1$, yet each basin has positive measure and is riddled with the other, so the asymptotic future of a point cannot be inferred from the futures of arbitrarily close neighbours.
- The parameter plane splits into three dynamical regimes, with $b=a/(1+a)$ and $b=a/(1-a)$ as the two blowout bifurcation curves; crossing one turns the corresponding attractor into a chaotic saddle while the other attractor retains a locally riddled basin.
- For $\nu_\psi$-almost every $x$, the stability index $\sigma_\nu(x,y)$ is constant on each side of the separating graph $\phi^*$, with explicit values $t^*_0\,(-\lambda_{\nu_\psi}(\phi_0))/\int\log|df|\,d\nu_\psi$ and $t^*_1\,(-\lambda_{\nu_\psi}(\phi_1))/\int\log|df|\,d\nu_\psi$.
- The Hausdorff dimension of the stability-index level sets equals $T(q)+q\sigma(q)$ for the $A_0$ side and $S(q)+q\sigma(q)$ for the $A_1$ side, so the full multifractal spectrum is described by the Legendre transform of a pressure function.
- In the SRB case $\psi=-\log|df|$, the stability-index theorems reduce to the previously known intermingled-basin results, confirming that the new pressure-based formulas contain the earlier theory as a special case.
Reading between the lines
- The same two-parameter mechanism should work for any piecewise expanding Markov base map with two branches whose fiber maps have negative Schwarzian derivative, since the paper's Remark 2.1 states the arguments extend and the only quantitative input is the normal Lyapunov exponents.
- If the stability index is measured numerically at points just above and below the numerically located boundary curve for a parameter inside $\Gamma$, the predicted constant values from Theorem 4 give a direct, quantitative test of the thermodynamic formalism.
- The same pressure-based Legendre-transform construction could be applied to other skew products with finitely many invariant graphs to produce dimension spectra for riddled basins, although the existence of a single separating graph would need to be established separately in each case.
- If condition (H1) fails for a potential other than $-\log|df|$, the dimension formulas may still hold for the SRB measure, but the claimed generality over Gibbs measures would require an additional argument.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies a two-parameter family of skew-product maps F_{a,b}(x,y)=(2x mod 1, g_{a,b,x}(y)) with two invariant lines. It identifies an open region Γ in the (a,b)-plane in which the two invariant chaotic attractors have intermingled basins separated by an invariant graph φ*, proves the complementary chaotic-saddle regime, and locates blowout bifurcations on the two boundary curves. The second half uses thermodynamic formalism to compute Loynes exponents and the stability index for Gibbs measures satisfying condition (H1), and claims multifractal spectra of the stability-index level sets as Legendre transforms of pressure functions.
Significance. Theorems 1 and 2 rest on explicit derivative computations and give a clean, checkable example of intermingled basins and blowout bifurcations. The paper is also honest that the existence of the separating graph is imported from Keller's work. No parameters are fitted: the stability index and the pressure functions are defined directly from the dynamics. If the multifractal dimension formulas were established, they would be a valuable concrete illustration of the Legendre-transform formalism for stability-index level sets. As it stands, however, the multifractal part contains a load-bearing gap, so the full significance claimed in the abstract is not yet realized.
major comments (3)
- [Section 5, Eq. (5.3)] The displayed identification of A_{i,νψ}(σ) with a Birkhoff-average set is asserted rather than proved, and as written it is also incorrect in sign. For νψ-a.e. x and y < φ*(x), Theorem 4(2) gives σνψ(x,y) > 0, whereas the quotient Q0(x) = lim t*0 S_n log dg(φ0) / (-S_n log|df|) converges to t*0 λνψ(φ0) / ∫ log|df| dνψ < 0. Thus A0,νψ(σ) would need to be identified with {Q0 = -σ}, not {Q0 = σ}. The same sign issue occurs for A1 because (2.17) defines A1,νψ(σ) through σνψ = -σ. Moreover, Theorem 4 only computes the stability index on a νψ-full-measure set via the sequence (4.14); for σ different from the ergodic value the level set is νψ-null, and no pointwise estimates on that null set are supplied. Therefore (5.3) is an additional load-bearing assertion, and Theorems 5(2) and 6(2) do not follow from the preceding estimates.
- [Section 5.1, proof of Theorem 6] In Case (1) the paper concludes from ∫ log dg(φ1) dνq < 0 that φ* is defined νq-a.e., citing [Kel17, Proposition 1.6]. In the earlier use of Proposition 1.6, namely Corollary 2.12, the condition is applied after verifying negative Lyapunov exponents for both bounding graphs φ0 and φ1 with respect to the relevant measure. For q ≠ 0 the measure νq differs from νψ, and the text does not verify λνq(φ0) < 0. Hence the use of Proposition 1.6 is missing a hypothesis, and the Moran-cover argument that relies on the νq-defined graph is unsupported. The same omission affects Theorem 5, whose proof is only said to be similar.
- [Section 5, proof of Theorems 5 and 6] The theorems are stated as dimension formulas for the stability-index level sets A0,νψ(σ) and A1,νψ(σ), but the proof only attempts to analyze Birkhoff-average level sets after replacing them by (5.3). Even after correcting the sign, the paper would need to prove that the stability-index condition 'for every y < φ*(x)' is equivalent to a pointwise Birkhoff-average condition. No such equivalence is shown: Theorem 4 controls σνψ only for νψ-a.e. x and all y in the respective half-neighborhood, and the exceptional null sets are not controlled. The final dimension claims therefore need a genuinely new argument or a restricted statement, e.g. for the Birkhoff-average sets defined in (5.3).
minor comments (6)
- [Section 4.2, Case (2)] The formula 'σ0ν(x,y) = t*1 · -log λ_{a,b,ν}(φ1) / ∫ log|df| dν' should read '-λ_{a,b,ν}(φ1)' instead of '-log λ_{a,b,ν}(φ1)'.
- [Section 4, Proposition 4.7] Hypothesis (i) is mis-stated: the integral should be written as ∫ e^{t log dg^n_{a,b,x}(φ0)} dν, and the claimed bound cannot follow from boundedness of dg alone; the text should specify the intended inequality and its domain in t.
- [Section 4, Lemma 4.5 and Corollary 4.6] There are minor typos: 'Lemm 12.1.4' appears twice, and 'P(psi)' should be 'P(ψ)'.
- [Theorem 4 statement] The chain of equalities in item (1) is confusing because σ0νψ is used with different signs in the statement and in the proof; the notation should be aligned so that σνψ = σ1νψ - σ0νψ is transparent in both places.
- [Section 2, Condition (H1)] Condition (H1) is verified only for the SRB potential ψ = -log|df|. The paper should state explicitly that Theorems 3-6 for general ψ are conditional on this unverified hypothesis, or provide an example of a non-SRB potential satisfying it.
- [Title and abstract] The word 'parametetric' in the title is a typo, and in the abstract 'these attractors either exhibits' should be 'these attractors either exhibit'.
Circularity Check
No significant circularity: the paper's derived quantities are computed from the dynamics and external results, not fitted or defined into the conclusions; the main gaps are unproved identifications, not circular reductions.
full rationale
The derivation chain is not circular in the load-bearing sense. The normal Lyapunov exponents, pressure functions p_{i,psi}(t), Loynes exponents t*_i, and stability indices are computed from the explicit skew-product dynamics (2.2)-(2.3) and thermodynamic formalism, not fitted to the quantities they are claimed to predict. The intermingled-basin and separating-graph facts are imported from Keller's independent work ([Kel17, Proposition 1.6 and 2.2]) rather than from the authors' own prior results; none of the authors' self-citations ([FKG18], [KG20], [Rab+22], [KS04], [KS07], [KS08], [JK11], [JKM21]) is invoked as the proof of a central theorem. Theorems 5 and 6 apply the standard Pesin-Weiss multifractal formalism to pressure-defined functions T and S, so the Legendre-transform dimension formulas are consequences of that formalism rather than restatements of the definitions. The skeptical concerns are genuine correctness gaps, not circularity: equation (5.3) identifies the stability-index level sets A_{i,nu_psi}(sigma) with Birkhoff-average sets only for nu_psi-a.e. points, and for non-ergodic values sigma the sets are nu_psi-null, so no pointwise estimate is supplied; Section 5.1 also verifies only one of the two normal Lyapunov exponents needed for [Kel17, Proposition 1.6]. These omissions mean Theorems 5 and 6 are not fully established, but they do not make the argument circular: the equality (5.3) is asserted, not built into the definition of the stability index, and condition (H1) is stated as a hypothesis rather than smuggled in as a conclusion. Accordingly the circularity score is low.
Assumptions & free parameters
assumptions (5)
- standard math Existence, uniqueness, analyticity and convexity of equilibrium states and pressure for Holder continuous potentials on expanding Markov maps (Lemma 4.1).
- domain assumption The family F_{a,b} with f = 2x mod 1 and fiber maps g_a(y) = y + a y(1-y), g_b(y) = y - b y(1-y), a,b in (0,1/2), has the stated invariant sets and SRB measure nu_ac equal to Lebesgue on each branch.
- ad hoc to paper Condition (H1): existence of an equilibrium state nu_psi with P(psi)=0 and lambda_{nu_psi,a,b}(phi_i) < 0 for i=0,1.
- standard math [Kel17, Proposition 1.6] and [Kel17, Proposition 2.2] guarantee that if the two bounding graphs have negative normal exponents and positive intermingling, then l_{a,b} = r_{a,b} nu-a.e. and the basins are intermingled.
- standard math Koebe-type distortion estimate (4.10) and the distortion bounds (4.1)-(4.2) hold for the family.
Cite this review
Pith. "Pith review of Multifractal analysis of intermingled basins and blowout bifurcations in a parametetric family of skew product maps." pith.science (2026). https://pith.science/paper/MIA2NGZ4
@misc{pith2026250204869,
author = {Pith},
title = {Pith review of: Multifractal analysis of intermingled basins and blowout bifurcations in a parametetric family of skew product maps},
year = {2026},
howpublished = {\url{https://pith.science/paper/MIA2NGZ4}},
note = {Machine review of arXiv:2502.04869}
}
read the original abstract
In this paper we study a two-parameter family of planar maps characterized by two distinct invariant subspaces. The model reveals the existence of two chaotic attractors within these subspaces. We identify parameter values at which these attractors either exhibits a locally riddled basin of attraction or transitions into a chaotic saddle. In particular, we demonstrate that, for an open region in the parameter plane, their basins are intermingled. It is shown that a fractal boundary curve separates the basins of attraction of these two chaotic attractors, providing a detailed characterization of the riddled basin structure. Additionally, we show that the model undergoes a blowout bifurcation. An estimation of the stability index is examined using thermodynamic formalism. We also perform a multifractal analysis of the level sets of the stability index.
Figures
Figures from the paper (2 more)
Reference graph
Works this paper leans on
-
[1]
[AAN98] P. Ashwin, P. J. Aston, and M. Nicol. On the unfolding of a blowout bifurcation.Physica D: Nonlinear Phenomena111.1-4 (1998), 81–95. [ABS94] P. Ashwin, J. Buescu, and I. Stewart. Bubbling of attractors and synchronisation of chaotic oscillators.Physics Letters A193.2 (1994), 126–139. [ABS96] P. Ashwin, J. Buescu, and I. Stewart. From attractor to ...
-
[10]
1142/S0219493721400074. [Kan94] I. Kan. Open sets of diffeomorphisms having two attractors, each with an everywhere dense basin.Bulletin of the American Mathematical Society31.1 (1994), 68–74. [Kel14] G. Keller. Stability index for chaotically driven concave maps. J. Lond. Math. Soc. (2) 89.2 (2014), 603–622.doi: 10.1112/jlms/jdt070. [Kel17] G. Keller. St...
-
[17]
doi: 10.1142/S021812742050248X. [Kim+18] H. Kim, S. H. Lee, J. Davidsen, and S.-W. Son. Multistability and variations in basin of attraction in power-grid systems.New Journal of Physics20.11 (2018), 113006. [KS04] M. Kesseböhmer and B. O. Stratmann. A multifractal formalism for growth rates and applications to geometrically finite Kleinian groups.Ergodic ...
-
[42]
Cambridge University Press, Cambridge, 1998, pp
London Mathematical Society Student Texts. Cambridge University Press, Cambridge, 1998, pp. x+178.doi: 10.1017/ CBO9781107359987. [KG20] S. Karimi and F. H. Ghane. Analysis of coexistence and extinction in a two-species competition model. Internat. J. Bifur. Chaos Appl. Sci. Engrg.30.16 (2020), 2050248,
work page 2020
-
[53]
Fields Inst. Commun. Amer. Math. Soc., Provi- dence, RI, 2008, pp. 1–21. [Bue12] J.Buescu. Exotic attractors: From Liapunov stability to riddled basins.Vol.153.Birkhäuser,
work page 2008
-
[79]
The stability index for dynamically defined Weierstrass functions
Graduate Texts in Mathematics. Springer-Verlag, New York-Berlin, 1982, pp. ix+250. [WW17] C. P. Walkden and T. Withers. The stability index for dynamically defined Weierstrass functions. arXiv preprint arXiv:1709.02451(2017). Department of Mathematics, Ferdowsi University of Mashhad, Iran Email address: htina@um.ac.ir Institute for Dynamical Systems, Facu...
work page Pith review arXiv 2017
-
[151]
Cambridge University Press, Cambridge, 2016, pp
Cambridge Studies in Advanced Mathematics. Cambridge University Press, Cambridge, 2016, pp. xvi+530. doi: 10.1017/CBO9781316422601. [Wal82] P. Walters. An introduction to ergodic theory. Vol
-
[279]
Pesin and H
[PW97] Y. Pesin and H. Weiss. A multifractal analysis of equilibrium measures for conformal expanding maps and Moran-like geometric constructions.Journal of Statistical Physics 86 (1997), 233–275. [RA16] U. A. M. Roslan and P. Ashwin. Local and global stability indices for a riddled basin attractor of a piecewise linear map.Dynamical Systems31.3 (2016), 3...
1997
Show all 14 references
-
[887]
Parry and M
32 REFERENCES [PP90] W. Parry and M. Pollicott. Zeta functions and the periodic orbit structure of hyperbolic dynamics. Astérisque 187.188 (1990), 1–268. [PS75] D. Plachky and J. Steinebach. A theorem about probabilities of large deviations with an application to queuing theor...
1990
-
[1239]
Jaerisch and M
[JK11] J. Jaerisch and M. Kesseböhmer. Regularity of multifractal spectra of conformal iterated function systems.Trans. Amer. Math. Soc.363.1 (2011), 313–330.doi: 10.1090/S0002- 9947-2010-05326-7. REFERENCES 31 [JKM21] J. Jaerisch, M. Kesseböhmer, and S. Munday. A multifractal...
2011 doi
-
[2012]
Cazelles
[Caz01] B. Cazelles. Dynamics with riddled basins of attraction in models of interacting popula- tions. Chaos, Solitons & Fractals12.2 (2001), 301–311. [Daz+16] A. Daza, A. Wagemakers, B. Georgeot, D. Guéry-Odelin, and M. A. Sanjuán. Basin entropy: a new tool to analyze uncert...
2001
-
[3528]
[Jäg03] T. H. Jäger. Quasiperiodically forced interval maps with negative Schwarzian derivative. Nonlinearity 16.4 (2003),
2003
-
[4134]
[Ott+94] E. Ott, J. Alexander, I. Kan, J. C. Sommerer, and J. A. Yorke. The transition to chaotic attractors with riddled basins.Physica D: Nonlinear Phenomena76.4 (1994), 384–410. [PA11] O. Podvigina and P. Ashwin. On local attraction properties and a stability index for hete...
1994
-
[5329]
[HCP94] J. F. Heagy, T. L. Carroll, and L. M. Pecora. Experimental and numerical evidence for riddled basins in coupled chaotic systems.Physical Review Letters73.26 (1994),
1994
Reviewed August 8, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.