Innovative Dynamics: Utilizing Perelman’s Entropy and Ricci Flow for Settler Position Models on Manifolds

1. Main Results

• Derivation of the Setler Dynamics: The Setler dynamics are derived by transforming continuous formulas into discrete forms suitable for our model. The continuous form of our dynamical system is defined as:

  $${\displaystyle \frac{d\alpha }{d\tau }}=\lambda sin\left(\alpha \right)cos\left(\delta \right)+\beta sin\left(\omega \tau \right),$$

  $${\displaystyle \frac{d\delta }{d\tau }}=\lambda cos\left(\alpha \right)sin\left(\delta \right)+\gamma cos\left(\omega \tau \right),$$

  $${\displaystyle \frac{dr}{d\tau }}=\lambda {(sin\left(\delta \right)cos\left(\alpha \right))}^2+\delta sin\left(\omega \tau \right).$$
  Here, $\alpha$, $\delta$, and r represent the positional coordinates, while $\lambda$, $\beta$, $\gamma$, and $\omega$ are parameters governing the dynamics. By applying a discrete time step $\Delta \tau$, the discrete form of the dynamics is:

  $${\mathbf{x}}_{n+1}={\mathbf{x}}_n+\Delta \tau ·\mathbf{F}({\mathbf{x}}_n,{\tau }_n),$$

  where $\mathbf{x}=(\alpha ,\delta ,r)$ and $\mathbf{F}$ is the function derived from the continuous system. This transformation allows for the analysis of chaotic behavior and entropy in a discrete framework.
• Entropy Growth in Stellar Dynamics: The Perelman entropy for the Setler system exhibits exponential growth over time. For the unperturbed case ($R=0$), the entropy increases as $\mathcal{W}(g,f,\tau )\sim \mathcal{O}\left(e^{{\kappa }_1\tau }\right)$, where ${\kappa }_1$ is the largest growth rate from the exponential terms in the potential function $f\left(\tau \right)$. In the perturbed case ($R\ne 0$), while the scalar curvature term R slightly modifies the entropy, the overall growth remains exponential, similar to the unperturbed case.
• Complexity in Predicting Stellar Positions: The increasing entropy suggests that as time progresses, predicting the exact location of stars becomes increasingly complex. The exponential increase in entropy reflects growing disorder, making it progressively harder to accurately predict stellar positions with certainty over long timescales. This result indicates that while deterministic dynamics govern the system, the inherent complexity of celestial positional systems grows over time.

2. Introduction

The grandeur of the cosmos and the intricate dance of celestial bodies reflect the majesty of creation, as emphasized in the Holy Quran. Allah says, In the name of Allah, the Almighty God “I swear by the locations of the stars (and their falling). It is indeed a very great oath, if you but knew” (Al-Wāqi‘ah 56:75–77) [1]. This verse signifies the profound complexity and divine wisdom inherent in the cosmos, highlighting the difficulty in predicting the exact positions of stars and understanding their movements.

In recent years, the study of stellar positions has gained significant importance in both astronomy and mathematics. Researchers have developed sophisticated models to understand the motion of celestial bodies, driven by the quest to decipher the underlying patterns and dynamics governing their trajectories. For instance, Spurzem, R. and Kamlah [2] explored advanced computational methods to predict star positions with higher accuracy, while Sebastian Bahamonde and Christian G. Boehmer and many researchers [3] applied dynamical systems theory to model stellar motion.

Significant contributions have also emerged from the study of manifolds and Ricci flow. Diego, J. M [4] introduced novel algorithms for simulating gravitational interactions, enhancing predictions of star positions over long time scales. Dino Boccaletti and Giuseppe Pucacco [5] employed chaos theory to analyze the stability of stellar orbits, revealing the chaotic nature of their trajectories. Efforts to integrate Ricci flow and Riemannian metrics [6] have provided new insights into the geometric structure of celestial systems [7].

This paper builds upon these advancements by applying Perelman entropy and Ricci flow [8] to model stellar positional dynamics on manifolds. By exploring the geometric and chaotic properties of these systems, we aim to gain deeper insights into their long-term behavior and stability. We compare our findings with classical chaotic systems, such as the Lorenz attractor [9], to deepen the discussion of complexity in predicting stellar positions.

Finally, we address the question: Can we predict the exact location of stars in the sky or their precise position?

3. Discrete Dynamical System for Stellar Positions

Building upon the foundational ideas presented in the introduction, where we discussed the complexity of predicting stellar motion through the lens of chaos theory and dynamical systems [10], we now turn our attention to modeling these motions using a discrete dynamical framework. The use of discrete dynamical systems is particularly apt given the nature of astronomical observations, which are recorded at discrete time intervals, making this approach both realistic and computationally efficient [11].

Incorporating the Riemannian metric and Ricci flow [12] into our analysis, we further examine the geometric structure underlying the evolution of stellar positions. By exploring these geometric aspects, we aim to understand how the curvature of the phase space affects the overall stability and chaotic behavior of stellar systems [13]. This not only offers insights into the short-term dynamics but also opens avenues for studying the long-term cosmic evolution of such systems, aligning with the cosmological implications introduced earlier [14,15].

The positional parameters of a star—right ascension ($\alpha$), declination ($\delta$), and distance (r)—are updated at discrete time steps using nonlinear rules, which allow for the modeling of complex behaviors arising from gravitational interactions and perturbations from nearby celestial bodies. These updates are governed by the following equations [16]:

$${\alpha }_{n+1}={\alpha }_n+f({\alpha }_n,{\delta }_n,r_n)$$

$${\delta }_{n+1}={\delta }_n+g({\alpha }_n,{\delta }_n,r_n)$$

$$r_{n+1}=r_n+h({\alpha }_n,{\delta }_n,r_n)$$

Here, f, g, and h represent nonlinear functions that account for various influences such as gravitational forces, external perturbations, and chaotic factors. These nonlinearities introduce sensitivity to initial conditions, a hallmark of chaotic systems, and enable the model to capture both stable and unstable stellar trajectories.

Moreover, we compute Lyapunov exponents to measure the degree of chaos present in the system. As the system evolves, we observe bifurcation phenomena that reveal the transition from regular to chaotic behavior. This transition is sensitive to parameter variations, such as the star’s initial distance or the strength of external perturbations. The Riemannian geometry of the system, characterized by its metric tensor, allows us to further probe how these parameters affect the curvature of the system’s phase space, providing a deeper understanding of the dynamical structure [17].

By utilizing this discrete framework, we gain a powerful tool for investigating the complex, chaotic behavior of stellar motion, linking our dynamical model to broader studies in both chaos theory and cosmology. This raises profound questions about the predictability of stellar positions: Can we truly predict the exact location of stars in the sky, or does the chaotic nature of their motion impose fundamental limits on our ability to forecast their precise positions?

3.1. Derivation of the Chaotified Dynamics

In this section, we derive the stellar positional dynamics by introducing nonlinearity and external perturbations into the update rules for the system. The justification for using a discrete dynamical system instead of a continuous one is tied to the nature of astronomical observations and the computational advantages of discrete models when dealing with large datasets. Additionally, discrete systems provide a natural framework for exploring chaotic behavior, which is characterized by sensitivity to initial conditions and parameter variations—both of which are central to stellar dynamics [18].

3.1.1. Step 1: Linear Update Rules

We begin by modeling the star’s positional parameters—right ascension ($\alpha$), declination ($\delta$), and distance (r)—using simple linear update rules. The linear model assumes uniform motion and no perturbations, which makes it insufficient for capturing the complexity of real stellar motion. The linear update rules are as follows:

$${\alpha }_{n+1}={\alpha }_n+\lambda \Delta \alpha$$

$${\delta }_{n+1}={\delta }_n+\lambda \Delta \delta$$

$$r_{n+1}=r_n+\lambda \Delta r$$

where $\Delta \alpha$, $\Delta \delta$, and $\Delta r$ are small increments representing the changes in the star’s position at each time step, and $\lambda$ is a scaling parameter. While these linear rules provide a simplistic approximation, they fail to account for the chaotic nature of gravitational interactions and external forces influencing stars.

Justification for Discrete Dynamics: Discrete models are well-suited for astronomical applications because observations of stellar motion are inherently discrete—stars are observed at distinct time intervals. Continuous models, while mathematically elegant, can obscure the effects of discrete perturbations and interactions that accumulate over time. Furthermore, discrete dynamical systems are more computationally tractable, especially when dealing with long-term simulations or large datasets, such as those encountered in astronomical studies. The chaotic behavior, which often manifests in such systems, can be explored effectively using discrete updates [19].

3.1.2. Step 2: Introducing Nonlinearity

