🌀 nonlinear sciences

Exploring Chaotic Motion of a Particle in the Centre of a Galaxy with a Prolate Halo

This study employs pseudo-Newtonian potentials and dynamical analysis tools like Poincaré sections and Maximum Lyapunov Exponents to demonstrate how the spin of a central supermassive black hole, modeled by the Artemova-Bjornsson-Novikov potential within a prolate shell-like halo, modulates the transition between order and chaos in the orbits of nearby particles.

Original authors: Uditi Nag, Yeasin Ali, Suparna Roychowdhury

Published 2026-02-18

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Uditi Nag, Yeasin Ali, Suparna Roychowdhury

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a galaxy not as a smooth, perfect whirlpool of stars, but as a busy, chaotic city center. At the very heart of this city sits a Supermassive Black Hole (SMBH). Think of this black hole as a giant, invisible whirlpool drain that pulls everything toward it with immense force.

Now, imagine this drain isn't sitting in empty space. It's surrounded by a giant, fuzzy cloud of invisible matter called a halo. In this specific study, the researchers are looking at a halo that is shaped like a rugby ball (elongated and stretched out), rather than a perfect sphere. They call this a "prolate" halo.

Here is the story of how a single star (or a "test particle") moves through this complex environment, explained simply:

1. The Setup: A Tug-of-War

The star is caught in a tug-of-war between two main forces:

The Black Hole: It wants to suck the star in. Because black holes are so heavy and spin, they twist space itself (like a spinning top twisting a blanket).
The Halo: The rugby-ball-shaped cloud of matter pushes and pulls on the star in weird ways because it's not perfectly round. It has "lopsided" parts (dipole) and "bulging" parts (quadrupole).

2. The Problem: Predicting the Path is Hard

If you throw a ball in a vacuum, you can easily predict where it will land. But in this galaxy, the rules are messy. The star's path is non-linear, meaning a tiny change in where it starts can send it careening in a completely different direction. This is called chaos.

To figure this out without needing a supercomputer to solve Einstein's complex equations (which describe how gravity warps space), the scientists used a "cheat code." They used a Pseudo-Newtonian Potential.

The Analogy: Imagine trying to simulate a rollercoaster. Instead of calculating every tiny vibration of the track and the wind resistance (General Relativity), you use a simplified set of rules that feels like the real thing but is much easier to calculate. This "ABN potential" is their cheat code to mimic the spinning black hole's gravity.

3. The Experiment: The "Strobe Light" View

How do you see if a path is chaotic? You can't just watch the star for a million years. Instead, the researchers used a technique called a Poincaré Section.

The Analogy: Imagine taking a photo of the star every single time it crosses a specific invisible line (the equator of the galaxy) moving upward.
- Regular Motion: If the photos form a perfect, smooth circle or line, the star is dancing in a predictable rhythm. It's like a clockwork toy.
- Chaotic Motion: If the photos are scattered randomly all over the page, the star is lost in a storm. It's like throwing confetti into a hurricane.
- Sticky Motion: Sometimes the star gets "stuck" in a messy zone, hovering between order and chaos, like a fly buzzing around a lightbulb before finally escaping.

4. The Big Discoveries

The researchers played with two main knobs to see how they changed the star's dance:

Knob A: The "Lopsidedness" of the Halo (The Ratio P)

They increased the difference between the "bulge" and the "side" of the halo.
Result: The more lopsided the halo, the more chaotic the star became. It's like if the rugby ball was actually a jagged, bumpy rock; the star would bounce around wildly.

Knob B: The Spin of the Black Hole (The Parameter a)

They made the black hole spin faster.
Result: Surprisingly, spinning faster made things more stable.
The Analogy: Think of a spinning top. When it spins slowly, it wobbles and falls over easily (chaos). When it spins incredibly fast, it stands up straight and stays stable. The black hole's spin acts like a stabilizer, organizing the star's path and reducing the chaos, even when the halo is very lopsided.

5. Newton vs. Einstein (The "Cheat Code" Check)

The scientists compared their "cheat code" (Newtonian physics with a twist) against the "real thing" (Relativistic physics).

When the black hole spins slowly, both methods gave similar results.
When the black hole spins fast, the "real" physics showed slightly more chaos than the "cheat code," but the general trend was the same: Spin = Stability.

The Takeaway

This paper tells us that in the messy center of a galaxy, chaos is the default, but the spin of the central black hole acts as a conductor. Even if the surrounding galaxy is lumpy and uneven, a fast-spinning black hole can help keep the stars in a somewhat orderly dance, preventing them from getting completely lost in the gravitational storm.

It's a reminder that in the universe, even in the most chaotic places, a little bit of spin can bring a sense of order.

1. Problem Statement

The central regions of galaxies host supermassive black holes (SMBHs) surrounded by complex matter distributions, such as nuclear star clusters and extended halos. While general relativity (GR) provides the most accurate description of SMBH environments, solving Einstein's equations for realistic galactic systems is computationally prohibitive. Furthermore, standard analytical potentials (e.g., Hernquist, NFW) often fail to capture specific structural asymmetries like prolate halos, shells, or rings.

