Explorerโ€บResearch Paperโ€บComputer Science
Research PaperResearchia:202610.04001

Astrodynamics Simulation and Orbital Analysis of the LHS 1903 Planetary System

Chirayu Gupta

Abstract

Astrodynamics Simulation and Orbital Analysis of the LHS 1903 Planetary System Researcher: Chirayu Gupta Age: 15 Current level: High-school student Research area: Astrodynamics, exoplanetary science, computational astrophysics Research Question To what extent is the LHS 1903 planetary system dynamically stable over short-term timescales when modeled using Newtonian N-body gravitational simulations? Overview I independently investigated the orbital dynamics of the LHS 1903 planetary system, a four-planet system orbiting an M-type red dwarf approximately 116 light-years from Earth. The system is particularly interesting because of its compact architecture and the different planetary types and orbital distances within the system. I used published planetary parameters from the NASA Exoplanet Archive and implemented them in NASA's General Mission Analysis Tool (GMAT) to investigate the system's orbital behaviour over time. Method I obtained published orbital parameters for LHS 1903b, c, d and e and converted the relevant quantities into formats compatible with GMAT. The simulation incorporated: Semi-major axes Orbital eccentricities Planetary masses Inclination Right ascension of the ascending node Argument of periapsis True anomaly Because GMAT does not allow the creation of an arbitrary star as a central celestial object in the configuration I used, I represented LHS 1903 using the default Sun object and assigned it the relevant properties of the LHS 1903 star. The planets were configured as satellites of this central body, and their trajectories were propagated using GMAT's numerical modelling capabilities. For orbital parameters that are not fully constrained observationally, standard reference values were assumed to construct a consistent initial configuration. Results The system was propagated for approximately 167 simulated years. During the simulation: All four planets maintained their orbital trajectories. No orbital crossings were observed. No obvious instability events occurred. The orbital elements remained largely unchanged over the simulated period. The planets maintained nested orbital paths around the host star. The results therefore indicate that the modeled system remained dynamically stable over the simulated short-term timescale. Limitations The simulation is not intended to establish that LHS 1903 is stable over geological or astronomical timescales. Important limitations include: The simulation duration was only approximately 167 years. Several orbital-orientation parameters were assumed because they are not fully constrained for the planets. The model used point-mass representations. Relativistic effects were not included. Stellar activity was not modeled. The model does not establish long-term stability over millions or billions of years. Therefore, the current result should be interpreted as evidence of stability within the specific simulated configuration and timescale, rather than as a definitive prediction of the system's long-term evolution. Future Research I would like to develop this project further by investigating: Much longer numerical integrations. The influence of different assumed orbital orientations. Possible orbital resonances between planets. The sensitivity of stability to variations in the published orbital parameters. Improved physical modelling of the host star and planetary interactions. Whether the apparent stability persists under substantially different initial conditions. Motivation for Further Research This project began as an independent investigation into whether a real exoplanetary system could be studied computationally using publicly available astronomical data. I am now interested in taking the work beyond an independent simulation and learning how professional researchers approach planetary-system dynamics, numerical modelling, uncertainty, and the interpretation of simulation results. I would particularly value the opportunity to work under the guidance of a researcher in exoplanetary science, planetary-system dynamics, or computational astrophysics.

Submitted: October 4, 2026Subjects: Computer Science; Research Paper

Description / Details