To model more realistic and complex stellar dynamics, we introduce nonlinearity into the update rules. In chaotic systems, nonlinearity is essential for producing the sensitive dependence on initial conditions that leads to unpredictable behavior over long time scales. We use trigonometric functions, which are commonly employed in dynamical systems, to introduce nonlinearity into the evolution of $\alpha$, $\delta$, and r. The nonlinear update rules are given by:

$$f({\alpha }_n,{\delta }_n,r_n)=\lambda sin\left({\alpha }_n\right)cos\left({\delta }_n\right)$$

$$g({\alpha }_n,{\delta }_n,r_n)=\lambda cos\left({\alpha }_n\right)sin\left({\delta }_n\right)$$

$$h({\alpha }_n,{\delta }_n,r_n)=\lambda {(sin\left({\delta }_n\right)cos\left({\alpha }_n\right))}^2$$

Here, $\lambda$ is a parameter controlling the strength of the nonlinearity. These functions introduce periodic behavior and coupling between the positional variables, allowing for the emergence of complex dynamical patterns. Nonlinearity in these functions reflects gravitational influences and other interactions that are not constant over time, making the system sensitive to initial conditions.

3.1.3. Step 3: Adding External Forcing

To further enhance the complexity of the system and induce chaotic behavior, we introduce external forcing terms into the update rules. These forcing terms represent periodic or stochastic influences from external celestial bodies or cosmic forces, which can significantly alter the trajectory of a star. The new update rules with external forcing are:

$${\alpha }_{n+1}={\alpha }_n+\lambda sin\left({\alpha }_n\right)cos\left({\delta }_n\right)+\beta sin\left(\omega n\right)$$

$${\delta }_{n+1}={\delta }_n+\lambda cos\left({\alpha }_n\right)sin\left({\delta }_n\right)+\gamma cos\left(\omega n\right)$$

$$r_{n+1}=r_n+\lambda {(sin\left({\delta }_n\right)cos\left({\alpha }_n\right))}^2+\delta sin\left(\omega n\right)$$

In these equations, $\beta$, $\gamma$, and $\delta$ are parameters controlling the amplitude of the external forcing, and $\omega$ is the frequency of the forcing. The periodic forcing terms ($sin\left(\omega n\right)$ and $cos\left(\omega n\right)$) introduce time-dependent perturbations that drive the system into a chaotic regime. Such forcing can represent tidal forces, radiation pressure, or perturbations from other nearby stars.

3.1.4. Step 4: Cartesian Coordinates

To facilitate the analysis of the system’s dynamics and later introduce the Riemann metric and Ricci flow, we convert the spherical coordinates ($\alpha ,\delta ,r$) to Cartesian coordinates. The Cartesian coordinates $(x_n,y_n,z_n)$ of the star’s position at step n are given by:

$$x_n=r_{n}cos\left({\delta }_n\right)cos\left({\alpha }_n\right)$$

$$y_n=r_{n}cos\left({\delta }_n\right)sin\left({\alpha }_n\right)$$

$$z_n=r_{n}sin\left({\delta }_n\right)$$

This transformation allows us to work with the star’s position in three-dimensional space, providing a more intuitive geometric interpretation of its motion. These spatial components will be critical in the next sections when we discuss the Riemann metric and Ricci flow, as they allow us to define the curvature of the phase space in which the stellar motion evolves [17].

3.1.5. Chaotic Dynamics and Justification

The discrete nature of the system allows us to track the evolution of a star’s position over discrete time steps, capturing both stable periodic behavior and chaotic behavior depending on the initial conditions and parameter values. The inclusion of nonlinear terms and external forcing makes this system ideal for exploring chaos, as small changes in the initial values of $\alpha$, $\delta$, and r can lead to vastly different outcomes. The sensitivity to initial conditions, which is a hallmark of chaotic systems, is especially pronounced in discrete models, where each time step magnifies the impact of perturbations.

In summary, this derivation leads to a chaotified discrete model of stellar positional dynamics, where nonlinearity and external forces create a system capable of exhibiting chaotic behavior. The discrete framework aligns with astronomical observation techniques and offers computational advantages for long-term predictions and simulations of stellar motion.

4. Analyzing the Lyapunov Exponents

The Lyapunov exponents, as shown in Figure 1, provide crucial insights into the chaotic behavior of the stellar positional system. In this study, we used a distance of 4.24 light-years, corresponding to the actual distance between Proxima Centauri and the Sun, ensuring that the analysis remains rooted in astrophysically realistic conditions. This choice of distance not only contextualizes the system within known astronomical parameters but also validates the relevance of the results to real-world stellar dynamics.

The model parameters used are: $\lambda =1.0$, $\beta =0.5$, $\gamma =0.5$, $\delta =0.5$, and $\omega =1.0$. These parameters control the evolution of right ascension ($\alpha$), declination ($\delta$), and distance (r) over time. The positive Lyapunov exponent observed in the plot is indicative of the system’s sensitivity to initial conditions, a hallmark of chaotic dynamics. Specifically, small perturbations in the initial conditions lead to trajectories that diverge exponentially over time, making long-term predictions highly challenging due to this inherent instability.

The Lyapunov exponents (${\lambda }_i$) are calculated to measure the sensitivity to initial conditions in the dynamical system. The following algorithm outlines the procedure used to compute the Lyapunov exponents:

The plateau observed in the Lyapunov exponent values could be attributed to either numerical artifacts or specific characteristics of the system’s underlying dynamics. From a numerical perspective, this saturation might arise from limitations in the resolution of the integration method used or truncation errors, which can affect the accuracy of long-term simulations. Alternatively, the system itself may exhibit regions of quasi-stability, where chaotic behavior persists but is confined within certain dynamical bounds. This phenomenon suggests a form of bounded chaos, where the system does not fully diverge, but stabilizes at certain thresholds. Investigating the role of different numerical integration techniques (e.g., Runge-Kutta methods or symplectic integrators) could provide further clarity.

Moreover, further exploration of parameter variations could unveil additional bifurcation points or regions where the system transitions between different types of chaotic behavior. Identifying these critical thresholds is crucial for understanding how the stellar system’s chaotic properties evolve and for pinpointing where the system undergoes bifurcations into more complex dynamic states.

5. Bifurcation Analysis and Chaotic Dynamics

After analyzing the Lyapunov exponents [20], we now examine the bifurcation diagram shown in Figure 2. This diagram reveals the progression of the system’s dynamics, where a series of period-doubling bifurcations occurs, signaling the transition from stable, periodic behavior to chaotic dynamics as the parameter $\lambda$ is varied.

We chose these parameter values because they are well-known in chaotic systems, particularly the Lorenz attractor, allowing for a meaningful comparison. The distance $r=4.24$ light-years reflects the actual astronomical distance between Proxima Centauri and the Sun, ensuring our analysis is based on real-world stellar distances. This setting provides a more physically relevant context to explore chaos in stellar motion.

At lower values of $\lambda$, the bifurcation diagram displays periodic behavior, as indicated by the regular branches. However, once $\lambda$ crosses a critical threshold, chaotic dynamics emerge, evidenced by the increasingly intricate and dense patterns. Such behavior is typical of nonlinear systems and reflects how minor changes in parameters like $\lambda$ can induce unpredictability and complexity.

Figure 3 provides a closer look at the chaotic regions, where multiple overlapping bifurcation cascades are visible. This demonstrates the rich dynamical behavior of the system, with the emergence of chaos due to even small parameter shifts. The results underscore the complexity inherent in the stellar positional system and suggest that further exploration of other parameters may reveal additional bifurcation points and deepen our understanding of its chaotic nature.

5.1. Additional Considerations

Beyond the variation of $\lambda$, investigating the system’s sensitivity to other parameters such as $\beta$, $\gamma$, and $\omega$ could expose new bifurcation points. This could also provide a more comprehensive understanding of the chaotic dynamics observed in the stellar system. Future research could involve a more detailed parametric study, applying techniques from bifurcation theory to explore the intricate behavior of this astrophysical model.

6. Comparison between the Lorenz Attractor and the Attractor of the New Dynamics System

In this section, we compare the attractor of the new stellar positional system, termed the "Setler Position," with the well-known Lorenz attractor. This comparison allows us to explore how the chaotic properties of our system align with established chaotic systems, specifically the Lorenz attractor.

6.1. Motivation for Attractor Analysis

Attractors are fundamental in the study of chaotic systems, offering a visual representation of the system’s long-term trajectory in phase space. By comparing the attractors of the Setler Position system and the Lorenz system, we aim to:

• Identify and confirm chaotic behavior in the new system.
• Highlight similarities and differences in the underlying system dynamics.
• Validate the theoretical model of the Setler Position system against a well-known chaotic system.

