Code for Nereid as a Regular Satellite of Neptune

Description

Dynamics component:

The files problem_isolation_mass.c and simulation_zenodo.c contain the code that was used to run two sets of simulations on the tidal evolution of Triton into a disk of satellites around Neptune. The first (problem_isolation_mass.c, and the simulation_runs.zip) set of simulations were pre-review of the paper,  and did not include the sun and used a smaller total disk than the mass of the Uranian satellites. The second set of runs, with simulation_zenodo.c and nereid_uranus_like_simulations.zip contain the final versions shown in the paper. These simulations were developed loosely using the tides_spin_kozai example from Lu et al. 2023.

Running these two .c 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 Neptune and Triton respectively in the next six (or, just the first , and then triplets of a/e/i for each successive satellite. The corresponding .dat files have the initial satellite system configuration, and for the second set of runs, the collisional history, diagram of the system configuration, as well as a .bin file that can be used to re-start the sim.

The Jupyter notebook also contains a cell explaining how to plot the dynamical simulations.

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 JWST 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 or permanently archived at https://zenodo.org/records/17394194). The method of the JWST reduction code is described in Wong+ 2024  (https://doi.org/10.3847/PSJ/ad2fc3). Alternatively, the spectra can be obtained from digitization of the data shown in https://www.nature.com/articles/s41550-024-02433-2/ and https://iopscience.iop.org/article/10.3847/PSJ/ade807. 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), or digitized from the same work. The nereid spectrum itself is available as nereid_ian.dat.

Files