-
Processing, evaluating and understanding FMRI data with afni_proc.py
Authors:
Richard C. Reynolds,
Daniel R. Glen,
Gang Chen,
Ziad S. Saad,
Robert W. Cox,
Paul A. Taylor
Abstract:
FMRI data are noisy, complicated to acquire, and typically go through many steps of processing before they are used in a study or clinical practice. Being able to visualize and understand the data from the start through the completion of processing, while being confident that each intermediate step was successful, is challenging. AFNI's afni_proc$.$py is a tool to create and run a processing pipel…
▽ More
FMRI data are noisy, complicated to acquire, and typically go through many steps of processing before they are used in a study or clinical practice. Being able to visualize and understand the data from the start through the completion of processing, while being confident that each intermediate step was successful, is challenging. AFNI's afni_proc$.$py is a tool to create and run a processing pipeline for FMRI data. With its flexible features, afni_proc$.$py allows users to both control and evaluate their processing at a detailed level. It has been designed to keep users informed about all processing steps: it does not just process the data, but first outputs a fully commented processing script that the users can read, query, interpret and refer back to. Having this full provenance is important for being able to understand each step of processing; it also promotes transparency and reproducibility by keeping the record of individual-level processing and modeling specifics in a single, shareable place. Additionally, afni_proc$.$py creates pipelines that contain several automatic self-checks for potential problems during runtime. The output directory contains a dictionary of relevant quantities that can be programmatically queried for potential issues and a systematic, interactive quality control (QC) HTML. All of these features help users evaluate and understand their data and processing in detail. We describe these and other aspects of afni_proc$.$py here using a set of task-based and resting state FMRI example commands.
△ Less
Submitted 22 August, 2024; v1 submitted 7 June, 2024;
originally announced June 2024.
-
The status and challenges for prostate SBRT treatments in United States proton therapy centers: An NRG Oncology practice survey
Authors:
Jiajian Shen,
Paige A. Taylor,
Carlos E. Vargas,
Minglei Kang,
Jatinder Saini,
Jun Zhou,
Peilong Wang,
Wei Liu,
Charles B. Simone II,
Ying Xiao,
Liyong Lin
Abstract:
A survey was designed to inquire about the practice of proton SBRT treatment for prostate cancer. The survey was distributed to all 30 proton therapy centers in the United States that participate in the National Clinical Trial Network in Feb. 2023. The survey focused on usage, patient selection criteria, prescriptions, target contours, dose constraints, treatment plan optimization and evaluation m…
▽ More
A survey was designed to inquire about the practice of proton SBRT treatment for prostate cancer. The survey was distributed to all 30 proton therapy centers in the United States that participate in the National Clinical Trial Network in Feb. 2023. The survey focused on usage, patient selection criteria, prescriptions, target contours, dose constraints, treatment plan optimization and evaluation methods, patient-specific QA, and IGRT methods. Results: We received responses from 25 centers (83% participation). Only 8 respondent proton centers (32%) reported performing SBRT of the prostate. The remaining 17 centers cited three primary reasons for not offering this treatment: no clinical need, lack of volumetric imaging, and/or lack of clinical evidence. Only 1 center cited the reduction in overall reimbursement as a concern for not offering prostate SBRT. Several common practices among the 8 centers offering SBRT for the prostate were noted, such as using Hydrogel spacers, fiducial markers, and MRI for target delineation. Most proton centers (87.5%) utilized pencil beam scanning (PBS) delivery and completed Imaging and Radiation Oncology Core (IROC) phantom credentialing. Treatment planning typically used parallel opposed lateral beams, and consistent parameters for setup and range uncertainties were used for plan optimization and robustness evaluation. Measurements-based patient-specific QA, beam delivery every other day, fiducial contours for IGRT, and total doses of 35-40 GyRBE were consistent across all centers. However, there was no consensus on the risk levels for patient selection. Conclusion: Prostate SBRT is used in about 1/3 of proton centers in the US. There was a significant consistency in practices among proton centers treating with proton SBRT. It is possible that the adoption of proton SBRT may become more common if proton SBRT is more commonly offered in clinical trials.
△ Less
Submitted 27 February, 2024;
originally announced February 2024.
-
Proton Pencil-Beam Scanning Stereotactic Body Radiation Therapy and Hypofractionated Radiation Therapy for Thoracic Malignancies: Patterns of Practice Survey and Recommendations for Future Development from NRG Oncology and PTCOG
Authors:
Wei Liu,
Hongying Feng,
Paige A. Taylor,
Minglei Kang,
Jiajian Shen,
Jatinder Saini,
Jun Zhou,
Huan B. Giap,
Nathan Y. Yu,
Terence S. Sio,
Pranshu Mohindra,
Joe Y. Chang,
Jeffrey D. Bradley,
Ying Xiao,
Charles B. Simone II,
Liyong Lin
Abstract:
Stereotactic body radiation therapy (SBRT) and hypofractionation using pencil-beam scanning (PBS) proton therapy (PBSPT) is an attractive option for thoracic malignancies. Combining the advantages of target coverage conformity and critical organ sparing from both PBSPT and SBRT, this new delivery technique has great potential to improve the therapeutic ratio, particularly for tumors near critical…
▽ More
Stereotactic body radiation therapy (SBRT) and hypofractionation using pencil-beam scanning (PBS) proton therapy (PBSPT) is an attractive option for thoracic malignancies. Combining the advantages of target coverage conformity and critical organ sparing from both PBSPT and SBRT, this new delivery technique has great potential to improve the therapeutic ratio, particularly for tumors near critical organs. Safe and effective implementation of PBSPT SBRT/hypofractionation to treat thoracic malignancies is more challenging than the conventionally-fractionated PBSPT due to concerns of amplified uncertainties at the larger dose per fraction. NRG Oncology and Particle Therapy Cooperative Group (PTCOG) Thoracic Subcommittee surveyed US proton centers to identify practice patterns of thoracic PBSPT SBRT/hypofractionation. From these patterns, we present recommendations for future technical development of proton SBRT/hypofractionation for thoracic treatment. Amongst other points, the recommendations highlight the need for volumetric image guidance and multiple CT-based robust optimization and robustness tools to minimize further the impact of uncertainties associated with respiratory motion. Advances in direct motion analysis techniques are urgently needed to supplement current motion management techniques.
△ Less
Submitted 1 February, 2024;
originally announced February 2024.
-
False Discovery Rate and Localizing Power
Authors:
Anderson M. Winkler,
Paul A. Taylor,
Thomas E. Nichols,
Chris Rorden
Abstract:
False discovery rate (FDR) is commonly used for correction for multiple testing in neuroimaging studies. However, when using two-tailed tests, making directional inferences about the results can lead to vastly inflated error rate, even approaching 100\% in some cases. This happens because FDR only provides weak control over the error rate, meaning that the proportion of error is guaranteed only gl…
▽ More
False discovery rate (FDR) is commonly used for correction for multiple testing in neuroimaging studies. However, when using two-tailed tests, making directional inferences about the results can lead to vastly inflated error rate, even approaching 100\% in some cases. This happens because FDR only provides weak control over the error rate, meaning that the proportion of error is guaranteed only globally over all tests, not within subsets, such as among those in only one or another direction. Here we consider and evaluate different strategies for FDR control with two-tailed tests, using both synthetic and real imaging data. Approaches that separate the tests by direction of the hypothesis test, or by the direction of the resulting test statistic, more properly control the directional error rate and preserve FDR benefits, albeit with a doubled risk of errors under complete absence of signal. Strategies that combine tests in both directions, or that use simple two-tailed p-values, can lead to invalid directional conclusions, even if these tests remain globally valid. To enable valid thresholding for directional inference, we suggest that imaging software should allow the possibility that the user sets asymmetrical thresholds for the two sides of the statistical map. While FDR continues to be a valid, powerful procedure for multiple testing correction, care is needed when making directional inferences for two-tailed tests, or more broadly, when making any localized inference.
△ Less
Submitted 7 January, 2024;
originally announced January 2024.
-
Orbital Period Change of Dimorphos Due to the DART Kinetic Impact
Authors:
Cristina A. Thomas,
Shantanu P. Naidu,
Peter Scheirich,
Nicholas A. Moskovitz,
Petr Pravec,
Steven R. Chesley,
Andrew S. Rivkin,
David J. Osip,
Tim A. Lister,
Lance A. M. Benner,
Marina Brozović,
Carlos Contreras,
Nidia Morrell,
Agata Rożek,
Peter Kušnirák,
Kamil Hornoch,
Declan Mages,
Patrick A. Taylor,
Andrew D. Seymour,
Colin Snodgrass,
Uffe G. Jørgensen,
Martin Dominik,
Brian Skiff,
Tom Polakis,
Matthew M. Knight
, et al. (24 additional authors not shown)
Abstract:
The Double Asteroid Redirection Test (DART) spacecraft successfully performed the first test of a kinetic impactor for asteroid deflection by impacting Dimorphos, the secondary of near-Earth binary asteroid (65803) Didymos, and changing the orbital period of Dimorphos. A change in orbital period of approximately 7 minutes was expected if the incident momentum from the DART spacecraft was directly…
▽ More
The Double Asteroid Redirection Test (DART) spacecraft successfully performed the first test of a kinetic impactor for asteroid deflection by impacting Dimorphos, the secondary of near-Earth binary asteroid (65803) Didymos, and changing the orbital period of Dimorphos. A change in orbital period of approximately 7 minutes was expected if the incident momentum from the DART spacecraft was directly transferred to the asteroid target in a perfectly inelastic collision, but studies of the probable impact conditions and asteroid properties indicated that a considerable momentum enhancement ($β$) was possible. In the years prior to impact, we used lightcurve observations to accurately determine the pre-impact orbit parameters of Dimorphos with respect to Didymos. Here we report the change in the orbital period of Dimorphos as a result of the DART kinetic impact to be -33.0 +/- 1.0 (3$σ$) minutes. Using new Earth-based lightcurve and radar observations, two independent approaches determined identical values for the change in the orbital period. This large orbit period change suggests that ejecta contributed a significant amount of momentum to the asteroid beyond what the DART spacecraft carried.
△ Less
Submitted 3 March, 2023;
originally announced March 2023.
-
Constraining the Limitations of NEATM-like Models: A Case Study with Near-Earth Asteroid (285263) 1998 QE2
Authors:
Samuel A. Myers,
Ellen S. Howell,
Christopher Magri,
Ronald J. Vervack Jr,
Yanga R. Fernández,
Sean E. Marshall,
Patrick A. Taylor
Abstract:
Near-Earth Asteroids (NEAs) are a key testbed for investigations into planet formation, asteroid dynamics, and planetary defense initiatives. These studies rely on understanding NEA sizes, albedo distributions, and regolith properties. Simple thermal models are a commonly used method for determining these properties, however they have inherent limitations due to the simplifying assumptions they ma…
▽ More
Near-Earth Asteroids (NEAs) are a key testbed for investigations into planet formation, asteroid dynamics, and planetary defense initiatives. These studies rely on understanding NEA sizes, albedo distributions, and regolith properties. Simple thermal models are a commonly used method for determining these properties, however they have inherent limitations due to the simplifying assumptions they make about asteroid shapes and properties. With the recent collapse of the Arecibo Telescope and a decrease of direct size measurements, as well as future facilities such as LSST and NEO Surveyor coming online soon, these models will play an increasingly important role in our knowledge of the NEA population. Therefore, it is key to understand the limits of these models. In this work we constrain the limitations of simple thermal models by comparing model results to more complex thermophysical models, radar data, and other existing analyses. Furthermore, we present a method for placing tighter constraints on inferred NEA properties using simple thermal models. These comparisons and constraints are explored using the NEA (285263) 1998 QE2 as a case study. We analyze QE2 with a simple thermal model and data from both the NASA IRTF SpeX instrument and NEOWISE mission. We determine an albedo between 0.05 and 0.10 and thermal inertia between 0 and 425 J m$^{-2}$ s$^{-1/2}$ K$^{-1}$. We find that overall the simple thermal model is able to well constrain the properties of QE2, however we find that model uncertainties can be influenced by topography, viewing geometry, and the wavelength range of data used.
△ Less
Submitted 12 January, 2023;
originally announced January 2023.
-
Physical properties of near-Earth asteroid (2102) Tantalus from multi-wavelength observations
Authors:
Agata Rożek,
Stephen C. Lowry,
Benjamin Rozitis,
Lord R. Dover,
Patrick A. Taylor,
Anne Virkki,
Simon F. Green,
Colin Snodgrass,
Alan Fitzsimmons,
Justyn Campbell-White,
Sedighe Sajadian,
Valerio Bozza,
Martin J. Burgdorf,
Martin Dominik,
R. Figuera Jaimes,
Tobias C. Hinse,
Markus Hundertmark,
Uffe G. Jørgensen,
Penélope Longa-Peña,
Markus Rabus,
Sohrab Rahvar,
Jesper Skottfelt,
John Southworth
Abstract:
Between 2010 and 2017 we have collected new optical and radar observations of the potentially hazardous asteroid (2102)~Tantalus from the ESO NTT and Danish telescopes at the La Silla Observatory and from the Arecibo planetary radar. The object appears to be nearly spherical, showing a low amplitude light-curve variation and limited large-scale features in the radar images. The spin-state is diffi…
▽ More
Between 2010 and 2017 we have collected new optical and radar observations of the potentially hazardous asteroid (2102)~Tantalus from the ESO NTT and Danish telescopes at the La Silla Observatory and from the Arecibo planetary radar. The object appears to be nearly spherical, showing a low amplitude light-curve variation and limited large-scale features in the radar images. The spin-state is difficult to constrain with the available data; including a certain light-curve subset significantly changes the spin-state estimates, and the uncertainties on period determination are significant. Constraining any change in rotation rate was not possible, despite decades of observations. The convex lightcurve-inversion model, with rotational pole at $λ=210{\pm}41$° and $β=-30{\pm}35$°, is more flattened than the two models reconstructed by including radar observations: with prograde ($λ=36{\pm}23$°, $β=30{\pm}15$°), and with retrograde rotation mode ($λ=180{\pm}24$°, $β=-30{\pm}16$°). Using data from WISE we were able to determine that the prograde model produces the best agreement in size determination between radar and thermophysical modelling. Radar measurements indicate possible variation in surface properties, suggesting one side might have lower radar albedo and be rougher at centimetre-to-decimetre scale than the other. However, further observations are needed to confirm this. Thermophysical analysis indicates a surface covered in fine-grained regolith, consistent with radar albedo and polarisation ratio measurements. Finally, geophysical investigation of the spin-stability of Tantalus shows that it could be exceeding its critical spin-rate via cohesive forces.
△ Less
Submitted 28 June, 2022;
originally announced June 2022.
-
Asteroid 16 Psyche: Shape, Features, and Global Map
Authors:
Michael K. Shepard,
Katherine de Kleer,
Saverio Cambioni,
Patrick A. Taylor,
Anne K. Virkki,
Edgard G. Rivera-Valentin,
Carolina Rodriguez Sanchez-Vahamonde,
Luisa Fernanda Zambrano-Marin,
Christopher Magri,
David Dunham,
John Moore,
Maria Camarca
Abstract:
We develop a shape model of asteroid 16 Psyche using observations acquired in a wide range of wavelengths: Arecibo S-band delay-Doppler imaging, Atacama Large Millimeter Array (ALMA) plane-of-sky imaging, adaptive optics (AO) images from Keck and the Very Large Telescope (VLT), and a recent stellar occultation. Our shape model has dimensions 278 (-4/+8) km x 238 (-4/+6) km x 171 (-1/+5) km, an eff…
▽ More
We develop a shape model of asteroid 16 Psyche using observations acquired in a wide range of wavelengths: Arecibo S-band delay-Doppler imaging, Atacama Large Millimeter Array (ALMA) plane-of-sky imaging, adaptive optics (AO) images from Keck and the Very Large Telescope (VLT), and a recent stellar occultation. Our shape model has dimensions 278 (-4/+8) km x 238 (-4/+6) km x 171 (-1/+5) km, an effective spherical diameter Deff = 222 -1/+4 km, and a spin axis (ecliptic lon, lat) of (36 deg, -8 deg) +/- 2 deg. We survey all the features previously reported to exist, tentatively identify several new features, and produce a global map of Psyche. Using 30 calibrated radar echoes, we find Psyche's overall radar albedo to be 0.34 +/- 0.08 suggesting that the upper meter of regolith has a significant metal (i.e., Fe-Ni) content. We find four regions of enhanced or complex radar albedo, one of which correlates well with a previously identified feature on Psyche, and all of which appear to correlate with patches of relatively high optical albedo. Based on these findings, we cannot rule out a model of Psyche as a remnant core, but our preferred interpretation is that Psyche is a differentiated world with a regolith composition analogous to enstatite or CH/CB chondrites and peppered with localized regions of high metal concentrations. The most credible formation mechanism for these regions is ferrovolcanism as proposed by Johnson et al. (Nature Astronomy vol 4, January 2020, 41-44).
△ Less
Submitted 7 October, 2021;
originally announced October 2021.
-
Astro2020 Activities and Projects White Paper: Arecibo Observatory in the Next Decade
Authors:
D. Anish Roshi,
L. D. Anderson,
E. Araya,
D. Balser,
W. Brisken,
C. Brum,
D. Campbell,
S. Chatterjee,
E. Churchwell,
J. Condon,
J. Cordes,
F. Cordova,
Y. Fernandez,
J. Gago,
T. Ghosh,
P. F. Goldsmith,
C. Heiles,
D. Hickson,
B. Jeffs,
K. M. Jones,
J. Lautenbach,
B. M. Lewis,
R. S. Lynch,
P. K. Manoharan,
S. Marshall
, et al. (20 additional authors not shown)
Abstract:
The white paper discusses Arecibo Observatory's plan for facility improvements and activities over the next decade. The facility improvements include: (a) improving the telescope surface, pointing and focusing to achieve superb performance up to ~12.5 GHz; (b) equip the telescope with ultrawide-band feeds; (c) upgrade the instrumentation with a 4 GHz bandwidth high dynamic range digital link and a…
▽ More
The white paper discusses Arecibo Observatory's plan for facility improvements and activities over the next decade. The facility improvements include: (a) improving the telescope surface, pointing and focusing to achieve superb performance up to ~12.5 GHz; (b) equip the telescope with ultrawide-band feeds; (c) upgrade the instrumentation with a 4 GHz bandwidth high dynamic range digital link and a universal backend and (d) augment the VLBI facility by integrating the 12m telescope for phase referencing. These upgrades to the Arecibo telescope are critical to keep the national facility in the forefront of research in radio astronomy while maintaining its dominance in radar studies of near-Earth asteroids, planets and satellites. In the next decade, the Arecibo telescope will play a synergistic role with the upcoming facilities such as ngVLA, SKA and the now commissioned FAST telescope. Further, the observatory will be actively engaged in mentoring and training programs for students from a diverse background.
△ Less
Submitted 13 July, 2019;
originally announced July 2019.
-
Local Limit Theorems for the Random Conductance Model and Applications to the Ginzburg-Landau $\nablaφ$ Interface Model
Authors:
Sebastian Andres,
Peter A. Taylor
Abstract:
We study a continuous-time random walk on $\mathbb{Z}^d$ in an environment of random conductances taking values in $(0,\infty)$. For a static environment, we extend the quenched local limit theorem to the case of a general speed measure, given suitable ergodicity and moment conditions on the conductances and on the speed measure. Under stronger moment conditions, an annealed local limit theorem is…
▽ More
We study a continuous-time random walk on $\mathbb{Z}^d$ in an environment of random conductances taking values in $(0,\infty)$. For a static environment, we extend the quenched local limit theorem to the case of a general speed measure, given suitable ergodicity and moment conditions on the conductances and on the speed measure. Under stronger moment conditions, an annealed local limit theorem is also derived. Furthermore, an annealed local limit theorem is exhibited in the case of time-dependent conductances, under analogous moment and ergodicity assumptions. This dynamic local limit theorem is then applied to prove a scaling limit result for the space-time covariances in the Ginzburg-Landau $\nablaφ$ model. We also show that the associated Gibbs distribution scales to a Gaussian free field. These results apply to convex potentials for which the second derivative may be unbounded.
△ Less
Submitted 4 February, 2021; v1 submitted 11 July, 2019;
originally announced July 2019.
-
Arecibo Radar Observations of Near-Earth Asteroid (3200) Phaethon During the 2017 Apparition
Authors:
Patrick A. Taylor,
Edgard G. Rivera-Valentin,
Lance A. M. Benner,
Sean E. Marshall,
Anne K. Virkki,
Flaviane C. F. Venditti,
Luisa F. Zambrano-Marin,
Sriram S. Bhiravarasu,
Betzaida Aponte-Hernandez,
Carolina Rodriguez Sanchez-Vahamonde,
Jon D. Giorgini
Abstract:
We report Arecibo S-band (2380 MHz; 12.6 cm) radar observations of near-Earth asteroid (3200) Phaethon during the December 2017 apparition when Phaethon passed within 0.07 au of Earth. Radar images with a resolution of 75 m per pixel reveal a roughly spheroidal shape more than 6 km in diameter at the equator with several discernible surface features hundreds of meters in extent. These include a po…
▽ More
We report Arecibo S-band (2380 MHz; 12.6 cm) radar observations of near-Earth asteroid (3200) Phaethon during the December 2017 apparition when Phaethon passed within 0.07 au of Earth. Radar images with a resolution of 75 m per pixel reveal a roughly spheroidal shape more than 6 km in diameter at the equator with several discernible surface features hundreds of meters in extent. These include a possible crater more than 1 km across located below 30$^{\circ}$ latitude and a roughly 600-m radar-dark region near one of the poles. Overall, the radar images of Phaethon are reminiscent of those of (101955) Bennu, target of the OSIRIS-REx mission. As such, the shape of Phaethon is suspected to have an equatorial ridge similar to the top-shaped models of several other radar-observed near-Earth asteroids as well as the optical images of (162173) Ryugu returned by the Hayabusa2 spacecraft. Preliminary analysis of the radar data finds no satellites and gives no indication of a dusty coma at the time of these observations.
△ Less
Submitted 3 February, 2019;
originally announced February 2019.
-
Stability of Spatial Smoothness and Cluster-Size Threshold Estimates in FMRI using AFNI
Authors:
Robert W Cox,
Paul A Taylor
Abstract:
In a recent analysis of FMRI datasets [K Mueller et al, Front Hum Neurosci 11:345], the estimated spatial smoothness parameters and the statistical significance of clusters were found to depend strongly on the resampled voxel size (for the same data, over a range of 1 to 3 mm) in one popular FMRI analysis software package (SPM12). High sensitivity of thresholding results on such an arbitrary param…
▽ More
In a recent analysis of FMRI datasets [K Mueller et al, Front Hum Neurosci 11:345], the estimated spatial smoothness parameters and the statistical significance of clusters were found to depend strongly on the resampled voxel size (for the same data, over a range of 1 to 3 mm) in one popular FMRI analysis software package (SPM12). High sensitivity of thresholding results on such an arbitrary parameter as final spatial grid size is an undesirable feature in a processing pipeline. Here, we examine the stability of spatial smoothness and cluster-volume threshold estimates with respect to voxel resampling size in the AFNI software package's pipeline. A publicly available collection of resting-state and task FMRI datasets from 78 subjects was analyzed using standard processing steps in AFNI. We found that the spatial smoothness and cluster-volume thresholds are fairly stable over the voxel resampling size range of 1 to 3 mm, in contradistinction to the reported results from SPM12.
△ Less
Submitted 21 September, 2017;
originally announced September 2017.
-
Yarkovsky Drift Detections for 247 Near-Earth Asteroids
Authors:
Adam H. Greenberg,
Jean-Luc Margot,
Ashok K. Verma,
Patrick A. Taylor,
Susan E. Hodge
Abstract:
The Yarkovsky effect is a thermal process acting upon the orbits of small celestial bodies, which can cause these orbits to slowly expand or contract with time. The effect is subtle (da/dt ~ 10^-4 au/My for a 1 km diameter object) and is thus generally difficult to measure. We analyzed both optical and radar astrometry for 600 near-Earth asteroids (NEAs) for the purpose of detecting and quantifyin…
▽ More
The Yarkovsky effect is a thermal process acting upon the orbits of small celestial bodies, which can cause these orbits to slowly expand or contract with time. The effect is subtle (da/dt ~ 10^-4 au/My for a 1 km diameter object) and is thus generally difficult to measure. We analyzed both optical and radar astrometry for 600 near-Earth asteroids (NEAs) for the purpose of detecting and quantifying the Yarkovsky effect. We present 247 NEAs with measured drift rates, which is the largest published set of Yarkovsky detections. This large sample size provides an opportunity to examine the Yarkovsky effect in a statistical manner. In particular, we describe two independent population-based tests that verify the measurement of Yarkovsky orbital drift. First, we provide observational confirmation for the Yarkovsky effect's theoretical size dependence of 1/D, where D is diameter. Second, we find that the observed ratio of negative to positive drift rates in our sample is 2.34, which, accounting for bias and sampling uncertainty, implies an actual ratio of $2.7^{+0.3}_{-0.7}$. This ratio has a vanishingly small probability of occurring due to chance or statistical noise. The observed ratio of retrograde to prograde rotators is two times lower than the ratio expected from numerical predictions from NEA population studies and traditional assumptions about the sense of rotation of NEAs originating from various main belt escape routes. We also examine the efficiency with which solar energy is converted into orbital energy and find a median efficiency in our sample of 12%. We interpret this efficiency in terms of NEA spin and thermal properties.
△ Less
Submitted 7 February, 2020; v1 submitted 18 August, 2017;
originally announced August 2017.
-
FMRI Clustering and False Positive Rates
Authors:
Robert W. Cox,
Gang Chen,
Daniel R. Glen,
Richard C. Reynolds,
Paul A. Taylor
Abstract:
Recently, Eklund et al. (2016) analyzed clustering methods in standard FMRI packages: AFNI (which we maintain), FSL, and SPM [1]. They claimed: 1) false positive rates (FPRs) in traditional approaches are greatly inflated, questioning the validity of "countless published fMRI studies"; 2) nonparametric methods produce valid, but slightly conservative, FPRs; 3) a common flawed assumption is that th…
▽ More
Recently, Eklund et al. (2016) analyzed clustering methods in standard FMRI packages: AFNI (which we maintain), FSL, and SPM [1]. They claimed: 1) false positive rates (FPRs) in traditional approaches are greatly inflated, questioning the validity of "countless published fMRI studies"; 2) nonparametric methods produce valid, but slightly conservative, FPRs; 3) a common flawed assumption is that the spatial autocorrelation function (ACF) of FMRI noise is Gaussian-shaped; and 4) a 15-year-old bug in AFNI's 3dClustSim significantly contributed to producing "particularly high" FPRs compared to other software. We repeated simulations from [1] (Beijing-Zang data [2], see [3]), and comment on each point briefly.
△ Less
Submitted 15 February, 2017;
originally announced February 2017.
-
FMRI Clustering in AFNI: False Positive Rates Redux
Authors:
Robert W. Cox,
Gang Chen,
Daniel R. Glen,
Richard C. Reynolds,
Paul A. Taylor
Abstract:
Recent reports of inflated false positive rates (FPRs) in FMRI group analysis tools by Eklund et al. (2016) have become a large topic within (and outside) neuroimaging. They concluded that: existing parametric methods for determining statistically significant clusters had greatly inflated FPRs ("up to 70%," mainly due to the faulty assumption that the noise spatial autocorrelation function is Gaus…
▽ More
Recent reports of inflated false positive rates (FPRs) in FMRI group analysis tools by Eklund et al. (2016) have become a large topic within (and outside) neuroimaging. They concluded that: existing parametric methods for determining statistically significant clusters had greatly inflated FPRs ("up to 70%," mainly due to the faulty assumption that the noise spatial autocorrelation function is Gaussian- shaped and stationary), calling into question potentially "countless" previous results; in contrast, nonparametric methods, such as their approach, accurately reflected nominal 5% FPRs. They also stated that AFNI showed "particularly high" FPRs compared to other software, largely due to a bug in 3dClustSim. We comment on these points using their own results and figures and by repeating some of their simulations. Briefly, while parametric methods show some FPR inflation in those tests (and assumptions of Gaussian-shaped spatial smoothness also appear to be generally incorrect), their emphasis on reporting the single worst result from thousands of simulation cases greatly exaggerated the scale of the problem. Importantly, FPR statistics depend on "task" paradigm and voxelwise p-value threshold; as such, we show how results of their study provide useful suggestions for FMRI study design and analysis, rather than simply a catastrophic downgrading of the field's earlier results. Regarding AFNI (which we maintain), 3dClustSim's bug-effect was greatly overstated - their own results show that AFNI results were not "particularly" worse than others. We describe further updates in AFNI for characterizing spatial smoothness more appropriately (greatly reducing FPRs, though some remain >5%); additionally, we outline two newly implemented permutation/randomization-based approaches producing FPRs clustered much more tightly about 5% for voxelwise p<=0.01.
△ Less
Submitted 15 February, 2017;
originally announced February 2017.
-
Asteroid 1566 Icarus's size, shape, orbit, and Yarkovsky drift from radar observations
Authors:
Adam H. Greenberg,
Jean-Luc Margot,
Ashok K. Verma,
Patrick A. Taylor,
Shantanu P. Naidu,
Marina. Brozovic,
Lance A. M. Benner
Abstract:
Near-Earth asteroid (NEA) 1566 Icarus ($a=1.08$ au, $e=0.83$, $i=22.8^{\circ}$) made a close approach to Earth in June 2015 at 22 lunar distances (LD). Its detection during the 1968 approach (16 LD) was the first in the history of asteroid radar astronomy. A subsequent approach in 1996 (40 LD) did not yield radar images. We describe analyses of our 2015 radar observations of Icarus obtained at the…
▽ More
Near-Earth asteroid (NEA) 1566 Icarus ($a=1.08$ au, $e=0.83$, $i=22.8^{\circ}$) made a close approach to Earth in June 2015 at 22 lunar distances (LD). Its detection during the 1968 approach (16 LD) was the first in the history of asteroid radar astronomy. A subsequent approach in 1996 (40 LD) did not yield radar images. We describe analyses of our 2015 radar observations of Icarus obtained at the Arecibo Observatory and the DSS-14 antenna at Goldstone. These data show that the asteroid is a moderately flattened spheroid with an equivalent diameter of 1.44 km with 18% uncertainties, resolving long-standing questions about the asteroid size. We also solve for Icarus' spin axis orientation ($λ=270^{\circ}\pm10^{\circ}, β=-81^{\circ}\pm10^{\circ}$), which is not consistent with the estimates based on the 1968 lightcurve observations. Icarus has a strongly specular scattering behavior, among the highest ever measured in asteroid radar observations, and a radar albedo of $\sim$2\%, among the lowest ever measured in asteroid radar observations. The low cross-section suggests a high-porosity surface, presumably related to Icarus' cratering, spin, and thermal histories. Finally, we present the first use of our orbit determination software for the generation of observational ephemerides, and we demonstrate its ability to determine subtle perturbations on NEA orbits by measuring Icarus' orbit-averaged drift in semi-major axis ($(-4.62\pm0.48) \times 10^{-4}$ au/My, or $\sim$60 m per revolution). Our Yarkovsky rate measurement resolves a discrepancy between two published rates that did not include the 2015 radar astrometry.
△ Less
Submitted 13 January, 2017; v1 submitted 21 December, 2016;
originally announced December 2016.
-
Physical Characterization of ~2-meter Diameter Near-Earth Asteroid 2015 TC25: A possible boulder from E-type Asteroid (44) Nysa
Authors:
Vishnu Reddy,
Juan A. Sanchez,
William F. Bottke,
Audrey Thirouin,
Edgard G. Rivera-Valentin,
Michael S. Kelley,
William Ryan,
Edward A. Cloutis,
Stephen C. Tegler,
Eileen V. Ryan,
Patrick A. Taylor,
James E. Richardson,
Nicholas Moskovitz,
Lucille Le Corre
Abstract:
Small near-Earth asteroids (>20 meters) are interesting because they are progenitors for meteorites in our terrestrial collection. Crucial to our understanding of the effectiveness of our atmosphere in filtering low-strength impactors is the physical characteristics of these small near-Earth asteroids (NEAs). In the past, characterization of small NEAs has been a challenge because of the difficult…
▽ More
Small near-Earth asteroids (>20 meters) are interesting because they are progenitors for meteorites in our terrestrial collection. Crucial to our understanding of the effectiveness of our atmosphere in filtering low-strength impactors is the physical characteristics of these small near-Earth asteroids (NEAs). In the past, characterization of small NEAs has been a challenge because of the difficulty in detecting them prior to close Earth flyby. In this study we physically characterized the 2-meter diameter near-Earth asteroid 2015 TC25 using ground-based optical, near-infrared and radar assets during a close flyby of the Earth (distance 69,000 miles) in Oct. 2015. Our observations suggest that its surface composition is similar to aubrites, a rare class of high albedo differentiated meteorites. Aubrites make up only 0.14 % of all know meteorites in our terrestrial meteorite collection. 2015 TC25 is also a very fast rotator with a rotation period of 133 seconds. We compared spectral and dynamical properties of 2015 TC25 and found the best candidate source body in the inner main belt to be the 70-km diameter E-type asteroid (44) Nysa. We attribute difference in spectral slope between the two objects to the lack of regolith on the surface of 2015 TC25. Using the albedo of E-type asteroids (50-60%) we refine the diameter of 2015 TC25 to 2-meters making it one of the smallest NEA ever to be characterized.
△ Less
Submitted 3 December, 2016; v1 submitted 30 November, 2016;
originally announced December 2016.
-
Thermal Properties and an Improved Shape Model for Near-Earth Asteroid (162421) 2000 ET70
Authors:
Sean E. Marshall,
Ellen S. Howell,
Christopher Magri,
Ronald J. Vervack Jr.,
Donald B. Campbell,
Yanga R. Fernández,
Michael C. Nolan,
Jenna L. Crowell,
Michael D. Hicks,
Kenneth J. Lawrence,
Patrick A. Taylor
Abstract:
We present thermal properties and an improved shape model for potentially hazardous asteroid (162421) 2000 ET70. In addition to the radar data from 2000 ET70's apparition in 2012, our model incorporates optical lightcurves and infrared spectra that were not included in the analysis of Naidu et al. (2013, Icarus 226, 323-335). We confirm the general "clenched fist" appearance of the Naidu et al. mo…
▽ More
We present thermal properties and an improved shape model for potentially hazardous asteroid (162421) 2000 ET70. In addition to the radar data from 2000 ET70's apparition in 2012, our model incorporates optical lightcurves and infrared spectra that were not included in the analysis of Naidu et al. (2013, Icarus 226, 323-335). We confirm the general "clenched fist" appearance of the Naidu et al. model, but compared to their model, our best-fit model is about 10% longer along its long principal axis, nearly identical along the intermediate axis, and about 25% shorter along the short axis. We find the asteroid's dimensions to be 2.9 km $\times$ 2.2 km $\times$ 1.5 km (with relative uncertainties of about 10%, 15%, and 25%, respectively). With the available data, 2000 ET70's period and pole position are degenerate with each other. The radar and lightcurve data together constrain the pole direction to fall along an arc that is about twenty-three degrees long and eight degrees wide. Infrared spectra from the NASA InfraRed Telescope Facility (IRTF) provide an additional constraint on the pole. Thermophysical modeling, using our SHERMAN software, shows that only a subset of the pole directions, about twelve degrees of that arc, are compatible with the infrared data. Using all of the available data, we find that 2000 ET70 has a sidereal rotation period of 8.944 hours ($\pm$ 0.009 h) and a north pole direction of ecliptic coordinates $(52^{\circ}, -60^{\circ}) \pm 6^{\circ}$. The infrared data, acquired over several dates, require that the thermal properties (albedo, thermal inertia, surface roughness) must change across the asteroid's surface. By incorporating the detailed shape model and spin state into our thermal modeling, the multiple ground-based observations at different viewing geometries have allowed us to constrain the levels of the variations in the surface properties.
△ Less
Submitted 10 March, 2017; v1 submitted 14 October, 2016;
originally announced October 2016.
-
Capabilities of Earth-based radar facilities for near-Earth asteroid observations
Authors:
Shantanu. P. Naidu,
Lance. A. M. Benner,
Jean-Luc Margot,
Michael. W. Busch,
Patrick. A. Taylor
Abstract:
We evaluated the planetary radar capabilities at Arecibo, the Goldstone 70-m DSS-14 and 34-m DSS-13 antennas, the 70-m DSS-43 antenna at Canberra, the Green Bank Telescope, and the Parkes Radio Telescope in terms of their relative sensitivities and the number of known near-Earth asteroids (NEAs) detectable per year in monostatic and bistatic configurations. In the 2015 calendar year, monostatic ob…
▽ More
We evaluated the planetary radar capabilities at Arecibo, the Goldstone 70-m DSS-14 and 34-m DSS-13 antennas, the 70-m DSS-43 antenna at Canberra, the Green Bank Telescope, and the Parkes Radio Telescope in terms of their relative sensitivities and the number of known near-Earth asteroids (NEAs) detectable per year in monostatic and bistatic configurations. In the 2015 calendar year, monostatic observations with Arecibo and DSS-14 were capable of detecting 253 and 131 NEAs respectively, with signal-to-noise ratios (SNRs) greater than 30/track. Combined, the two observatories were capable of detecting 276 NEAs. Of these, Arecibo detected 77 and Goldstone detected 32, or 30% and 24% the numbers that were possible. The two observatories detected an additional 18 and 7 NEAs respectively, with SNRs of less than 30/track. This indicates that a substantial number of potential targets are not being observed. The bistatic configuration with DSS-14 transmitting and the Green Bank Telescope receiving was capable of detecting about 195 NEAs, or ~50% more than with monostatic observations at DSS-14. Most of the detectable asteroids were targets of opportunity that were discovered less than 15 days before the end of their observing windows. About 50% of the detectable asteroids have absolute magnitudes > 25, which corresponds diameters < ~30 m.
△ Less
Submitted 27 June, 2016; v1 submitted 4 April, 2016;
originally announced April 2016.
-
Radar Imaging and Characterization of Binary Near-Earth Asteroid (185851) 2000 DP107
Authors:
Shantanu P. Naidu,
Jean-Luc Margot,
Patrick A. Taylor,
Michael C. Nolan,
Michael W. Busch,
Lance A. M. Benner,
Marina Brozovic,
Jon D. Giorgini,
Joseph S. Jao,
Chris Magri
Abstract:
Potentially hazardous asteroid (185851) 2000 DP107 was the first binary near-Earth asteroid to be imaged. Radar observations in 2000 provided images at 75 m resolution that revealed the shape, orbit, and spin-up formation mechanism of the binary. The asteroid made a more favorable flyby of the Earth in 2008, yielding images at 30 m resolution. We used these data to obtain shape models for the two…
▽ More
Potentially hazardous asteroid (185851) 2000 DP107 was the first binary near-Earth asteroid to be imaged. Radar observations in 2000 provided images at 75 m resolution that revealed the shape, orbit, and spin-up formation mechanism of the binary. The asteroid made a more favorable flyby of the Earth in 2008, yielding images at 30 m resolution. We used these data to obtain shape models for the two components and to improve the estimates of the mutual orbit, component masses, and spin periods. The primary has a sidereal spin period of 2.7745 +/- 0.0007 h and is roughly spheroidal with an equivalent diameter of 863 m +/- 5 %. It has a mass of 4.656 +/- 0.43 x 10^11 kg and a density of 1381 +/- 244 kg m^{-3}. It exhibits an equatorial ridge similar to the (66391) 1999 KW4 primary, however the equatorial ridge in this case is not as regular and has a ~300 m diameter concavity on one side. The secondary has a sidereal spin period of 1.77 +/- 0.02 days commensurate with the orbital period. The secondary is slightly elongated and has overall dimensions of 377 x 314 x 268 m (6 % uncertainties). Its mass is 0.178 +/- 0.021 x 10^{11} kg and its density is 1047 +/- 230 kg m^{-3}. The mutual orbit has a semi-major axis of 2.659 +/- 0.08 km, an eccentricity of 0.019 +/- 0.01, and a period of 1.7556 +/- 0.0015 days. The normalized total angular momentum of this system exceeds the amount required for the expected spin-up formation mechanism. An increase of angular momentum from non-gravitational forces after binary formation is a possible explanation. The two components have similar radar reflectivity, suggesting a similar composition consistent with formation by spin-up. The secondary appears to exhibit a larger circular polarization ratio than the primary, suggesting a rougher surface or subsurface at radar wavelength scales.
△ Less
Submitted 29 June, 2015; v1 submitted 5 March, 2015;
originally announced March 2015.
-
Orbit and Bulk Density of the OSIRIS-REx Target Asteroid (101955) Bennu
Authors:
Steven R. Chesley,
Davide Farnocchia,
Michael C. Nolan,
David Vokrouhlicky,
Paul W. Chodas,
Andrea Milani,
Federica Spoto,
Benjamin Rozitis,
Lance A. M. Benner,
William F. Bottke,
Michael W. Busch,
Joshua P. Emery,
Ellen S. Howell,
Dante S. Lauretta,
Jean-Luc Margot,
Patrick A. Taylor
Abstract:
The target asteroid of the OSIRIS-REx asteroid sample return mission, (101955) Bennu (formerly 1999 RQ$_{36}$), is a half-kilometer near-Earth asteroid with an extraordinarily well constrained orbit. An extensive data set of optical astrometry from 1999--2013 and high-quality radar delay measurements to Bennu in 1999, 2005, and 2011 reveal the action of the Yarkovsky effect, with a mean semimajor…
▽ More
The target asteroid of the OSIRIS-REx asteroid sample return mission, (101955) Bennu (formerly 1999 RQ$_{36}$), is a half-kilometer near-Earth asteroid with an extraordinarily well constrained orbit. An extensive data set of optical astrometry from 1999--2013 and high-quality radar delay measurements to Bennu in 1999, 2005, and 2011 reveal the action of the Yarkovsky effect, with a mean semimajor axis drift rate $da/dt = (-19.0 \pm 0.1)\times 10^{-4}$ au/Myr or $284\pm 1.5\;\rm{m/yr}$. The accuracy of this result depends critically on the fidelity of the observational and dynamical model. As an example, neglecting the relativistic perturbations of the Earth during close approaches affects the orbit with $3σ$ significance in $da/dt$.
The orbital deviations from purely gravitational dynamics allow us to deduce the acceleration of the Yarkovsky effect, while the known physical characterization of Bennu allows us to independently model the force due to thermal emissions. The combination of these two analyses yields a bulk density of $ρ= 1260\pm70\,\rm{kg/m^3}$, which indicates a macroporosity in the range $40\pm10$% for the bulk densities of likely analog meteorites, suggesting a rubble-pile internal structure. The associated mass estimate is $(7.8\pm0.9)\times 10^{10}\, \rm{kg}$ and $GM = 5.2\pm0.6\,\rm{m^3/s^2}$.
Bennu's Earth close approaches are deterministic over the interval 1654--2135, beyond which the predictions are statistical in nature. In particular, the 2135 close approach is likely within the lunar distance and leads to strong scattering and therefore numerous potential impacts in subsequent years, from 2175--2196. The highest individual impact probability is $9.5\times 10^{-5}$ in 2196, and the cumulative impact probability is $3.7\times 10^{-4}$, leading to a cumulative Palermo Scale of -1.70.
△ Less
Submitted 22 February, 2014;
originally announced February 2014.
-
Tidal End States of Binary Asteroid Systems with a Nonspherical Component
Authors:
Patrick A. Taylor,
Jean-Luc Margot
Abstract:
We derive the locations of the fully synchronous end states of tidal evolution for binary asteroid systems having one spherical component and one oblate- or prolate-spheroid component. Departures from a spherical shape, at levels observed among binary asteroids, can result in the lack of a stable tidal end state for particular combinations of the system mass fraction and angular momentum, in which…
▽ More
We derive the locations of the fully synchronous end states of tidal evolution for binary asteroid systems having one spherical component and one oblate- or prolate-spheroid component. Departures from a spherical shape, at levels observed among binary asteroids, can result in the lack of a stable tidal end state for particular combinations of the system mass fraction and angular momentum, in which case the binary must collapse to contact. We illustrate our analytical results with near-Earth asteroids (8567) 1996 HW1, (66391) 1999 KW4, and 69230 Hermes.
△ Less
Submitted 4 February, 2014; v1 submitted 2 December, 2013;
originally announced December 2013.
-
Radar Imaging and Physical Characterization of Near-Earth Asteroid (162421) 2000 ET70
Authors:
Shantanu P. Naidu,
Jean-Luc Margot,
Michael W. Busch,
Patrick A. Taylor,
Michael C. Nolan,
Marina Brozovic,
Lance A. M. Benner,
Jon D. Giorgini,
Christopher Magri
Abstract:
We observed near-Earth asteroid (162421) 2000 ET70 using the Arecibo and Goldstone radar systems over a period of 12 days during its close approach to the Earth in February 2012. We obtained continuous wave spectra and range-Doppler images with range resolutions as fine as 15 m. Inversion of the radar images yields a detailed shape model with an effective spatial resolution of 100 m. The asteroid…
▽ More
We observed near-Earth asteroid (162421) 2000 ET70 using the Arecibo and Goldstone radar systems over a period of 12 days during its close approach to the Earth in February 2012. We obtained continuous wave spectra and range-Doppler images with range resolutions as fine as 15 m. Inversion of the radar images yields a detailed shape model with an effective spatial resolution of 100 m. The asteroid has overall dimensions of 2.6 km X 2.2 km X 2.1 km (5% uncertainties) and a surface rich with kilometer-scale ridges and concavities. This size, combined with absolute magnitude measurements, implies an extremely low albedo (~2%). It is a principal axis rotator and spins in a retrograde manner with a sidereal spin period of 8.96 +/- 0.01 hours. In terms of gravitational slopes evaluated at scales of 100 m, the surface seems mostly relaxed with over 99% of the surface having slopes less than 30 degrees, but there are some outcrops at the north pole that may have steeper slopes. Our precise measurements of the range and velocity of the asteroid, combined with optical astrometry, enables reliable trajectory predictions for this potentially hazardous asteroid in the interval 460-2813.
△ Less
Submitted 29 May, 2013; v1 submitted 28 January, 2013;
originally announced January 2013.
-
An analytic toy model for relativistic accretion in Kerr spacetime
Authors:
Emilio Tejeda,
Paul A. Taylor,
John C. Miller
Abstract:
We present a relativistic model for the stationary axisymmetric accretion flow of a rotating cloud of non-interacting particles falling onto a Kerr black hole. Based on a ballistic approximation, streamlines are described analytically in terms of timelike geodesics, while a simple numerical scheme is introduced for calculating the density field. A novel approach is presented for describing all of…
▽ More
We present a relativistic model for the stationary axisymmetric accretion flow of a rotating cloud of non-interacting particles falling onto a Kerr black hole. Based on a ballistic approximation, streamlines are described analytically in terms of timelike geodesics, while a simple numerical scheme is introduced for calculating the density field. A novel approach is presented for describing all of the possible types of orbit by means of a single analytic expression. This model is a useful tool for highlighting purely relativistic signatures in the accretion flow dynamics coming from a strong gravitational field with frame-dragging. In particular, we explore the coupling due to this between the spin of the black hole and the angular momentum of the infalling matter. Moreover, we demonstrate how this analytic solution may be used for benchmarking general relativistic numerical hydrodynamics codes by comparing it against results of smoothed particle hydrodynamics simulations for a collapsar-like setup. These simulations are performed first for a ballistic flow (with zero pressure) and then for a hydrodynamical one where we measure the effects of pressure gradients on the infall, thus exploring the extent of applicability of the ballistic approximation.
△ Less
Submitted 19 December, 2012; v1 submitted 3 July, 2012;
originally announced July 2012.
-
Measuring the Effects of Artificial Viscosity in SPH Simulations of Rotating Fluid Flows
Authors:
Paul A. Taylor,
John C. Miller
Abstract:
A commonly cited drawback of SPH is the introduction of spurious shear viscosity by the artificial viscosity term in situations involving rotation. Existing approaches for quantifying its effect include approximate analytic formulae and disc-averaged be- haviour in specific ring-spreading simulations, based on the kinematic effects produced by the artificial viscosity. These methods have disadvant…
▽ More
A commonly cited drawback of SPH is the introduction of spurious shear viscosity by the artificial viscosity term in situations involving rotation. Existing approaches for quantifying its effect include approximate analytic formulae and disc-averaged be- haviour in specific ring-spreading simulations, based on the kinematic effects produced by the artificial viscosity. These methods have disadvantages, in that they typically are applicable to a very small range of physical scenarios, have a large number of simplifying assumptions, and often are tied to specific SPH formulations which do not include corrective (e.g., Balsara) or time-dependent viscosity terms. In this study we have developed a simple, generally applicable and practical technique for evaluating the local effect of artificial viscosity directly from the creation of specific entropy for each SPH particle. This local approach is simple and quick to implement, and it al- lows a detailed characterization of viscous effects as a function of position. Several advantages of this method are discussed, including its ease in evaluation, its greater accuracy and its broad applicability. In order to compare this new method with ex- isting ones, simple disc flow examples are used. Even in these basic cases, the very roughly approximate nature of the previous methods is shown. Our local method pro- vides a detailed description of the effects of the artificial viscosity throughout the disc, even for extended examples which implement Balsara corrections. As a further use of this approach, explicit dependencies of the effective viscosity in terms of SPH and flow parameters are estimated from the example cases. In an appendix, a method for the initial placement of SPH particles is discussed which is very effective in reducing numerical fluctuations.
△ Less
Submitted 24 June, 2012; v1 submitted 21 December, 2011;
originally announced December 2011.
-
Binary Asteroid Systems: Tidal End States and Estimates of Material Properties
Authors:
Patrick A. Taylor,
Jean-Luc Margot
Abstract:
The locations of the fully despun, double synchronous end states of tidal evolution are derived for spherical components. With the exception of nearly equal-mass binaries, binary asteroid systems are in the midst of lengthy tidal evolutions, far from their fully synchronous tidal end states. Calculations of material strength indicate that binaries in the main belt with 100-km-scale primary compone…
▽ More
The locations of the fully despun, double synchronous end states of tidal evolution are derived for spherical components. With the exception of nearly equal-mass binaries, binary asteroid systems are in the midst of lengthy tidal evolutions, far from their fully synchronous tidal end states. Calculations of material strength indicate that binaries in the main belt with 100-km-scale primary components are consistent with being made of monolithic or fractured rock as expected for binaries likely formed from sub-catastrophic impacts in the early solar system. To tidally evolve in their dynamical lifetime, near-Earth binaries with km-scale primaries or smaller must be much weaker mechanically than their main-belt counterparts even if formed in the main belt prior to injection into the near-Earth region. Small main-belt binaries with primary components less than 10 km in diameter, depending on their ages, could either be as strong as large main-belt binaries or as weak as near-Earth binaries because the inherent uncertainty in the age of a binary system can affect the calculation of material strength by orders of magnitude. Several other issues are considered, though these typically affect the calculation of material strength by no more than a factor of two. We also find indirect evidence within all three groups of binary asteroids that the inter-component separation may evolve via another mechanism(s) with the binary YORP effect being a likely candidate.
△ Less
Submitted 25 March, 2011; v1 submitted 24 January, 2011;
originally announced January 2011.
-
Radar Observations and the Shape of Near-Earth Asteroid 2008 EV5
Authors:
Michael W. Busch,
Steven J. Ostro,
Lance A. M. Benner,
Marina Brozovic,
Jon D. Giorgini,
Joseph S. Jao,
Daniel J. Scheeres,
Christopher Magri,
Michael C. Nolan,
Ellen S. Howell,
Patrick A. Taylor,
Jean-Luc Margot,
Walter Brisken
Abstract:
We observed the near-Earth asteroid 2008 EV5 with the Arecibo and Goldstone planetary radars and the Very Long Baseline Array during December 2008. EV5 rotates retrograde and its overall shape is a 400 /pm 50 m oblate spheroid. The most prominent surface feature is a ridge parallel to the asteroid's equator that is broken by a concavity 150 m in diameter. Otherwise the asteroid's surface is notabl…
▽ More
We observed the near-Earth asteroid 2008 EV5 with the Arecibo and Goldstone planetary radars and the Very Long Baseline Array during December 2008. EV5 rotates retrograde and its overall shape is a 400 /pm 50 m oblate spheroid. The most prominent surface feature is a ridge parallel to the asteroid's equator that is broken by a concavity 150 m in diameter. Otherwise the asteroid's surface is notably smooth on decameter scales. EV5's radar and optical albedos are consistent with either rocky or stony-iron composition. The equatorial ridge is similar to structure seen on the rubble-pile near-Earth asteroid (66391) 1999 KW4 and is consistent with YORP spin-up reconfiguring the asteroid in the past. We interpret the concavity as an impact crater. Shaking during the impact and later regolith redistribution may have erased smaller features, explaining the general lack of decameter-scale surface structure.
△ Less
Submitted 19 January, 2011;
originally announced January 2011.
-
Tidal Evolution of Close Binary Asteroid Systems
Authors:
Patrick A. Taylor,
Jean-Luc Margot
Abstract:
We provide a generalized discussion of tidal evolution to arbitrary order in the expansion of the gravitational potential between two spherical bodies of any mass ratio. To accurately reproduce the tidal evolution of a system at separations less than five times the radius of the larger primary component, the tidal potential due to the presence of a smaller secondary component is expanded in terms…
▽ More
We provide a generalized discussion of tidal evolution to arbitrary order in the expansion of the gravitational potential between two spherical bodies of any mass ratio. To accurately reproduce the tidal evolution of a system at separations less than five times the radius of the larger primary component, the tidal potential due to the presence of a smaller secondary component is expanded in terms of Legendre polynomials to arbitrary order rather than truncated at leading order as is typically done in studies of well-separated system like the Earth and Moon. The equations of tidal evolution including tidal torques, the changes in spin rates of the components, and the change in semimajor axis (orbital separation) are then derived for binary asteroid systems with circular and equatorial mutual orbits. Accounting for higher-order terms in the tidal potential serves to speed up the tidal evolution of the system leading to underestimates in the time rates of change of the spin rates, semimajor axis, and mean motion in the mutual orbit if such corrections are ignored. Special attention is given to the effect of close orbits on the calculation of material properties of the components, in terms of the rigidity and tidal dissipation function, based on the tidal evolution of the system. It is found that accurate determinations of the physical parameters of the system, e.g., densities, sizes, and current separation, are typically more important than accounting for higher-order terms in the potential when calculating material properties. In the scope of the long-term tidal evolution of the semimajor axis and the component spin rates, correcting for close orbits is a small effect, but for an instantaneous rate of change in spin rate, semimajor axis, or mean motion, the close-orbit correction can be on the order of tens of percent.
△ Less
Submitted 7 January, 2011;
originally announced January 2011.
-
Orbits of Near-Earth Asteroid Triples 2001 SN263 and 1994 CC: Properties, Origin, and Evolution
Authors:
Julia Fang,
Jean-Luc Margot,
Marina Brozovic,
Michael C. Nolan,
Lance A. M. Benner,
Patrick A. Taylor
Abstract:
Three-body model fits to Arecibo and Goldstone radar data reveal the nature of two near-Earth asteroid triples. Triple-asteroid system 2001 SN263 is characterized by a primary of ~10^13 kg, an inner satellite ~1% as massive orbiting at ~3 primary radii in ~0.7 days, and an outer satellite ~2.5% as massive orbiting at ~13 primary radii in ~6.2 days. 1994 CC is a smaller system with a primary of mas…
▽ More
Three-body model fits to Arecibo and Goldstone radar data reveal the nature of two near-Earth asteroid triples. Triple-asteroid system 2001 SN263 is characterized by a primary of ~10^13 kg, an inner satellite ~1% as massive orbiting at ~3 primary radii in ~0.7 days, and an outer satellite ~2.5% as massive orbiting at ~13 primary radii in ~6.2 days. 1994 CC is a smaller system with a primary of mass ~2.6 \times 10^11 kg and two satellites ~2% and ~1% as massive orbiting at distances of ~5.5 and ~19.5 primary radii. Their orbital periods are ~1.2 and ~8.4 days. Examination of resonant arguments shows that the satellites are not currently in a mean-motion resonance. Precession of the apses and nodes are detected in both systems (2001 SN263 inner body: d{\varpi}/dt ~1.1 deg/day, 1994 CC inner body: d{\varpi}/dt ~ -0.2 deg/day), which is in agreement with analytical predictions of the secular evolution due to mutually interacting orbits and primary oblateness. Nonzero mutual inclinations between the orbital planes of the satellites provide the best fits to the data in both systems (2001 SN263: ~14 degrees, 1994 CC: ~16 degrees). Our best-fit orbits are consistent with nearly circular motion, except for 1994 CC's outer satellite which has an eccentric orbit of e ~ 0.19. We examine several processes that can generate the observed eccentricity and inclinations, including the Kozai and evection resonances, past mean-motion resonance crossings, and close encounters with terrestrial planets. In particular, we find that close planetary encounters can easily excite the eccentricities and mutual inclinations of the satellites' orbits to the currently observed values.
△ Less
Submitted 28 February, 2011; v1 submitted 9 December, 2010;
originally announced December 2010.
-
Long duration gamma-ray bursts: hydrodynamic instabilities in collapsar disks
Authors:
Paul A. Taylor,
John C. Miller,
Philipp Podsiadlowski
Abstract:
We present 3D numerical simulations of the early evolution of long-duration gamma-ray bursts in the collapsar scenario. Starting from the core-collapse of a realistic progenitor model, we follow the formation and evolution of a central black hole and centrifugally balanced disk. The dense, hot accretion disk produces freely-escaping neutrinos and is hydrodynamically unstable to clumping and to for…
▽ More
We present 3D numerical simulations of the early evolution of long-duration gamma-ray bursts in the collapsar scenario. Starting from the core-collapse of a realistic progenitor model, we follow the formation and evolution of a central black hole and centrifugally balanced disk. The dense, hot accretion disk produces freely-escaping neutrinos and is hydrodynamically unstable to clumping and to forming non-axisymmetric (m=1, 2) modes. We show that these spiral structures, which form on dynamical timescales, can efficiently transfer angular momentum outward and can drive the high required accretion rates (>=0.1-1 M_sun) for producing a jet. We utilise the smoothed particle hydrodynamics code, Gadget-2, modified to implement relevant microphysics, such as cooling by neutrinos, a plausible treatment approximating the central object and relativistic effects. Finally, we discuss implications of this scenario as a source of energy to produce relativistically beamed gamma-ray jets.
△ Less
Submitted 23 October, 2010; v1 submitted 23 June, 2010;
originally announced June 2010.
-
Migration of Interplanetary Dust
Authors:
S. I. Ipatov,
J. C. Mather,
P. A. Taylor
Abstract:
We numerically investigate the migration of dust particles with initial orbits close to those of the numbered asteroids, observed trans-Neptunian objects, and Comet Encke. The fraction of silicate asteroidal particles that collided with the Earth during their lifetime varied from 1.1% for 100 micron particles to 0.008% for 1 micron particles. Almost all asteroidal particles with diameter d>4 mic…
▽ More
We numerically investigate the migration of dust particles with initial orbits close to those of the numbered asteroids, observed trans-Neptunian objects, and Comet Encke. The fraction of silicate asteroidal particles that collided with the Earth during their lifetime varied from 1.1% for 100 micron particles to 0.008% for 1 micron particles. Almost all asteroidal particles with diameter d>4 microns collided with the Sun. The peaks in the migrating asteroidal dust particles' semi-major axis distribution at the n:(n+1) resonances with Earth and Venus and the gaps associated with the 1:1 resonances with these planets are more pronounced for larger particles. The probability of collisions of cometary particles with the Earth is smaller than for asteroidal particles, and this difference is greater for larger particles.
△ Less
Submitted 25 August, 2003;
originally announced August 2003.
-
Migration of Asteroidal Dust
Authors:
S. I. Ipatov,
J. C. Mather,
P. A. Taylor
Abstract:
We numerically investigate the migration of dust particles with initial orbits close to those of the numbered asteroids. The fraction of silicate particles that collided with the Earth during their lifetimes varied from 0.2% for 40 micron particles to 0.008% for 1 micron particles. Almost all particles with diameter d>4 microns collided with the Sun. The peaks in the migrating asteroidal dust pa…
▽ More
We numerically investigate the migration of dust particles with initial orbits close to those of the numbered asteroids. The fraction of silicate particles that collided with the Earth during their lifetimes varied from 0.2% for 40 micron particles to 0.008% for 1 micron particles. Almost all particles with diameter d>4 microns collided with the Sun. The peaks in the migrating asteroidal dust particles' semi-major axis distribution at the n/(n+1) resonances with Earth and Venus and the gaps associated with the 1:1 resonances with these planets are more pronounced for larger particles.
△ Less
Submitted 17 March, 2003;
originally announced March 2003.