Showing 1–39 of 39 results for author: Jackson, A P

RZ Piscium Hosts a Compact and Highly Perturbed Debris Disk

Authors: Kate Y. L. Su, Grant M. Kennedy, George H. Rieke, A. Meredith Hughes, Yu-Chia Lin, Jamar Kittling, Alan P. Jackson, Ramya M. Anche, Hauyu Baobab Liu

Abstract: RZ Piscium (RZ Psc) is well-known in the variable star field because of its numerous, irregular optical dips in the past five decades, but the nature of the system is heavily debated in the literature. We present multiyear infrared monitoring data from Spitzer and WISE to track the activities of the inner debris production, revealing stochastic infrared variability as short as weekly timescales th… ▽ More RZ Piscium (RZ Psc) is well-known in the variable star field because of its numerous, irregular optical dips in the past five decades, but the nature of the system is heavily debated in the literature. We present multiyear infrared monitoring data from Spitzer and WISE to track the activities of the inner debris production, revealing stochastic infrared variability as short as weekly timescales that is consistent with destroying a 90-km-size asteroid every year. ALMA 1.3 mm data combined with spectral energy distribution modeling show that the disk is compact ($\sim$0.1--13 au radially) and lacks cold gas. The disk is found to be highly inclined and has a significant vertical scale height. These observations confirm that RZ Psc hosts a close to edge-on, highly perturbed debris disk possibly due to migration of recently formed giant planets which might be triggered by the low-mass companion RZ Psc B if the planets formed well beyond the snowlines. △ Less

Submitted 18 October, 2023; originally announced October 2023.

Comments: 16 pages, 5 figures, accepted for publication in ApJ

Mercury's formation within the Early Instability Scenario

Authors: Matthew S. Clement, John E. Chambers, Nathan A. Kaib, Sean N. Raymond, Alan P. Jackson

Abstract: The inner solar system's modern orbital architecture provides inferences into the epoch of terrestrial planet formation; a ~100 Myr time period of planet growth via collisions with planetesimals and other proto-planets. While classic numerical simulations of this scenario adequately reproduced the correct number of terrestrial worlds, their semi-major axes and approximate formation timescales, the… ▽ More The inner solar system's modern orbital architecture provides inferences into the epoch of terrestrial planet formation; a ~100 Myr time period of planet growth via collisions with planetesimals and other proto-planets. While classic numerical simulations of this scenario adequately reproduced the correct number of terrestrial worlds, their semi-major axes and approximate formation timescales, they struggled to replicate the Earth-Mars and Venus-Mercury mass ratios. In a series of past independent investigations, we demonstrated that Mars' mass is possibly the result of Jupiter and Saturn's early orbital evolution, while Mercury's diminutive size might be the consequence of a primordial mass deficit in the region. Here, we combine these ideas in a single modeled scenario designed to simultaneously reproduce the formation of all four terrestrial planets and the modern orbits of the giant planets in broad strokes. By evaluating our Mercury analogs' core mass fractions, masses, and orbital offsets from Venus, we favor a scenario where Mercury forms through a series of violent erosive collisions between a number of ~Mercury-mass embryos in the inner part of the terrestrial disk. We also compare cases where the gas giants begin the simulation locked in a compact 3:2 resonant configuration to a more relaxed 2:1 orientation and find the former to be more successful. In 2:1 cases, the entire Mercury-forming region is often depleted due to strong sweeping secular resonances that also tend to overly excite the orbits of Earth and Venus as they grow. While our model is quite successful at replicating Mercury's massive core and dynamically isolated orbit, the planets' low mass remains extremely challenging to match. Finally, we discuss the merits and drawbacks of alternative evolutionary scenarios and initial disk conditions. △ Less

Submitted 23 January, 2023; originally announced January 2023.

Comments: 23 pages, 9 figures, 3 tables, accepted for publication in Icarus

A Star-sized Impact-Produced Dust Clump in the Terrestrial Zone of the HD 166191 System

Authors: Kate Y. L. Su, Grant M. Kennedy, Everett Schlawin, Alan P. Jackson, George H. Rieke

Abstract: We report on five years of 3-5 micron photometry measurements obtained by warm Spitzer to track the dust debris emission in the terrestrial zone of HD 166191 in combination with simultaneous optical data. We show that the debris production in this young (~10 Myr) system increased significantly in early 2018 and reached a record high level (almost double by mid 2019) by the end of the Spitzer missi… ▽ More We report on five years of 3-5 micron photometry measurements obtained by warm Spitzer to track the dust debris emission in the terrestrial zone of HD 166191 in combination with simultaneous optical data. We show that the debris production in this young (~10 Myr) system increased significantly in early 2018 and reached a record high level (almost double by mid 2019) by the end of the Spitzer mission (early 2020), suggesting intense collisional activity in its terrestrial zone likely due to either initial assembling of terrestrial planets through giant impacts or dynamical shake-up from unseen planet-mass objects or recent planet migration. This intense activity is further highlighted by detecting a star-size dust clump, passing in front of the star, in the midst of its infrared brightening. We constrain the minimum size and mass of the clump using multiwavelength transit profiles and conclude that the dust clump is most likely created by a large impact involving objects of several hundred km in size with an apparent period of 142 days (i.e., 0.62 au assuming a circular orbit). The system's evolutionary state (right after the dispersal of its gas-rich disk) makes it extremely valuable to learn about the process of terrestrial planet formation and planetary architecture through future observations. △ Less

Submitted 4 March, 2022; originally announced March 2022.

Comments: 17 pages, 5 figures, accepted by ApJ

RW Aur A : SpeX Spectral Evidence for Differentiated Planetesimal Formation, Migration and Destruction in an 3 Myr Old Excited CTTS System

Authors: C. M. Lisse, M. L. Sitko, S. J. Wolk, H. M. Günther, S. Brittain, J. D. Green, J. Steckloff, B. Johnson, C. C. Espaillat, M. Koutoukali, S. Y. Moorman, A. P. Jackson