<div><div><span style="font-size: 17.6px;">Astrodynamics Simulation and Orbital Analysis of the LHS 1903 Planetary System</span></div><div><span style="font-size: 17.6px;">Researcher: Chirayu Gupta</span></div><div><span style="font-size: 17.6px;">Age: 15</span></div><div><span style="font-size: 17.6px;">Current level: High-school student</span></div><div><span style="font-size: 17.6px;">Research area: Astrodynamics, exoplanetary science, computational astrophysics</span></div><div><span style="font-size: 17.6px;">Research Question</span></div><div><span style="font-size: 17.6px;">To what extent is the LHS 1903 planetary system dynamically stable over short-term timescales when modeled using Newtonian N-body gravitational simulations?</span></div><div><span style="font-size: 17.6px;">Overview</span></div><div><span style="font-size: 17.6px;">I independently investigated the orbital dynamics of the LHS 1903 planetary system, a four-planet system orbiting an M-type red dwarf approximately 116 light-years from Earth.</span></div><div><span style="font-size: 17.6px;">The system is particularly interesting because of its compact architecture and the different planetary types and orbital distances within the system. I used published planetary parameters from the NASA Exoplanet Archive and implemented them in NASA's General Mission Analysis Tool (GMAT) to investigate the system's orbital behaviour over time.</span></div><div><span style="font-size: 17.6px;">Method</span></div><div><span style="font-size: 17.6px;">I obtained published orbital parameters for LHS 1903b, c, d and e and converted the relevant quantities into formats compatible with GMAT.</span></div><div><span style="font-size: 17.6px;">The simulation incorporated:</span></div><div><span style="font-size: 17.6px;">Semi-major axes</span></div><div><span style="font-size: 17.6px;">Orbital eccentricities</span></div><div><span style="font-size: 17.6px;">Planetary masses</span></div><div><span style="font-size: 17.6px;">Inclination</span></div><div><span style="font-size: 17.6px;">Right ascension of the ascending node</span></div><div><span style="font-size: 17.6px;">Argument of periapsis</span></div><div><span style="font-size: 17.6px;">True anomaly</span></div><div><span style="font-size: 17.6px;">Because GMAT does not allow the creation of an arbitrary star as a central celestial object in the configuration I used, I represented LHS 1903 using the default Sun object and assigned it the relevant properties of the LHS 1903 star.</span></div><div><span style="font-size: 17.6px;">The planets were configured as satellites of this central body, and their trajectories were propagated using GMAT's numerical modelling capabilities.</span></div><div><span style="font-size: 17.6px;">For orbital parameters that are not fully constrained observationally, standard reference values were assumed to construct a consistent initial configuration.</span></div><div><span style="font-size: 17.6px;">Results</span></div><div><span style="font-size: 17.6px;">The system was propagated for approximately 167 simulated years.</span></div><div><span style="font-size: 17.6px;">During the simulation:</span></div><div><span style="font-size: 17.6px;">All four planets maintained their orbital trajectories.</span></div><div><span style="font-size: 17.6px;">No orbital crossings were observed.</span></div><div><span style="font-size: 17.6px;">No obvious instability events occurred.</span></div><div><span style="font-size: 17.6px;">The orbital elements remained largely unchanged over the simulated period.</span></div><div><span style="font-size: 17.6px;">The planets maintained nested orbital paths around the host star.</span></div><div><span style="font-size: 17.6px;">The results therefore indicate that the modeled system remained dynamically stable over the simulated short-term timescale.</span></div><div><span style="font-size: 17.6px;">Limitations</span></div><div><span style="font-size: 17.6px;">The simulation is not intended to establish that LHS 1903 is stable over geological or astronomical timescales.</span></div><div><span style="font-size: 17.6px;">Important limitations include:</span></div><div><span style="font-size: 17.6px;">The simulation duration was only approximately 167 years.</span></div><div><span style="font-size: 17.6px;">Several orbital-orientation parameters were assumed because they are not fully constrained for the planets.</span></div><div><span style="font-size: 17.6px;">The model used point-mass representations.</span></div><div><span style="font-size: 17.6px;">Relativistic effects were not included.</span></div><div><span style="font-size: 17.6px;">Stellar activity was not modeled.</span></div><div><span style="font-size: 17.6px;">The model does not establish long-term stability over millions or billions of years.</span></div><div><span style="font-size: 17.6px;">Therefore, the current result should be interpreted as evidence of stability within the specific simulated configuration and timescale, rather than as a definitive prediction of the system's long-term evolution.</span></div><div><span style="font-size: 17.6px;">Future Research</span></div><div><span style="font-size: 17.6px;">I would like to develop this project further by investigating:</span></div><div><span style="font-size: 17.6px;">Much longer numerical integrations.</span></div><div><span style="font-size: 17.6px;">The influence of different assumed orbital orientations.</span></div><div><span style="font-size: 17.6px;">Possible orbital resonances between planets.</span></div><div><span style="font-size: 17.6px;">The sensitivity of stability to variations in the published orbital parameters.</span></div><div><span style="font-size: 17.6px;">Improved physical modelling of the host star and planetary interactions.</span></div><div><span style="font-size: 17.6px;">Whether the apparent stability persists under substantially different initial conditions.</span></div><div><span style="font-size: 17.6px;">Motivation for Further Research</span></div><div><span style="font-size: 17.6px;">This project began as an independent investigation into whether a real exoplanetary system could be studied computationally using publicly available astronomical data.</span></div><div><span style="font-size: 17.6px;">I am now interested in taking the work beyond an independent simulation and learning how professional researchers approach planetary-system dynamics, numerical modelling, uncertainty, and the interpretation of simulation results.</span></div><div><span style="font-size: 17.6px;">I would particularly value the opportunity to work under the guidance of a researcher in exoplanetary science, planetary-system dynamics, or computational astrophysics.</span></div><div><br></div></div>

Please sign in to join the discussion.

No comments yet. Be the first to share your thoughts!

Access Paper
Submission Info
Date:
Oct 4, 2026
Topic:
Research Paper
Area:
Computer Science
Comments:
0
Bookmark