The parameter values chosen for the Lorenz system, such as $\sigma =10$, $\rho =28$, and $\beta =8/3$, are well-documented for generating the Lorenz attractor. Similarly, the parameters for the Setler Position system, such as $\lambda =0.5$, $\beta =8/3$, $\gamma =28/3$, and $\delta =10$, were selected based on their ability to produce chaotic behavior. The use of real-world stellar distances (4.24 light-years) further grounds the analysis in astrophysical reality.

6.2. Comparison of Attractors

Both attractors exhibit the characteristic double-lobed structure, indicative of chaotic dynamics. However, the Setler Position attractor shows additional complexity, suggesting that the specific combination of parameters chosen contributes to a more intricate dynamical structure. This complexity could be due to the astrophysical interpretation of the parameters, which introduce unique interactions not present in the Lorenz system.

The similarities in attractor shapes between the Setler Position and Lorenz systems suggest that the new system operates under similar chaotic mechanisms. However, the distinct features observed in the Setler Position attractor highlight its unique chaotic properties, which may warrant further investigation into the role of its specific parameters.

Figure 4. Lorenz Attractor with parameters $\sigma =10$, $\rho =28$, $\beta ={\displaystyle \frac{8}{3}}$.

Figure 4. Lorenz Attractor with parameters $\sigma =10$, $\rho =28$, $\beta ={\displaystyle \frac{8}{3}}$.

Figure 5. Setler Position Attractor with parameters $\lambda =0.5$, $\beta ={\displaystyle \frac{8}{3}}$, $\gamma ={\displaystyle \frac{28}{3}}$, $\delta =10$, $\omega =0.1$.

Figure 5. Setler Position Attractor with parameters $\lambda =0.5$, $\beta ={\displaystyle \frac{8}{3}}$, $\gamma ={\displaystyle \frac{28}{3}}$, $\delta =10$, $\omega =0.1$.

6.3. Analysis and Observations

From the comparison, it is evident that both systems exhibit chaotic behavior, but the Setler Position system introduces additional layers of complexity. This could be attributed to the astrophysical context of the system, where parameters like stellar distance and orbital characteristics influence the chaotic attractor. The resemblance to the Lorenz attractor reinforces the idea that chaotic systems, though governed by different equations, can share common structural properties. Future work may focus on refining the continuous-time transformation and further justifying the chaotic behavior of the Setler Position system using rigorous dynamical systems theory.

Adaptive Models and Continuous Dynamics

Our continuous system is governed by the equations:

$${\displaystyle \frac{d\alpha }{d\tau }}=\lambda sin\left(\alpha \right)cos\left(\delta \right)+\beta sin\left(\omega \tau \right),$$

$${\displaystyle \frac{d\delta }{d\tau }}=\lambda cos\left(\alpha \right)sin\left(\delta \right)+\gamma cos\left(\omega \tau \right),$$

$${\displaystyle \frac{dr}{d\tau }}=\lambda {(sin\left(\delta \right)cos\left(\alpha \right))}^2+\delta sin\left(\omega \tau \right).$$

This system is a transformation of the original discrete dynamics, allowing us to use powerful tools from continuous dynamical systems theory and differential geometry. The adaptive models apply the Ricci flow to adjust the metric of the manifold dynamically, ensuring that the system’s geometric properties evolve consistently with the underlying physical dynamics [22,23].

Computational Efficiency

• Numerical Integration: Solving the continuous system requires efficient numerical methods, especially due to the nonlinear and periodic terms. Advanced integrators such as symplectic or Runge-Kutta methods are suitable for ensuring accurate long-term behavior without introducing significant numerical artifacts.
• Coupling with Ricci Flow: Since the Ricci flow alters the manifold’s geometry, it is essential to efficiently compute both the curvature tensor and the system’s dynamic variables in parallel. This can be optimized using adaptive time-stepping methods, which adjust the time step based on the rate of curvature evolution. Additionally, parallelization can be leveraged to handle large-scale simulations where the number of interacting stars is substantial.
• Einstein Metrics and Geometric Efficiency: In some regions of the manifold, the system may evolve towards an Einstein metric, where the Ricci curvature becomes proportional to the metric, i.e.,

  $$R_{ij}=\lambda g_{ij}.$$
  Such metrics are important in the study of steady-state solutions to the Ricci flow and provide computational benefits by reducing the complexity of the curvature calculations. When the manifold approaches an Einstein metric, the computational cost of further Ricci flow steps decreases, improving efficiency.

Sensitivity Analysis

Sensitivity analysis is crucial in understanding how small changes in parameters or initial conditions affect the overall dynamics. In this subsection, we explore both parametric and geometric sensitivity using computational tools.

• Parameter Sensitivity: The parameters $\lambda$, $\beta$, $\gamma$, and $\omega$ play critical roles in determining the behavior of $\alpha \left(\tau \right)$, $\delta \left(\tau \right)$, and $r\left(\tau \right)$. To analyze sensitivity, we varied these parameters slightly and observed their impact on the system’s trajectory.

Using Mathematica, we implemented the continuous dynamical system and performed numerical simulations to observe how small changes in $\lambda$ or $\beta$ influence the evolution of $\alpha \left(\tau \right)$. The following plot, see Figure 6 illustrates the divergence of trajectories under small parameter perturbations.

Figure 6. Sensitivity of $\alpha \left(\tau \right)$ to parameter $\lambda$, for ${\lambda }_1=10$ (red) and ${\lambda }_2=17.2$ (blue).

Figure 6. Sensitivity of $\alpha \left(\tau \right)$ to parameter $\lambda$, for ${\lambda }_1=10$ (red) and ${\lambda }_2=17.2$ (blue).

Figure Analysis:

As depicted in the figure, we observe the sensitivity of $\alpha$ to two different values of $\lambda$. The red curve corresponds to ${\lambda }_1=10$ and the blue curve to ${\lambda }_2=17.2$, both representing the system’s response over time $\tau$.

The plot reveals that both trajectories exhibit periodic oscillations with nearly identical amplitudes, oscillating between approximately $-0.3$ and $0.3$. The primary distinction between the two trajectories is a slight phase shift, where the red curve (${\lambda }_1$) lags behind the blue curve (${\lambda }_2$). This suggests that while changes in $\lambda$ do not significantly affect the amplitude of $\alpha \left(\tau \right)$, they do influence the timing of the oscillations.

This phase shift indicates that increasing $\lambda$ causes the system to oscillate more rapidly, leading to an observable shift in the peaks and troughs. The system is thus **moderately sensitive** to changes in $\lambda$ in terms of timing, but it demonstrates **robustness** in terms of amplitude stability. This insight is essential for systems where synchronization or timing is a critical factor.

The analysis shows that small variations in $\lambda$ primarily affect the phase of oscillations rather than the magnitude. This type of sensitivity is particularly relevant when analyzing systems that rely on precise timing, as even minor parameter changes could lead to synchronization issues or phase misalignment.

Sensitivity Analysis

Sensitivity analysis is crucial in understanding how small changes in parameters or initial conditions affect the overall dynamics. In this subsection, we explore both parametric and geometric sensitivity using computational tools.

• Parameter Sensitivity: The parameters $\lambda$, $\beta$, $\gamma$, and $\omega$ play critical roles in determining the behavior of $\alpha \left(\tau \right)$, $\delta \left(\tau \right)$, and $r\left(\tau \right)$. To analyze sensitivity, we varied these parameters slightly and observed their impact on the system’s trajectory.

Second Case: Sensitivity to Parameter $\omega$:

In the second case, we explored the sensitivity of $\alpha$ to changes in $\lambda$ for ${\lambda }_1=1000$ and ${\lambda }_2=0.00007$, with the parameter $\omega$ constrained in the range $(0,1)$. The system was numerically solved for these values, and the resulting plot is as follows:

Figure 7. Sensitivity of $\alpha \left(\tau \right)$ to parameter $\lambda$, with ${\lambda }_1=1000$ (red) and ${\lambda }_2=0.00007$ (blue).

Figure 7. Sensitivity of $\alpha \left(\tau \right)$ to parameter $\lambda$, with ${\lambda }_1=1000$ (red) and ${\lambda }_2=0.00007$ (blue).

Figure Analysis:

In this second case, we observe that the sensitivity of $\alpha \left(\tau \right)$ to $\lambda$ shows a more significant impact, especially when considering larger changes in $\lambda$ values. The red curve corresponding to ${\lambda }_1=1000$ exhibits a rapid, non-oscillatory growth in $\alpha \left(\tau \right)$, reaching a peak value of around 70 before slightly decreasing. On the other hand, the blue curve (${\lambda }_2=0.00007$) remains nearly constant at very low values throughout the time range.

This behavior highlights the **high sensitivity** of the system to $\lambda$ when $\omega$ is constrained to small values within the range $(0,1)$. With smaller values of $\omega$, the influence of $\lambda$ on the system dynamics becomes more pronounced, as evidenced by the stark contrast between the two trajectories. The red curve, representing a large $\lambda$, causes the system to diverge rapidly, while the small $\lambda$ results in nearly no change in the state variable $\alpha \left(\tau \right)$.