The core challenge addressed in this study is understanding how the interplay between a central rotating SMBH and a non-spherical (prolate) galactic halo influences the orbital dynamics of test particles. Specifically, the authors investigate the transition from regular to chaotic motion driven by:

Halo Asymmetry: Modeled via dipole and quadrupole multipole moments.
SMBH Spin: The rotational effects of the black hole (Kerr-like).
Relativistic Effects: Comparing Newtonian dynamics against pseudo-Newtonian approximations of GR.

2. Methodology

The authors employ a hybrid modeling approach combining pseudo-Newtonian potentials with multipole expansions to simulate the gravitational field.

Potential Formulation:
- Central SMBH: Instead of a simple Newtonian monopole, the study uses the Artemova–Björnsson–Novikov (ABN) potential. This pseudo-Newtonian potential approximates the gravitational field of a Kerr-type (rotating) black hole, capturing relativistic effects like frame-dragging within a Newtonian framework.
- Galactic Halo: The surrounding prolate halo is modeled using a multipole expansion truncated at the quadrupole order.
  - Monopole: Represents the central SMBH mass.
  - Dipole ( $D$ ): Accounts for halo asymmetries.
  - Quadrupole ( $Q < 0$ ): Models the prolate (elongated) shape of the halo shell.
- Total Potential ( $\Phi_g$ ): The sum of the ABN potential, the dipole term ($Dz$), and the quadrupole term ( $-\frac{Q}{2}(2z^2 - \rho^2)$ ).
Dynamical Frameworks:
- Newtonian Dynamics: Solved using standard equations of motion derived from the effective potential $V(\rho, z) = \Phi_g + \frac{L^2}{2\rho^2}$ .
- Relativistic Dynamics: Uses modified equations of motion that incorporate the ABN potential's relativistic corrections while maintaining non-relativistic definitions for momentum to allow direct comparison.
Numerical Analysis:
- Poincaré Sections: Trajectories are integrated using MATLAB (ode45) in the equatorial plane ( $z=0$ ) with fixed angular momentum ( $L=4.2$ ) and energy ( $E=0.976$ ). The system is visualized via 2D slices $(\rho, p_\rho)$ at crossings where $z=0$ and $p_z > 0$ .
- Chaos Indicators:
  - Visual Inspection: Distinguishing between smooth invariant curves (regular), KAM-like tori (sticky/quasi-periodic), and scattered points (chaotic).
  - Maximum Lyapunov Exponent (MLE): Calculated to quantify dynamical instability and measure the exponential divergence of nearby trajectories, providing a coordinate-independent metric for chaos.

3. Key Contributions

Novel Potential Combination: The study uniquely combines the ABN pseudo-Newtonian potential (for SMBH spin) with a multipole expansion (dipole + quadrupole) to model a prolate halo, a configuration not extensively explored in previous literature.
Comparative Analysis: It provides a direct, systematic comparison between Newtonian and relativistic (pseudo-Newtonian) regimes under varying spin and halo parameters.
Quantification of Spin Effects: The work rigorously demonstrates how the spin parameter of the SMBH acts as a stabilizing mechanism against chaos in the presence of non-spherical halos.

4. Results

The study yields several critical findings regarding the stability of orbits:

Impact of Halo Asymmetry ( $|P|$ ):
- The ratio $P = D/Q$ (dipole-to-quadrupole) is a primary driver of chaos.
- As the magnitude $|P|$ increases (indicating stronger dipole asymmetry), the system transitions from regular to chaotic motion. This is confirmed by both Poincaré sections (scattering of points) and increasing Maximum Lyapunov Exponents (MLE).
Impact of SMBH Spin ( $a$ ):
- There is a distinct inverse relationship between spin and chaos.
- Increasing the spin parameter $a$ systematically suppresses chaotic behavior.
- At high spin values ( $a = 0.9$ ), the chaotic regions in the phase space shrink significantly, and the system exhibits more regular behavior compared to low spin ( $a = 0.1$ ).
Newtonian vs. Relativistic Dynamics:
- Low Spin/Low Asymmetry: Newtonian and relativistic models show nearly identical chaotic behavior.
- High Asymmetry ( $|P|=300$ ): Relativistic dynamics exhibit stronger chaos than Newtonian dynamics.
- High Spin ( $a=0.9$ ): The difference between the two frameworks diminishes; high spin reduces chaos in both models to comparable levels.
Sticky Orbits: The study identifies "sticky" orbits (quasi-periodic) that align with KAM tori, which act as transitional states between regular and fully chaotic motion.

5. Significance

This research advances the understanding of galactic nucleus dynamics by:

Refining Stability Models: It demonstrates that SMBH spin is a crucial factor in stabilizing stellar orbits against the disruptive effects of non-spherical galactic halos.
Methodological Utility: It validates the use of pseudo-Newtonian potentials (like ABN) as efficient and accurate tools for studying complex relativistic dynamics without the computational cost of full GR simulations.
Astrophysical Implications: The findings suggest that in galaxies with high-spin SMBHs, the chaotic mixing of stars in the nuclear region may be less pronounced than in low-spin systems, potentially affecting the long-term evolution of nuclear star clusters and the rate of tidal disruption events.
Diagnostic Tools: The correlation between the MLE and Poincaré sections provides a robust framework for identifying subtle dynamical instabilities in complex gravitational potentials.