Abstract: We present 2007 - 2020 SpeX VISNIR spectral monitoring of the highly variable RW Aur A CTTS. We find direct evidence for a highly excited, IR bright, asymmetric, and time variable system. Comparison of the spectral and temporal trends found determines 5 different components: (1) a stable continuum from 0.7 - 1.3 um, with approx color temperature 4000K, produced by the CTTS photospheric surface; (2… ▽ More We present 2007 - 2020 SpeX VISNIR spectral monitoring of the highly variable RW Aur A CTTS. We find direct evidence for a highly excited, IR bright, asymmetric, and time variable system. Comparison of the spectral and temporal trends found determines 5 different components: (1) a stable continuum from 0.7 - 1.3 um, with approx color temperature 4000K, produced by the CTTS photospheric surface; (2) variable hydrogen emission lines emitted from hot excited hydrogen in the CTTSs protostellar atmosphere/accretion envelope; (3) hot CO gas in the CTTSs protostellar atmosphere/accretion envelope; (4) highly variable 1.8-5.0 um thermal continuum emission with color temperature ranging from 1130 to 1650K, due to a surrounding accretion disk that is spatially variable and has an inner wall at r = 0.04 AU and T = 1650K, and outer edges at approx 1200K; and (5) transient, bifurcated signatures of abundant Fe II + associated SI, SiI, and SrI in the systems jet structures. The bifuracted signatures first appeared in 2015, but these collapsed and disappeared into a small single peak protostellar atmosphere feature by late 2020. The temporal evolution of RW Aur As spectral signatures is consistent with a dynamically excited CTTS system forming differentiated Vesta-sized planetesimals in an asymmetric accretion disk and migrating them inward to be destructively accreted. By contrast, nearby, coeval binary companion RW Aur B evinces only (1) a stable WTTS photospheric continuum from 0.7 - 1.3 um + (3) cold CO gas in absorption + (4) stable 1.8-5.0 um thermal disk continuum emission with color temperature approx 1650K. △ Less

Submitted 6 March, 2022; v1 submitted 25 January, 2022; originally announced January 2022.

Comments: 38 Pages, 16 Figures, 2 Tables, accepted for publication in Astrophys J. 24-Jan-2022

Carbon monoxide gas produced by a giant impact in the inner region of a young system

Authors: Tajana Schneiderman, Luca Matrà, Alan P. Jackson, Grant M. Kennedy, Quentin Kral, Sebastián Marino, Karin I. Öberg, Kate Y. L. Su, David J. Wilner, Mark C. Wyatt

Abstract: Models of terrestrial planet formation predict that the final stages of planetary assembly, lasting tens of millions of years beyond the dispersal of young protoplanetary disks, are dominated by planetary collisions. It is through these giant impacts that planets like the young Earth grow to their final mass and achieve long-term stable orbital configurations. A key prediction is that these impact… ▽ More Models of terrestrial planet formation predict that the final stages of planetary assembly, lasting tens of millions of years beyond the dispersal of young protoplanetary disks, are dominated by planetary collisions. It is through these giant impacts that planets like the young Earth grow to their final mass and achieve long-term stable orbital configurations. A key prediction is that these impacts produce debris. To date, the most compelling observational evidence for post-impact debris comes from the planetary system around the nearby 23 Myr-old A star HD 172555. This system shows large amounts of fine dust with an unusually steep size distribution and atypical dust composition, previously attributed to either a hypervelocity impact or a massive asteroid belt. Here, we report the spectrally resolved detection of a CO gas ring co-orbiting with dusty debris between ~6-9 au - a region analogous to the outer terrestrial planet region of our Solar System. Taken together, the dust and CO detections favor a giant impact between large, volatile-rich bodies. This suggests that planetary-scale collisions, analogous to the Moon-forming impact, can release large amounts of gas as well as debris, and that this gas is observable, providing a window into the composition of young planets. △ Less

Submitted 28 October, 2021; originally announced October 2021.

Comments: 21 pages (double spaced, including references and figures), 3 figures, 1 table, published in Nature 20/10/2021, published version available online through SharedIt initiative at https://rdcu.be/cAlEy

Journal ref: Nature 598, 425-428 (2021)

Matters Arising on "Breakup of a long-period comet as the origin of the dinosaur extinction" by Siraj & Loeb

Authors: Steven J. Desch, Alan P Jackson, Jessica L. Noviello, Ariel D. Anbar

Abstract: The recent publication by Siraj & Loeb (2021; Nature Scientific Reports 11, 3803) attempts to revive the debate over whether the Chicxulub impactor was a comet or an asteroid. They calculate that ~20% of long-period comets impacting Earth will have first been disrupted by passage inside the Sun's Roche limit, generating thousands of fragments, each the needed size of the Chicxulub impactor. This w… ▽ More The recent publication by Siraj & Loeb (2021; Nature Scientific Reports 11, 3803) attempts to revive the debate over whether the Chicxulub impactor was a comet or an asteroid. They calculate that ~20% of long-period comets impacting Earth will have first been disrupted by passage inside the Sun's Roche limit, generating thousands of fragments, each the needed size of the Chicxulub impactor. This would increase the impact rate of comets by a factor ~15, making them as likely to hit the Earth as an asteroid. They also argue that a comet would be a factor of 10 more likely to match the geochemical constraints, which indicate the Chicxulub impactor was carbonaceous chondrite-like. These conclusions are based on misinterpretations of the literature. Siraj & Loeb [1] overestimate the number of fragments produced during tidal disruption of a comet: tens of fragments are produced, not thousands. They also conflate 'carbonaceous chondrite' with specific types of carbonaceous chondrite, and ignore the evidence of iridium, making comets seem more likely than asteroids to match the Chicxulub impactor, when in fact they likely can be ruled out. Rather than a comet, an asteroidal impactor similar to CM or CR carbonaceous chondrites is strongly favored. △ Less

Submitted 2 June, 2021; originally announced June 2021.

Comments: 4 pages. Submitted to Nature Scientific Reports

Dynamical avenues for Mercury's origin I: The lone survivor of a primordial generation of short-period proto-planets

Authors: Matthew S. Clement, John E. Chambers, Alan P. Jackson

Abstract: The absence of planets interior to Mercury continues to puzzle terrestrial planet formation models, particularly when contrasted with the relatively high derived occurrence rates of short-period planets around Sun-like stars. Recent work proposed that the majority of systems hosting hot super-Earths attain their orbital architectures through an epoch of dynamical instability after forming in quasi… ▽ More The absence of planets interior to Mercury continues to puzzle terrestrial planet formation models, particularly when contrasted with the relatively high derived occurrence rates of short-period planets around Sun-like stars. Recent work proposed that the majority of systems hosting hot super-Earths attain their orbital architectures through an epoch of dynamical instability after forming in quasi-stable, tightly packed configurations. Isotopic evidence seems to suggest that the formation of objects in the super-Earth mass regime is unlikely to have occurred in the solar system as the terrestrial-forming disk is thought to have been significantly mass-deprived starting around 2 Myr after CAI; a consequence of either Jupiter's growth or an intrinsic disk feature. Nevertheless, terrestrial planet formation models and high-resolution investigations of planetesimal dynamics in the gas disk phase occasionally find that quasi-stable proto-planets with masses comparable to that of Mars emerge in the vicinity of Mercury's modern orbit. In this paper, we investigate whether it is possible for a primordial configuration of such objects to be cataclysmically destroyed in a manner that leaves Mercury behind as the sole survivor without disturbing the other terrestrial worlds. We use numerical simulations to show that this scenario is plausible. In many cases, the surviving Mercury analog experiences a series of erosive impacts; thereby boosting its Fe/Si ratio. A caveat of our proposed genesis scenario for Mercury is that Venus typically experiences at least one late giant impact. △ Less

Submitted 22 April, 2021; originally announced April 2021.

Comments: 18 pages, 15 figures, 1 table, accepted for publication in AJ. Summary available at https://youtu.be/Mu762MjRxtE (video)

1I/'Oumuamua as an N2 ice fragment of an exo-pluto surface. II. Generation of N2 ice fragments and the origin of 'Oumuamua

Authors: Steven J. Desch, Alan P. Jackson

Abstract: The origin of the interstellar object 1I/'Oumuamua, has defied explanation. In a companion paper (Jackson & Desch, 2021), we show that a body of N2 ice with axes 45 m x 44 m x 7.5 m at the time of observation would be consistent with its albedo, non-gravitational acceleration, and lack of observed CO or CO2 or dust. Here we demonstrate that impacts on the surfaces of Pluto-like Kuiper belt objects… ▽ More The origin of the interstellar object 1I/'Oumuamua, has defied explanation. In a companion paper (Jackson & Desch, 2021), we show that a body of N2 ice with axes 45 m x 44 m x 7.5 m at the time of observation would be consistent with its albedo, non-gravitational acceleration, and lack of observed CO or CO2 or dust. Here we demonstrate that impacts on the surfaces of Pluto-like Kuiper belt objects (KBOs) would have generated and ejected ~10^14 collisional fragments--roughly half of them H2O ice fragments and half of them N2 ice fragments--due to the dynamical instability that depleted the primordial Kuiper belt. We show consistency between these numbers and the frequency with which we would observe interstellar objects like 1I/'Oumuamua, and more comet-like objects like 2I/Borisov, if other stellar systems eject such objects with efficiency like that of the Sun; we infer that differentiated KBOs and dynamical instabilities that eject impact-generated fragments may be near-universal among extrasolar systems. Galactic cosmic rays would erode such fragments over 4.5 Gyr, so that fragments are a small fraction (~0.1%) of long-period Oort comets, but C/2016 R2 may be an example. We estimate 'Oumuamua was ejected about 0.4-0.5 Gyr ago, from a young (~10^8 yr) stellar system, which we speculate was in the Perseus arm. Objects like 'Oumuamua may directly probe the surface compositions of a hitherto-unobserved type of exoplanet: "exo-plutos". 'Oumuamua may be the first sample of an exoplanet brought to us. △ Less

Submitted 15 March, 2021; originally announced March 2021.

Comments: 29 pages, 1 figure. Companion to Jackson & Desch (2021), "1I/'Oumuamua as an N2 ice fragment of an exo-pluto surface. I. Size and compositional constraints"

1I/'Oumuamua as an N2 ice fragment of an exo-Pluto surface: I. size and compositional constraints

Authors: Alan P. Jackson, Steven J. Desch

Abstract: The origin of the interstellar object 1I/'Oumuamua has defied explanation. We perform calculations of the non-gravitational acceleration that would be experienced by bodies composed of a range of different ices and demonstrate that a body composed of N2 ice would satisfy the available constraints on the non-gravitational acceleration, size and albedo, and lack of detectable emission of CO or CO2 o… ▽ More The origin of the interstellar object 1I/'Oumuamua has defied explanation. We perform calculations of the non-gravitational acceleration that would be experienced by bodies composed of a range of different ices and demonstrate that a body composed of N2 ice would satisfy the available constraints on the non-gravitational acceleration, size and albedo, and lack of detectable emission of CO or CO2 or dust. We find that 'Oumuamua was small, with dimensions 45 m x 44 m x 7.5 m at the time of observation at 1.42 au from the Sun, with a high albedo of 0.64. This albedo is consistent with the N2 surfaces of bodies like Pluto and Triton. We estimate 'Oumuamua was ejected about 0.4-0.5 Gyr ago from a young stellar system, possibly in the Perseus arm. Objects like 'Oumuamua may directly probe the surface compositions of a hitherto-unobserved type of exoplanet: "exo-plutos". In a companion paper (Desch & Jackson, 2021) we demonstrate that dynamical instabilities like the one experienced by the Kuiper belt, in other stellar systems, plausibly could generate and eject large numbers of N2 ice fragments. 'Oumuamua may be the first sample of an exoplanet brought to us. △ Less

Submitted 15 March, 2021; originally announced March 2021.

Comments: 24 pages, 3 figures, 2 tables. Accepted for publication in JGR: Planets. Code available at https://doi.org/10.5281/zenodo.45584

Mid-infrared Studies of HD 113766 and HD 172555: Assessing Variability in the Terrestrial Zone of Young Exoplanetary Systems

Authors: Kate Y. L. Su, George H. Rieke, Carl Melis, Alan P. Jackson, Paul S. Smith, Huan Y. A. Meng, Andras Gaspar

Abstract: We present multi-epoch infrared photometry and spectroscopy obtained with warm Spitzer, Subaru and SOFIA to assess variability for the young ($\sim$20 Myr) and dusty debris systems around HD 172555 and HD 113766A. No variations (within 0.5%) were found for the former at either 3.6 or 4.5 $μ$m, while significant non-periodic variations (peak-to-peak of $\sim$10-15% relative to the primary star) wer… ▽ More We present multi-epoch infrared photometry and spectroscopy obtained with warm Spitzer, Subaru and SOFIA to assess variability for the young ($\sim$20 Myr) and dusty debris systems around HD 172555 and HD 113766A. No variations (within 0.5%) were found for the former at either 3.6 or 4.5 $μ$m, while significant non-periodic variations (peak-to-peak of $\sim$10-15% relative to the primary star) were detected for the latter. Relative to the Spitzer IRS spectra taken in 2004, multi-epoch mid-infrared spectra reveal no change in either the shape of the prominent 10 $μ$m solid-state features or the overall flux levels (no more than 20%) for both systems, corroborating that the population of sub-$μ$m-sized grains that produce the pronounced solid-state features is stable over a decadal timescale. We suggest that these sub-$μ$m-sized grains were initially generated in an optically thick clump of debris of mm-sized vapor condensates resulting from a recent violent impact between large asteroidal or planetary bodies. Because of the shielding from the stellar photons provided by this clump, intense collisions led to an over-production of fine grains that would otherwise be ejected from the system by radiation pressure. As the clump is sheared by its orbital motion and becomes optically thin, a population of very fine grains could remain in stable orbits until Poynting-Robertson drag slowly spirals them into the star. We further suggest that the 3-5 $μ$m disk variation around HD 113766A is consistent with a clump/arc of such fine grains on a modestly eccentric orbit in its terrestrial zone. △ Less

Submitted 11 June, 2020; originally announced June 2020.

Comments: to be published in ApJ

HD 145263: Spectral Observations of Silica Debris Disk Formation via Extreme Space Weathering?

Authors: C. M. Lisse, H. Y. A. Meng, M. L. Sitko, A. Morlok, B. C. Johnson, A. P. Jackson, R. J. Vervack Jr., C. H. Chen, S. J. Wolk, M. D. Lucas, M. Marengo, D. T. Britt

Abstract: We report here time domain infrared spectroscopy and optical photometry of the HD145263 silica-rich circumstellar disk system taken from 2003 through 2014. We find an F4V host star surrounded by a stable, massive 1e22 - 1e23 kg (M_Moon to M_Mars) dust disk. No disk gas was detected, and the primary star was seen rotating with a rapid ~1.75 day period. After resolving a problem with previously repo… ▽ More We report here time domain infrared spectroscopy and optical photometry of the HD145263 silica-rich circumstellar disk system taken from 2003 through 2014. We find an F4V host star surrounded by a stable, massive 1e22 - 1e23 kg (M_Moon to M_Mars) dust disk. No disk gas was detected, and the primary star was seen rotating with a rapid ~1.75 day period. After resolving a problem with previously reported observations, we find the silica, Mg-olivine, and Fe-pyroxene mineralogy of the dust disk to be stable throughout, and very unusual compared to the ferromagnesian silicates typically found in primordial and debris disks. By comparison with mid-infrared spectral features of primitive solar system dust, we explore the possibility that HD 145263's circumstellar dust mineralogy occurred with preferential destruction of Fe-bearing olivines, metal sulfides, and water ice in an initially comet-like mineral mix and their replacement by Fe-bearing pyroxenes, amorphous pyroxene, and silica. We reject models based on vaporizing optical stellar megaflares, aqueous alteration, or giant hypervelocity impacts as unable to produce the observed mineralogy. Scenarios involving unusually high Si abundances are at odds with the normal stellar absorption near-infrared feature strengths for Mg, Fe, and Si. Models involving intense space weathering of a thin surface patina via moderate (T < 1300 K) heating and energetic ion sputtering due to a stellar superflare from the F4V primary are consistent with the observations. The space weathered patina should be reddened, contain copious amounts of nanophase Fe, and should be transient on timescales of decades unless replenished. △ Less

Submitted 15 March, 2020; originally announced March 2020.

Comments: 41 Pages, 5 Figures, 5 Tables, Accepted for publication in the Astrophysical Journal

Automated crater shape retrieval using weakly-supervised deep learning

Authors: Mohamad Ali-Dib, Kristen Menou, Alan P. Jackson, Chenchong Zhu, Noah Hammond

Abstract: Crater ellipticity determination is a complex and time consuming task that so far has evaded successful automation. We train a state of the art computer vision algorithm to identify craters in Lunar digital elevation maps and retrieve their sizes and 2D shapes. The computational backbone of the model is MaskRCNN, an "instance segmentation" general framework that detects craters in an image while s… ▽ More Crater ellipticity determination is a complex and time consuming task that so far has evaded successful automation. We train a state of the art computer vision algorithm to identify craters in Lunar digital elevation maps and retrieve their sizes and 2D shapes. The computational backbone of the model is MaskRCNN, an "instance segmentation" general framework that detects craters in an image while simultaneously producing a mask for each crater that traces its outer rim. Our post-processing pipeline then finds the closest fitting ellipse to these masks, allowing us to retrieve the crater ellipticities. Our model is able to correctly identify 87\% of known craters in the longitude range we hid from the network during training and validation (test set), while predicting thousands of additional craters not present in our training data. Manual validation of a subset of these "new" craters indicates that a majority of them are real, which we take as an indicator of the strength of our model in learning to identify craters, despite incomplete training data. The crater size, ellipticity, and depth distributions predicted by our model are consistent with human-generated results. The model allows us to perform a large scale search for differences in crater diameter and shape distributions between the lunar highlands and maria, and we exclude any such differences with a high statistical significance. The predicted test set catalogue and trained model are available here: https://github.com/malidib/Craters_MaskRCNN/. △ Less

Submitted 11 March, 2020; v1 submitted 20 June, 2019; originally announced June 2019.

Comments: 59 pages, 13 figures, Accepted for publication in Icarus

Extreme Debris Disk Variability -- Exploring the Diverse Outcomes of Large Asteroid Impacts During the Era of Terrestrial Planet Formation

Authors: Kate Y. L. Su, Alan P. Jackson, Andras Gaspar, George H. Rieke, Ruobing Dong, Johan Olofsson, G. M. Kennedy, Zoe M. Leinhardt, Renu Malhotra, Michael Hammer, Huan Y. A. Meng, W. Rujopakarn, Joseph E. Rodriguez, Joshua Pepper, D. E. Reichart, David James, Keivan G. Stassun

Abstract: The most dramatic phases of terrestrial planet formation are thought to be oligarchic and chaotic growth, on timescales of up to 100-200 Myr, when violent impacts occur between large planetesimals of sizes up to proto-planets. Such events are marked by the production of large amounts of debris as has been observed in some exceptionally bright and young debris disks (termed extreme debris disks). H… ▽ More The most dramatic phases of terrestrial planet formation are thought to be oligarchic and chaotic growth, on timescales of up to 100-200 Myr, when violent impacts occur between large planetesimals of sizes up to proto-planets. Such events are marked by the production of large amounts of debris as has been observed in some exceptionally bright and young debris disks (termed extreme debris disks). Here we report five years of Spitzer measurements of such systems around two young solar-type stars: ID8 and P1121. The short-term (weekly to monthly) and long-term (yearly) disk variability is consistent with the aftermaths of large impacts involving large asteroid-size bodies. We demonstrate that an impact-produced clump of optically thick dust, under the influence of the dynamical and viewing geometry effects, can produce short-term modulation in the disk light curves. The long-term disk flux variation is related to the collisional evolution within the impact-produced fragments once released into a circumstellar orbit. The time-variable behavior observed in the P1121 system is consistent with a hypervelocity impact prior to 2012 that produced vapor condensates as the dominant impact product. Two distinct short-term modulations in the ID8 system argue for two violent impacts at different times and locations. Its long-term variation is consistent with the collisional evolution of two different populations of impact-produced debris dominated by either vapor condensates or escaping boulders. The bright, variable emission from the dust produced in large impacts from extreme debris disks provides a unique opportunity to study violent events during the era of terrestrial planet formation. △ Less

Submitted 19 April, 2019; v1 submitted 25 March, 2019; originally announced March 2019.

Comments: accepted for publication in AJ. typos and latex referencing fixed

Oort cloud asteroids: Collisional evolution, the Nice Model, and the Grand Tack

Authors: Andrew Shannon, Alan P. Jackson, Mark C. Wyatt

Abstract: If the Solar system had a history of planet migration, the signature of that migration may be imprinted on the populations of asteroids and comets that were scattered in the planets' wake. Here, we consider the dynamical and collisional evolution of inner Solar system asteroids which join the Oort cloud. We compare the Oort cloud asteroid populations produced by migration scenarios based on the `N… ▽ More If the Solar system had a history of planet migration, the signature of that migration may be imprinted on the populations of asteroids and comets that were scattered in the planets' wake. Here, we consider the dynamical and collisional evolution of inner Solar system asteroids which join the Oort cloud. We compare the Oort cloud asteroid populations produced by migration scenarios based on the `Nice' and `Grand Tack' scenarios, as well as a null hypothesis where the planets have not migrated, to the detection of one such object, C/2014 S3 (PANSTARRS). Our simulations find that the discovery of C/2014 S3 (PANSTARRS) only has a greater than one percent chance of occurring if the Oort cloud asteroids evolved on to Oort cloud orbits when the Solar system was not more than about one million years old, as this early transfer to the Oort cloud is necessary to keep the amount of collisional evolution low. We argue this only occurs when a giant (greater than thirty Earth masses) planet orbits at 1 ~ 2 au, and thus our results strongly favour a `Grand Tack'-like migration having occurred early in the Solar system's history. △ Less

Submitted 7 March, 2019; originally announced March 2019.

Comments: Resubmitted to MNRAS after revision from referee's report

Can a machine learn the outcome of planetary collisions?

Authors: Diana Valencia, Emaad Paracha, Alan P. Jackson

Abstract: Planetary-scale collisions are common during the last stages of formation of solid planets, including the Solar system terrestrial planets. The problem of growing planets has been divided into studying the gravitational interaction of embryos relevant in million year timescales and treated with N-body codes and the collision between objects with a timescale of hours to days and treated with smooth… ▽ More Planetary-scale collisions are common during the last stages of formation of solid planets, including the Solar system terrestrial planets. The problem of growing planets has been divided into studying the gravitational interaction of embryos relevant in million year timescales and treated with N-body codes and the collision between objects with a timescale of hours to days and treated with smoothed-particle hydrodynamics. These are now being coupled with simple parameterized models. We set out to investigate if machine learning techniques can offer a better solution by predicting the outcome of collisions which can then be used in N-body simulations. We considered three different supervised machine learning approaches: gradient boosting regression trees, nested models, and gaussian processes. We found that the former produced the best results, and that it was slightly surpassed by ensembling different algorithms. With gaussian processes, we found the regions of parameter space that may yield the most information to machine learning algorithms. Thus, we suggest SPH calculations to focus first on mass ratios above 0.5. △ Less

Submitted 11 February, 2019; originally announced February 2019.

Comments: 13 pages, 11 figure, 2 tables. Submitted to ApJ

Effect of Re-impacting Debris on the Solidification of the Lunar Magma Ocean

Authors: Viranga Perera, Alan P. Jackson, Linda T. Elkins-Tanton, Erik Asphaug

Abstract: Anorthosites that comprise the bulk of the lunar crust are believed to have formed during solidification of a Lunar Magma Ocean (LMO) in which these rocks would have floated to the surface. This early flotation crust would have formed a thermal blanket over the remaining LMO, prolonging solidification. Geochronology of lunar anorthosites indicates a long timescale of LMO cooling, or re-melting and… ▽ More Anorthosites that comprise the bulk of the lunar crust are believed to have formed during solidification of a Lunar Magma Ocean (LMO) in which these rocks would have floated to the surface. This early flotation crust would have formed a thermal blanket over the remaining LMO, prolonging solidification. Geochronology of lunar anorthosites indicates a long timescale of LMO cooling, or re-melting and re-crystallization in one or more late events. To better interpret this geochronology, we model LMO solidification in a scenario where the Moon is being continuously bombarded by returning projectiles released from the Moon-forming giant impact. More than one lunar mass of material escaped the Earth-Moon system onto heliocentric orbits following the giant impact, much of it to come back on returning orbits for a period of 100 Myr. If large enough, these projectiles would have punctured holes in the nascent floatation crust of the Moon, exposing the LMO to space and causing more rapid cooling. We model these scenarios using a thermal evolution model of the Moon that allows for production (by cratering) and evolution (solidification and infill) of holes in the flotation crust that insulates the LMO. For effective hole production, solidification of the magma ocean can be significantly expedited, decreasing the cooling time by more than a factor of 5. If hole production is inefficient, but shock conversion of projectile kinetic energy to thermal energy is efficient, then LMO solidification can be somewhat prolonged, lengthening the cooling time by 50% or more. △ Less

Submitted 12 April, 2018; originally announced April 2018.

Dynamical and biological panspermia constraints within multi-planet exosystems

Authors: Dimitri Veras, David J. Armstrong, James A. Blake, Jose F. Gutierrez-Marcos, Alan P. Jackson, Hendrik Schaefer

Abstract: As discoveries of multiple planets in the habitable zone of their parent star mount, developing analytical techniques to quantify extrasolar intra-system panspermia will become increasingly important. Here, we provide user-friendly prescriptions that describe the asteroid impact characteristics which would be necessary to transport life both inwards and outwards within these systems within a singl… ▽ More As discoveries of multiple planets in the habitable zone of their parent star mount, developing analytical techniques to quantify extrasolar intra-system panspermia will become increasingly important. Here, we provide user-friendly prescriptions that describe the asteroid impact characteristics which would be necessary to transport life both inwards and outwards within these systems within a single framework. Our focus is on projectile generation and delivery and our expressions are algebraic, eliminating the need for the solution of differential equations. We derive a probability distribution function for life-bearing debris to reach a planetary orbit, and describe the survival of micro-organisms during planetary ejection, their journey through interplanetary space, and atmospheric entry. △ Less

Submitted 11 November, 2018; v1 submitted 12 February, 2018; originally announced February 2018.

Comments: Published in Astrobiology, but colour plots are in the arXiv version only

The Taurus Boundary of Stellar/Substellar (TBOSS) Survey II. Disk Masses from ALMA Continuum Observations

Authors: K. Ward-Duong, J. Patience, J. Bulger, G. van der Plas, F. Ménard, C. Pinte, A. P. Jackson, G. Bryden, N. J. Turner, P. Harvey, A. Hales, R. J. De Rosa

Abstract: We report 885$μ$m ALMA continuum flux densities for 24 Taurus members spanning the stellar/substellar boundary, with spectral types from M4 to M7.75. Of the 24 systems, 22 are detected at levels ranging from 1.0-55.6 mJy. The two non-detections are transition disks, though other transition disks in the sample are detected. Converting ALMA continuum measurements to masses using standard scaling law… ▽ More We report 885$μ$m ALMA continuum flux densities for 24 Taurus members spanning the stellar/substellar boundary, with spectral types from M4 to M7.75. Of the 24 systems, 22 are detected at levels ranging from 1.0-55.6 mJy. The two non-detections are transition disks, though other transition disks in the sample are detected. Converting ALMA continuum measurements to masses using standard scaling laws and radiative transfer modeling yields dust mass estimates ranging from $\sim$0.3-20M$_{\oplus}$. The dust mass shows a declining trend with central object mass when combined with results from submillimeter surveys of more massive Taurus members. The substellar disks appear as part of a continuous sequence and not a distinct population. Compared to older Upper Sco members with similar masses across the substellar limit, the Taurus disks are brighter and more massive. Both Taurus and Upper Sco populations are consistent with an approximately linear relationship in $M_{dust}$ to $M_{star}$, although derived power-law slopes depend strongly upon choices of stellar evolutionary model and dust temperature relation. The median disk around early M-stars in Taurus contains a comparable amount of mass in small solids as the average amount of heavy elements in Kepler planetary systems on short-period orbits around M-dwarf stars, with an order of magnitude spread in disk dust mass about the median value. Assuming a gas:dust ratio of 100:1, only a small number of low-mass stars and brown dwarfs have a total disk mass amenable to giant planet formation, consistent with the low frequency of giant planets orbiting M-dwarfs. △ Less

Submitted 20 December, 2017; originally announced December 2017.

Comments: 41 pages and 32 figures, with all tables and appendices presented here in their entirety. Accepted for publication in AJ (November 26, 2017)

Ejection of rocky and icy material from binary star systems: Implications for the origin and composition of 1I/`Oumuamua

Authors: Alan P. Jackson, Daniel Tamayo, Noah Hammond, Mohamad Ali-Dib, Hanno Rein

Abstract: In single star systems like our own Solar system, comets dominate the mass budget of bodies that are ejected into interstellar space, since they form further away and are less tightly bound. However 1I/`Oumuamua, the first interstellar object detected, appears asteroidal in its spectra and in its lack of detectable activity. We argue that the galactic budget of interstellar objects like 1I/`Oumuam… ▽ More In single star systems like our own Solar system, comets dominate the mass budget of bodies that are ejected into interstellar space, since they form further away and are less tightly bound. However 1I/`Oumuamua, the first interstellar object detected, appears asteroidal in its spectra and in its lack of detectable activity. We argue that the galactic budget of interstellar objects like 1I/`Oumuamua should be dominated by planetesimal material ejected during planet formation in circumbinary systems, rather than in single star systems or widely separated binaries. We further show that in circumbinary systems, rocky bodies should be ejected in comparable numbers to icy ones. This suggests that a substantial fraction of additional interstellar objects discovered in the future should display an active coma. We find that the rocky population, of which 1I/`Oumuamua seems to be a member, should be predominantly sourced from A-type and late B-star binaries. △ Less

Submitted 28 March, 2018; v1 submitted 12 December, 2017; originally announced December 2017.

Comments: 5 pages, 2 figures. accepted in MNRAS Letters, updated with reviewer comments

Constraints on the pre-impact orbits of Solar System giant impactors

Authors: Alan P. Jackson, Travis S. J. Gabriel, Erik I. Asphaug

Abstract: We provide a fast method for computing constraints on impactor pre-impact orbits, applying this to the late giant impacts in the Solar System. These constraints can be used to make quick, broad comparisons of different collision scenarios, identifying some immediately as low-probability events, and narrowing the parameter space in which to target follow-up studies with expensive N-body simulations… ▽ More We provide a fast method for computing constraints on impactor pre-impact orbits, applying this to the late giant impacts in the Solar System. These constraints can be used to make quick, broad comparisons of different collision scenarios, identifying some immediately as low-probability events, and narrowing the parameter space in which to target follow-up studies with expensive N-body simulations. We benchmark our parameter space predictions, finding good agreement with existing N-body studies for the Moon. We suggest that high-velocity impact scenarios in the inner Solar System, including all currently proposed single impact scenarios for the formation of Mercury, should be disfavoured. This leaves a multiple hit-and-run scenario as the most probable currently proposed for the formation of Mercury. △ Less

Submitted 14 November, 2017; originally announced November 2017.

Comments: 13 pages, 9 figures, accepted in MNRAS

Cosmic Chandlery with Thermonuclear Supernovae

Authors: Alan C. Calder, Brendan K. Krueger, Aaron P. Jackson, Don E. Willcox, Broxton J. Miles, Dean M. Townsley

Abstract: Thermonuclear (Type Ia) supernovae are bright stellar explosions, the light curves of which can be calibrated to allow for use as "standard candles" for measuring cosmological distances. Contemporary research investigates how the brightness of an event may be influenced by properties of the progenitor system that follow from properties of the host galaxy such as composition and age. The goals are… ▽ More Thermonuclear (Type Ia) supernovae are bright stellar explosions, the light curves of which can be calibrated to allow for use as "standard candles" for measuring cosmological distances. Contemporary research investigates how the brightness of an event may be influenced by properties of the progenitor system that follow from properties of the host galaxy such as composition and age. The goals are to better understand systematic effects and to assess the intrinsic scatter in the brightness, thereby reducing uncertainties in cosmological studies. We present the results from ensembles of simulations in the single-degenerate paradigm addressing the influence of age and metallicity on the brightness of an event and compare our results to observed variations of brightness that correlate with properties of the host galaxy. We also present results from "hybrid" progenitor models that incorporate recent advances in stellar evolution. △ Less

Submitted 6 December, 2016; originally announced December 2016.

Comments: 14 pages, accepted to the proceedings of Astronum2016 - 11th International Conference on Numerical Modeling of Space Plasma Flows, June 6-10, 2016

How to design a planetary system for different scattering outcomes: giant impact sweet spot, maximising exocomets, scattered disks

Authors: M. C. Wyatt, A. Bonsor, A. P. Jackson, S. Marino, A. Shannon

Abstract: This paper considers the dynamics of scattering of planetesimals or planetary embryos by a planet on a circumstellar orbit. We classify six regions in the planet's mass versus semimajor axis parameter space according to the dominant outcome for scattered objects: ejected, accreted, remaining, escaping, Oort Cloud, depleted Oort Cloud. We use these outcomes to consider which planetary system archit… ▽ More This paper considers the dynamics of scattering of planetesimals or planetary embryos by a planet on a circumstellar orbit. We classify six regions in the planet's mass versus semimajor axis parameter space according to the dominant outcome for scattered objects: ejected, accreted, remaining, escaping, Oort Cloud, depleted Oort Cloud. We use these outcomes to consider which planetary system architectures maximise the observability of specific signatures, given that signatures should be detected first around systems with optimal architectures (if such systems exist in nature). Giant impact debris is most readily detectable for 0.1-10Mearth planets at 1-5au, depending on detection method and spectral type. While A stars have putative giant impact debris at 4-6au consistent with this sweet spot, that of FGK stars is typically <<1au contrary to expectations; an absence of 1-3au giant impact debris could indicate a low frequency of terrestrial planets there. Three principles maximise cometary influx from exo-Kuiper belts: a chain of closely separated planets interior to the belt, none of which is a Jupiter-like ejector; planet masses not increasing strongly with distance (for a net inward torque on comets); ongoing replenishment of comets, possibly by embedded low-mass planets. A high Oort Cloud comet influx requires no ejector and architectures that maximise the Oort Cloud population. Cold debris disks are usually considered classical Kuiper belt analogues. Here we consider the possibility of detecting scattered disk analogues, which could be betrayed by a broad radial profile and lack of small grains, as well as spherical 100-1000au mini-Oort Clouds. Some implications for escaping planets around young stars, detached planets akin to Sedna, and the formation of super-Earths, are also discussed. △ Less

Submitted 3 October, 2016; originally announced October 2016.

Comments: 25 pages, MNRAS in press

The spherical Brazil Nut Effect and its significance to asteroids

Authors: Viranga Perera, Alan P. Jackson, Erik Asphaug, Ronald-Louis Ballouz

Abstract: Many asteroids are likely rubble-piles that are a collection of smaller objects held together by gravity and possibly cohesion. These asteroids are seismically shaken by impacts, which leads to excitation of their constituent particles. As a result it has been suggested that their surfaces and sub-surface interiors may be governed by a size sorting mechanism known as the Brazil Nut Effect. We stud… ▽ More Many asteroids are likely rubble-piles that are a collection of smaller objects held together by gravity and possibly cohesion. These asteroids are seismically shaken by impacts, which leads to excitation of their constituent particles. As a result it has been suggested that their surfaces and sub-surface interiors may be governed by a size sorting mechanism known as the Brazil Nut Effect. We study the behavior of a model asteroid that is a spherical, self-gravitating aggregate with a binary size-distribution of particles under the action of applied seismic shaking. We find that above a seismic threshold, larger particles rise to the surface when friction is present, in agreement with previous studies that focussed on cylindrical and rectangular box configurations. Unlike previous works we also find that size sorting takes place even with zero friction, though the presence of friction does aid the sorting process above the seismic threshold. Additionally we find that while strong size sorting can take place near the surface, the innermost regions remain unsorted under even the most vigorous shaking. △ Less

Submitted 18 August, 2016; originally announced August 2016.

Journal ref: Icarus, Volume 278, 1 November 2016, Pages 194-203

Insights into planet formation from debris disks: II. Giant impacts in extrasolar planetary systems

Authors: Mark C. Wyatt, Alan P. Jackson

Abstract: Giant impacts refer to collisions between two objects each of which is massive enough to be considered at least a planetary embryo. The putative collision suffered by the proto-Earth that created the Moon is a prime example, though most Solar System bodies bear signatures of such collisions. Current planet formation models predict that an epoch of giant impacts may be inevitable, and observations… ▽ More Giant impacts refer to collisions between two objects each of which is massive enough to be considered at least a planetary embryo. The putative collision suffered by the proto-Earth that created the Moon is a prime example, though most Solar System bodies bear signatures of such collisions. Current planet formation models predict that an epoch of giant impacts may be inevitable, and observations of debris around other stars are providing mounting evidence that giant impacts feature in the evolution of many planetary systems. This chapter reviews giant impacts, focussing on what we can learn about planet formation by studying debris around other stars. Giant impact debris evolves through mutual collisions and dynamical interactions with planets. General aspects of this evolution are outlined, noting the importance of the collision-point geometry. The detectability of the debris is discussed using the example of the Moon-forming impact. Such debris could be detectable around another star up to 10Myr post-impact, but model uncertainties could reduce detectability to a few 100yr window. Nevertheless the 3% of young stars with debris at levels expected during terrestrial planet formation provide valuable constraints on formation models; implications for super-Earth formation are also discussed. Variability recently observed in some bright disks promises to illuminate the evolution during the earliest phases when vapour condensates may be optically thick and acutely affected by the collision-point geometry. The outer reaches of planetary systems may also exhibit signatures of giant impacts, such as the clumpy debris structures seen around some stars. △ Less

Submitted 15 March, 2016; originally announced March 2016.

Comments: 43 pages, accepted. Chapter in International Space Science Institute (ISSI) Book on "The Disk in Relation to the Formation of Planets and their Proto-atmospheres", to be published in Space Science Reviews by Springer

On Measuring the Metallicity of a Type Ia Supernova's Progenitor

Authors: Broxton J. Miles, Daniel R. van Rossum, Dean M. Townsley, F. X. Timmes, Aaron P. Jackson, Alan C. Calder, Edward F. Brown

Abstract: In Type Ia Supernovae (\sneia), the relative abundances of chemical elements are affected by the neutron excess in the composition of the progenitor white dwarf. Since these products leave signatures in the spectra near maximum light, spectral features may be used to constrain the composition of the progenitor. We calculate the nucleosynthetic yields for three \snia simulations, assuming single de… ▽ More In Type Ia Supernovae (\sneia), the relative abundances of chemical elements are affected by the neutron excess in the composition of the progenitor white dwarf. Since these products leave signatures in the spectra near maximum light, spectral features may be used to constrain the composition of the progenitor. We calculate the nucleosynthetic yields for three \snia simulations, assuming single degenerate, Chandrasekhar mass progenitors, for a wide range of progenitor metallicities, and calculate synthetic light curves and spectra to explore correlations between progenitor metallicity and the strength of spectral features. We use two 2D simulations of the deflagration-detonation-transition scenario with different $^{56}$Ni yields and the W7 simulation to control for differences between explosion models and total yields. While the overall yields of intermediate mass elements (16 $<$ A $\leq$ 40) differ between the three cases, trends in the yields are similar. With increasing metallicity, $^{28}$Si yields remain nearly constant, $^{40}$Ca yields decline, and Ti and $^{54}$Fe yields increase. In the synthetic spectra, we identify two features at 30 days post explosion that appear to deepen with progenitor metallicity: a Ti feature around 4200\,Å and a Fe feature around 5200\,Å\@. In all three simulations, their pseudo equivalent widths show a systematic trend with progenitor metallicity. This suggests that these two features may allow differentiation among progenitor metallicities of observed \sneia and potentially help reduce the intrinsic Hubble scatter. △ Less

Submitted 3 May, 2016; v1 submitted 24 August, 2015; originally announced August 2015.

Comments: 20 pages, 18 figures, accepted by the Astrophysical Journal

Eight billion asteroids in the Oort cloud

Authors: Andrew Shannon, Alan P. Jackson, Dimitri Veras, Mark Wyatt

Abstract: The Oort cloud is usually thought of as a collection of icy comets inhabiting the outer reaches of the Solar system, but this picture is incomplete. We use simulations of the formation of the Oort cloud to show that ~4% of the small bodies in the Oort cloud should have formed within 2.5 au of the Sun, and hence be ice-free rock-iron bodies. If we assume these Oort cloud asteroids have the same siz… ▽ More The Oort cloud is usually thought of as a collection of icy comets inhabiting the outer reaches of the Solar system, but this picture is incomplete. We use simulations of the formation of the Oort cloud to show that ~4% of the small bodies in the Oort cloud should have formed within 2.5 au of the Sun, and hence be ice-free rock-iron bodies. If we assume these Oort cloud asteroids have the same size distribution as their cometary counterparts, the Large Synoptic Survey Telescope should find roughly a dozen Oort cloud asteroids during ten years of operations. Measurement of the asteroid fraction within the Oort cloud can serve as an excellent test of the Solar system's formation and dynamical history. Oort cloud asteroids could be of particular concern as impact hazards as their high mass density, high impact velocity, and low visibility make them both hard to detect and hard to divert or destroy. However, they should be a rare class of object, and we estimate globally catastrophic collisions should only occur about once per billion years. △ Less

Submitted 27 October, 2014; originally announced October 2014.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal

Molecular Gas Clumps from the Destruction of Icy Bodies in the $β$ Pictoris Debris Disk

Authors: W. R. F. Dent, M. C. Wyatt, A. Roberge, J. -C. Augereau, S. Casassus, S. Corder, J. S. Greaves, I. de Gregorio-Monsalvo, A. Hales, A. P. Jackson, A. Meredith Hughes, A. -M. Lagrange, B. Matthews, D. Wilner

Abstract: Many stars are surrounded by disks of dusty debris formed in the collisions of asteroids, comets and dwarf planets. But is gas also released in such events? Observations at submm wavelengths of the archetypal debris disk around $β$ Pictoris show that 0.3% of a Moon mass of carbon monoxide orbits in its debris belt. The gas distribution is highly asymmetric, with 30% found in a single clump 85AU fr… ▽ More Many stars are surrounded by disks of dusty debris formed in the collisions of asteroids, comets and dwarf planets. But is gas also released in such events? Observations at submm wavelengths of the archetypal debris disk around $β$ Pictoris show that 0.3% of a Moon mass of carbon monoxide orbits in its debris belt. The gas distribution is highly asymmetric, with 30% found in a single clump 85AU from the star, in a plane closely aligned with the orbit of the inner planet, $β$ Pic b. This gas clump delineates a region of enhanced collisions, either from a mean motion resonance with an unseen giant planet, or from the remnants of a collision of Mars-mass planets. △ Less

Submitted 4 April, 2014; originally announced April 2014.

Journal ref: Science (2014) 343, 1490

Debris from giant impacts between planetary embryos at large orbital radii

Authors: Alan P. Jackson, Mark C. Wyatt, Amy Bonsor, Dimitri Veras

Abstract: We consider the observational signatures of giant impacts between planetary embryos. While the debris released in the impact remains in a clump for only a single orbit, there is a much longer lasting asymmetry caused by the fact that all debris must pass through the collision-point. The resulting asymmetry is stationary, it does not orbit the star. The debris is concentrated in a clump at the coll… ▽ More We consider the observational signatures of giant impacts between planetary embryos. While the debris released in the impact remains in a clump for only a single orbit, there is a much longer lasting asymmetry caused by the fact that all debris must pass through the collision-point. The resulting asymmetry is stationary, it does not orbit the star. The debris is concentrated in a clump at the collision-point, with a more diffuse structure on the opposite side. The asymmetry lasts for typically around 1000 orbital periods of the progenitor, which can be several Myr at distances of ~50 AU. We describe how the appearance of the asymmetric disc depends on the mass and eccentricity of the progenitor, as well as viewing orientation. The wavelength of observation, which determines the grain sizes probed, is also important. Notably, the increased collision rate of the debris at the collision-point makes this the dominant production site for any secondary dust and gas created. For dust small enough to be removed by radiation pressure, and gas with a short lifetime, this causes their distribution to resemble a jet emanating from the (stationary) collision-point. We suggest that the asymmetries seen at large separations in some debris discs, like Beta Pictoris, could be the result of giant impacts. If so this would indicate that planetary embryos are present and continuing to grow at several tens of AU at ages of up to tens of Myr. △ Less

Submitted 7 March, 2014; originally announced March 2014.

Comments: 20 pages, 17 figures, 2 appendices, accepted in MNRAS

Power-Law Wrinkling Turbulence-Flame Interaction Model for Astrophysical Flames

Authors: Aaron P. Jackson, Dean M. Townsley, Alan C. Calder

Abstract: We extend a model for turbulence-flame interactions (TFI) to consider astrophysical flames with a particular focus on combustion in type Ia supernovae. The inertial range of the turbulent cascade is nearly always under-resolved in simulations of astrophysical flows, requiring the use of a model in order to quantify the effects of subgrid-scale wrinkling of the flame surface. We provide implementat… ▽ More We extend a model for turbulence-flame interactions (TFI) to consider astrophysical flames with a particular focus on combustion in type Ia supernovae. The inertial range of the turbulent cascade is nearly always under-resolved in simulations of astrophysical flows, requiring the use of a model in order to quantify the effects of subgrid-scale wrinkling of the flame surface. We provide implementation details to extend a well-tested TFI model to low-Prandtl number flames for use in the compressible hydrodynamics code FLASH. A local, instantaneous measure of the turbulent velocity is calibrated for FLASH and verification tests are performed. Particular care is taken to consider the relation between the subgrid rms turbulent velocity and the turbulent flame speed, especially for high-intensity turbulence where the turbulent flame speed is not expected to scale with the turbulent velocity. Finally, we explore the impact of different TFI models in full-star, three-dimensional simulations of type Ia supernovae. △ Less

Submitted 18 February, 2014; originally announced February 2014.

Comments: 20 pages, 12 figures, accepted to the Astrophysical Journal

The influence of chemical composition on models of Type Ia supernovae

Authors: Alan C. Calder, Brendan K. Krueger, Aaron P. Jackson, Dean M. Townsley

Abstract: Type Ia supernovae are bright stellar explosions distinguished by standardizable light curves that allow for their use as distance indicators for cosmological studies. Despite the highly successful use of these events in this capacity, many fundamental questions remain. Contemporary research investigates how properties of the progenitor system that follow from the host galaxy such as composition a… ▽ More Type Ia supernovae are bright stellar explosions distinguished by standardizable light curves that allow for their use as distance indicators for cosmological studies. Despite the highly successful use of these events in this capacity, many fundamental questions remain. Contemporary research investigates how properties of the progenitor system that follow from the host galaxy such as composition and age influence the brightness of an event with the goal of better understanding and assessing the intrinsic scatter in the brightness. We provide an overview of these supernovae and proposed progenitor systems, all of which involve one or more compact stars known as white dwarfs. We describe contemporary research investigating how the composition and structure of the progenitor white dwarf systematically influences the explosion outcome assuming the progenitor is a single white dwarf that has gained mass from a companion. We present results illustrating some of these systematic effects from our research. △ Less

Submitted 9 March, 2013; originally announced March 2013.

Comments: 22 pages, 3 figures, Frontiers of Physics, in press

Evaluating Systematic Dependencies of Type Ia Supernovae: The Influence of Central Density

Authors: Brendan K. Krueger, Aaron P. Jackson, Alan C. Calder, Dean M. Townsley, Edward F. Brown, Francis X. Timmes

Abstract: We present a study exploring a systematic effect on the brightness of type Ia supernovae using numerical models that assume the single-degenerate paradigm. Our investigation varied the central density of the progenitor white dwarf at flame ignition, and considered its impact on the explosion yield, particularly the production and distribution of radioactive Ni-56, which powers the light curve. We… ▽ More We present a study exploring a systematic effect on the brightness of type Ia supernovae using numerical models that assume the single-degenerate paradigm. Our investigation varied the central density of the progenitor white dwarf at flame ignition, and considered its impact on the explosion yield, particularly the production and distribution of radioactive Ni-56, which powers the light curve. We performed a suite of two-dimensional simulations with randomized initial conditions, allowing us to characterize the statistical trends that we present. The simulations indicate that production of Fe-group material is statistically independent of progenitor central density, but the mass of stable Fe-group isotopes is tightly correlated with central density, with a decrease in the production of Ni-56 at higher central densities. These results imply progenitors with higher central densities produce dimmer events. We provide details of the post-explosion distribution of Ni-56 in the models, including the lack of a consistent centrally-located deficit of Ni-56, which may be compared to observed remnants. By performing a self-consistent extrapolation of our model yields and considering the main-sequence lifetime of the progenitor star and the elapsed time between the formation of the white dwarf and the onset of accretion, we develop a brightness-age relation that improves our prediction of the expected trend for single degenerates and we compare this relation with observations. △ Less

Submitted 9 August, 2012; originally announced August 2012.

Comments: 25 pages, 18 figures, accepted to ApJ on July 30, 2012

Journal ref: Astrophysical Journal 757 (2012) 175

Debris from terrestrial planet formation: the Moon-forming collision

Authors: Alan P. Jackson, Mark C. Wyatt

Abstract: We study the evolution of debris created in the giant impacts expected during the final stages of terrestrial planet formation. The starting point is the debris created in a simulation of the Moon-forming impact. The dynamical evolution is followed for 10 Myr including the effects of Earth, Venus, Mars and Jupiter. The spatial distribution evolves from a clump in the first few months to an asymmet… ▽ More We study the evolution of debris created in the giant impacts expected during the final stages of terrestrial planet formation. The starting point is the debris created in a simulation of the Moon-forming impact. The dynamical evolution is followed for 10 Myr including the effects of Earth, Venus, Mars and Jupiter. The spatial distribution evolves from a clump in the first few months to an asymmetric ring for the first 10 kyr and finally becoming an axisymmetric ring by about 1 Myr after the impact. By 10 Myr after the impact 20% of the particles have been accreted onto Earth and 17% onto Venus, with 8% ejected by Jupiter and other bodies playing minor roles. However, the fate of the debris also depends strongly on how fast it is collisionally depleted, which depends on the poorly constrained size distribution of the impact debris. Assuming that the debris is made up of 30% by mass mm-cm-sized vapour condensates and 70% boulders up to 500 km, we find that the condensates deplete rapidly on ~1000 yr timescales, whereas the boulders deplete predominantly dynamically. By considering the luminosity of dust produced in collisions within the boulder-debris distribution we find that the Moon-forming impact would have been readily detectable around other stars in Spitzer 24 micron surveys for around 25 Myr after the impact, with levels of emission comparable to many known hot dust systems. The vapour condensates meanwhile produce a short-lived, optically thick, spike of emission. We use these surveys to make an estimate of the fraction of stars that form terrestrial planets, F_TPF. Since current terrestrial planet formation models invoke multiple giant impacts, the low fraction of 10-100 Myr stars found to have warm (~150 K) dust implies that F_TPF ~<10%. △ Less

Submitted 19 June, 2012; originally announced June 2012.

Comments: 23 pages, 19 figures (14 colour), 1 table. Accepted in MNRAS

Planetary evaporation by UV & X-ray radiation: basic hydrodynamics

Authors: James E. Owen, Alan P. Jackson

Abstract: We consider the evaporation of close in planets by the star's intrinsic EUV and X-ray radiation. We calculate evaporation rates by solving the hydrodynamical problem for planetary evaporation including heating from both X-ray and EUV radiation. We show that most close-in planets (a<0.1 AU) are evaporating hydrodynamically, with the evaporation occurring in two distinct regimes: X-ray driven, in wh… ▽ More We consider the evaporation of close in planets by the star's intrinsic EUV and X-ray radiation. We calculate evaporation rates by solving the hydrodynamical problem for planetary evaporation including heating from both X-ray and EUV radiation. We show that most close-in planets (a<0.1 AU) are evaporating hydrodynamically, with the evaporation occurring in two distinct regimes: X-ray driven, in which the X-ray heated flow contains a sonic point, and EUV driven, in which the X-ray region is entirely sub-sonic. The mass-loss rates scale as L_X/a^2 for X-ray driven evaporation, and as Phi_*^{1/2}/a for EUV driven evaporation at early times, with mass-loss rates of order 10e10-10e14 g/s. No exact scaling exists for the mass-loss rate with planet mass and planet radius, however, in general evaporation proceeds more rapidly for planets with lower densities and higher masses. Furthermore, we find that in general the transition from X-ray driven to EUV driven evaporation occurs at lower X-ray luminosities for planets closer to their parent stars and for planets with lower densities. Coupling our evaporation models to the evolution of the high energy radiation - which falls with time - we are able to follow the evolution of evaporating planets. We find that most planets start off evaporating in the X-ray driven regime, but switch to EUV driven once the X-ray luminosity falls below a critical value. The evolution models suggest that while `hot Jupiters' are evaporating, they are not evaporating at a rate sufficient to remove the entire gaseous envelope on Gyr time-scales. However, we do find that close in Neptune mass planets are more susceptible to complete evaporation of their envelopes. Thus we conclude that planetary evaporation is more important for lower mass planets, particularly those in the `hot Neptune'/`super Earth' regime. △ Less

Submitted 11 June, 2012; originally announced June 2012.

Comments: 18 Pages, Accepted for publication in MNRAS

On Simulating Type Ia Supernovae

Authors: Alan C. Calder, Brendan K. Krueger, Aaron P. Jackson, Dean M. Townsley, Edward F. Brown, Francis X. Timmes

Abstract: Type Ia supernovae are bright stellar explosions distinguished by standardizable light curves that allow for their use as distance indicators for cosmological studies. Despite their highly successful use in this capacity, the progenitors of these events are incompletely understood. We describe simulating type Ia supernovae in the paradigm of a thermonuclear runaway occurring in a massive white dwa… ▽ More Type Ia supernovae are bright stellar explosions distinguished by standardizable light curves that allow for their use as distance indicators for cosmological studies. Despite their highly successful use in this capacity, the progenitors of these events are incompletely understood. We describe simulating type Ia supernovae in the paradigm of a thermonuclear runaway occurring in a massive white dwarf star. We describe the multi-scale physical processes that realistic models must incorporate and the numerical models for these that we employ. In particular, we describe a flame-capturing scheme that addresses the problem of turbulent thermonuclear combustion on unresolved scales. We present the results of our study of the systematics of type Ia supernovae including trends in brightness following from properties of the host galaxy that agree with observations. We also present performance results from simulations on leadership-class architectures. △ Less

Submitted 4 May, 2012; originally announced May 2012.

Comments: 13 pages, 3 figures, accepted to proceedings of the Conference on Computational Physics, Oct. 30 - Nov. 3, 2011

The coronal X-ray - age relation and its implications for the evaporation of exoplanets

Authors: Alan P. Jackson, Timothy A. Davis, Peter J. Wheatley

Abstract: We study the relationship between coronal X-ray emission and stellar age for late-type stars, and the variation of this relationship with spectral type. We select 717 stars from 13 open clusters and find that the ratio of X-ray to bolometric luminosity during the saturated phase of coronal emission decreases from 10^-3.1 for late K-dwarfs to 10^-4.3 for early F-type stars (across the range 0.29<(B… ▽ More We study the relationship between coronal X-ray emission and stellar age for late-type stars, and the variation of this relationship with spectral type. We select 717 stars from 13 open clusters and find that the ratio of X-ray to bolometric luminosity during the saturated phase of coronal emission decreases from 10^-3.1 for late K-dwarfs to 10^-4.3 for early F-type stars (across the range 0.29<(B-V)_0<1.41). Our determined saturation timescales vary between 10^7.6 and 10^8.3 years, though with no clear trend across the whole FGK range. We apply our X-ray emission - age relations to the investigation of the evaporation history of 121 known transiting exoplanets using a simple energy-limited model of evaporation and taking into consideration Roche lobe effects and different heating/evaporation efficiencies. We confirm that a linear cut-off of the planet distribution in the M^2/R^3 versus a^-2 plane is an expected result of population modification by evaporation and that the known transiting exoplanets display such a cut-off. We find that for an evaporation efficiency of 25% we expect around 1 in 10 of the known transiting exoplanets to have lost > 5% of their mass since formation. In addition we provide estimates of the minimum formation mass for which a planet could be expected to survive for 4 Gyrs for a range of stellar and planetary parameters. We emphasise the importance of the earliest periods of a planet's life for its evaporation history with 75% expected to occur within the first Gyr. This raises the possibility of using evaporation histories to distinguish between different migration scenarios. For planets with available spin-orbit angles no difference is found between the distributions of planets with misaligned orbits and those with aligned orbits. This suggests that misalignment occurs early in the life of the planetary system, though more data is needed. △ Less

Submitted 10 April, 2012; v1 submitted 31 October, 2011; originally announced November 2011.

Comments: 20 pages, 14 figures, accepted by MNRAS. Full versions of Tables 3, 4 and 5 available with the online journal version

Evaluating Systematic Dependencies of Type Ia Supernovae

Authors: Alan C. Calder, Brendan K. Krueger, Aaron P. Jackson, Dean M. Townsley, Francis X. Timmes, Edward F. Brown, David A. Chamulak

Abstract: Type Ia supernovae are bright stellar explosions thought to occur when a thermonuclear runaway consumes roughly a solar mass of degenerate stellar material. These events produce and disseminate iron-peak elements, and properties of their light curves allow for standardization and subsequent use as cosmological distance indicators. The explosion mechanism of these events remains, however, only part… ▽ More Type Ia supernovae are bright stellar explosions thought to occur when a thermonuclear runaway consumes roughly a solar mass of degenerate stellar material. These events produce and disseminate iron-peak elements, and properties of their light curves allow for standardization and subsequent use as cosmological distance indicators. The explosion mechanism of these events remains, however, only partially understood. Many models posit the explosion beginning with a deflagration born near the center of a white dwarf that has gained mass from a stellar companion. In order to match observations, models of this single-degenerate scenario typically invoke a subsequent transition of the (subsonic) deflagration to a (supersonic) detonation that rapidly consumes the star. We present an investigation into the systematics of thermonuclear supernovae assuming this paradigm. We utilize a statistical framework for a controlled study of two-dimensional simulations of these events from randomized initial conditions. We investigate the effect of the composition and thermal history of the progenitor on the radioactive yield, and thus brightness, of an event. Our results offer an explanation for some observed trends of mean brightness with properties of the host galaxy. △ Less

Submitted 17 November, 2010; originally announced November 2010.

Comments: 14 pages, to appear in the Proceedings of the SciDAC 2010 meeting

Evaluating Systematic Dependencies of Type Ia Supernovae: The Influence of Deflagration to Detonation Density

Authors: Aaron P. Jackson, Alan C. Calder, Dean M. Townsley, David A. Chamulak, Edward F. Brown, F. X. Timmes

Abstract: We explore the effects of the deflagration to detonation transition (DDT) density on the production of Ni-56 in thermonuclear supernova explosions (type Ia supernovae). Within the DDT paradigm, the transition density sets the amount of expansion during the deflagration phase of the explosion and therefore the amount of nuclear statistical equilibrium (NSE) material produced. We employ a theoretica… ▽ More We explore the effects of the deflagration to detonation transition (DDT) density on the production of Ni-56 in thermonuclear supernova explosions (type Ia supernovae). Within the DDT paradigm, the transition density sets the amount of expansion during the deflagration phase of the explosion and therefore the amount of nuclear statistical equilibrium (NSE) material produced. We employ a theoretical framework for a well-controlled statistical study of two-dimensional simulations of thermonuclear supernovae with randomized initial conditions that can, with a particular choice of transition density, produce a similar average and range of Ni-56 masses to those inferred from observations. Within this framework, we utilize a more realistic "simmered" white dwarf progenitor model with a flame model and energetics scheme to calculate the amount of Ni-56 and NSE material synthesized for a suite of simulated explosions in which the transition density is varied in the range 1-3x10^7 g/cc. We find a quadratic dependence of the NSE yield on the log of the transition density, which is determined by the competition between plume rise and stellar expansion. By considering the effect of metallicity on the transition density, we find the NSE yield decreases by 0.055 +/- 0.004 solar masses for a 1 solar metallicity increase evaluated about solar metallicity. For the same change in metallicity, this result translates to a 0.067 +/- 0.004 solar mass decrease in the Ni-56 yield, slightly stronger than that due to the variation in electron fraction from the initial composition. Observations testing the dependence of the yield on metallicity remain somewhat ambiguous, but the dependence we find is comparable to that inferred from some studies. △ Less

Submitted 7 July, 2010; originally announced July 2010.

Comments: 15 pages, 13 figures, accepted to ApJ on July 6, 2010

Journal ref: Astrophys.J. 720 (2010) 99-113

On variations of the brightness of type Ia supernovae with the age of the host stellar population

Authors: Brendan K. Krueger, Aaron P. Jackson, Dean M. Townsley, Alan C. Calder, Edward F. Brown, F. X. Timmes

Abstract: Recent observational studies of type Ia supernovae (SNeIa) suggest correlations between the peak brightness of an event and the age of the progenitor stellar population. This trend likely follows from properties of the progenitor white dwarf (WD), such as central density, that follow from properties of the host stellar population. We present a statistically well-controlled, systematic study utiliz… ▽ More Recent observational studies of type Ia supernovae (SNeIa) suggest correlations between the peak brightness of an event and the age of the progenitor stellar population. This trend likely follows from properties of the progenitor white dwarf (WD), such as central density, that follow from properties of the host stellar population. We present a statistically well-controlled, systematic study utilizing a suite of multi-dimensional SNeIa simulations investigating the influence of central density of the progenitor WD on the production of Fe-group material, particularly radioactive Ni-56, which powers the light curve. We find that on average, as the progenitor's central density increases, production of Fe-group material does not change but production of Ni-56 decreases. We attribute this result to a higher rate of neutronization at higher density. The central density of the progenitor is determined by the mass of the WD and the cooling time prior to the onset of mass transfer from the companion, as well as the subsequent accretion heating and neutrino losses. The dependence of this density on cooling time, combined with the result of our central density study, offers an explanation for the observed age-luminosity correlation: a longer cooling time raises the central density at ignition thereby producing less Ni-56 and thus a dimmer event. While our ensemble of results demonstrates a significant trend, we find considerable variation between realizations, indicating the necessity for averaging over an ensemble of simulations to demonstrate a statistically significant result. △ Less

Submitted 6 July, 2010; originally announced July 2010.

Comments: 5 pages, 4 figures, 1 table, accepted to ApJL

Journal ref: Astrophysical Journal Letters 719 (2010) L5-L9

Evaluating Systematic Dependencies of Type Ia Supernovae: The Influence of Progenitor Ne22 Content on Dynamics

Authors: Dean M. Townsley, Aaron P. Jackson, Alan C. Calder, David A. Chamulak, Edward F. Brown, F. X. Timmes

Abstract: We present a theoretical framework for formal study of systematic effects in Supernovae Type Ia (SN Ia) that utilizes 2-d simulations to implement a form of the deflagration-detonation transition (DDT) explosion scenario. The framework is developed from a randomized initial condition that leads to a sample of simulated SN Ia whose Ni56 masses have a similar average and range to those observed, a… ▽ More We present a theoretical framework for formal study of systematic effects in Supernovae Type Ia (SN Ia) that utilizes 2-d simulations to implement a form of the deflagration-detonation transition (DDT) explosion scenario. The framework is developed from a randomized initial condition that leads to a sample of simulated SN Ia whose Ni56 masses have a similar average and range to those observed, and have many other modestly realistic features such as the velocity extent of intermediate mass elements. The intended purpose is to enable statistically well-defined studies of both physical and theoretical parameters of the SN Ia explosion simulation. We present here a thorough description of the outcome of the SN Ia explosions produced by our current simulations. A first application of this framework is utilized to study the dependence of the SN Ia on the Ne22 content, which is known to be directly influenced by the progenitor stellar population's metallicity. Our study is very specifically tailored to measure how the Ne22 content influences the competition between the rise of plumes of burned material and the expansion of the star before these plumes reach DDT conditions. This competition controls the amount of material in nuclear statistical equilibrium (NSE) and therefore Ni56 produced by setting the density at which nucleosynthesis takes place during the detonation phase of the explosion. Although the outcome following from any particular ignition condition can change dramatically with Ne22 content, with a sample of 20 ignition conditions we find that the systematic change in the expansion of the star prior to detonation is not large enough to compete with the dependence on initial neutron excess discussed by Timmes, Brown & Truran (2003). (Abridged) △ Less

Submitted 23 June, 2009; originally announced June 2009.

Comments: 25 pages 15 figures, accepted to ApJ

Journal ref: Astrophys.J.701:1582-1604,2009