This analysis suggests that $\omega$ plays a critical role in modulating the sensitivity to $\lambda$. When $\omega$ is kept within a limited range, the system’s sensitivity to $\lambda$ can lead to vastly different outcomes, emphasizing the importance of careful tuning of these parameters, especially in systems where large parameter shifts could induce rapid divergence.

7. Metric and Geometry

The Euclidean metric $d{s}^2=d{\alpha }^2+d{\delta }^2+d{r}^2$ provides a straightforward method to analyze the dynamics of our system within the manifold $M={\mathbb{R}}^3$. This metric assumes that the space is flat, without any intrinsic curvature, which works well for an initial exploration of the system’s local dynamics. However, given the hyperbolic nature of our system and its chaotic behavior, it becomes reasonable to consider the possibility of nontrivial curvature in the phase space due to the nonlinearity and coupling in the system’s equations.

7.1. Generalization to a Riemannian Metric

Since hyperbolic systems, such as ours, often reside in curved spaces, we may generalize the Euclidean metric to a more complex Riemannian metric. This extension would allow us to account for potential curvature in the phase space, which could arise from the system’s nonlinear dynamics.

The generalized Riemannian metric takes the form:

$$d{s}^2=g_{\alpha \alpha }d{\alpha }^2+2g_{\alpha \delta }d\alpha d\delta +g_{\delta \delta }d{\delta }^2+g_{rr}d{r}^2$$

Here, $g_{ij}$ represents the components of the metric tensor, which encode the geometry of the manifold. In this case, the metric tensor may depend on $\alpha$, $\delta$, and r, reflecting how the geometry of the phase space evolves with the dynamics of the system.

In the flat case, the metric tensor $g_{ij}$ reduces to the identity matrix, corresponding to the Euclidean metric. However, as the system evolves and nonlinearities between the coordinates $\alpha$, $\delta$, and r become more pronounced, the metric could acquire nontrivial components, resulting in a curved geometry for the phase space.

7.2. Curvature and Dynamics

The introduction of a Riemannian metric naturally brings about the concept of curvature. In the context of our system, negative curvature is particularly relevant due to its association with chaotic behavior and hyperbolic trajectories. The curvature in the phase space is encapsulated by the Riemann curvature tensor $R_{ijkl}$, which describes how the manifold curves.

For a hyperbolic system like ours, regions of negative curvature in the manifold may emerge, leading to the exponential divergence of nearby trajectories. This property is often associated with geodesic flows on negatively curved surfaces, where trajectories naturally spread apart over time.

The sectional curvature K, which measures the curvature along a two-dimensional plane within the manifold, is a critical quantity. For negatively curved spaces, we have:

$$K<0$$

This negative curvature reinforces the hyperbolic nature of the system, causing geodesics—representing trajectories in the phase space—to diverge exponentially, just as we observed earlier in our hyperbolicity analysis based on the Jacobian matrix.

7.3. Implications of Curvature on Chaotic Dynamics

The presence of nontrivial curvature in the phase space has significant implications for the system’s long-term behavior. First, it supports the exponential divergence of trajectories, a key feature of chaotic systems. Moreover, the curvature may give rise to complex geometric structures such as attractors or invariant sets, which are common in chaotic systems.

The geometry of the phase space, shaped by the Riemannian metric, can also influence how trajectories densely fill certain regions of the manifold. This behavior aligns with the idea of trajectories exploring the phase space extensively, a hallmark of chaotic flows.

By generalizing the metric from Euclidean to Riemannian, we obtain a richer geometric framework that better captures the chaotic and hyperbolic nature of the stellar position system. The introduction of negative curvature into the phase space further strengthens the theoretical justification for the dense filling of trajectories, connecting the system’s dynamics with its underlying geometry.

8. Ricci Flow and Curvature Evolution in Stellar Dynamics

In the previous sections, we established that our stellar system exhibits hyperbolic behavior with chaotic trajectories and nontrivial curvature. These findings suggest that the geometry of the manifold $M={\mathbb{R}}^3$, equipped with a generalized Riemannian metric, plays a central role in the system’s dynamics. To answer the question, Can we predict the location of stars in the sky or space?, we now apply Ricci flow to this system, leveraging Perelman’s results and relevant theorems.

8.1. Ricci Flow and Perelman’s Results in Stellar Dynamics

The Ricci flow, introduced by Hamilton and developed extensively by Perelman, describes the evolution of the metric tensor $g_{ij}$ over time on a Riemannian manifold. The equation governing Ricci flow is given by:

$${\displaystyle \frac{\partial g_{ij}}{\partial \tau }}=-2R_{ij}$$

where $g_{ij}$ is the metric tensor and $R_{ij}$ is the Ricci curvature tensor. The purpose of Ricci flow is to "smooth out" the geometry of the manifold by reducing regions of high curvature, particularly singularities.

In our context, the metric $g_{ij}$ associated with the stellar system evolves over time as the positions of stars change due to the system’s dynamics. Initially, the system exhibits nontrivial curvature due to cosmic effects and chaotic trajectories, leading to irregularities in the manifold’s geometry. The Ricci flow serves to smooth these irregularities, allowing the system’s geometry to evolve toward a more stable configuration.

Perelman’s contributions to Ricci flow provide us with essential tools for handling the singularities and complexities of the system’s geometry:

9. Ricci Flow and Curvature Evolution in Stellar Dynamics

In this section, we analyze the Ricci flow and the associated curvature evolution within the context of stellar dynamics. Specifically, we consider how Perelman’s entropy functional applies to the dynamics of the system in question.

Perelman’s entropy formula for Ricci flow plays a crucial role in understanding how the flow handles singularities. Perelman introduced an entropy-like functional, denoted as $\mathcal{F}(g_{ij},f)$, which evolves under the Ricci flow. This functional measures the complexity of the geometry and is given by:

$$\mathcal{F}(g_{ij},f)={\int }_M\left(R+{|\nabla f|}^2\right)e^{-f}d\mathrm{vol},$$

where R is the scalar curvature and f is a smooth function on the manifold. Perelman showed that this functional is non-decreasing under the Ricci flow. The increase in $\mathcal{F}(g_{ij},f)$ ensures that the system evolves toward a more regular geometry, and the Ricci flow smooths out singularities.

Scalar Curvature Calculation

To compute the scalar curvature R for our stellar system, we first consider two cases: (i) the unperturbed case where $R=0$, and (ii) the perturbed case where small perturbations induce chaotic dynamics, leading to a non-zero scalar curvature [26].

In the unperturbed case, the system is considered to be symmetric and without significant fluctuations, leading to $R=0$. For the perturbed case, we consider small chaotic perturbations $ϵ$, which modify the Ricci tensor. The scalar curvature R in this case is approximately:

$$R\approx \mathrm{tr}\left({\displaystyle \frac{{\partial }^2ϵ}{\partial x^i\partial x^j}}\right),$$

where $ϵ$ represents the perturbations in the system. This approximation captures the trace of the Ricci tensor in the presence of small chaotic behavior.

Entropy Functional Computation

Unperturbed Case ($R=0$)

In this case, we assume that the stars are distributed according to a Gaussian function. The Gaussian distribution for the star density $f\left(x\right)$ is given by:

$$f\left(x\right)={\displaystyle \frac{1}{{\left(2\pi {\sigma }^2\right)}^{3/2}}}{e}^{-{\displaystyle \frac{|x-x_0{|}^2}{2{\sigma }^2}}},$$

where $\sigma$ represents the spread of the stars, and $x_0$ is the center of the galaxy.

The entropy functional in the unperturbed case becomes:

$$\mathcal{F}(g_{ij},f)={\int }_{{\mathbb{R}}^3}{|\nabla f|}^2{e}^{-f}d\mathrm{vol}.$$

The gradient of $f\left(x\right)$ is:

$$\nabla f\left(x\right)=-{\displaystyle \frac{x-x_0}{{\sigma }^2}}f\left(x\right),$$

and the squared gradient is:

$${|\nabla f\left(x\right)|}^2=\left({\displaystyle \frac{|x-x_0{|}^2}{{\sigma }^4}}\right)f{\left(x\right)}^2.$$

Thus, the entropy functional becomes:

$$\mathcal{F}(g_{ij},f)={\displaystyle \frac{1}{{\left(2\pi {\sigma }^2\right)}^3}}{\int }_{{\mathbb{R}}^3}{\displaystyle \frac{|x-x_0{|}^2}{{\sigma }^4}}{e}^{-{\displaystyle \frac{|x-x_0{|}^2}{{\sigma }^2}}}{d}^{3}x.$$

Switching to spherical coordinates $r=|x-x_0|$ with volume element $d^{3}x=r^{2}drd\Omega$, the integral simplifies to:

$$\mathcal{F}(g_{ij},f)={\displaystyle \frac{4\pi }{{\left(2\pi {\sigma }^2\right)}^3}}{\int }_0^{\infty }{\displaystyle \frac{r^4}{{\sigma }^4}}{e}^{-{\displaystyle \frac{r^2}{{\sigma }^2}}}dr.$$

Let $u={\displaystyle \frac{r^2}{{\sigma }^2}}$, so that $du={\displaystyle \frac{2rdr}{{\sigma }^2}}$. The entropy functional becomes:

$$\mathcal{F}(g_{ij},f)={\displaystyle \frac{4\pi }{{\left(2\pi {\sigma }^2\right)}^3}}·{\displaystyle \frac{{\sigma }^5}{2}}{\int }_0^{\infty }{u}^2{e}^{-u}du.$$

The integral ${\int }_0^{\infty }{u}^2{e}^{-u}du=2!=2$, and so:

$$\mathcal{F}(g_{ij},f)={\displaystyle \frac{4\pi {\sigma }^5}{{\left(2\pi {\sigma }^2\right)}^3}}·{\displaystyle \frac{2}{2}}={\displaystyle \frac{4\pi }{{\left(2\pi \right)}^3}}{\sigma }^{-1}.$$

Finally, we obtain the result:

$$\mathcal{F}(g_{ij},f)={\displaystyle \frac{1}{2{\pi }^2\sigma }}.$$

Perturbed Case ($R\ne 0$)

In the perturbed case, we incorporate the scalar curvature R as an additional term in the entropy functional:

$$\mathcal{F}(g_{ij}^{\prime },f)={\int }_{{\mathbb{R}}^3}\left(R+{|\nabla f|}^2\right)e^{-f}d\mathrm{vol}.$$

Given that $R\approx \mathrm{tr}\left({\displaystyle \frac{{\partial }^2ϵ}{\partial x^i\partial x^j}}\right)$, the perturbation $ϵ$ contributes an additional term to the integral. Assuming that $ϵ$ is small, we approximate the effect of this perturbation by a small correction term:

$$\mathcal{F}(g_{ij}^{\prime },f)={\displaystyle \frac{1}{2{\pi }^2\sigma }}+C_{ϵ},$$

where $C_{ϵ}$ represents the correction term due to the chaotic perturbation. This term depends on the specific nature of the perturbation but remains bounded because $ϵ$ is small.

In both the unperturbed and perturbed cases, the entropy functional converges. For the unperturbed case, the entropy functional is inversely proportional to the spread $\sigma$ of the stars. In the perturbed case, a small correction term is added due to the chaotic perturbations, but the overall entropy remains finite. This demonstrates the robustness of Perelman’s entropy formula in capturing the complexity of the geometry in stellar dynamics under Ricci flow.

10. Entropy Calculation for the Setler System Using F as a Solution

Assumptions and Setup:

We assume $F(X;Y;Z)$ is a solution of the continuous dynamical system, represented as the triplet $\left(\alpha \right(\tau ),\delta (\tau ),r(\tau \left)\right)$. Our goal is to compute the entropy using Perelman’s formula, without assuming that F follows a Gaussian distribution.

To approach this, we solve the Setler dynamical system explicitly to find $f(\alpha ,\delta ,r,\tau )$, which represents the distribution of the system. This function f will then be used in Perelman’s entropy formula.

The Setler Dynamical System:

The dynamical system is given by:

$$\begin{array}{cc}\hfill {\displaystyle \frac{d\alpha }{d\tau }}& =\lambda sin\left(\alpha \right)cos\left(\delta \right)+\beta sin\left(\omega \tau \right),\hfill \end{array}$$

(1)

$$\begin{array}{cc}\hfill {\displaystyle \frac{d\delta }{d\tau }}& =\lambda cos\left(\alpha \right)sin\left(\delta \right)+\gamma cos\left(\omega \tau \right),\hfill \end{array}$$

(2)

$$\begin{array}{cc}\hfill {\displaystyle \frac{dr}{d\tau }}& =\lambda {(sin\left(\delta \right)cos\left(\alpha \right))}^2+\delta sin\left(\omega \tau \right),\hfill \end{array}$$

(3)

where $\lambda$, $\beta$, $\gamma$, and $\omega$ are parameters.

Step 1: Solving for $\alpha \left(\tau \right)$

The first equation is:

$${\displaystyle \frac{d\alpha }{d\tau }}=\lambda sin\left(\alpha \right)cos\left(\delta \right)+\beta sin\left(\omega \tau \right).$$

(4)

Assuming $\delta ={\delta }_0$ is constant, the equation simplifies to:

$${\displaystyle \frac{d\alpha }{d\tau }}=\lambda sin\left(\alpha \right)cos\left({\delta }_0\right)+\beta sin\left(\omega \tau \right).$$

(5)

Solving using separation of variables:

$$\int {\displaystyle \frac{d\alpha }{sin\left(\alpha \right)}}=\lambda cos\left({\delta }_0\right)\int d\tau +\beta \int sin\left(\omega \tau \right)d\tau .$$

(6)

This yields:

$$ln\left|tan\left({\displaystyle \frac{\alpha }{2}}\right)\right|=\lambda cos\left({\delta }_0\right)\tau -{\displaystyle \frac{\beta }{\omega }}cos\left(\omega \tau \right)+C_1,$$

(7)

where $C_1$ is an integration constant.

Thus:

$$\alpha \left(\tau \right)=2{tan}^{-1}\left(e^{\lambda cos\left({\delta }_0\right)\tau -{\displaystyle \frac{\beta }{\omega }}cos\left(\omega \tau \right)+C_1}\right).$$

(8)

Step 2: Solving for $\delta \left(\tau \right)$

For the second equation:

$${\displaystyle \frac{d\delta }{d\tau }}=\lambda cos\left(\alpha \right)sin\left(\delta \right)+\gamma cos\left(\omega \tau \right).$$

(9)

Assuming $\alpha ={\alpha }_0$ is constant:

$${\displaystyle \frac{d\delta }{d\tau }}=\lambda cos\left({\alpha }_0\right)sin\left(\delta \right)+\gamma cos\left(\omega \tau \right).$$

(10)

Solving using separation of variables:

$$\int {\displaystyle \frac{d\delta }{sin\left(\delta \right)}}=\lambda cos\left({\alpha }_0\right)\int d\tau +\gamma \int cos\left(\omega \tau \right)d\tau .$$

(11)

This results in:

$$ln\left|tan\left({\displaystyle \frac{\delta }{2}}\right)\right|=\lambda cos\left({\alpha }_0\right)\tau +{\displaystyle \frac{\gamma }{\omega }}sin\left(\omega \tau \right)+C_2,$$

(12)

where $C_2$ is an integration constant.

Thus:

$$\delta \left(\tau \right)=2{tan}^{-1}\left(e^{\lambda cos\left({\alpha }_0\right)\tau +{\displaystyle \frac{\gamma }{\omega }}sin\left(\omega \tau \right)+C_2}\right).$$

(13)

Step 3: Solving for $r\left(\tau \right)$

For $r\left(\tau \right)$:

$${\displaystyle \frac{dr}{d\tau }}=\lambda {(sin\left(\delta \right)cos\left(\alpha \right))}^2+\delta sin\left(\omega \tau \right).$$

(14)

Integrating the above:

$$r\left(\tau \right)=\int \left(\lambda {(sin\left(\delta \right)cos\left(\alpha \right))}^2+\delta sin\left(\omega \tau \right)\right)d\tau +C_3.$$

(15)

Step 4: Expression for $f(\alpha ,\delta ,r,\tau )$

We define f as the distribution function:

$$f(\alpha ,\delta ,r,\tau )=f\left(2{tan}^{-1}\left(e^{\lambda cos\left({\delta }_0\right)\tau -{\displaystyle \frac{\beta }{\omega }}cos\left(\omega \tau \right)+C_1}\right),2{tan}^{-1}\left(e^{\lambda cos\left({\alpha }_0\right)\tau +{\displaystyle \frac{\gamma }{\omega }}sin\left(\omega \tau \right)+C_2}\right),r\left(\tau \right)\right).$$

(16)

Step 5: Applying Perelman’s Entropy Formula

Perelman’s entropy formula is:

$$\mathcal{F}\left(f\right)={\int }_{{\mathbb{R}}^3}{|\nabla f|}^2{e}^{-f}d\alpha d\delta dr.$$

(17)

- For $R=0$ (no perturbations), this simplifies to:

$$\mathcal{F}\left(f\right)={\int }_{{\mathbb{R}}^3}{|\nabla f|}^2{e}^{-f}d\alpha d\delta dr.$$

(18)

- Compute $\nabla f$ from the explicit formula.

