A Kinematic Analysis of S-Cluster Trajectories

1. Introduction

For over three centuries, classical mechanics has accurately described celestial motions in weak-field regimes. However, this paradigm breaks down in the dense Galactic Center, dominated by the supermassive black hole Sagittarius A* (Sgr A*) [1,2]. Recent breakthroughs, particularly the imaging of Sgr A* [3] and the precise tracking of the dense stellar cluster (S-Cluster) in its surroundings by the GRAVITY collaboration [4,5,6,7], have opened a new window into the strong-field gravitational regime. In this extreme environment, the gravitational potential is so intense that Newtonian approximations fail: space becomes curved, time dilates, and closed Keplerian ellipses are replaced by complex precessing rosettes.

State-of-the-art studies primarily investigate phenomena such as strong-field gravitational lensing and black hole shadows [24], the Bayesian estimation of compact object parameters through near-infrared astrometry [25], the analytical modeling of post-Newtonian orbital effects and relativistic propagation delays for natural bodies like the S-stars or hypothetical binary pulsars [26,27]. While the observational community has successfully utilized stars like S2 as test candidates to verify General Relativity [14,31], the literature lacks a comprehensive framework dedicated to relativistic astrodynamics and active mission analysis in the strong-field regime. Indeed, the transition from passive observation to the kinematic evaluation of a spacecraft navigating near a supermassive black hole remains a largely unexplored gap. In this context, the classical Keplerian framework is no longer a valid approximation but a source of structural divergence: if an active probe were to operate in the vicinity of Sgr A*, relying on classical orbital mechanics, such as traditional Lambert solvers, would not merely result in minor kinematic errors, but lead to macroscopic trajectory divergences.

Therefore, the primary objective of this study is to investigate the orbital dynamics in the strong-field regime of Sgr A* and to quantify the physical limitations of classical models.

The study defines and measures the "Relativistic Miss Distance": the macroscopic discrepancy between the position predicted by classical mechanics and the actual location in curved spacetime. Ultimately, this work aims to mathematically isolate the kinematic phenomena caused by spacetime curvature, such as the attractive cubic term in the effective potential, the step-wise Schwarzschild precession, and the chronometric delay, demonstrating that rigorous, purely relativistic modeling is an absolute requirement for the accurate description of any body’s kinematic evolution in the core of the Milky Way.

2. Mathematical Model and Numerical Strategy

This section outlines the mathematical and computational framework developed to simulate spacecraft kinematics within the strong-field regime of Sgr A*.

2.1. the Schwarzschild Metric Assumption

In this study, the supermassive black hole Sgr A* has been modeled using the static, spherically symmetric Schwarzschild metric. While recent observational constraints suggest Sgr A* possesses a non-zero spin, the static approximation is rigorously justified for the spatial regime of the target S-Cluster stars (e.g., S62 and S4714; S22 and S91). At their pericenter distances (r ∼${10}^2-{10}^3{r}_s$1), the gravitational potential is heavily dominated by the mass monopole moment [17]. Consequently, the primary relativistic deviation from Newtonian dynamics is the prograde Schwarzschild precession, scaling as $\propto 1/r^2$ in the effective potential [8]. In contrast, spin-induced frame-dragging effects (Lense-Thirring precession) scale as $\propto J/r^3$ (J is the intrinsic angular momentum of the black hole). and remain sub-dominant, acting as negligible background noise rather than a primary dynamical driver at these orbital distances [17,32]. Therefore, the black hole’s mass is considered the sole relativistic driver for this kinematic analysis.

2.2. Equations of Motion and Proper Time

Diverging from classical mechanics, where coordinate time t serves as an absolute independent variable, General Relativity requires the formulation of kinematics using proper time $\tau$, the time measured by a clock co-moving with the spacecraft. Consequently, the state vector is expanded into a nine-dimensional array comprising the spacetime coordinates ($t,r,\theta ,\varphi$), the contravariant components of the 4-velocity $u^{\mu }={\displaystyle \frac{d{x}^{\mu }}{d\tau }}$, and the instantaneous vehicle mass m.

A fundamental geometric constraint of the Pseudo-Riemannian manifold is the 4-velocity normalization

$$g_{\mu \nu }{u}^{\mu }{u}^{\nu }=-c^2,$$

(1)

This invariant condition enforces a strictly timelike trajectory and is algorithmically utilized to continuously determine the temporal component $u^t={\displaystyle \frac{dt}{d\tau }}$ rigorously accounting for gravitational time dilation near Sgr A*.

