Fractional Logistic Map Dynamics and Bifurcations

ISEF Category: Mathematics

Subcategory: Analysis  ·  Difficulty: Advanced  ·  Setup: University Lab  ·  Time: Full Year

The Hook

A tiny change in an equation can flip a system from calm growth to chaos. That is the weird power of nonlinear maps. With this project, you can study how a simple rule for population-like growth changes when you adjust the exponents. You will turn math into pictures, patterns, and a chaos detector called a Lyapunov exponent.

What Is It?

The logistic map is a classic equation that updates a number over and over again. Think of it like a rule for a tiny economy or a population. Each new value depends on the last one, but growth gets squeezed as the value gets closer to one. That tension can create steady behavior, repeating cycles, or chaos.

This project changes the usual logistic map by adding two exponents, α and β. Those exponents reshape how the system grows and how it slows down. You can think of them like two knobs that bend the curve before the next step happens. Small changes in those knobs can move the system from smooth behavior to wild swings.

A bifurcation diagram shows where those behavior shifts happen. A Lyapunov exponent measures whether nearby starting points move apart or stay close. If the exponent is positive, the system tends to amplify tiny differences. If it is negative, the system tends to settle down.

Why This Is a Good Topic

This is a strong science fair topic because you can test it with clear inputs and clear outputs. You change α, β, or r, then record what the map does after many iterations. That gives you a clean path to graphs, pattern finding, and mathematical proof. The topic also connects to real systems that show feedback, like populations, finance, and some biological processes. You can learn how to model chaos, build computational experiments, and compare numerical results with theory.

Research Questions

• How does changing α affect the onset of period-doubling in the fractional-order logistic map?
• How does changing β alter the shape of the bifurcation diagram for fixed r?
• What is the effect of varying the starting value x_0 on the long-term attractor when α and β are fixed?
• Does the Lyapunov exponent become positive in the same parameter region where the bifurcation diagram becomes chaotic?
• To what extent can a closed-form Lyapunov exponent predict numerical chaos in a one-parameter sub-family?
• Which values of α and β produce the widest stable interval before the map becomes chaotic?

Basic Materials

• Laptop or desktop computer with Python installed.
• Python notebook environment, such as Jupyter or Google Colab.
• Spreadsheet software for organizing parameter sweeps.
• Graphing tool for plotting bifurcation diagrams and time series.
• Calculator or simple scripting environment for checking formulas by hand.
• Notebook for recording parameter choices, outputs, and observations.

Advanced Materials

• University or school workstation with higher RAM for large parameter sweeps.
• Python with NumPy, SciPy, and Matplotlib for simulation and plotting.
• Symbolic math software such as Mathematica, Maple, or SymPy for algebraic checks.
• Version control tool such as Git for tracking code changes.
• High-resolution plotting software for publication-quality figures.
• Access to a mentor or advisor familiar with nonlinear dynamical systems.

Software & Tools

• Python: Simulates the map, sweeps parameters, and computes Lyapunov estimates.
• Jupyter Notebook: Keeps code, notes, and plots in one place.
• NumPy: Handles fast array-based iteration across many parameter values.
• Matplotlib: Draws bifurcation diagrams and time-series plots.
• SymPy: Helps check symbolic steps when you try to derive a closed form.

Experiment Steps

1. Define the exact one-parameter sub-family you will analyze, and decide which variables stay fixed.
2. Choose the numerical range for r, α, β, and x_0, and justify why those ranges matter.
3. Build an iteration script that records both transient behavior and long-term behavior.
4. Plan how you will detect fixed points, cycles, and chaotic regions from your output.
5. Derive or verify the Lyapunov exponent formula for the sub-family, then compare it with simulation.
6. Organize your results into bifurcation plots, stability tables, and a short proof or derivation outline.

Common Pitfalls

• Treating the first few iterates as data, which can hide the real long-term pattern.
• Using too few parameter values, which makes the bifurcation diagram look blurry or misleading.
• Ignoring domain restrictions on x, α, or β, which can produce invalid values or complex-number issues.
• Comparing theory and simulation without matching the same sub-family, which makes the Lyapunov check meaningless.
• Assuming a positive Lyapunov estimate always means random-looking noise, when it can also come from numerical error or poor resolution.

What Makes This Competitive

A stronger version of this project goes beyond making a pretty bifurcation plot. You would define the parameter space carefully, prove the formula for a special case, and test whether the formula matches numerical behavior across many starting values. You could also compare your map with the standard logistic map and explain exactly what the fractional exponents change. Careful error checks, clean visuals, and a clear stability argument would make the work much stronger.

Project Variations

• Study how the same map behaves when you hold β fixed and vary α across a wider range.
• Compare the fractional-order logistic map with the classic logistic map and measure how the bifurcation structure changes.
• Analyze how different initial values x_0 affect the time needed to reach the long-term pattern for the same parameter set.

Learn More

• MIT OpenCourseWare, Nonlinear Dynamics and Chaos: Search MIT OpenCourseWare for lecture materials on iterative maps, stability, and bifurcations.
• NIST Digital Library of Mathematical Functions: Use this for background on special functions and careful notation when your derivation gets technical.
• MathWorld, Logistic Map entry: A quick reference for standard logistic map behavior and related terms.
• arXiv: Search for preprints on fractional logistic maps, bifurcation diagrams, and Lyapunov exponents.
• SpringerOpen or peer-reviewed journal articles on nonlinear dynamics: Search for review papers and recent studies on generalized logistic maps.