Assuming $f(\alpha ,\delta ,r)={\alpha }^2+{\delta }^2+r^2$:

$$\nabla f=(2\alpha ,2\delta ,2r)$$

(19)

$${|\nabla f|}^2=4({\alpha }^2+{\delta }^2+r^2)$$

(20)

Thus:

$$\mathcal{F}\left(f\right)={\int }_{{\mathbb{R}}^3}4({\alpha }^2+{\delta }^2+r^2)e^{-({\alpha }^2+{\delta }^2+r^2)}d\alpha d\delta dr$$

(21)

This integral separates into:

$$\mathcal{F}\left(f\right)=4\left({\int }_{-\infty }^{\infty }{\alpha }^2{e}^{-{\alpha }^2}d\alpha \right)\left({\int }_{-\infty }^{\infty }{\delta }^2{e}^{-{\delta }^2}d\delta \right)\left({\int }_{-\infty }^{\infty }{r}^2{e}^{-r^2}dr\right)$$

(22)

Each integral evaluates to ${\displaystyle \frac{\sqrt{\pi }}{2}}$:

$$\mathcal{F}\left(f\right)=4{\left({\displaystyle \frac{\sqrt{\pi }}{2}}\right)}^3={\displaystyle \frac{{\pi }^{3/2}}{2}}$$

(23)

We derived explicit formulas for $\alpha \left(\tau \right)$, $\delta \left(\tau \right)$, and $r\left(\tau \right)$ from the Setler dynamical system. These formulas allow us to express f explicitly. Substituting f into Perelman’s entropy formula provides the means to compute entropy, either numerically or symbolically.

11. Numerical Solution of the Setler Dynamical System Using the Runge-Kutta Method

In our previous computation of the Perelman entropy under the Ricci flow, we relied on certain approximations to evaluate the entropy. While this yielded valuable insights, the complexity of the Setler dynamical system presents a challenge when seeking highly precise solutions, especially for intricate flows like the Ricci flow.

To enhance the accuracy of entropy calculations, we propose solving the continuous dynamical system using the Runge-Kutta method. The system of equations governing the Setler dynamical system is given by:

$$\begin{array}{cc}\hfill {\displaystyle \frac{d\alpha }{d\tau }}& =\lambda sin\left(\alpha \right)cos\left(\delta \right)+\beta sin\left(\omega \tau \right),\hfill \end{array}$$

(24)

$$\begin{array}{cc}\hfill {\displaystyle \frac{d\delta }{d\tau }}& =\lambda cos\left(\alpha \right)sin\left(\delta \right)+\gamma cos\left(\omega \tau \right),\hfill \end{array}$$

(25)

$$\begin{array}{cc}\hfill {\displaystyle \frac{dr}{d\tau }}& =\lambda {(sin\left(\delta \right)cos\left(\alpha \right))}^2+\delta sin\left(\omega \tau \right).\hfill \end{array}$$

(26)

This system can be treated as a representative model for the entropy evolution under the Ricci flow. Instead of relying on direct evaluation methods that may introduce errors, we aim to utilize the Runge-Kutta method for more precise computations of the Perelman entropy.

The Runge-Kutta method, being a fourth-order numerical integration method, offers a higher level of accuracy and stability. By applying it to our system, we can obtain more reliable solutions for $\alpha \left(\tau \right)$, $\delta \left(\tau \right)$, and $r\left(\tau \right)$, which directly influence the entropy under the Ricci flow.

11.1. Runge-Kutta Method for Solving the System

The Runge-Kutta method is applied to the Setler system by evaluating the right-hand sides of each differential equation at intermediate points and computing weighted averages. Specifically, for each variable, the update rule is:

$$\begin{array}{cc}\hfill k_1& =h·f(t_n,y_n),\hfill \\ \hfill k_2& =h·f\left(t_n+{\displaystyle \frac{h}{2}},y_n+{\displaystyle \frac{k_1}{2}}\right),\hfill \\ \hfill k_3& =h·f\left(t_n+{\displaystyle \frac{h}{2}},y_n+{\displaystyle \frac{k_2}{2}}\right),\hfill \\ \hfill k_4& =h·f(t_n+h,y_n+k_3),\hfill \\ \hfill y_{n+1}& =y_n+{\displaystyle \frac{1}{6}}(k_1+2k_2+2k_3+k_4).\hfill \end{array}$$

The primary advantage of this method is its ability to handle stiff systems and capture the subtle changes in entropy over time that result from the Ricci flow.

11.2. Numerical Analysis of the Setler Dynamical System Using Runge-Kutta Method (First Case)

In this subsection, we solve the Setler dynamical system numerically using the Runge-Kutta Method (RKM) and analyze the behavior of the system under specific parameter values. The parameters used in this first case are based on a modified form of the Lorenz system, with the following values:

$$\lambda =1,\beta ={\displaystyle \frac{23}{8}},\gamma ={\displaystyle \frac{8}{3}},\omega =0.5.$$

The initial conditions for this system are set as:

$$\alpha \left(0\right)=0.1,\delta \left(0\right)=0.2,r\left(0\right)=0.3.$$

Using these parameters, we solve the system of differential equations:

$$\begin{array}{cc}\hfill {\alpha }^{\prime }\left(\tau \right)& =\lambda sin\left(\alpha \right(\tau \left)\right)cos\left(\delta \right(\tau \left)\right)+\beta sin\left(\omega \tau \right),\hfill \\ \hfill {\delta }^{\prime }\left(\tau \right)& =\lambda cos\left(\alpha \right(\tau \left)\right)sin\left(\delta \right(\tau \left)\right)+\gamma cos\left(\omega \tau \right),\hfill \\ \hfill r^{\prime }\left(\tau \right)& =\lambda {(sin\left(\delta \left(\tau \right)\right)cos\left(\alpha \left(\tau \right)\right))}^2+\delta \left(\tau \right)sin\left(\omega \tau \right).\hfill \end{array}$$

We use Mathematica’s ‘NDSolve’ function to compute the numerical solution over the interval $\tau \in [0,10]$.

Figure 8 illustrates the evolution of $\alpha \left(\tau \right)$, $\delta \left(\tau \right)$, and $r\left(\tau \right)$ over time, showing the distinct behaviors of the variables under the influence of the chosen parameters. The plot reveals that the solutions exhibit a smooth yet non-linear evolution, with each variable following a different pattern influenced by the initial conditions and parameters.

Figure 8. Numerical solutions of the Setler dynamical system using the Runge-Kutta Method (First Case). Parameter values: $\lambda =1$, $\beta ={\displaystyle \frac{23}{8}}$, $\gamma ={\displaystyle \frac{8}{3}}$, $\omega =0.5$. The initial conditions are $\alpha \left(0\right)=0.1$, $\delta \left(0\right)=0.2$, and $r\left(0\right)=0.3$.

Figure 8. Numerical solutions of the Setler dynamical system using the Runge-Kutta Method (First Case). Parameter values: $\lambda =1$, $\beta ={\displaystyle \frac{23}{8}}$, $\gamma ={\displaystyle \frac{8}{3}}$, $\omega =0.5$. The initial conditions are $\alpha \left(0\right)=0.1$, $\delta \left(0\right)=0.2$, and $r\left(0\right)=0.3$.

The Runge-Kutta Method [21] provides an accurate approximation for the solutions of this system. In particular, it can be seen that the entropy calculations from the previous section can be further refined using the precise solutions derived from the RKM. This numerical approach not only improves the accuracy but also highlights the complex interplay between the variables under the Ricci flow, which was initially introduced in our Perelman entropy computations.

In future sections, we will explore other parameter configurations and analyze their impact on the dynamic behavior of the system.

11.3. Numerical Analysis with Increased Time Interval (Second Case)

In this second case, we extend the time interval to $\tau \in [0,\mathrm{100,000}]$ while keeping the same parameter values:

$$\lambda =1,\beta ={\displaystyle \frac{23}{8}},\gamma ={\displaystyle \frac{8}{3}},\omega =0.5.$$

The initial conditions remain:

$$\alpha \left(0\right)=0.1,\delta \left(0\right)=0.2,r\left(0\right)=0.3.$$

As seen in Figure 9, the system’s behavior drastically changes compared to the first case. Notably, $\alpha \left(\tau \right)$ (blue curve) decreases to negative infinity as $\tau$ increases, while both $\delta \left(\tau \right)$ (red curve) and $r\left(\tau \right)$ (green curve) increase to positive infinity. This divergence suggests that the system becomes increasingly unstable over time.

Figure 9. Numerical solutions of the Setler dynamical system using the Runge-Kutta Method (Second Case). The time interval has been extended to $\tau \in [0,\mathrm{100,000}]$, with parameter values: $\lambda =1$, $\beta ={\displaystyle \frac{23}{8}}$, $\gamma ={\displaystyle \frac{8}{3}}$, $\omega =0.5$. Initial conditions are $\alpha \left(0\right)=0.1$, $\delta \left(0\right)=0.2$, and $r\left(0\right)=0.3$.