For a propelled spacecraft, the free-fall geodesic path is altered by the engine thrust, leading to the forced geodesics equation:

$${\displaystyle \frac{d{u}^{\mu }}{d\tau }}=-{\mathsf{\Gamma }}_{\alpha \beta }^{\mu }{u}^{\alpha }{u}^{\beta }+{\displaystyle \frac{f^{\mu }}{m}},$$

(2)

where ${\mathsf{\Gamma }}_{\alpha \beta }^{\mu }$ are the Christoffel symbols encoding the spacetime curvature, and $f^{\mu }$ is the propulsive 4-force [9]. To preserve the invariant rest mass of the vehicle, the thrust vector maintains strict orthogonality to the 4-velocity [10]

$$g_{\mu \nu }{f}^{\mu }{u}^{\nu }=0.$$

(3)

Because the resulting system of coupled Ordinary Differential Equations (ODEs) is highly non-linear, analytical solutions are unavailable; furthermore, the extreme eccentricities of the S-Cluster stars induce severe timescale variations. To prevent truncation errors and strictly preserve the 4-velocity norm, the relativistic propagation is executed utilizing an 8th-order explicit Runge-Kutta integration scheme (DOP853) [11]. This high-fidelity numerical strategy ensures that the macroscopic trajectory divergences observed are genuine physical manifestations of curved spacetime rather than numerical artifacts.

2.3. S-Cluster Experimental Scenarios

Prior to deployment in the Galactic Center, the numerical framework was rigorously validated within the Solar System via an Earth-Mars transfer simulation. This weak-field test confirmed that, in environments with negligible spacetime curvature, the relativistic propagator naturally converges to classical Newtonian dynamics, yielding a negligible integration drift ($<5.5$ km) when benchmarked against NASA Mars Design Reference Architecture 5.0 [12].

Following this validation, the kinematic analysis was applied to the S-Cluster, dividing the study into two distinct dynamical regimes to systematically probe the gravitational potential of Sgr A*.

2.3.1. Deep Gravity Regime (S62 → S4714)

The first experimental scenario targets the innermost region of the S-Cluster, focusing on the highly eccentric stars S62 and S4714 ($e\approx 0.976$, $e\approx 0.985$; for the full list of orbital parameters, see Table A1 in Appendix A)2. Because their pericenters penetrate deeply into the relativistic potential well of Sgr A* (see Table A2 in Appendix A), these targets maximize strong-field signatures, such as the macroscopic Schwarzschild precession of the pericenter and severe gravitational time dilation. This regime serves as the primary stress-test to quantify the structural breakdown of classical Keplerian mechanics and the resulting spatial divergence [6,13].

2.3.2. Outer Cluster Regime (S22 → S91)

To validate the methodology across the wider galactic potential, the second scenario investigates a transfer between the stars S22 and S91. Situated at significantly larger radial distances ($a>{10}^4$ AU, see Table A1 in Appendix A)3 and exhibiting moderate eccentricities, these stars reside in a weak-field relativistic regime. Here, General Relativity transitions into a perturbative force, and the orbital geometry remains quasi-Keplerian. This scenario acts as a critical control test, demonstrating the numerical stability of the relativistic solver over multi-century flight times and confirming that the accumulated trajectory drift is a secular physical effect rather than a numerical artifact [6,14].

3. Results

This section presents the numerical results of the trajectory analysis within the Galactic Center. By comparing classical Newtonian predictions with fully relativistic propagation, we quantify the structural divergence of the spacecraft’s path and isolate the specific physical mechanisms responsible for the breakdown of classical astrodynamics in the strong-field regime.

3.1. Optimal Transfer Strategies and Mission Profiles

Before detailing the microscopic dynamical effects of General Relativity, it is necessary to establish the macroscopic mission profiles. To identify the optimal transfer trajectories, a classical Newtonian Lambert solver was employed across extensive departure and time-of-flight search grids, later coupled with a relativistic propagator. Table 1 summarizes the absolute minimum-energy solutions identified for the selected S-cluster regimes.

In the Deep Gravity Regime (S62 → S4714 transfer), the extreme orbital velocities of the targets dictate prohibitive energy requirements. A direct transfer (0 revolutions) demands an extreme $\Delta V$ of 1436.9 km/s to intercept the target in 5.41 years. Relaxing the geometric constraints to permit a multi-revolution phasing (3 full revolutions around Sgr A*), allows the spacecraft to execute massive orbital plane changes at more favorable true anomalies, reducing the energetic cost to 1335.8 km/s, albeit extending the flight time to 56.39 years.

