Code for Nereid as a Regular Satellite of Neptune
Description
Dynamics component:
The files problem_isolation_mass.c and problem_only_nereid.c contain the code that was used to run the two sets of simulations, the first with a disk of satellites with which Triton interacts, and a second with tidal circularization of Triton into a system with only Nereid and Neptune. These were developed based on the tides_spin_kozai example from Lu et al. 2023. In the problem_only_nereid.c file, lines 90 and 127 change the tidal Q (quality factor). To change the tidal Q, the file must be recompiled.
Running these files requires the installation of rebound and reboundx. The versions used are 4.4.6 and 4.4.1 respectively. First, place the two .c files into the rebound/reboundx/examples/tides_spin_kozai/ directory. Then, compile either of the files by editing the filename of the .c file in the makefile. The simulation runs shown in the paper are attached in a zip file. It contains three .txt files with time in the first column, the a/e/i for Triton in the next three, and then a/e/i for each successive satellite. The corresponding .dat files have the initial satellite system configuration.
Spectroscopy component:
Attached is the Jupyter notebook (.ipynb file) containing the routines used to generate plots 1 and 2 as well as perform the band measurements. The spectra can be found on a MAST doi (10.17909/2bj1-h545), and were all reduced using the jwstspec code which can be found at (https://github.com/ianyuwong/jwstspec). The method of the JWST reduction code is described in Wong+ 2024 (https://doi.org/10.3847/PSJ/ad2fc3). The spectrum of Hyperion is digitized from Clark+ 2012 (https://doi.org/10.1016/j.icarus.2012.01.008) and the VIMS spectrum of Phoebe by request of author from Fraser and Brown 2018 (https://doi.org/10.3847/1538-3881/aac213).