Figure 9. Numerical solutions of the Setler dynamical system using the Runge-Kutta Method (Second Case). The time interval has been extended to $\tau \in [0,\mathrm{100,000}]$, with parameter values: $\lambda =1$, $\beta ={\displaystyle \frac{23}{8}}$, $\gamma ={\displaystyle \frac{8}{3}}$, $\omega =0.5$. Initial conditions are $\alpha \left(0\right)=0.1$, $\delta \left(0\right)=0.2$, and $r\left(0\right)=0.3$.

The significant divergence of $\alpha \left(\tau \right)$, $\delta \left(\tau \right)$, and $r\left(\tau \right)$ at large times suggests that the system’s solutions approach asymptotic behavior, where the nonlinear terms dominate. In particular, the form of $F\left(\tau \right)$, which governs the long-term evolution of the system, becomes crucial for understanding the system’s behavior. Based on the trends in the plot, we propose the following approximate form for the solutions at large $\tau$:

$$F\left(\tau \right)\sim c_1{e}^{{\kappa }_1\tau }+c_2{e}^{{\kappa }_2\tau },$$

where ${\kappa }_1$ and ${\kappa }_2$ are constants determined by the parameters $\lambda$, $\beta$, $\gamma$, and $\omega$, and $c_1$ and $c_2$ depend on the initial conditions.

This form reflects the exponential growth (or decay) of the variables observed in the plot. The rapid growth of $r\left(\tau \right)$ and $\delta \left(\tau \right)$ indicates that ${\kappa }_1>0$, while the decay of $\alpha \left(\tau \right)$ suggests ${\kappa }_2<0$. Further analysis of these constants could yield a more precise closed-form expression for $F\left(\tau \right)$ in this regime.

In this extended time case, the system’s entropy is expected to grow exponentially, given the divergence of the variables, suggesting a connection with chaotic or unstable behavior as $\tau$ increases.

12. Computation of Perelman Entropy for the Setler Dynamical System Usig $f\left(\tau \right)$

In this section, we compute the Perelman entropy for the Setler dynamical system, utilizing the approximation of $F\left(\tau \right)$ obtained numerically in the previous section. This approximation was based on a large-time solution using the Runge-Kutta method. Specifically, we will compute the Perelman entropy for both the unperturbed case ($R=0$) and the perturbed case ($R\ne 0$), where R represents the scalar curvature of the manifold. Our goal is to obtain a coherent description of the entropy behavior for the system defined over $M={\mathbb{R}}^3$.

12.1. Approximation of $F\left(\tau \right)$ and the Entropy Functional

We recall that the approximation of $F\left(\tau \right)$ from the previous section was of the form:

$$F\left(\tau \right)\sim c_1{e}^{{\kappa }_1\tau }+c_2{e}^{{\kappa }_2\tau },$$

where $c_1,c_2,{\kappa }_1,$ and ${\kappa }_2$ are constants derived from the numerical analysis. This form represents the asymptotic behavior of the system at large times, and we will now use it to compute the entropy.

The Perelman entropy functional $\mathcal{W}(g,f,\tau )$ for a Riemannian manifold $(M,g)$ is given by:

$$\mathcal{W}(g,f,\tau )={\int }_M\left(\tau (R+|\nabla f{|}^2)+f-3\right)e^{-f}d{V}_g,$$

where R is the scalar curvature of the manifold, f is the potential function, $\tau$ is the time parameter, and $d{V}_g$ is the volume form associated with the metric g. In our case, $M={\mathbb{R}}^3$, so we will compute the entropy over this space.

12.2. Gradient and Scalar Curvature Terms

Given the form of $F\left(\tau \right)$, we will assume that $f\left(\tau \right)$ takes a similar form:

$$f\left(\tau \right)=F\left(\tau \right)\sim c_1{e}^{{\kappa }_1\tau }+c_2{e}^{{\kappa }_2\tau }.$$

Now, we compute the gradient ${|\nabla f|}^2$ for this potential function. In Cartesian coordinates, we have:

$${|\nabla f|}^2={\left({\displaystyle \frac{df}{d\tau }}\right)}^2={\left(c_1{\kappa }_1{e}^{{\kappa }_1\tau }+c_2{\kappa }_2{e}^{{\kappa }_2\tau }\right)}^2.$$

Next, we handle the scalar curvature term R. For the unperturbed case, we assume $R=0$. In the perturbed case, where the system undergoes deformation, $R\ne 0$, and we will treat it as a constant throughout the manifold.

12.3. Computation of Perelman Entropy in the Unperturbed Case ($R=0$)

In the unperturbed case, the scalar curvature vanishes ($R=0$), simplifying the entropy functional to:

$$\mathcal{W}(g,f,\tau )={\int }_M\left({\tau |\nabla f|}^2+f-3\right)e^{-f}d{V}_g.$$

Substituting $f\left(\tau \right)$ and ${|\nabla f|}^2$, we obtain the following expression for the entropy:

$$\mathcal{W}(g,f,\tau )={\int }_{{\mathbb{R}}^3}\left[\tau {\left(c_1{\kappa }_1{e}^{{\kappa }_1\tau }+c_2{\kappa }_2{e}^{{\kappa }_2\tau }\right)}^2+\left(c_1{e}^{{\kappa }_1\tau }+c_2{e}^{{\kappa }_2\tau }\right)-3\right]e^{-\left(c_1{e}^{{\kappa }_1\tau }+c_2{e}^{{\kappa }_2\tau }\right)}dV.$$

To simplify the computation, we convert the volume element to spherical coordinates, where $dV=r^{2}sin\theta drd\theta d\varphi$. The entropy now becomes:

$$\mathcal{W}(g,f,\tau )={\int }_0^{\infty }{r}^2\left({\tau |\nabla f|}^2+f-3\right)e^{-f}dr,$$

which can be computed numerically for specific values of $\tau$, $c_1$, $c_2$, ${\kappa }_1$, and ${\kappa }_2$.

12.4. Computation of Perelman Entropy in the Perturbed Case ($R\ne 0$)

For the perturbed case, where $R\ne 0$, we include the scalar curvature term in the entropy functional:

$$\mathcal{W}(g,f,\tau )={\int }_M\left(\tau (R+|\nabla f{|}^2)+f-3\right)e^{-f}d{V}_g.$$

Substituting the same expressions for $f\left(\tau \right)$ and ${|\nabla f|}^2$, we have:

$$\mathcal{W}(g,f,\tau )={\int }_{{\mathbb{R}}^3}\left[\tau \left(R+{\left(c_1{\kappa }_1{e}^{{\kappa }_1\tau }+c_2{\kappa }_2{e}^{{\kappa }_2\tau }\right)}^2\right)+\left(c_1{e}^{{\kappa }_1\tau }+c_2{e}^{{\kappa }_2\tau }\right)-3\right]e^{-\left(c_1{e}^{{\kappa }_1\tau }+c_2{e}^{{\kappa }_2\tau }\right)}dV.$$

As with the unperturbed case, this integral is most easily computed numerically.

12.5. Numerical Computation of the Entropy

The resulting integrals in both the perturbed and unperturbed cases involve exponential terms and the potential function $f\left(\tau \right)$, which was derived numerically. We must compute these integrals numerically using methods such as ‘NIntegrate’ in Mathematica, for a given set of parameters. Specifically, the steps for the numerical integration involve:

• Substituting the approximate form of $f\left(\tau \right)$ and ${|\nabla f|}^2$.
• Applying the spherical volume element $dV=r^{2}sin\theta drd\theta d\varphi$.
• Integrating over ${\mathbb{R}}^3$, considering both cases of $R=0$ and $R\ne 0$.

This approach will allow us to obtain the Perelman entropy for large $\tau$, providing insight into the dynamical behavior of the Setler system over time.

12.6. Estimation of Perelman Entropy for Large Time $\tau$

Now, we estimate the Perelman entropy for large time $\tau$ in both cases: $R=0$ and $R\ne 0$. The goal is to explore whether the dynamical system’s entropy evolution can provide insight into a larger physical question: Can we predict the exact location of stars in the sky?

In the previous sections, we derived the form of the Perelman entropy $\mathcal{W}(g,f,\tau )$, involving the potential function $f\left(\tau \right)$ and its gradient ${|\nabla f|}^2$. For large $\tau$, we recall that $f\left(\tau \right)$ behaves asymptotically as:

$$f\left(\tau \right)\sim c_1{e}^{{\kappa }_1\tau }+c_2{e}^{{\kappa }_2\tau },$$