Conversely, the Outer Cluster transfer (S22 → S91) requires 912.6 km/s over a vast 537-year duration. This indicates that while the local spacetime curvature is significantly weaker in this region, the global kinematics of the S-cluster still demand propulsion systems fundamentally beyond contemporary chemical or nuclear-thermal limits (as detailed in Appendix B).

Crucially, the "Final Error" metric in Table 1 quantifies the spatial divergence at the arrival epoch when propagating these classical initial conditions through the full Schwarzschild metric. These macroscopic drifts (up to 14.22 AU) demonstrate that while the classical solver identifies the correct energetic basin of attraction, it fails to provide a precise interception. This discrepancy necessitates the detailed relativistic error analysis and dynamical characterization presented in the subsequent sections.

3.2. Relativistic Dynamics at Pericenter

To fully interpret the macroscopic navigation errors that accumulate over the transfer arc, it is necessary to investigate the kinematics at the point of maximum gravitational stress: the periapsis. The following analysis characterizes the effective potential landscape and the velocity profile to determine the dynamical root causes of the trajectory divergence.

3.2.1. The Effective Potential Landscape

The spacecraft’s energetic boundaries and orbital geometry are governed by the Effective Potential $V_{eff}\left(r\right)$. The fundamental divergence from classical mechanics stems from the General Relativistic potential’s additional attractive cubic term [9,15,16]:

$$V_{GR}\left(r\right)=-{\displaystyle \frac{\mu }{r}}+{\displaystyle \frac{h^2}{2r^2}}-{\displaystyle \frac{\mu h^2}{c^2{r}^3}}$$

(4)

where h is the specific angular momentum. Unlike the classical Newtonian centrifugal barrier $(+h^2/2r^2)$, which prevents collapse, the relativistic correction dominates at small radii, creating a purely attractive potential well.

In the Schwarzschild metric, the ISCO (Innermost Stable Circular Orbit) represents the absolute theoretical limit for stable circular motion around a black hole, located at $r_{ISCO}=6GM/c^2=3r_s$. Beneath this critical radius, the stable minimum of the effective potential disappears, transforming into an inflection point. This boundary defines the so-called "plunge region": within this altitude, no stable circular orbits can exist, and any infinitesimal perturbation causes a particle to spiral inevitably into the central singularity, regardless of its angular momentum [15]. While the selected targets S62 and S4714 do not cross the ISCO, their deep penetration into the potential well and their extreme proximity to this boundary during periapsis passages further emphasize the structural instability and the highly non-linear nature of the dynamical environment the spacecraft must navigate.

In the Outer Cluster Regime (Figure 1b), the correction at pericenter ($r\approx 11921$ AU) is seven orders of magnitude weaker ($\sim -7·{10}^2$ J/kg) than the primary potential. This negligible disturbance cannot reshape the closed elliptical orbit and acts merely as a constant, weak drag, confirming a quasi-Newtonian local environment.

Conversely, in the Deep Gravity Regime (Figure 1a), the relativistic correction reaches a massive $-1.2·{10}^9$ J/kg at the transfer pericenter ($r\approx 187$ AU). This profound energy mismatch implies that the spacecraft is subjected to an intense, unmodeled specific force that sharply bends the trajectory. Ultimately, this demonstrates that the classical Lambert solver fails because it computes the injection velocity based on a significantly "shallower" potential well, completely ignoring the structural morphing of the orbital geometry induced by spacetime curvature.

3.2.2. Velocity Profile near Pericenter Passage

To isolate the instantaneous dynamical cause of the cumulative trajectory drift, a kinematic analysis is conducted at the periapsis. Because gravitational time dilation prevents a direct temporal comparison between the two models, the radial distance r is adopted as the common independent variable. Both Relativistic and Newtonian velocity profiles are computed and interpolated onto a dense logarithmic radial grid to extract the exact velocity deviation $\Delta v=v_{GR}-v_{Newt}$ at every point along the trajectory.

In the Deep Gravity Regime (Figure 2a), the spacecraft reaches periapsis velocities exceeding ${10}^4$ km/s. At the closest approach ($r\approx 187$ AU), the deviation peaks at a massive $1231.27$ m/s. This $\sim 1.2$ km/s discrepancy proves that General Relativity acts as a dominant force rather than a perturbative correction. If an onboard guidance system were to calculate an impulsive burn using Newtonian assumptions at this distance, the resulting $\Delta V$ would be fundamentally flawed, confirming the structural breakdown of classical mechanics in this environment.

