-
The DECam Ecliptic Exploration Project (DEEP) VI: first multi-year observations of trans-Neptunian objects
Authors:
Hayden Smotherman,
Pedro H. Bernardinelli,
Stephen K. N. Portillo,
Andrew J. Connolly,
J. Bryce Kalmbach,
Steven Stetzler,
Mario Juric,
Dino Bektesvic,
Zachary Langford,
Fred C. Adams,
William J. Oldroyd,
Matthew J. Holman,
Colin Orion Chandler,
Cesar Fuentes,
David W. Gerdes,
Hsing Wen Lin,
Larissa Markwardt,
Andrew McNeill,
Michael Mommert,
Kevin J. Napier,
Matthew J. Payne,
Darin Ragozzine,
Andrew S. Rivkin,
Hilke Schlichting,
Scott S. Sheppard
, et al. (3 additional authors not shown)
Abstract:
We present the first set of trans-Neptunian objects (TNOs) observed on multiple nights in data taken from the DECam Ecliptic Exploration Project (DEEP). Of these 110 TNOs, 105 do not coincide with previously known TNOs and appear to be new discoveries. Each individual detection for our objects resulted from a digital tracking search at TNO rates of motion, using two to four hour exposure sets, and…
▽ More
We present the first set of trans-Neptunian objects (TNOs) observed on multiple nights in data taken from the DECam Ecliptic Exploration Project (DEEP). Of these 110 TNOs, 105 do not coincide with previously known TNOs and appear to be new discoveries. Each individual detection for our objects resulted from a digital tracking search at TNO rates of motion, using two to four hour exposure sets, and the detections were subsequently linked across multiple observing seasons. This procedure allows us to find objects with magnitudes $m_{VR} \approx 26$. The object discovery processing also included a comprehensive population of objects injected into the images, with a recovery and linking rate of at least $94\%$. The final orbits were obtained using a specialized orbit fitting procedure that accounts for the positional errors derived from the digital tracking procedure. Our results include robust orbits and magnitudes for classical TNOs with absolute magnitudes $H \sim 10$, as well as a dynamically detached object found at 76 au (semi-major axis $a\approx 77 \, \mathrm{au}$). We find a disagreement between our population of classical TNOs and the CFEPS-L7 three component model for the Kuiper belt.
△ Less
Submitted 5 October, 2023;
originally announced October 2023.
-
The DECam Ecliptic Exploration Project (DEEP) III: Survey characterization and simulation methods
Authors:
Pedro H. Bernardinelli,
Hayden Smotherman,
Zachary Langford,
Stephen K. N. Portillo,
Andrew J. Connolly,
J. Bryce Kalmbach,
Steven Stetzler,
Mario Juric,
William J. Oldroyd,
Hsing Wen Lin,
Fred C. Adams,
Colin Orion Chandler,
Cesar Fuentes,
David W. Gerdes,
Matthew J. Holman,
Larissa Markwardt,
Andrew McNeill,
Michael Mommert,
Kevin J. Napier,
Matthew J. Payne,
Darin Ragozzine,
Andrew S. Rivkin,
Hilke Schlichting,
Scott S. Sheppard,
Ryder Strauss
, et al. (2 additional authors not shown)
Abstract:
We present a detailed study of the observational biases of the DECam Ecliptic Exploration Project's (DEEP) B1 data release and survey simulation software that enables direct statistical comparisons between models and our data. We inject a synthetic population of objects into the images, and then subsequently recover them in the same processing as our real detections. This enables us to characteriz…
▽ More
We present a detailed study of the observational biases of the DECam Ecliptic Exploration Project's (DEEP) B1 data release and survey simulation software that enables direct statistical comparisons between models and our data. We inject a synthetic population of objects into the images, and then subsequently recover them in the same processing as our real detections. This enables us to characterize the survey's completeness as a function of apparent magnitudes and on-sky rates of motion. We study the statistically optimal functional form for the magnitude, and develop a methodology that can estimate the magnitude and rate efficiencies for all survey's pointing groups simultaneously. We have determined that our peak completeness is on average 80\% in each pointing group, and our magnitude drops to $25\%$ of this value at $m_{25} = 26.22$. We describe the freely available survey simulation software and its methodology. We conclude by using it to infer that our effective search area for objects at 40 au is $14.8°^2$, and that our lack of dynamically cold distant objects means that there at most $8\times 10^3$ objects with $60 < a < 80$ au and absolute magnitudes $H \leq 8$.
△ Less
Submitted 5 October, 2023;
originally announced October 2023.
-
Ensembles of unified crust and core equations of state in a nuclear-multimessenger astrophysics environment
Authors:
William G. Newton,
Lauren Balliet,
Srdan Budimir,
Gabriel Crocombe,
Brianna Douglas,
Thomas Blake Head,
Zach Langford,
Luis Rivera,
Josh Sanford
Abstract:
We present an ensemble of unified neutron star crust and core equations of state, constructed using an extended Skyrme energy density functional through the crust and outer core, and appended by two piecewise polytropes at higher densities. The equations of state are parameterized by the first three coefficients in the density expansion of the symmetry energy $J,L$ and $K_{\rm sym}$, the moment of…
▽ More
We present an ensemble of unified neutron star crust and core equations of state, constructed using an extended Skyrme energy density functional through the crust and outer core, and appended by two piecewise polytropes at higher densities. The equations of state are parameterized by the first three coefficients in the density expansion of the symmetry energy $J,L$ and $K_{\rm sym}$, the moment of inertia of a 1.338 M$_{\odot}$ star $I_{1.338}$ and the maximum neutron star mass $M_{\rm max}$. We construct an ensemble with uniform priors on all five parameters, and then apply data filters to the ensemble to explore the effect of combining neutron skin data from PREX with astrophysical measurements of radii and tidal deformabilities from NICER and LIGO/VIRGO. Neutron skins are calculated directly using the EDFs. We demonstrate that both the nuclear data and astrophysical data play a role in constraining crust properties such as the mass, thickness and moment of inertia of the crust and the nuclear pasta layers therein, and that astrophysical data better constrains $K_{\rm sym}$ than PREX data.
△ Less
Submitted 9 March, 2022; v1 submitted 22 December, 2021;
originally announced December 2021.
-
A differentiable N-body code for transit timing and dynamical modeling. I. Algorithm and derivatives
Authors:
Eric Agol,
David M. Hernandez,
Zachary Langford
Abstract:
When fitting N-body models to astronomical data - including transit times, radial velocity, and astrometric positions at observed times - the derivatives of the model outputs with respect to the initial conditions can help with model optimization and posterior sampling. Here we describe a general-purpose symplectic integrator for arbitrary orbital architectures, including those with close encounte…
▽ More
When fitting N-body models to astronomical data - including transit times, radial velocity, and astrometric positions at observed times - the derivatives of the model outputs with respect to the initial conditions can help with model optimization and posterior sampling. Here we describe a general-purpose symplectic integrator for arbitrary orbital architectures, including those with close encounters, which we have recast to maintain numerical stability and precision for small step sizes. We compute the derivatives of the N-body coordinates and velocities as a function of time with respect to the initial conditions and masses by propagating the Jacobian along with the N-body integration. For the first time we obtain the derivatives of the transit times with respect to the initial conditions and masses using the chain rule, which is quicker and more accurate than using finite differences or automatic differentiation. We implement this algorithm in an open source package, NbodyGradient.jl, written in the Julia language, which has been used in the optimization and error analysis of transit-timing variations in the TRAPPIST-1 system. We present tests of the accuracy and precision of the code, and show that it compares favorably in speed to other integrators which are written in C.
△ Less
Submitted 3 June, 2021;
originally announced June 2021.
-
Refining the transit timing and photometric analysis of TRAPPIST-1: Masses, radii, densities, dynamics, and ephemerides
Authors:
Eric Agol,
Caroline Dorn,
Simon L. Grimm,
Martin Turbet,
Elsa Ducrot,
Laetitia Delrez,
Michael Gillon,
Brice-Olivier Demory,
Artem Burdanov,
Khalid Barkaoui,
Zouhair Benkhaldoun,
Emeline Bolmont,
Adam Burgasser,
Sean Carey,
Julien de Wit,
Daniel Fabrycky,
Daniel Foreman-Mackey,
Jonas Haldemann,
David M. Hernandez,
James Ingalls,
Emmanuel Jehin,
Zachary Langford,
Jeremy Leconte,
Susan M. Lederer,
Rodrigo Luger
, et al. (10 additional authors not shown)
Abstract:
We have collected transit times for the TRAPPIST-1 system with the Spitzer Space Telescope over four years. We add to these ground-based, HST and K2 transit time measurements, and revisit an N-body dynamical analysis of the seven-planet system using our complete set of times from which we refine the mass ratios of the planets to the star. We next carry out a photodynamical analysis of the Spitzer…
▽ More
We have collected transit times for the TRAPPIST-1 system with the Spitzer Space Telescope over four years. We add to these ground-based, HST and K2 transit time measurements, and revisit an N-body dynamical analysis of the seven-planet system using our complete set of times from which we refine the mass ratios of the planets to the star. We next carry out a photodynamical analysis of the Spitzer light curves to derive the density of the host star and the planet densities. We find that all seven planets' densities may be described with a single rocky mass-radius relation which is depleted in iron relative to Earth, with Fe 21 wt% versus 32 wt% for Earth, and otherwise Earth-like in composition. Alternatively, the planets may have an Earth-like composition, but enhanced in light elements, such as a surface water layer or a core-free structure with oxidized iron in the mantle. We measure planet masses to a precision of 3-5%, equivalent to a radial-velocity (RV) precision of 2.5 cm/sec, or two orders of magnitude more precise than current RV capabilities. We find the eccentricities of the planets are very small; the orbits are extremely coplanar; and the system is stable on 10 Myr timescales. We find evidence of infrequent timing outliers which we cannot explain with an eighth planet; we instead account for the outliers using a robust likelihood function. We forecast JWST timing observations, and speculate on possible implications of the planet densities for the formation, migration and evolution of the planet system.
△ Less
Submitted 14 January, 2021; v1 submitted 2 October, 2020;
originally announced October 2020.