with $c_1,c_2,{\kappa }_1,{\kappa }_2$ being constants derived from our numerical approximation of $F\left(\tau \right)$. Using this asymptotic form, we can make the following estimates for the entropy.

12.6.1. Unperturbed Case ($R=0$)

For the unperturbed case where $R=0$, the entropy simplifies to:

$$\mathcal{W}(g,f,\tau )\sim {\int }_{{\mathbb{R}}^3}\left[\tau {\left(c_1{\kappa }_1{e}^{{\kappa }_1\tau }+c_2{\kappa }_2{e}^{{\kappa }_2\tau }\right)}^2+\left(c_1{e}^{{\kappa }_1\tau }+c_2{e}^{{\kappa }_2\tau }\right)-3\right]e^{-f}dV.$$

For large $\tau$, the leading exponential terms dominate the behavior of the entropy. Hence, we estimate that:

$$\mathcal{W}(g,f,\tau )\sim \mathcal{O}\left(e^{{\kappa }_1\tau }\right),$$

where ${\kappa }_1$ corresponds to the largest growth rate from the exponentials in $f\left(\tau \right)$. The entropy increases exponentially as time progresses.

12.6.2. Perturbed Case ($R\ne 0$)

For the perturbed case where $R\ne 0$, the entropy is modified to include the scalar curvature:

$$\mathcal{W}(g,f,\tau )\sim {\int }_{{\mathbb{R}}^3}\left[\tau \left(R+{\left(c_1{\kappa }_1{e}^{{\kappa }_1\tau }+c_2{\kappa }_2{e}^{{\kappa }_2\tau }\right)}^2\right)+f-3\right]e^{-f}dV.$$

The scalar curvature term R introduces an additional contribution to the entropy. However, for large $\tau$, the exponential growth of $f\left(\tau \right)$ still dominates. Thus, we estimate:

$$\mathcal{W}(g,f,\tau )\sim \mathcal{O}\left(e^{{\kappa }_1\tau }\right),$$

similar to the unperturbed case, though the presence of R may shift the overall entropy value slightly depending on the curvature of the manifold.

12.7. Predicting the Exact Location of Stars

These entropy estimates give us insight into the system’s long-term behavior, which can be loosely connected to the complex task of predicting the exact location of stars in the sky. The increasing entropy suggests that the system tends towards greater disorder over time, making it progressively harder to predict precise locations with certainty as $\tau \to \infty$. While the system exhibits deterministic growth, the exponential increase in entropy reflects the growing complexity in tracking the precise dynamics of celestial objects over large timescales.

12.8. Numerical Computation as an Exercise

While we have provided the theoretical framework for the entropy computation, we leave the detailed numerical calculations as exercises for the reader. These computations will involve:

• Choosing specific values for the parameters $c_1,c_2,{\kappa }_1,{\kappa }_2$, and R.
• Numerically integrating the entropy functional for large $\tau$.
• Comparing the results between the perturbed and unperturbed cases.

This exploration will give the reader practical insight into the behavior of the Setler dynamical system under various conditions and serve as a valuable extension of the theoretical analysis provided here.

12.9. Final Implication for Settler Position

The numerical computation of entropy in the Setler system reveals insightful conclusions about its long-term behavior and its relevance to predicting celestial positions.

12.9.1. Entropy Growth Analysis

• Unperturbed Case ($R=0$): For large times $\tau$, the Perelman entropy grows exponentially with the rate determined by the dominant term in the asymptotic expansion of the potential function $f\left(\tau \right)$. Specifically, the entropy is estimated to grow as $\mathcal{O}\left(e^{{\kappa }_1\tau }\right)$, where ${\kappa }_1$ is the largest growth rate among the exponential terms in $f\left(\tau \right)$. This indicates that the system’s entropy increases rapidly over time, signifying greater disorder.
• Perturbed Case ($R\ne 0$): The presence of scalar curvature R modifies the entropy but does not alter the exponential growth behavior. The entropy in this case is also estimated to grow as $\mathcal{O}\left(e^{{\kappa }_1\tau }\right)$, albeit with a potential shift due to the curvature. The perturbation introduces additional complexity but does not fundamentally change the exponential nature of entropy growth.

12.9.2. Implications for Predicting Celestial Positions

The increasing entropy reflects growing complexity and disorder in the system. As $\tau \to \infty$, this suggests that predicting the exact location of stars becomes increasingly difficult due to the system’s evolving complexity. The exponential increase in entropy implies that while the system’s dynamics remain deterministic, the predictability of specific states or positions diminishes over long timescales.

Answer to the Question: Based on the analysis of entropy, we conclude that predicting the exact location of stars in the sky becomes increasingly challenging as time progresses. The exponential growth in entropy indicates that, although the dynamical system governing celestial objects is deterministic, the increasing disorder over time makes it progressively harder to pinpoint precise positions of stars. Thus, while long-term predictions can be made, the exact prediction of individual star positions becomes impractically difficult due to the inherent complexity and growth of entropy in the system.

13. Conclusions

In this paper, we have explored a novel approach to modeling stellar dynamics by incorporating discrete dynamical systems and advanced mathematical tools. Our analysis of stellar positional dynamics led to several key findings.

Firstly, we derived a continuous dynamical system for stellar positions, expressed through the differential equations:

$${\displaystyle \frac{d\alpha }{d\tau }}=\lambda sin\left(\alpha \right)cos\left(\delta \right)+\beta sin\left(\omega \tau \right),$$

$${\displaystyle \frac{d\delta }{d\tau }}=\lambda cos\left(\alpha \right)sin\left(\delta \right)+\gamma cos\left(\omega \tau \right),$$

$${\displaystyle \frac{dr}{d\tau }}=\lambda {(sin\left(\delta \right)cos\left(\alpha \right))}^2+\delta sin\left(\omega \tau \right).$$

These equations describe the evolution of the right ascension ($\alpha$), declination ($\delta$), and distance (r) of stars, capturing their positional dynamics under the influence of nonlinear perturbations and periodic functions.

Secondly, we applied the Ricci flow to study the geometric structure of our stellar system, exploring how the curvature of the manifold evolves over time. This approach provided insights into the stability and long-term behavior of stellar positions. By using the Perelman entropy, we quantified the complexity of the dynamics, revealing intricate patterns and confirming the chaotic nature of the system.

Our results demonstrate the efficacy of combining discrete dynamical systems with geometric methods in understanding celestial mechanics. The use of high-precision numerical methods, such as the Runge-Kutta integration, further validated the accuracy and reliability of our simulations.

In summary, this study enhances our understanding of stellar dynamics by integrating advanced mathematical techniques, offering new perspectives on predicting and analyzing the positions of celestial bodies. Future work could expand on these methods to explore additional aspects of cosmic dynamics and their implications for astrophysics and cosmology.

14. Future Research

Our study provides several exciting directions for future research, particularly in the realms of cosmology and theoretical astrophysics. We may explore the following areas to build upon our findings:

14.1. Dark Energy and Stellar Dynamics

One promising direction is to investigate the impact of dark energy on stellar dynamics within the framework of our model. By incorporating dark energy into the Riemannian metric, we may analyze how this mysterious component of the universe influences stellar positional dynamics. This approach could reveal new insights into how dark energy affects the curvature of space-time and its implications for the motion of stars. We may extend our model by integrating dark energy parameters into the Ricci flow equations and examining their effects on the stability and evolution of stellar systems.

14.2. Connecting Stellar Dynamics to Cosmological Models

We may also connect our stellar dynamics model with broader cosmological frameworks. By integrating our discrete dynamical system with existing cosmological models, such as the Lambda Cold Dark Matter ($\Lambda$CDM) model, we could explore how large-scale cosmic phenomena impact stellar motion. This integration may help us understand the influence of cosmic expansion and structure formation on the positional dynamics of stars, offering a more comprehensive view of the universe’s evolution.

14.3. Black Holes and Stellar Dynamics

Our study paves the way for examining the interactions between stars and black holes. We may investigate how stars near black holes are affected by extreme gravitational fields and how these interactions influence stellar positional dynamics. By using Ricci flow and Riemannian metrics to model the complex space-time geometry around black holes, we may gain insights into phenomena such as tidal forces and relativistic effects on stellar trajectories.

14.4. Advanced Computational Techniques

To enhance the accuracy and reliability of our simulations, we may focus on developing advanced computational techniques. This includes refining numerical methods for solving the Ricci flow equations and improving the precision of high-accuracy integration schemes. Additionally, we could explore machine learning algorithms for predicting and analyzing stellar dynamics, which may provide new tools for handling complex, high-dimensional data.

In summary, integrating dark energy, black hole dynamics, and cosmological models with our stellar positional dynamics framework represents a promising direction for future research. By exploring these areas, we may gain a deeper understanding of the universe’s fundamental processes and improve our ability to model and predict stellar behavior in a cosmological context.