Conversely, the analysis of the Outer Cluster (Figure 2b) offers a stark contrast. Despite the high orbital velocities of the target stars ($\sim 600$ km/s), the local spacetime curvature is significantly flatter. At periapsis ($r\approx 11921$ AU), the maximum deviation is merely $3.50$ m/s. This confirms that the local dynamics remain quasi-Newtonian and that the macroscopic spatial drift observed over the mission is strictly caused by the secular integration of this minute error over the 500-year flight time.

Figure 2. Velocity Profile at Pericenter Passage. Top: the absolute orbital velocity magnitude v as a function of the radial distance. Bottom: the instantaneous velocity deviation $\Delta v=v_{GR}-v_{Newt}$.

Figure 2. Velocity Profile at Pericenter Passage. Top: the absolute orbital velocity magnitude v as a function of the radial distance. Bottom: the instantaneous velocity deviation $\Delta v=v_{GR}-v_{Newt}$.

3.3. Secular Evolution and Orbital Reshaping

While the pericenter analysis isolates the instantaneous relativistic perturbation, its cumulative effect profoundly alters the global orbital geometry over time. To fully characterize the transition from quasi-Keplerian stability to the strong-field regime, it is necessary to investigate the secular evolution of the trajectories, demonstrating how spacetime curvature physically reshapes the orbital elements during the flight.

3.3.1. Schwarzschild Precession of Pericenter

The most prominent strong-field signature is the Schwarzschild Precession, causing a secular advance of the argument of pericenter ($\omega$) and generating characteristic unclosed "rosette" orbits. For an orbit with semi-major axis a and eccentricity e, the weak-field approximation predicts an angular shift per orbit equal to [14,17,18]:

$$\Delta \omega \approx {\displaystyle \frac{6\pi GM}{c^{2}a(1-e^2)}}.$$

(5)

To visualize the magnitude of this effect, the trajectories of the four reference stars were propagated over multiple periods.

In the Deep Gravity Regime (Figure 3a,b), the breakdown of the closed-orbit approximation is visually immediate. Due to the extreme eccentricities ($e>0.97$), the evolution of $\omega$ is not continuous but step-wise: it remains nearly constant during the long apocenter coasting phase and undergoes a sharp, impulsive jump exactly at periapsis. The resulting precession rates are massive, peaking at numerically integrated values of $\approx 0.{1361}^{°}$/yr for S62 and $\approx 0.{1578}^{°}$/yr for S4714. This behavior constitutes a structural cause of the Newtonian guidance failure: while the classical Lambert solver assumes a fixed inertial target, the real target’s orbital plane physically rotates by nearly a degree during the transfer arc, causing a severe geometric mismatch at the predicted intercept point.

Conversely, the Outer Cluster anchors (Figure 4a,b) exhibit a distinct, quasi-Newtonian morphology. Because of their lower eccentricities ($e\approx 0.3-0.4$), the relativistic perturbation is less impulsive and more distributed, manifesting as a continuous sinusoidal oscillation superimposed on a negligible linear trend ($\approx 0.{0000}^{°}$/yr). The trajectories appear as perfectly closed elliptical traces, confirming that, for these distant targets, the geometric rotation is dynamically suppressed. This validates the use of fixed Keplerian ephemerids for mission planning in the weak-field limit, while highlighting that in the inner core, the target’s relativistic rotation is a dominant factor that cannot be ignored.

A more extensive compilation of results is present in Table A3 in Appendix A.

3.3.2. Oscillation of Osculating Orbital Elements

In classical Keplerian mechanics, the semi-major axis a and eccentricity e are constants of motion, fixed by the initial energy and angular momentum. However, in the General Relativistic framework, these quantities lose their status as strict invariants. The non-Newtonian $1/r^3$ term in the effective potential induces continuous variations in the orbital geometry, meaning that the physical trajectory continuously "morphs" mid-flight. To quantify this structural divergence from the Lambert solver’s static assumptions, the osculating Keplerian elements were computed at every integration step across the three reference scenarios; the result is shown in Figure 5.

In the Deep Gravity Regime, the spacecraft experiences a dramatic, localized violation of Keplerian mechanics. As the vehicle approaches the periapsis, the orbital elements undergo a sharp, impulsive V-shaped excursion: the semi-major axis contracts by approximately $0.39\%$, and the eccentricity decreases by $0.08\%$. This confirms that the interaction is impulsive, effectively injecting the spacecraft into a different orbit after the strong-field scattering event.

