This KenyonBromley-psc1_readme.txt file was generated on 20181214 by Scott Kenyon Links to Publication Field updated. 2021-12-09, SES ------------------- GENERAL INFORMATION ------------------- 1. Title of Dataset: A Pluto-Charon Sonata: The Dynamical Architecture of the Circumbinary Satellite System 2. Author Information Principal Investigator Contact Information Name: Scott J. Kenyon Institution: Smithsonian Astrophysical Observatory Address: 60 Garden Street, Cambridge, MA 02138 USA Email: skenyon@cfa.harvard.edu Associate or Co-investigator Contact Information Name: Benjamin C. Bromley Institution: Department of Physics & Astronomy, University of Utah Address: 201 JFB, Salt Lake City, UT 84112 USA Email: bromley@physics.utah.edu Alternate Contact Information Name: n/a Institution: Address: Email: 3. Date of data collection: 20180515 to 20181214 4. Geographic location of data collection: NASA discover computer (computer simulations) transferred by ftp to Smithsonian Astrophysical Observatory (Cambridge, MA) and to Dept of Physics & Astronomy, University of Utah (Salt Lake City, UT). 5. Information about funding sources that supported the collection of the data: NASA Emerging World Program: grant NNX17AE24G -------------------------- SHARING/ACCESS INFORMATION -------------------------- 1. Licenses/restrictions placed on the data: none 2. Links to publications that cite or use the data: Kenyon, S. J., & Bromley, B. C. (2019). A Pluto–Charon Sonata: The Dynamical Architecture of the Circumbinary Satellite System. The Astronomical Journal, 157(2), 79. https://doi.org/10.3847/1538-3881/aafa72 3. Links to other publicly accessible locations of the data: none 4. Links/relationships to ancillary data sets: none 5. Was data derived from another source? If yes, list source(s): none 6. Recommended citation for the data: ApJ, TBD --------------------- DATA & FILE OVERVIEW --------------------- 1. File list for directory psc1: Computer files for the paper ``A Pluto-Charon Sonata: The Dynamical Architecture of the Circumbinary Satellite System' (Kenyon & Bromley 2019, ApJ, TBD). The basic data files are binary; two C programs are included to read data from the binary files. (a) Binary files pcs1-0001[a-n]: PC binary only, e0 = 1.e-4, N = 20 (0001a) to 150 (0001n) symplectic steps (in steps of 10) pcs1-0002[a-n]: PC binary only, e0 = 1.e-5, N = 20 (0001a) to 150 (0001n) symplectic steps (in steps of 10) pcs1-0003[a-n]: PC binary only, e0 = 1.e-6, N = 20 (0001a) to 150 (0001n) symplectic steps (in steps of 10) pcs1-0004[a-n]: PC binary only, e0 = 1.e-7, N = 20 (0001a) to 150 (0001n) symplectic steps (in steps of 10) pcs1-0011a: PC binary + tracers, e0 = 1.e-5, retrograde orbits, a0 = 1.00-1.50 a_pc pcs1-0011b: PC binary + tracers, e0 = 1.e-5, retrograde orbits, a0 = 1.45-2.10 a_pc pcs1-0011c: PC binary + tracers, e0 = 1.e-5, retrograde orbits, a0 = 1.95-2.65 a_pc pcs1-0011d: PC binary + tracers, e0 = 1.e-5, retrograde orbits, a0 = 2.60-3.25 a_pc pcs1-0012a: PC binary + tracers, e0 = 1.e-5, prograde orbits, a0 = 1.00-1.50 a_pc pcs1-0012b: PC binary + tracers, e0 = 1.e-5, prograde orbits, a0 = 1.45-2.10 a_pc pcs1-0012c: PC binary + tracers, e0 = 1.e-5, prograde orbits, a0 = 1.95-2.65 a_pc pcs1-0012d: PC binary + tracers, e0 = 1.e-5, prograde orbits, a0 = 2.60-3.25 a_pc pcs1-0013a: PC binary + tracers, e0 = 1.e-5, polar orbits, a0 = 1.00-1.50 a_pc pcs1-0013b: PC binary + tracers, e0 = 1.e-5, polar orbits, a0 = 1.45-2.10 a_pc pcs1-0013c: PC binary + tracers, e0 = 1.e-5, polar orbits, a0 = 1.95-2.65 a_pc pcs1-0013d: PC binary + tracers, e0 = 1.e-5, polar orbits, a0 = 2.60-3.25 a_pc pcs1-0021a: PC binary + SNKH + tracers, e0 = 3.e-3, prograde orbits, a0 = 0.975 a_s - 1.025 a_n pcs1-0021b: PC binary + SNKH + tracers, e0 = 3.e-3, prograde orbits, a0 = 0.975 a_n - 1.025 a_k pcs1-0021c: PC binary + SNKH + tracers, e0 = 3.e-3, prograde orbits, a0 = 0.975 a_k - 1.025 a_h pcs1-0021d: PC binary + SNKH + tracers, e0 = 3.e-3, prograde orbits, a0 = 0.975-1.125 a_h pcs1-0021e: PC binary + SNKH + tracers, e0 = 3.e-3, prograde orbits, a0 = 0.700-1.025 a_s pcs1-0022a: PC binary + SNKH + tracers, e0 = 3.e-3, polar orbits, a0 = 0.975 a_s - 1.025 a_n pcs1-0022b: PC binary + SNKH + tracers, e0 = 3.e-3, polar orbits, a0 = 0.975 a_n - 1.025 a_k pcs1-0022c: PC binary + SNKH + tracers, e0 = 3.e-3, polar orbits, a0 = 0.975 a_k - 1.025 a_h pcs1-0022d: PC binary + SNKH + tracers, e0 = 3.e-3, polar orbits, a0 = 0.975-1.125 a_h pcs1-0022e: PC binary + SNKH + tracers, e0 = 3.e-3, polar orbits, a0 = 0.700-1.025 a_s pcs1-01nn: PC binary + SNKH + 2km sat, prograde orbits, a0 < a_s pcs1-02nn: PC binary + SNKH + 2km sat, prograde orbits, a0 ~= a_n pcs1-03nn: PC binary + SNKH + 2km sat, prograde orbits, a0 ~= a_h pcs1-04nn: PC binary + SNKH + 2km sat, prograde orbits, a0 > a_h pcs1-05nn: PC binary + SNKH + 2km sat, polar orbits, a0 = a_n - a_k pcs1-06nn: PC binary + SNKH + 2km sat, polar orbits, a0 ~= a_k - a_h pcs1-07nn: PC binary + SNKH + 2km sat, polar orbits, a0 > a_h (b) C programs summary.c generates basic summary of timesteps in a binary file usage example: "summary pcs1-0013d" extr6d.c extracts (a,x,y,x,vx,vy,vz,e,i) for a set of time steps from a binary file usage example: "extr6d pcs1-0013d 1068 10" will extract every 10th time step makes files psc1-0013d.0 pcs1-0013d.10 ... pcs1-0013d.1060 pcs1-0013d.1068 usage example: "extr6d pcs1-0013d 100 1" will extract the first 101 time steps 2. Relationship between files: each file contains results for an independent calculation 3. Additional related data collected that was not included in the current data package: none 4. Are there multiple versions of the dataset? no -------------------------- METHODOLOGICAL INFORMATION -------------------------- 1. Description of methods used for collection/generation of data: The data were generated by computer simulations using the C++ code "Orchestra", a proprietary hybrid code that follows the dynamical evolution of solids and gas orbiting a central object. Algorithms in the code are described in the following papers (author names abbreviated to B for Bromley, K for Kenyon, and L for Jane X Luu along with a year for publication date, AJ = Astronomical Journal, ApJ = Astrophysical Journal, S=Supplement): KL1998, AJ 115:2136; KL1999, AJ 118:1101; KB2001,AJ 121:538; KB2002,AJ 123:1757; KB2004, AJ 127:513; BK2006, AJ 131:2737; KB2006, AJ 131:1837; KB2008, ApJS 179:451; KB2010, ApJS 188:242; BK2011, ApJ 731:101; KB2012, AJ 143:63; KB2014, AJ 147:8. Initial conditions for these simulations described in the published paper. 2. Methods for processing the data: Various C and fortran programs are used to analyze the data for the calculations. Several C programs needed to extract information from the computer generated binary output files are included with the dataset. The C programs include basic summaries of the structure of the data files and the usage to extract data from each binary file. 3. Instrument- or software-specific information needed to interpret the data: Appropriate software is included in directory. 4. Standards and calibration information, if appropriate: none 5. Environmental/experimental conditions: all calculations were run on the NASA discover cluster 6. Describe any quality-assurance procedures performed on the data: Aside from tests summarized in the papers described in item 1, test calculations are summarized in the Appendix of each paper and compared to an appropriate benchmark. 7. People involved with sample collection, processing, analysis and/or submission: Scott Kenyon and Ben Bromley ------------------------------------------------------------------ DATA-SPECIFIC INFORMATION FOR: binary data files ------------------------------------------------------------------ 1. Binary filenames have the format "pcs1-nnnn[a-z]" where nnnn is a numerical sequence with an appended letter a-z. Each binary file contains the full set of records for a complete calculation. The relationship between the "nnnn" in the filename and the parameters for each calculation in the journal paper is summarized in Data & File Overview above. 2. Binary files have a sequence of records, where one record corresponds to one time step in a complete calculation. The records contain the time (t), the time step number (it), the number of particles with mass (Nb), the number of bytes allocated to the physical information for each massive particle (= 16 * Nb), the number of massless particles (Nt), and the number of bytes allocated for physical information for each massless particle (=16 *Nt). Each record consists of (i) a short header (which contains the double precision variable t and the five integers it, Nb, 16*Nb, Nt, and 16*Nt), (ii) a set of Nb integers pointing to the radial annulus of each massive object, (iii) a set of Nb integers pointing to the mass bin within an annulus for each massive object, (iv) a set of 16 * Nb double precision numbers for the physical parameters of each massive object, (v) a set of Nt integers pointing to the radial annulus of each massless object, (vi) a set of Nt integers pointing to the mass bin within an annulus for each massless object, and (vii) a set of 16 * Nb double precision numbers for the physical parameters of each massless object. This sequence repeats nt times, where nt is the number of recorded time steps. The included C program "summary.c" extracts basic information for each time step from a binary file and reports this information to the user. The included C program "extr6d.c" extracts position and velocity information for each object in each record of a binary file and creates a set of ascii files (one for each time step) with this information. 3. Physical variables for massive and massless objects: each record for a massive or massless object consists of 16 double precision numbers with indices 0 to 15 inclusive. index 0: M = total mass (zero for massless particles) index 1: Mg = mass in gas index 2: M*rho = product of total mass and mass density in solids index 3: R = physical radius index 4: L = luminosity index 5: A = semimajor axis of orbit around central object index 6: X = X coordinate index 7: Y = Y coordinate index 8: Z = Z coordinate index 9: Vx = X velocity index 10: Vy = Y velocity index 11: Vz = Z velocity index 12: A = semimajor axis (also used for dA/dt in main program) index 13: E = eccentricity index 14: I = eccentricity index 15: M = mass (also used for dM/dt in main program) 4. C programs used to extract data from binary files For calculations with massive particles but no massless particles, programs extract physical information on the massive particles for all time steps in a single calculation. From the supplied routines, it is possible to construct code to return information on other physical variables for each particle. (i) summary.c returns an ascii file with t, it, Nb, Nt and other basic information for all time steps in a single calculation (ii) extr6d.c returns an ascii file with (x,y,z) and (vx,vy,vz) for each particle