Conversely, the Outer Cluster transfer exhibits a qualitative shift in behavior. The evolution of a and e transitions from sharp spikes to smooth, long-period variations. The relative variance is minute (∼$0.0008\%$), bordering on numerical precision limits. Nevertheless, the distinct curvilinear profile confirms that these are genuine physical oscillations caused by the weak but non-zero gradient of the relativistic correction over vast radial distances rather than numerical artifacts.

Finally, the Earth-Mars control transfer yields a perfectly flat variance ($0.0000\%$, within machine precision). This result serves as a definitive validation of the numerical framework, rigorously confirming that the orbital variations observed in the Galactic Center are authentic manifestations of curved spacetime rather than numerical artifacts or energy leakage in the solver.

3.3.3. Gravitational Time Dilation

In General Relativity, time is not an absolute background parameter: its flow depends on the local gravitational potential and the velocity of the observer. This creates a physical discrepancy between the coordinate time (t), measured by a distant observer and used to define planetary ephemerides, and the proper time ($\tau$), measured by the atomic clocks on board the spacecraft [15]:

$$d\tau ={\left(1-{\displaystyle \frac{2GM}{r}}\right)}^{1/2}dt.$$

(6)

From an engineering perspective, this phenomenon manifests as a chronometric desynchronization drift, quantified as ${\delta }_{lag}\left(t\right)=t-\tau$, as shown in Figure 6. This clock drift critically impacts navigation (ranging), communications (Doppler tracking), and on-board scheduling.

In the Deep Gravity Regime, the transfer exhibits a distinct non-linear behavior. As the spacecraft approaches the pericenter passage, the curve steepens significantly due to the extreme orbital velocities (∼$10\%c$) and the depth of the gravitational potential ($1-r_s/r$). The total accumulated lag over the ∼2000-day mission reaches a critical 5.81 hours. An uncorrected desynchronization of nearly 6 hours implies that the spacecraft’s internal timeline is completely decoupled from the mission control’s coordinate timeline, causing any pre-scheduled maneuver to fail entirely.

In the Outer Cluster scenario, the time dilation accumulates almost linearly. The total lag reaches 14.16 hours, but this drift is accumulated over a vast flight time of more than 500 years. Here, the effect is not driven by a sudden plunge into the potential well, but by the relentless, secular integration of a minute potential difference over centuries.

Finally, the Earth-Mars control transfer serves as a baseline: the total accumulated lag over the 200-day transfer is merely 0.21 seconds. Considering the spacecraft’s orbital velocity (∼30 km/s), this temporal discrepancy translates to an along-track position error of approximately 6 km. While this validates the physical consistency of the model in weak fields, it highlights a fundamental scale difference: in the Solar System, neglecting time dilation results in a kilometers-level error correctable via trajectory tuning, whereas in the Galactic Center, it results in a macroscopic failure of the entire navigation timeline.

Figure 6. General Relativistic Time Dilation Accumulation. Top: Deep Gravity regime, highlighting the sharp acceleration of the drift at pericenter; Middle: Outer Cluster regime, exhibiting a linear, secular accumulation over centuries; Bottom: Solar System validation (Earth-Mars transfer), showing a negligible lag of 0.21 seconds.

Figure 6. General Relativistic Time Dilation Accumulation. Top: Deep Gravity regime, highlighting the sharp acceleration of the drift at pericenter; Middle: Outer Cluster regime, exhibiting a linear, secular accumulation over centuries; Bottom: Solar System validation (Earth-Mars transfer), showing a negligible lag of 0.21 seconds.

3.4. the Structural Failure of Classical Mechanics

Having characterized the fundamental dynamical roots of the trajectory divergence, this section quantifies the macroscopic operational consequences for mission design. To systematically assess the breakdown of the classical approximation, a global error analysis is performed across the entire temporal search domain.

3.4.1. Relativistic Miss Distance Maps

To evaluate the impact accuracy of the Newtonian Lambert solver, Relativistic Miss Distance Maps were generated. For every node in the departure-duration search grid, the algorithm computes the vector norm of the spatial discrepancy between the expected Newtonian arrival point and the actual relativistic position at the final time $t_f$:

$$ϵ=\left|\right|{\mathbf{r}}_{GR}\left(t_f\right)-{\mathbf{r}}_{Newt}\left(t_f\right)\left|\right|.$$

(7)

In the Deep Gravity Regime (Figure 7a), the error magnitude exhibits a structural and unavoidable failure of classical mechanics. Even in the absolute best-case geometric configuration, the spacecraft misses the target star by a minimum of $1.4·{10}^7$ km. The majority of the solution space presents errors in the order of ${10}^{10}-{10}^{11}$ km. This macroscopic divergence definitively proves that there are no optimal launch windows where Newtonian approximations accidentally cancel out relativistic effects; the classical solution functions strictly as a coarse initialization guess.

In the Outer Cluster Regime (Figure 7b), the miss distance reaches a maximum of ∼9 AU. While this represents a minute relative error (≈$0.1\%$) compared to the colossal orbital scale, confirming that General Relativity acts locally as a perturbative force, the absolute drift still exceeds the target’s Sphere of Influence. The classical approximation successfully guides the spacecraft to the general vicinity, but an autonomous Relativistic Differential Corrector remains necessary for a precision intercept.

Finally, the Solar System validation (Figure 7c) confirms the physical consistency of the model. The maximum integration error drops to negligible values ($\le 1.5·{10}^5$ km), ensuring that 100% of classical trajectories successfully intersect the Martian Sphere of Influence. Furthermore, the error field’s strict 2.13-year periodicity perfectly matches the Earth-Mars synodic cycle, proving that the accumulated drift is highly deterministic, geometry-dependent, and an authentic product of spacetime curvature rather than a numerical artifact.

4. Discussion

The findings presented in this study demonstrate that navigating the extreme gravitational environment of Sagittarius A* requires a fundamental paradigm shift in astrodynamics. While classical Lambert solvers successfully identify the topologically optimal transfer strategies and estimate the general energetic requirements (effectively placing the spacecraft within the correct basin of attraction) their open-loop application leads to macroscopic targeting errors. The minimum Miss Distance of ${10}^7$ km observed in the Deep Gravity regime definitively proves that Newtonian mechanics is not merely an approximation, but a structurally inadequate framework for precision interception in the Galactic Center. In this domain, General Relativity acts as a dominant dynamical driver, and the chaotic accumulation of spacetime curvature effects renders ballistic propagation completely unfeasible.

Overcoming these limitations necessitates several theoretical and operational advancements. First, the field must transition from hybrid estimators to fully Relativistic Two-Point Boundary Value Problem solvers [28]. Unlike classical algorithms relying on Keplerian conics, the exact solution to the forced geodesic equations in curved spacetime would require the inversion of elliptic integrals (e.g., utilizing the Weierstrass ℘-function formalism) [29,30]. Developing such semi-analytical solvers would allow for the determination of transfer arcs without relying on weak-field assumptions or computationally expensive shooting methods.

Secondly, the assumption of static Keplerian ephemerides for the target stars introduces a systematic limit to precision. As demonstrated by the step-wise Schwarzschild precession of S62 and S4714, the target physically rotates away from the predicted intercept point. Future models must integrate the relativistic equations of motion for the destination stars simultaneously with the spacecraft. By evolving the target dynamically within the same spacetime structure, the guidance algorithm can aim for the physical position of the star rather than a flawed geometric prediction.

Furthermore, to achieve the highest astrophysical fidelity, the mathematical architecture must eventually be extended from the static Schwarzschild approximation to the rotating Kerr metric.

Ultimately, the macroscopic divergence of the classical trajectories dictates a shift from pure trajectory design to autonomous navigation. The engineering culmination of this feasibility study highlights that precision interception in the Galactic Center strictly mandates the implementation of autonomous, closed-loop Guidance, Navigation, and Control (GNC) systems, capable of measuring relativistic observables in real-time to continuously compensate for gravitational anomalies.

5. Conclusions

This work establishes a high-fidelity computational framework to quantify the breakdown of classical astrodynamics in the strong-field regime of Sagittarius A*. Through the analysis of optimal transfers within the S-Cluster, it was demonstrated that the specific forces induced by the relativistic potential, the unclosed precession of the orbits, and the severe chronometric decoupling caused by gravitational time dilation systematically destroy the accuracy of Newtonian guidance. While advanced propulsion technologies (such as nuclear fusion or antimatter, detailed in Appendix B) make these maneuvers energetically sustainable, the topological distortion of spacetime demands a definitive departure from classical orbital mechanics. Future deep-space exploration within galactic cores will inextricably rely on the direct integration of General Relativity into real-time flight software, transforming theoretical geodesics into executable mission profiles.