Showing 1–50 of 57 results for author: Boley, A C

Satellite Visibility During the April 2024 Total Eclipse

Authors: Samantha M. Lawler, Hanno Rein, Aaron C. Boley

Abstract: On 8 April 2024, tens of millions of people across North America will be able to view a total solar eclipse. Such astronomical events have been important throughout history, but with nearly 10,000 satellites in orbit, we question whether total eclipses will now reveal a sky full of satellites, fundamentally changing this experience for humanity. Using the current population of Starlink satellites,… ▽ More On 8 April 2024, tens of millions of people across North America will be able to view a total solar eclipse. Such astronomical events have been important throughout history, but with nearly 10,000 satellites in orbit, we question whether total eclipses will now reveal a sky full of satellites, fundamentally changing this experience for humanity. Using the current population of Starlink satellites, we find that the brightest satellites would be naked-eye visible in dark skies, but the high sky brightness during totality will make them undetectable to the unaided eye. Our model does not take into account the effects of chance reflections from large, flat surfaces like solar panels, which we expect will cause glints and flares that could be visible from large satellites and abandoned rocket bodies. Time-lapse all-sky imaging might reveal satellites during the eclipse. △ Less

Submitted 27 March, 2024; originally announced March 2024.

Comments: submitted to RNAAS

Mass and Angular Momentum Transport in a Gravitationally Unstable Protoplanetary Disk with Improved 3D Radiative Hydrodynamics

Authors: Thomas Y. Steiman-Cameron, Richard H. Durisen, Aaron C. Boley, Scott Michael, Karna Desai, Caitlin R. McConnell

Abstract: During early phases of a protoplanetary disks's life, gravitational instabilities can produce significant mass transport, can dramatically alter disk structure, can mix and shock-process gas and solids, and may be instrumental in planet formation. We present a 3D grid-based radiative hydrodynamics study with varied resolutions of a 0.07 M$_\odot$ disk orbiting a 0.5 M$_\odot$ star as it settles ov… ▽ More During early phases of a protoplanetary disks's life, gravitational instabilities can produce significant mass transport, can dramatically alter disk structure, can mix and shock-process gas and solids, and may be instrumental in planet formation. We present a 3D grid-based radiative hydrodynamics study with varied resolutions of a 0.07 M$_\odot$ disk orbiting a 0.5 M$_\odot$ star as it settles over most of its radial extent into a quasi-steady asymptotic state that maintains approximate balance between heating produced by GIs and radiative cooling governed by realistic dust opacities. We assess disk stability criteria, thermodynamic properties, strengths of GIs, characteristics of density waves and torques produced by GIs, radial mass transport arising from these torques, and the level to which transport can be represented as local or nonlocal processes. Physical and thermal processes display distinct differences between inner optically thick and outer optically thin regions of the disk. In the inner region, gravitational torques are dominated by low-order Fourier components of the azimuthal mass distribution. These torques are strongly variable on the local dynamical time and are subject to rapid flaring presumably driven by recurrent swing amplification. In the outer region, m=1 torques dominate. Ring-like structures exhibiting strong noncircular motions, and vortices develop near the inner edge between 8 and 14 au. We find that GI-induced spiral modes erupt in a chaotic manner over the whole low-Q part of the disk, with many spiral modes appearing and disappearing, producing gravitoturbulence, but dominated by fluctuating large-scale modes, very different from a simple $α$-disk. △ Less

Submitted 31 October, 2023; originally announced November 2023.

Comments: 44 pages, 24 figures, to appear in The Astrophysical Journal

Evidence of Long-Term Period Variations in the Exoplanet Transit Database (ETD)

Authors: Simone R. Hagey, Billy Edwards, Aaron C. Boley

Abstract: We analyze a large number of citizen science data and identify eight Hot Jupiter systems that show evidence for deviations from a constant orbital period: HAT-P-19 b, HAT-P-32 b, TrES-1 b, TrES-2 b, TrES-5 b, WASP-4 b, WASP-10 b, and WASP-12 b. The latter system is already well known to exhibit strong evidence for tidal orbital decay and serves as an important control for this study. Several other… ▽ More We analyze a large number of citizen science data and identify eight Hot Jupiter systems that show evidence for deviations from a constant orbital period: HAT-P-19 b, HAT-P-32 b, TrES-1 b, TrES-2 b, TrES-5 b, WASP-4 b, WASP-10 b, and WASP-12 b. The latter system is already well known to exhibit strong evidence for tidal orbital decay and serves as an important control for this study. Several other systems we identify have disputed period drifts in the literature, allowing the results here to serve as an independent analysis. The citizen science data are from the Exoplanet Transit Database (ETD), which is a global project established in 2008 by the Variable Star and Exoplanet Section of the Czech Astronomical Society. With over 400 planets and 12,000 contributed observations spanning 15 years, the ETD is brimming with potential for studying the long-term orbital evolution of close-in Hot Jupiters. We use our results to discuss prioritization of targets for follow up investigations, which will be necessary to confirm the period drifts and their causes. △ Less

Submitted 21 September, 2022; originally announced September 2022.

Comments: 65 pages, 10 figures, accepted to AJ

Journal ref: AJ 164 (2022)

Investigating the risks of debris-generating ASAT tests in the presence of megaconstellations

Authors: Sarah Thiele, Aaron C. Boley

Abstract: The development of large constellations of satellites (i.e., so-called megaconstellations or satcons) is poised to increase the number of LEO satellites by more than an order of magnitude in the coming decades. Such a rapid growth of satellite numbers makes the consequences of major fragmentation events ever more problematic. In this study, we investigate the collisional risk to on-orbit infrastru… ▽ More The development of large constellations of satellites (i.e., so-called megaconstellations or satcons) is poised to increase the number of LEO satellites by more than an order of magnitude in the coming decades. Such a rapid growth of satellite numbers makes the consequences of major fragmentation events ever more problematic. In this study, we investigate the collisional risk to on-orbit infrastructure from kinetic anti-satellite (ASAT) weapon tests, using the 2019 Indian test as a model. We find that the probability of one or more collisions occurring over the lifetime of ASAT fragments increases significantly in a satcon environment compared with the orbital environment in 2019. For the case of 65,000 satellites in LEO, we find that the chance of one or more satellites being struck by ASAT fragments of size 1 cm or larger is more than 25% for a single test. Including sizes down to 3 mm in our models suggests that impacts will occur for any such event. Finally, we apply our methods to examine the November 2021 Russian ASAT test, also finding a significant collision probability over the lifetime of the fragments. The heavy commercialization of LEO demands a commitment to avoiding debris-generating ASAT tests. △ Less

Submitted 22 December, 2022; v1 submitted 23 November, 2021; originally announced November 2021.

Comments: 27 pages, 7 figures, 6 tables. Accepted for publication in the Journal of Astronautical Sciences

Plaskett 1.8 metre Observations of Starlink Satellites

Authors: Aaron C. Boley, Ewan Wright, Samantha Lawler, Paul Hickson, Dave Balam

Abstract: We present observations of 23 Starlink satellites in the $g'$ bandpass, obtained from the Dominion Astrophysical Observatory's Plaskett 1.8 metre telescope. The targets include a mixture of satellites with and without brightness mitigation measures (i.e., visors). At the time of the observations (16 July 2021), Starlink satellites were sunlight throughout the night, and even with strict elevation… ▽ More We present observations of 23 Starlink satellites in the $g'$ bandpass, obtained from the Dominion Astrophysical Observatory's Plaskett 1.8 metre telescope. The targets include a mixture of satellites with and without brightness mitigation measures (i.e., visors). At the time of the observations (16 July 2021), Starlink satellites were sunlight throughout the night, and even with strict elevation and azimuth limits, there were over 800 candidate Starlink arcs. The satellites altogether have a median absolute brightness (550 km) of $\overline{H}_g^{550} =5.3$ mag. Dividing the targets into those without and with visors, their median absolute magnitudes are $\overline{H}_g^{550}(no~visor)=5.1$ and $\overline{H}^{550}_g(visor)=5.7$ mag, respectively. While the visor sample is dimmer in aggregate, the absolute brightness distribution ranged from $H_g^{550}=4.3$ mag to 9.4 mag, with the brightest being a visored satellite and the dimmest a satellite with no mitigation. The intrinsic brightness dispersion among the full sample is $σ_g = 0.5$ mag. △ Less

Submitted 26 September, 2021; originally announced September 2021.

Comments: Submitted to AJ

Visibility Predictions for Near-Future Satellite Megaconstellations: Latitudes near 50 Degrees will Experience the Worst Light Pollution

Authors: S. M. Lawler, A. C. Boley, H. Rein

Abstract: Megaconstellations of thousands to tens of thousands of artificial satellites (satcons) are rapidly being developed and launched. These satcons will have negative consequences for observational astronomy research, and are poised to drastically interfere with naked-eye stargazing worldwide should mitigation efforts be unsuccessful. Here we provide predictions for the optical brightnesses and on-sky… ▽ More Megaconstellations of thousands to tens of thousands of artificial satellites (satcons) are rapidly being developed and launched. These satcons will have negative consequences for observational astronomy research, and are poised to drastically interfere with naked-eye stargazing worldwide should mitigation efforts be unsuccessful. Here we provide predictions for the optical brightnesses and on-sky distributions of several satcons, including Starlink, OneWeb, Kuiper, and StarNet/GW, for a total of 65,000 satellites on their filed or predicted orbits. We develop a simple model of satellite reflectivity, which is calibrated using published Starlink observations. We use this model to estimate the visible magnitudes and on-sky distributions for these satellites as seen from different places on Earth, in different seasons, and different times of night. For latitudes near 50 degrees North and South, satcon satellites make up a few percent of all visible point sources all night long near the summer solstice, as well as near sunrise and sunset on the equinoxes. Altering the satellites' altitudes only changes the specific impacts of the problem. Without drastic reduction of the reflectivities, or significantly fewer total satellites in orbit, satcons will significantly change the night sky worldwide. △ Less

Submitted 23 November, 2021; v1 submitted 9 September, 2021; originally announced September 2021.

Comments: accepted to AJ

Unlocking the Origins of Ultracool Dwarf Radio Emission

Authors: A. G. Hughes, A. C. Boley, R. A. Osten, J. A. White, M. Leacock

Abstract: Empirical trends in stellar X-ray and radio luminosities suggest that low mass ultracool dwarfs (UCDs) should not produce significant radio emission. Defying these expectations, strong non-thermal emission has been observed in a few UCDs in the 1-10 GHz range, with a variable component often attributed to global aurorae and a steady component attributed to other processes such as gyrosynchrotron e… ▽ More Empirical trends in stellar X-ray and radio luminosities suggest that low mass ultracool dwarfs (UCDs) should not produce significant radio emission. Defying these expectations, strong non-thermal emission has been observed in a few UCDs in the 1-10 GHz range, with a variable component often attributed to global aurorae and a steady component attributed to other processes such as gyrosynchrotron emission. While both auroral and gyrosynchrotron emission peak near the critical frequency, only the latter radiation is expected to extend into millimeter wavelengths. We present ALMA 97.5 GHz and VLA 33 GHz observations of a small survey of 5 UCDs. LP 349-25, LSR J1835+3259, and NLTT 33370 were detected at 97.5 GHz, while LP 423-31 and LP 415-20 resulted in non-detections at 33 GHz. A significant flare was observed in NLTT 33370 that reached a peak flux of 4880 +/- 360 microJy, exceeding the quiescent flux by nearly an order of magnitude, and lasting 20 seconds. These ALMA observations show bright 97.5 GHz emission with spectral indices ranging from alpha = -0.76 to alpha = -0.29, suggestive of optically thin gyrosynchrotron emission. If such emission traces magnetic reconnection events, then this could have consequences for both UCD magnetic models and the atmospheric stability of planets in orbit around them. Overall, our results provide confirmation that gyrosynchrotron radiation in radio loud UCDs can remain detectable into the millimeter regime. △ Less

Submitted 14 July, 2021; originally announced July 2021.

Comments: 10 pages, 4 figures, accepted in AJ

Precautionary Planetary Defence

Authors: Aaron C. Boley, Michael Byers

Abstract: The question of whether to attempt deflections during planetary defence emergencies has been subject to considerable decision-making analysis (Schmidt 2018; SMPAG Ad-Hoc Working Group on Legal Issues 2020). Hypothetical situations usually involve a newly discovered asteroid with a high impact probability on a set timescale. This paper addresses two further complexities: (1) limiting missions to an… ▽ More The question of whether to attempt deflections during planetary defence emergencies has been subject to considerable decision-making analysis (Schmidt 2018; SMPAG Ad-Hoc Working Group on Legal Issues 2020). Hypothetical situations usually involve a newly discovered asteroid with a high impact probability on a set timescale. This paper addresses two further complexities: (1) limiting missions to an asteroid due to the risk of a human-caused Earth impact; and (2) active management of asteroids to place them in "safe harbours", even when impact risks are otherwise below "decision to act" thresholds. We use Apophis as a case study, and address the two complexities in turn. △ Less

Submitted 21 June, 2021; originally announced June 2021.

Comments: Conference paper presented at the 7th IAA Planetary Defense Conference, 26-30 April 2021, Vienna, Austria. IAA-PDC-21-11-37

Collision rates of planetesimals near mean-motion resonances

Authors: Spencer C. Wallace, Thomas. R. Quinn, Aaron C. Boley

Abstract: In circumstellar discs, collisional grinding of planetesimals produces second-generation dust. While it remains unclear whether this ever becomes a major component of the total dust content, the presence of such dust, and potentially the substructure within, it can be used to explore a disc's physical conditions. A perturbing planet produces nonaxisymmetric structures and gaps in the dust, regardl… ▽ More In circumstellar discs, collisional grinding of planetesimals produces second-generation dust. While it remains unclear whether this ever becomes a major component of the total dust content, the presence of such dust, and potentially the substructure within, it can be used to explore a disc's physical conditions. A perturbing planet produces nonaxisymmetric structures and gaps in the dust, regardless of its origin. The dynamics of planetesimals, however, will be very different than that of small dust grains due to weaker gas interactions. Therefore, planetesimal collisions could create dusty disc structures that would not exist otherwise. In this work, we use N-body simulations to investigate the collision rate profile of planetesimals near mean-motion resonances. We find that a distinct bump or dip feature is produced in the collision profile, the presence of which depends on the libration width of the resonance and the separation between the peri- and apocenter distances of the edges of the resonance. The presence of one of these two features depends on the mass and eccentricity of the planet. Assuming that the radial dust emission traces the planetesimal collision profile, the presence of a bump or dip feature in the dust emission at the 2:1 mean-motion resonance can constrain the orbital properties of the perturbing planet. This assumption is valid, so long as radial drift does not play a significant role during the collisional cascade process. Under this assumption, these features in the dust emission should be marginally observable in nearby protoplanetary disks with ALMA. △ Less

Submitted 24 February, 2021; originally announced February 2021.

Comments: 17 pages, 12 figures, accepted for publication in MNRAS

Transit Duration Variations in Multi-Planet Systems

Authors: Aaron C. Boley, Christa Van Laerhoven, Agueda Paula Granados Contreras

Abstract: A planet's orbital orientation relative to an observer's line of sight determines the chord length for a transiting planet, i.e., the projected distance a transiting planet travels across the stellar disc. For a given circular orbit, the chord length determines the transit duration. Changes in the orbital inclination, the direction of the ascending node, or both, can alter this chord length and th… ▽ More A planet's orbital orientation relative to an observer's line of sight determines the chord length for a transiting planet, i.e., the projected distance a transiting planet travels across the stellar disc. For a given circular orbit, the chord length determines the transit duration. Changes in the orbital inclination, the direction of the ascending node, or both, can alter this chord length and thus result in transit duration variations (TDVs). Variation of the full orbital inclination vector can even lead to de-transiting or newly transiting planets for a system. We use Laplace-Lagrange secular theory to estimate the fastest nodal eigenfrequencies for over 100 short-period planetary systems. The highest eigenfrequency is an indicator of which systems should show the strongest TDVs. We further explore five cases (TRAPPIST-1, Kepler-11, K2-138, Kepler-445, and Kepler-334) using direct N-body simulations to characterize possible TDVs and to explore whether de-transiting planets could be possible for these systems. A range of initial conditions are explored, with each realization being consistent with the observed transits. We find that tens of percent of multiplanet systems have fast enough eigenfrequencies to expect large TDVs on decade timescales. Among the directly integrated cases, we find that de-transiting planets could occur on decade timescales and TDVs of 10 minutes per decade should be common. △ Less

Submitted 25 March, 2020; originally announced March 2020.

Comments: To appear in AJ

Meteoroid Stream Formation Due to the Extraction of Space Resources from Asteroids

Authors: Logan Fladeland, Aaron C. Boley, Michael Byers

Abstract: [Abridged] Asteroid mining is not necessarily a distant prospect. Hayabusa2 and OSIRIS-REx have recently rendezvoused with near-Earth asteroids and will return samples to Earth. While there is significant science motivation for these missions, there are also resource interests. Space agencies and commercial entities are particularly interested in ices and water-bearing minerals that could be used… ▽ More [Abridged] Asteroid mining is not necessarily a distant prospect. Hayabusa2 and OSIRIS-REx have recently rendezvoused with near-Earth asteroids and will return samples to Earth. While there is significant science motivation for these missions, there are also resource interests. Space agencies and commercial entities are particularly interested in ices and water-bearing minerals that could be used to produce rocket fuel in space. The internationally coordinated roadmaps of major space agencies depend on utilizing the natural resources of such celestial bodies. Several companies have already created plans for intercepting and extracting water and minerals from near-Earth objects, as even a small asteroid could have high economic worth. However, the low surface gravity of asteroids could make the release of mining waste and the subsequent formation of debris streams a consequence of asteroid mining. Strategies to contain material during extraction could still eventually require the purposeful jettison of waste to avoid managing unwanted mass. Using simulations, we explore the formation of mining debris streams by integrating particles released from four select asteroids. Radiation effects are included, and a range of debris sizes are explored. The simulation results are used to investigate the timescales for debris stream formation, the sizes of the streams, and the meteoroid fluxes compared with sporadic meteoroids. We find that for prodigious mining activities resulting in the loss of a few percent of the asteroid's mass or more, it is possible to produce streams that exceed the sporadic flux during stream crossing for some meteoroid sizes. The result of these simulations are intended to highlight potential unintended consequences that could result from NewSpace activity, which could help to inform efforts to develop international space resource guidelines. △ Less

Submitted 28 November, 2019; originally announced November 2019.

Comments: Conference paper for the 1st International Orbital Debris Conference, December 9-12, 2019, Sugar Land, TX, USA

Planetary Astronomy-Understanding the Origin of the Solar System

Authors: S. M. Lawler, A. C. Boley, M. Connors, W. Fraser, B. Gladman, C. L. Johnson, J. J. Kavelaars, G. Osinski, L. Philpott, J. Rowe, P. Wiegert, R. Winslow

Abstract: There is a vibrant and effective planetary science community in Canada. We do research in the areas of meteoritics, asteroid and trans-Neptunian object orbits and compositions, and space weather, and are involved in space probe missions to study planetary surfaces and interiors. For Canadian planetary scientists to deliver the highest scientific impact possible, we have several recommendations. Ou… ▽ More There is a vibrant and effective planetary science community in Canada. We do research in the areas of meteoritics, asteroid and trans-Neptunian object orbits and compositions, and space weather, and are involved in space probe missions to study planetary surfaces and interiors. For Canadian planetary scientists to deliver the highest scientific impact possible, we have several recommendations. Our top recommendation is to join LSST and gain access to the full data releases by hosting a data centre, which could be done by adding to the CADC, which is already highly involved in hosting planetary data and supporting computational modelling for orbital studies. We also support MSE, which can provide spectroscopy and thus compositional information for thousands of small bodies. We support a Canadian-led microsatellite, POEP, which will provide small body sizes by measuring occultations. We support the idea of piggybacking space weather instruments on other astronomical space probes to provide data for the space weather community. Many Canadian planetary scientists are involved in space probe missions, but through haphazard and temporary arrangements like co-appointments at US institutions, so we would like the community to support Canadian researchers to participate in these large, international missions. △ Less

Submitted 29 October, 2019; originally announced October 2019.

Comments: White paper submitted for the Canadian Long Range Plan LRP2020

Small and Moderate Aperture Telescopes for Research and Education

Authors: Aaron C. Boley, Terry Bridges, Paul Hickson, Harvey Richer, Brett Gladman, Jeremy Heyl, J. J. Kavelaars, Ingrid Stairs

Abstract: In this white paper (WP), we highlight several examples of small and moderate aperture telescopes that are being used for education and/or research. We further discuss potential costs for establishing new, small observatories, as well as joining existing international consortia. The WP includes a brief overview of select observing sites, with a discussion on how small telescopes at exceptional obs… ▽ More In this white paper (WP), we highlight several examples of small and moderate aperture telescopes that are being used for education and/or research. We further discuss potential costs for establishing new, small observatories, as well as joining existing international consortia. The WP includes a brief overview of select observing sites, with a discussion on how small telescopes at exceptional observing locations can be competitive, under certain circumstances, with larger and more expensive facilities located at poorer sites. In addition to research, these facilities enable many different types of educational experiences for wide range of people, from high school students to undergraduates to graduate students to postdocs. Canada should remain committed to partnering with large, international observatories such as CFHT, Gemini, and TMT, but it should also negotiate international agreements and commit funding to expand the use of small and moderate research observatories at domestic and international sites through coordination with the NRC, the Tri-Council, and the Canadian Foundation for Innovation. Both capital and operational costs (with site rental costs allowed) need to be included in support possibilities. CASCA should establish and maintain a small to moderate telescope expression of interest database that would help to facilitate Canadian institutions in organizing consortia, particularly for smaller institutions. The astronomical community should work with the NRC to make existing facilities more accessible to the astronomical community for research. This could involve, for example, automating the Plaskett and/or providing travel funds for supporting classical observing modes. Finally, a small to moderate aperture facility in the Arctic would be a world-class observatory and should be advanced over the next decade. △ Less

Submitted 1 October, 2019; originally announced October 2019.

Comments: White paper submitted to the Canadian Long Range Plan LRP2020

The Role of NewSpace in Furthering Canadian Astronomy

Authors: Aaron C. Boley, David Kendall, Michael Byers, Frederic J. Grandmont, Cameron Byers, Jennifer Busler, William MacDonald Evans, Brett Gladman, Tanya Harrison, Catherine Johnson

Abstract: [Highly abridged, from executive summary] As much as NewSpace presents opportunities, there are significant challenges that must be overcome, requiring engagement with policy makers to influence domestic and international space governance. Failure to do so could result in a range of long-lasting negative outcomes for science and space stewardship. How will the Canadian astronomical community engag… ▽ More [Highly abridged, from executive summary] As much as NewSpace presents opportunities, there are significant challenges that must be overcome, requiring engagement with policy makers to influence domestic and international space governance. Failure to do so could result in a range of long-lasting negative outcomes for science and space stewardship. How will the Canadian astronomical community engage with NewSpace? What are the implications for NewSpace on the astro-environment, including Earth orbits, lunar and cis-lunar orbits, and surfaces of celestial bodies? This white paper analyzes the rapid changes in space use and what those changes could mean for Canadian astronomers. Our recommendations are as follows: Greater cooperation between the astronomical and the Space Situational Awareness communities is needed. Build closer ties between the astronomical community and Global Affairs Canada (GAC). Establish a committee for evaluating the astro-environmental impacts of human space use, including on and around the Moon and other bodies. CASCA and the Tri-Council should coordinate to identify programs that would enable Canadian astronomers to participate in pay-for-use services at appropriate funding levels. CASCA should continue to foster a relationship with CSA, but also build close ties to the private space industry. Canadian-led deep space missions are within Canada's capabilities, and should be pursued. △ Less

Submitted 1 October, 2019; originally announced October 2019.

Comments: White paper submitted for the Canadian Long Range Plan LRP2020

The MESAS Project: Long wavelength follow-up observations of Sirius A

Authors: Jacob Aaron White, Jason Aufdenberg, Aaron C. Boley, Mark Devlin, Peter Hauschildt, Anna G. Hughes, A. Meredith Hughes, Brian Mason, Brenda Matthews, Attila Moór, Tony Mroczkowski, Charles Romero, Jonathan Sievers, Sara Stanchfield, Francisco Tapia, David J. Wilner

Abstract: Modeling the submillimeter to centimeter emission of stars is challenging due to a lack of sensitive observations at these long wavelengths. We launched an ongoing campaign to obtain new observations entitled Measuring the Emission of Stellar Atmospheres at Submillimeter/millimeter wavelengths (MESAS). Here we present ALMA, GBT, and VLA observations of Sirius A, the closest main-sequence A-type st… ▽ More Modeling the submillimeter to centimeter emission of stars is challenging due to a lack of sensitive observations at these long wavelengths. We launched an ongoing campaign to obtain new observations entitled Measuring the Emission of Stellar Atmospheres at Submillimeter/millimeter wavelengths (MESAS). Here we present ALMA, GBT, and VLA observations of Sirius A, the closest main-sequence A-type star, that span from 1.4 to 9.0 millimeters. These observations complement our previous millimeter data on Sirius A and are entirely consistent with the PHOENIX stellar atmosphere models constructed to explain them. We note that accurate models of long wavelength emission from stars are essential not only to understand fundamental stellar processes, but also to determine the presence of dusty debris in spatially unresolved observations of circumstellar disks. △ Less

Submitted 8 March, 2019; originally announced March 2019.

Comments: 7 pages, 1 figure, accepted for publication in ApJ

Science with an ngVLA: Stellar Emission as a Source of Flux Bias in Debris Disks

Authors: Jacob Aaron White, Jason Aufdenberg, Aaron C. Boley, Peter Hauschildt, A. Meredith Hughes, Brenda Matthews, Attila Moór, David J. Wilner

Abstract: Our understanding of stellar atmospheres and our ability to infer architectures of extrasolar planetary systems rely on understanding the emission of stars at submillimeter to centimeter wavelengths. In this chapter we describe how unconstrained stellar emission can interfere with the accurate characterization of circumstellar debris. The ngVLA is the only facility with the sensitivity that allows… ▽ More Our understanding of stellar atmospheres and our ability to infer architectures of extrasolar planetary systems rely on understanding the emission of stars at submillimeter to centimeter wavelengths. In this chapter we describe how unconstrained stellar emission can interfere with the accurate characterization of circumstellar debris. The ngVLA is the only facility with the sensitivity that allows for the observations of a broad range of stellar spectral types in a feasible amount of time. The observations will enable the building and testing of accurate models of stellar emission, which in turn are required for evaluating both the occurrence and abundance of debris over the proposed wavelength range of the ngVLA △ Less

Submitted 15 October, 2018; originally announced October 2018.

Comments: 6 pages, 2 figures, To be published in the ASP Monograph Series, "Science with a Next-Generation VLA", ed. E. J. Murphy (ASP, San Francisco, CA)

Gravitational instabilities in a protosolar-like disc III: molecular line detection and sensitivities

Authors: M. G. Evans, T. W. Hartquist, P. Caselli, A. C. Boley, J. D. Ilee, J. M. C. Rawlings

Abstract: At the time of formation, protoplanetary discs likely contain a comparable mass to their host protostars. As a result, gravitational instabilities (GIs) are expected to play a pivotal role in the early phases of disc evolution. However, as these young objects are heavily embedded, confirmation of GIs has remained elusive. Therefore, we use the radiative transfer code LIME to produce line images of… ▽ More At the time of formation, protoplanetary discs likely contain a comparable mass to their host protostars. As a result, gravitational instabilities (GIs) are expected to play a pivotal role in the early phases of disc evolution. However, as these young objects are heavily embedded, confirmation of GIs has remained elusive. Therefore, we use the radiative transfer code LIME to produce line images of a $0.17\,\mathrm{M}_{\odot}$ self-gravitating protosolar-like disc. We note the limitations of using LIME and explore methods to improve upon the default gridding. We synthesise noiseless observations to determine the sensitivities required to detect the spiral flux, and find that the line flux distribution does not necessarily correlate to the abundance density distribution; hence performing radiative transfer calculations is imperative. Moreover, the spiral features are seen in absorption, due to the GI-heated midplane and high extinction, which could be indicative of GI activity. If a small beamsize and appropriate molecular line are used then spatially resolving spirals in a protosolar-like disc should be possible with ALMA for an on-source time of 30 hr. Spectrally resolving non-axisymmetric structure takes only a tenth as long for a reasonable noise level, but attributing this structure to GI-induced activity would be tentative. Finally, we find that identifying finger-like features in PV diagrams of nearly edge-on discs, which are a direct indicator of spirals, is feasible with an on-source time of 19 hr, and hence likely offers the most promising means of confirming GI-driven spiral structure in young, embedded protoplanetary discs. △ Less

Submitted 10 October, 2018; originally announced October 2018.

Comments: Accepted for publication in MNRAS; 23 pages, 25 figures and 2 tables

Extended Millimeter Emission in the HD 141569 Circumstellar Disk Detected with ALMA

Authors: Jacob Aaron White, Aaron C. Boley

Abstract: We present archival ALMA observations of the HD 141569 circumstellar disk at 345, 230, and 100 GHz. These data detect extended millimeter emission that is exterior to the inner disk. We find through simultaneous visibility modeling of all three data sets that the system's morphology is described well by a two-component disk model. The inner disk ranges from approximately 16 to 45 au with a spectra… ▽ More We present archival ALMA observations of the HD 141569 circumstellar disk at 345, 230, and 100 GHz. These data detect extended millimeter emission that is exterior to the inner disk. We find through simultaneous visibility modeling of all three data sets that the system's morphology is described well by a two-component disk model. The inner disk ranges from approximately 16 to 45 au with a spectral index of 1.81 (q = 2.95) and the outer disk ranges from 95 to 300 au with a spectral index of 2.28 (q = 3.21). Azimuthally averaged radial emission profiles derived from the continuum images at each frequency show potential emission that is consistent with the visibility modeling. The analysis presented here shows that at ~5 Myr HD 141569's grain size distribution is steeper, and therefore more evolved, in the outer disk than in the inner disk. △ Less

Submitted 25 April, 2018; originally announced April 2018.

Comments: 11 pages, 4 figures, accepted to ApJ on April 23rd 2018

Magnetic Fields Recorded by Chondrules Formed in Nebular Shocks

Authors: Chuhong Mai, Steven J. Desch, Aaron C. Boley, Benjamin P. Weiss

Abstract: Recent laboratory efforts (Fu et al., 2014) have constrained the remanent magnetizations of chondrules and the magnetic field strengths at which the chondrules were exposed to as they cooled below their Curie points. An outstanding question is whether the inferred paleofields represent the background magnetic field of the solar nebula or were unique to the chondrule-forming environment. We investi… ▽ More Recent laboratory efforts (Fu et al., 2014) have constrained the remanent magnetizations of chondrules and the magnetic field strengths at which the chondrules were exposed to as they cooled below their Curie points. An outstanding question is whether the inferred paleofields represent the background magnetic field of the solar nebula or were unique to the chondrule-forming environment. We investigate the amplification of the magnetic field above background values for two proposed chondrule formation mechanisms, large-scale nebular shocks and planetary bow shocks. Behind large-scale shocks, the magnetic field parallel to the shock front is amplified by factors $\sim 10-30$, regardless of the magnetic diffusivity. Therefore, chondrules melted in these shocks probably recorded an amplified magnetic field. Behind planetary bow shocks, the field amplification is sensitive to the magnetic diffusivity. We compute the gas properties behind a bow shock around a 3000 km-radius planetary embryo, with and without atmospheres, using hydrodynamics models. We calculate the ionization state of the hot, shocked gas, including thermionic emission from dust, and thermal ionization of gas-phase potassium atoms, and the magnetic diffusivity due to Ohmic dissipation and ambipolar diffusion. We find that the diffusivity is sufficiently large that magnetic fields have already relaxed to background values in the shock downstream where chondrules acquire magnetizations, and that these locations are sufficiently far from the planetary embryos that chondrules should not have recorded a significant putative dynamo field generated on these bodies. We conclude that, if melted in planetary bow shocks, chondrules probably recorded the background nebular field. △ Less

Submitted 13 March, 2018; originally announced March 2018.

Comments: 17 pages, 11 figures, accepted for publication in ApJ

The Dynamics of Tightly-packed Planetary Systems in the Presence of an Outer Planet: case studies using Kepler-11 and Kepler-90

Authors: A. P. Granados Contreras, A. C. Boley

Abstract: We explore the effects of an undetected outer giant planet on the dynamics, observability, and stability of Systems with Tightly-packed Inner Planets (STIPs). We use direct numerical simulations along with secular theory and synthetic secular frequency spectra to analyze how analogues of Kepler-11 and Kepler-90 behave in the presence of a nearly co-planar, Jupiter-like outer perturber with semi-ma… ▽ More We explore the effects of an undetected outer giant planet on the dynamics, observability, and stability of Systems with Tightly-packed Inner Planets (STIPs). We use direct numerical simulations along with secular theory and synthetic secular frequency spectra to analyze how analogues of Kepler-11 and Kepler-90 behave in the presence of a nearly co-planar, Jupiter-like outer perturber with semi-major axes between 1 and 5.2 au. Most locations of the outer perturber do not affect the evolution of the inner planetary systems, apart from altering precession frequencies. However, there are locations at which an outer planet causes system instability due to, in part, secular eccentricity resonances. In Kepler-90, there is a range of orbital distances for which the outer perturber drives planets b and c, through secular interactions, onto orbits with inclinations that are $\sim16^\circ$ away from the rest of the planets. Kepler-90 is stable in this configuration. Such secular resonances can thus affect the observed multiplicity of transiting systems. We also compare the synthetic apsidal and nodal precession frequencies with the secular theory and find some misalignment between principal frequencies, indicative of strong interactions between the planets (consistent with the system showing TTVs). First-order libration angles are calculated to identify MMRs in the systems, for which two near-MMRs are shown in Kepler-90, with a 5:4 between b and c, as well as a 3:2 between g and h. △ Less

Submitted 5 February, 2018; originally announced February 2018.

Comments: 17 pages, accepted for publication in AJ

ALMA and VLA Observations of the HD 141569 System

Authors: Jacob Aaron White, A. C. Boley, M. A. MacGregor, A. M. Hughes, D. J. Wilner

Abstract: We present VLA 9 mm (33 GHz) observations of the HD 141569 system from semester 16A. The observations achieve a resolution of 0.25 arcsec ($\sim28$ au) and a sensitivity of $4.7~μ\rm Jy~beam^{-1}$. We find (1) a $52\pm 5~μ$Jy point source at the location of HD 141569A that shows potential variability, (2) the detected flux is contained within the SED-inferred central clearing of the disc meaning t… ▽ More We present VLA 9 mm (33 GHz) observations of the HD 141569 system from semester 16A. The observations achieve a resolution of 0.25 arcsec ($\sim28$ au) and a sensitivity of $4.7~μ\rm Jy~beam^{-1}$. We find (1) a $52\pm 5~μ$Jy point source at the location of HD 141569A that shows potential variability, (2) the detected flux is contained within the SED-inferred central clearing of the disc meaning the spectral index of the dust disc is steeper than previously inferred, and (3) the M dwarf companions are also detected and variable. Previous lower-resolution VLA observations (semester 14A) found a higher flux density, interpreted as solely dust emission. When combined with ALMA observations, the VLA 14A observations suggested the spectral index and grain size distribution of HD 141569's disc was shallow and an outlier among debris systems. Using archival ALMA observations of HD 141569 at 0.87 mm and 2.9 mm we find a dust spectral index of $α_{\rm mm} = 1.81\pm 0.20$. The VLA 16A flux corresponds to a brightness temperature of $\sim5\times10^{6}$ K, suggesting strong non-disc emission is affecting the inferred grain properties. The VLA 16A flux density of the M2V companion HD 141569B is $149\pm9~μ$Jy, corresponding to a brightness temperature of $\sim2\times10^{8}$ K and suggesting significant stellar variability when compared to the VLA14A observations, which are smaller by a factor of $\sim6$. △ Less

Submitted 20 November, 2017; originally announced November 2017.

Comments: Accepted for publication in MNRAS, 8 pages, 6 figures

On the Origin of Banded Structure in Dusty Protoplanetary Discs: HL Tau and TW Hya

Authors: Aaron C. Boley

Abstract: Recent observations of HL Tau revealed remarkably detailed structure within the system's circumstellar disc. A range of hypotheses have been proposed to explain the morphology, including, e.g., planet-disc interactions, condensation fronts, and secular gravitational instabilities. While embedded planets seem to be able to explain some of the major structure in the disc through interactions with ga… ▽ More Recent observations of HL Tau revealed remarkably detailed structure within the system's circumstellar disc. A range of hypotheses have been proposed to explain the morphology, including, e.g., planet-disc interactions, condensation fronts, and secular gravitational instabilities. While embedded planets seem to be able to explain some of the major structure in the disc through interactions with gas and dust, the substructure, such as low-contrast rings and bands, are not so easily reproduced. Here, we show that dynamical interactions between three planets (only two of which are modelled) and an initial population of large planetesimals can potentially explain both the major and minor banded features within the system. In this context, the small grains, which are coupled to the gas and reveal the disc morphology, are produced by the collisional evolution of the newly-formed planetesimals, which are ubiquitous in the system and are decoupled from the gas. △ Less

Submitted 14 October, 2017; originally announced October 2017.

Comments: To appear in ApJ

The chemistry of protoplanetary fragments formed via gravitational instabilities

Authors: J. D. Ilee, D. H. Forgan, M. G. Evans, C. Hall, R. Booth, C. J. Clarke, W. K. M. Rice, A. C. Boley, P. Caselli, T. W. Hartquist, J. M. C. Rawlings

Abstract: In this paper, we model the chemical evolution of a 0.25 M$_{\odot}$ protoplanetary disc surrounding a 1 M$_{\odot}$ star that undergoes fragmentation due to self-gravity. We use Smoothed Particle Hydrodynamics including a radiative transfer scheme, along with time-dependent chemical evolution code to follow the composition of the disc and resulting fragments over approximately 4000 yrs. Initially… ▽ More In this paper, we model the chemical evolution of a 0.25 M$_{\odot}$ protoplanetary disc surrounding a 1 M$_{\odot}$ star that undergoes fragmentation due to self-gravity. We use Smoothed Particle Hydrodynamics including a radiative transfer scheme, along with time-dependent chemical evolution code to follow the composition of the disc and resulting fragments over approximately 4000 yrs. Initially, four quasi-stable fragments are formed, of which two are eventually disrupted by tidal torques in the disc. From the results of our chemical modelling, we identify species that are abundant in the fragments (e.g. H$_{\rm 2}$O, H$_{\rm 2}$S, HNO, N$_{\rm 2}$, NH$_{\rm 3}$, OCS, SO), species that are abundant in the spiral shocks within the disc (e.g. CO, CH$_{\rm 4}$, CN, CS, H$_{\rm 2}$CO), and species which are abundant in the circumfragmentary material (e.g. HCO$^{\rm +}$). Our models suggest that in some fragments it is plausible for grains to sediment to the core before releasing their volatiles into the planetary envelope, leading to changes in, e.g., the C/O ratio of the gas and ice components. We would therefore predict that the atmospheric composition of planets generated by gravitational instability should not necessarily follow the bulk chemical composition of the local disc material. △ Less

Submitted 5 August, 2017; originally announced August 2017.

Comments: 17 pages, 7 figures, accepted for publication in MNRAS

Simulations of Small Solid Accretion onto Planetesimals in the Presence of Gas

Authors: A. Hughes, A. C. Boley

Abstract: The growth and migration of planetesimals in a young protoplanetary disc are fundamental to planet formation. In all models of early growth, there are several processes that can inhibit grains from reaching larger sizes. Nevertheless, observations suggest that growth of planetesimals must be rapid. If a small number of 100 km sized planetesimals do manage to form in the disc, then gas drag effects… ▽ More The growth and migration of planetesimals in a young protoplanetary disc are fundamental to planet formation. In all models of early growth, there are several processes that can inhibit grains from reaching larger sizes. Nevertheless, observations suggest that growth of planetesimals must be rapid. If a small number of 100 km sized planetesimals do manage to form in the disc, then gas drag effects could enable them to efficiently accrete small solids from beyond their gravitationally focused cross-section. This gas drag-enhanced accretion can allow planetesimals to grow at rapid rates, in principle. We present self-consistent hydrodynamics simulations with direct particle integration and gas drag coupling to estimate the rate of planetesimal growth due to pebble accretion. Wind tunnel simulations are used to explore a range of particle sizes and disc conditions. We also explore analytic estimates of planetesimal growth and numerically integrate planetesimal drift due to the accretion of small solids. Our results show that, for almost every case that we consider, there is a clearly preferred particle size for accretion that depends on the properties of the accreting planetesimal and the local disc conditions. For solids much smaller than the preferred particle size, accretion rates are significantly reduced as the particles are entrained in the gas and flow around the planetesimal. Solids much larger than the preferred size accrete at rates consistent with gravitational focusing. Our analytic estimates for pebble accretion highlight the timescales that are needed for the growth of large objects under different disc conditions and initial planetesimal sizes. △ Less

Submitted 1 August, 2017; originally announced August 2017.

Comments: To appear in MNRAS

Gravitational instabilities in a protosolar-like disc II: continuum emission and mass estimates

Authors: M. G. Evans, J. D. Ilee, T. W. Hartquist, P. Caselli, L. Szucs, S. J. D. Purser, A. C. Boley, R. H. Durisen, J. M. C. Rawlings

Abstract: Gravitational instabilities (GIs) are most likely a fundamental process during the early stages of protoplanetary disc formation. Recently, there have been detections of spiral features in young, embedded objects that appear consistent with GI-driven structure. It is crucial to perform hydrodynamic and radiative transfer simulations of gravitationally unstable discs in order to assess the validity… ▽ More Gravitational instabilities (GIs) are most likely a fundamental process during the early stages of protoplanetary disc formation. Recently, there have been detections of spiral features in young, embedded objects that appear consistent with GI-driven structure. It is crucial to perform hydrodynamic and radiative transfer simulations of gravitationally unstable discs in order to assess the validity of GIs in such objects, and constrain optimal targets for future observations. We utilise the radiative transfer code LIME to produce continuum emission maps of a $0.17\,\mathrm{M}_{\odot}$ self-gravitating protosolar-like disc. We note the limitations of using LIME as is and explore methods to improve upon the default gridding. We use CASA to produce synthetic observations of 270 continuum emission maps generated across different frequencies, inclinations and dust opacities. We find that the spiral structure of our protosolar-like disc model is distinguishable across the majority of our parameter space after 1 hour of observation, and is especially prominent at 230$\,$GHz due to the favourable combination of angular resolution and sensitivity. Disc mass derived from the observations is sensitive to the assumed dust opacities and temperatures, and therefore can be underestimated by a factor of at least 30 at 850$\,$GHz and 2.5 at 90$\,$GHz. As a result, this effect could retrospectively validate GIs in discs previously thought not massive enough to be gravitationally unstable, which could have a significant impact on the understanding of the formation and evolution of protoplanetary discs. △ Less

Submitted 1 June, 2017; originally announced June 2017.

Comments: Accepted for publication in MNRAS; 22 pages, 24 figures and 1 table

1.3 mm ALMA Observations of the Fomalhaut Debris System

Authors: J. A. White, A. C. Boley, W. R. F. Dent, E. B. Ford, S. Corder

Abstract: We present ALMA Band 6 observations (1.3 mm/233 GHz) of Fomalhaut and its debris disc. The observations achieve a sensitivity of 17 $μ$Jy and a resolution of 0.28 arcsec (2.1 au at a distance of 7.66 pc), which are the highest resolution observations to date of the millimetre grains in Fomalhaut's main debris ring. The ring is tightly constrained to $139^{+2}_{-3}$ au with a FWHM of $13\pm3$ au, f… ▽ More We present ALMA Band 6 observations (1.3 mm/233 GHz) of Fomalhaut and its debris disc. The observations achieve a sensitivity of 17 $μ$Jy and a resolution of 0.28 arcsec (2.1 au at a distance of 7.66 pc), which are the highest resolution observations to date of the millimetre grains in Fomalhaut's main debris ring. The ring is tightly constrained to $139^{+2}_{-3}$ au with a FWHM of $13\pm3$ au, following a Gaussian profile. The millimetre spectral index is constrained to $α_{mm} = -2.62\pm0.12$. We explore fitting debris disc models in the image plane, as well as fitting models using visibility data directly. The results are compared and the potential advantages/disadvantages of each approach are discussed. The detected central emission is indistinguishable from a point source, with a most probable flux of $0.90\pm 0.12$ mJy (including calibration uncertainties). This implies that any inner debris structure, as was inferred from far-Infrared observations, must contribute little to the total central emission. Moreover, the stellar flux is less than 70\% of that predicted by extrapolating a black body from the constrained stellar photosphere temperature. This result emphasizes that unresolved inner debris components cannot be fully characterized until the behaviour of the host star's intrinsic stellar emission at millimetre wavelengths is properly understood. △ Less

Submitted 5 December, 2016; originally announced December 2016.

Comments: 10 pages, 7 figures, Submitted to MNRAS after addressing referee's comments

ALMA Observations of HD141569's Circumstellar Disk

Authors: J. A. White, A. C. Boley, A. M. Hughes, K. M. Flaherty, E. Ford, D. Wilner, S. Corder, M. Payne

Abstract: We present ALMA band 7 (345 GHz) continuum and $^{12}$CO(J = 3-2) observations of the circumstellar disk surrounding HD141569. At an age of about 5 Myr, the disk has a complex morphology that may be best interpreted as a nascent debris system with gas. Our $870\rm~μm$ ALMA continuum observations resolve a dust disk out to approximately $ 56 ~\rm au$ from the star (assuming a distance of 116 pc) wi… ▽ More We present ALMA band 7 (345 GHz) continuum and $^{12}$CO(J = 3-2) observations of the circumstellar disk surrounding HD141569. At an age of about 5 Myr, the disk has a complex morphology that may be best interpreted as a nascent debris system with gas. Our $870\rm~μm$ ALMA continuum observations resolve a dust disk out to approximately $ 56 ~\rm au$ from the star (assuming a distance of 116 pc) with $0."38$ resolution and $0.07 ~ \rm mJy~beam^{-1}$ sensitivity. We measure a continuum flux density for this inner material of $3.8 \pm 0.4 ~ \rm mJy$ (including calibration uncertainties). The $^{12}$CO(3-2) gas is resolved kinematically and spatially from about 30 to 210 au. The integrated $^{12}$CO(3-2) line flux density is $15.7 \pm 1.6~\rm Jy~km~s^{-1}$. We estimate the mass of the millimeter debris and $^{12}$CO(3-2) gas to be $\gtrsim0.04~\rm M_{\oplus}$ and $\sim2\times 10^{-3}~\rm M_{\oplus}$, respectively. If the millimeter grains are part of a collisional cascade, then we infer that the inner disk ($<50$ au) has $\sim 160~\rm M_{\oplus}$ contained within objects less than 50 km in radius, depending on the planetesimal size distribution and density assumptions. MCMC modeling of the system reveals a disk morphology with an inclination of $53.4^{\circ}$ centered around a $\rm M=2.39~ M_{\odot}$ host star ($\rm Msin(i)=1.92~ M_{\odot}$). We discuss whether the gas in HD141569's disk may be second generation. If it is, the system can be used to study the clearing stages of planet formation. △ Less

Submitted 1 June, 2016; originally announced June 2016.

Comments: 35 pages, 7 figures, submitted to ApJ

Journal ref: The Astrophysical Journal, 829:6 (11pp), 2016 September 20

Planetary Embryo Bow Shocks as a Mechanism for Chondrule Formation

Authors: Christopher R. Mann, Aaron C. Boley, Melissa M. Morris

Abstract: We use radiation hydrodynamics with direct particle integration to explore the feasibility of chondrule formation in planetary embryo bow shocks. The calculations presented here are used to explore the consequences of a Mars-size planetary embryo traveling on a moderately excited orbit through the dusty, early environment of the solar system. The embryo's eccentric orbit produces a range of supers… ▽ More We use radiation hydrodynamics with direct particle integration to explore the feasibility of chondrule formation in planetary embryo bow shocks. The calculations presented here are used to explore the consequences of a Mars-size planetary embryo traveling on a moderately excited orbit through the dusty, early environment of the solar system. The embryo's eccentric orbit produces a range of supersonic relative velocities between the embryo and the circularly orbiting gas and dust, prompting the formation of bow shocks. Temporary atmospheres around these embryos, which can be created via volatile outgassing and gas capture from the surrounding nebula, can non-trivially affect thermal profiles of solids entering the shock. We explore the thermal environment of solids that traverse the bow shock at different impact radii, the effects that planetoid atmospheres have on shock morphologies, and the stripping efficiency of planetoidal atmospheres in the presence of high relative winds. Simulations are run using adiabatic and radiative conditions, with multiple treatments for the local opacities. Shock speeds of 5, 6, and 7 km/s are explored. We find that a high-mass atmosphere and inefficient radiative conditions can produce peak temperatures and cooling rates that are consistent with the constraints set by chondrule furnace studies. For most conditions, the derived cooling rates are potentially too high to be consistent with chondrule formation. △ Less

Submitted 16 May, 2016; originally announced May 2016.

Comments: 28 pages, 19 figures

MSC Class: 85-04

Journal ref: The Astrophysical Journal, Volume 818, Issue 2, article id. 103, 20 pp. (2016)

Resolved Gas Interior to the Dust Rings of the HD 141569 Disk

Authors: Kevin M. Flaherty, A. Meredith Hughes, Sean M. Andrews, Chunhua Qi, David J Wilner, Aaron C. Boley, Jacob A. White, Will Harney, Julia Zachary

Abstract: The disk around HD 141569 is one of a handful of systems whose weak infrared emission is consistent with a debris disk, but still has a significant reservoir of gas. Here we report spatially resolved mm observations of the CO(3-2) and CO(1-0) emission as seen with the SMA and CARMA. We find that the excitation temperature for CO is lower than expected from cospatial blackbody grains, similar to pr… ▽ More The disk around HD 141569 is one of a handful of systems whose weak infrared emission is consistent with a debris disk, but still has a significant reservoir of gas. Here we report spatially resolved mm observations of the CO(3-2) and CO(1-0) emission as seen with the SMA and CARMA. We find that the excitation temperature for CO is lower than expected from cospatial blackbody grains, similar to previous observations of analogous systems, and derive a gas mass that lies between that of gas-rich primordial disks and gas-poor debris disks. The data also indicate a large inner hole in the CO gas distribution and an outer radius that lies interior to the outer scattered light rings. This spatial distribution, with the dust rings just outside the gaseous disk, is consistent with the expected interactions between gas and dust in an optically thin disk. This indicates that gas can have a significant effect on the location of the dust within debris disks. △ Less

Submitted 11 January, 2016; originally announced January 2016.

Comments: Accepted to ApJ

The In Situ Formation of Giant Planets at Short Orbital Periods

Authors: A. C. Boley, A. P. Granados Contreras, B. Gladman

Abstract: We propose that two of the most surprising results so far among exoplanet discoveries are related: the existences of both hot Jupiters and the high frequency of multi-planet systems with periods $P\lesssim200$~days. In this paradigm, the vast majority of stars rapidly form along with multiple close-in planets in the mass range of Mars to super-Earths/mini-Neptunes. Such systems of tightly-packed i… ▽ More We propose that two of the most surprising results so far among exoplanet discoveries are related: the existences of both hot Jupiters and the high frequency of multi-planet systems with periods $P\lesssim200$~days. In this paradigm, the vast majority of stars rapidly form along with multiple close-in planets in the mass range of Mars to super-Earths/mini-Neptunes. Such systems of tightly-packed inner planets (STIPs) are metastable, with the time scale of the dynamical instability having a major influence on final planet types. In most cases, the planets consolidate into a system of fewer, more massive planets, but long after the circumstellar gas disk has dissipated. This can yield planets with masses above the traditional critical core of $\sim$10 $M_\oplus$, yielding short-period giants that lack abundant gas. A rich variety of physical states are also possible given the range of collisional outcomes and formation time of the close-in planets. However, when dynamical consolidation occurs before gas dispersal, a critical core can form that then grows via gas capture into a short-period gas giant. In this picture the majority of Hot and Warm Jupiters formed locally, rather than migrating down from larger distances. △ Less

Submitted 8 December, 2015; v1 submitted 14 October, 2015; originally announced October 2015.

Comments: Accepted for publication in ApJL

Gravitational instabilities in a protosolar-like disc I: dynamics and chemistry

Authors: M. G. Evans, J. D. Ilee, A. C. Boley, P. Caselli, R. H. Durisen, T. W. Hartquist, J. M. C. Rawlings

Abstract: To date, most simulations of the chemistry in protoplanetary discs have used 1+1D or 2D axisymmetric $α$-disc models to determine chemical compositions within young systems. This assumption is inappropriate for non-axisymmetric, gravitationally unstable discs, which may be a significant stage in early protoplanetary disc evolution. Using 3D radiative hydrodynamics, we have modelled the physical an… ▽ More To date, most simulations of the chemistry in protoplanetary discs have used 1+1D or 2D axisymmetric $α$-disc models to determine chemical compositions within young systems. This assumption is inappropriate for non-axisymmetric, gravitationally unstable discs, which may be a significant stage in early protoplanetary disc evolution. Using 3D radiative hydrodynamics, we have modelled the physical and chemical evolution of a 0.17 M$_{\odot}$ self-gravitating disc over a period of 2000 yr. The 0.8 M$_{\odot}$ central protostar is likely to evolve into a solar-like star, and hence this Class 0 or early Class I young stellar object may be analogous to our early Solar System. Shocks driven by gravitational instabilities enhance the desorption rates, which dominate the changes in gas-phase fractional abundances for most species. We find that at the end of the simulation, a number of species distinctly trace the spiral structure of our relatively low-mass disc, particularly CN. We compare our simulation to that of a more massive disc, and conclude that mass differences between gravitationally unstable discs may not have a strong impact on the chemical composition. We find that over the duration of our simulation, successive shock heating has a permanent effect on the abundances of HNO, CN and NH$_3$, which may have significant implications for both simulations and observations. We also find that HCO$^+$ may be a useful tracer of disc mass. We conclude that gravitational instabilities induced in lower mass discs can significantly, and permanently, affect the chemical evolution, and that observations with high-resolution instruments such as ALMA offer a promising means of characterising gravitational instabilities in protosolar discs. △ Less

Submitted 4 August, 2015; originally announced August 2015.

Comments: Accepted for publication in MNRAS; 19 pages, 17 figures and 5 tables

Overcoming the Meter Barrier and The Formation of Systems with Tightly-packed Inner Planets (STIPs)

Authors: Aaron C. Boley, Melissa A. Morris, Eric B. Ford

Abstract: We present a solution to the long outstanding meter barrier problem in planet formation theory. As solids spiral inward due to aerodynamic drag, they will enter disk regions that are characterized by high temperatures, densities, and pressures. High partial pressures of rock vapor can suppress solid evaporation, and promote collisions between partially molten solids, allowing rapid growth. This pr… ▽ More We present a solution to the long outstanding meter barrier problem in planet formation theory. As solids spiral inward due to aerodynamic drag, they will enter disk regions that are characterized by high temperatures, densities, and pressures. High partial pressures of rock vapor can suppress solid evaporation, and promote collisions between partially molten solids, allowing rapid growth. This process should be ubiquitous in planet-forming disks, which may be evidenced by the abundant class of Systems with Tightly-packed Inner Planets (STIPs) discovered by the NASA Kepler mission. △ Less

Submitted 24 July, 2014; originally announced July 2014.

Comments: To appear in ApJL

Core-assisted gas capture instability: a new mode of giant planet formation by gravitationally unstable discs

Authors: Sergei Nayakshin, Ravit Helled, Aaron C. Boley

Abstract: Giant planet formation in the core accretion (CA) paradigm is predicated by the formation of a core, assembled by the coagulation of grains and later by planetesimals within a protoplanetary disc. In contrast, in the disc instability paradigm, giant planet formation is believed to be independent of core formation: massive self-gravitating gas fragments cool radiatively and collapse as a whole. We… ▽ More Giant planet formation in the core accretion (CA) paradigm is predicated by the formation of a core, assembled by the coagulation of grains and later by planetesimals within a protoplanetary disc. In contrast, in the disc instability paradigm, giant planet formation is believed to be independent of core formation: massive self-gravitating gas fragments cool radiatively and collapse as a whole. We show that giant planet formation in the disc instability model may be also enhanced by core formation for reasons physically very similar to the CA paradigm. In the model explored here, efficient grain sedimentation within an initial fragment (rather than the disc) leads to the formation of a core composed of heavy elements. We find that massive atmospheres form around cores and undergo collapse as a critical core mass is exceeded, analogous to CA theory. The critical mass of the core to initiate such a collapse depends on the fragment mass and metallicity, as well as core luminosity, but ranges from less than 1 to as much as $\sim80$ Earth masses. We therefore suggest that there are two channels for the collapse of a gaseous fragment to planetary scales within the disc instability model: (i) H$_2$ dissociative collapse of the entire gaseous clump, and (ii) core-assisted gas capture, as presented here. We suggest that the first of these two is favoured in metal-poor environments and for fragments more massive than $\sim 5-10$ Jupiter masses, whereas the second is favored in metal-rich environments and fragments of lower mass. [Abridged] △ Less

Submitted 7 March, 2014; originally announced March 2014.

Comments: Accepted to MNRAS. 13 pages

High-Temperature Processing of Solids Through Solar Nebular Bow Shocks: 3D Radiation Hydrodynamics Simulations with Particles

Authors: A. C. Boley, M. A. Morris, S. J. Desch

Abstract: A fundamental, unsolved problem in Solar System formation is explaining the melting and crystallization of chondrules found in chondritic meteorites. Theoretical models of chondrule melting in nebular shocks has been shown to be consistent with many aspects of thermal histories inferred for chondrules from laboratory experiments; but, the mechanism driving these shocks is unknown. Planetesimals an… ▽ More A fundamental, unsolved problem in Solar System formation is explaining the melting and crystallization of chondrules found in chondritic meteorites. Theoretical models of chondrule melting in nebular shocks has been shown to be consistent with many aspects of thermal histories inferred for chondrules from laboratory experiments; but, the mechanism driving these shocks is unknown. Planetesimals and planetary embryos on eccentric orbits can produce bow shocks as they move supersonically through the disk gas, and are one possible source of chondrule-melting shocks. We investigate chondrule formation in bow shocks around planetoids through 3D radiation hydrodynamics simulations. A new radiation transport algorithm that combines elements of flux-limited diffusion and Monte Carlo methods is used to capture the complexity of radiative transport around bow shocks. An equation of state that includes the rotational, vibrational, and dissociation modes of H$_2$ is also used. Solids are followed directly in the simulations and their thermal histories are recorded. Adiabatic expansion creates rapid cooling of the gas, and tail shocks behind the embryo can cause secondary heating events. Radiative transport is efficient, and bow shocks around planetoids can have luminosities $\sim$few$\times10^{-8}$ L$_{\odot}$. While barred and radial chondrule textures could be produced in the radiative shocks explored here, porphyritic chondrules may only be possible in the adiabatic limit. We present a series of predicted cooling curves that merit investigation in laboratory experiments to determine whether the solids produced by bow shocks are represented in the meteoritic record by chondrules or other solids. △ Less

Submitted 27 August, 2013; v1 submitted 13 August, 2013; originally announced August 2013.

Comments: Accepted for publication in ApJ. Images have been resized to conform to arXiv limits, but are all readable upon adjusting the zoom. Changes from v1: Corrected typos discovered in proofs. Most changes are in the appendix

The Formation of Systems with Tightly-packed Inner Planets (STIPs) via Aerodynamic Drift

Authors: Aaron C. Boley, Eric B. Ford

Abstract: The NASA Kepler mission has revealed an abundant class of Systems with Tightly-packed Inner Planets (STIPs). The current paradigm for planet formation suggests that small planetesimals will quickly spiral into the host star due to aerodynamic drag, preventing rocky planet formation. In contrast, we find that aerodynamic drift, when acting on an ensemble of solids, can concentrate mass at short orb… ▽ More The NASA Kepler mission has revealed an abundant class of Systems with Tightly-packed Inner Planets (STIPs). The current paradigm for planet formation suggests that small planetesimals will quickly spiral into the host star due to aerodynamic drag, preventing rocky planet formation. In contrast, we find that aerodynamic drift, when acting on an ensemble of solids, can concentrate mass at short orbital periods in gaseous disks. Sublimation fronts may further aid this process. Kepler data suggest that the innermost known planets are found near the silicate sublimation zone. STIP planets should have a wide range of volatile fractions due to aerodynamic drift and H2 dissociation-driven gas accretion. We further propose that the low mass of Mars is evidence that the Solar System was once a proto-STIP. △ Less

Submitted 3 June, 2013; originally announced June 2013.

Comments: Submitted. Comments welcome

Simulated Observations of Young Gravitationally Unstable Protoplanetary Discs

Authors: T. A. Douglas, P. Caselli, J. D. Ilee, A. C. Boley, T. W. Hartquist, R. H. Durisen, J. M. C. Rawlings

Abstract: The formation and earliest stages of protoplanetary discs remain poorly constrained by observations. ALMA will soon revolutionise this field. Therefore, it is important to provide predictions which will be valuable for the interpretation of future high sensitivity and high angular resolution observations. Here we present simulated ALMA observations based on radiative transfer modelling of a relati… ▽ More The formation and earliest stages of protoplanetary discs remain poorly constrained by observations. ALMA will soon revolutionise this field. Therefore, it is important to provide predictions which will be valuable for the interpretation of future high sensitivity and high angular resolution observations. Here we present simulated ALMA observations based on radiative transfer modelling of a relatively massive (0.39 M_solar) self-gravitating disc embedded in a 10 M_solar dense core, with structure similar to the pre-stellar core L1544. We focus on simple species and conclude that C17O 3-2, HCO+ 3-2, OCS 26-25 and H2CO 404-303 lines can be used to probe the disc structure and kinematics at all scales. △ Less

Submitted 16 May, 2013; originally announced May 2013.

Comments: 12 pages, 15 figures, Accepted by MNRAS

The collapse of protoplanetary clumps formed through disc instability: 3D simulations of the pre-dissociation phase

Authors: M. Galvagni, T. Hayfield, A. C. Boley, L. Mayer, R. Roskar, P. Saha

Abstract: We present 3D smoothed particle hydrodynamics simulations of the collapse of clumps formed through gravitational instability in the outer part of a protoplanetary disc. The initial conditions are taken directly from a global disc simulation, and a realistic equation of state is used to follow the clumps as they contract over several orders of magnitude in density, approaching the molecular hydroge… ▽ More We present 3D smoothed particle hydrodynamics simulations of the collapse of clumps formed through gravitational instability in the outer part of a protoplanetary disc. The initial conditions are taken directly from a global disc simulation, and a realistic equation of state is used to follow the clumps as they contract over several orders of magnitude in density, approaching the molecular hydrogen dissociation stage. The effects of clump rotation, asymmetries, and radiative cooling are studied. Rotation provides support against fast collapse, but non-axisymmetric modes develop and efficiently transport angular momentum outward, forming a circumplanetary disc. This transport helps the clump reach the dynamical collapse phase, resulting from molecular hydrogen dissociation, on a thousand-year timescale, which is smaller than timescales predicted by some previous spherical 1D collapse models. Extrapolation to the threshold of the runaway hydrogen dissociation indicates that the collapse timescales can be shorter than inward migration timescales, suggesting that clumps could survive tidal disruption and deliver a proto-gas giant to distances of even a few AU from the central star. △ Less

Submitted 10 September, 2012; originally announced September 2012.

Comments: Accepted for publication in MNRAS

Swarm-NG: a CUDA Library for Parallel n-body Integrations with focus on Simulations of Planetary Systems

Authors: Saleh Dindar, Eric B. Ford, Mario Juric, Young In Yeo, Jianwei Gao, Aaron C. Boley, Benjamin Nelson, Jorg Peters

Abstract: We present Swarm-NG, a C++ library for the efficient direct integration of many n-body systems using highly-parallel Graphics Processing Unit (GPU), such as NVIDIA's Tesla T10 and M2070 GPUs. While previous studies have demonstrated the benefit of GPUs for n-body simulations with thousands to millions of bodies, Swarm-NG focuses on many few-body systems, e.g., thousands of systems with 3...15 bodi… ▽ More We present Swarm-NG, a C++ library for the efficient direct integration of many n-body systems using highly-parallel Graphics Processing Unit (GPU), such as NVIDIA's Tesla T10 and M2070 GPUs. While previous studies have demonstrated the benefit of GPUs for n-body simulations with thousands to millions of bodies, Swarm-NG focuses on many few-body systems, e.g., thousands of systems with 3...15 bodies each, as is typical for the study of planetary systems. Swarm-NG parallelizes the simulation, including both the numerical integration of the equations of motion and the evaluation of forces using NVIDIA's "Compute Unified Device Architecture" (CUDA) on the GPU. Swarm-NG includes optimized implementations of 4th order time-symmetrized Hermite integration and mixed variable symplectic integration, as well as several sample codes for other algorithms to illustrate how non-CUDA-savvy users may themselves introduce customized integrators into the Swarm-NG framework. To optimize performance, we analyze the effect of GPU-specific parameters on performance under double precision. Applications of Swarm-NG include studying the late stages of planet formation, testing the stability of planetary systems and evaluating the goodness-of-fit between many planetary system models and observations of extrasolar planet host stars (e.g., radial velocity, astrometry, transit timing). While Swarm-NG focuses on the parallel integration of many planetary systems,the underlying integrators could be applied to a wide variety of problems that require repeatedly integrating a set of ordinary differential equations many times using different initial conditions and/or parameter values. △ Less

Submitted 24 September, 2012; v1 submitted 6 August, 2012; originally announced August 2012.

Comments: Submitted to New Astronomy

Interactions Between Moderate- and Long-Period Giant Planets: Scattering Experiments for Systems in Isolation and with Stellar Flybys

Authors: Aaron C. Boley, Matthew J. Payne, Eric B. Ford

Abstract: The chance that a planetary system will interact with another member of its host star's nascent cluster would be greatly increased if gas giant planets form in situ on wide orbits. In this paper, we explore the outcomes of planet-planet scattering for a distribution of multiplanet systems that all have one of the planets on an initial orbit of 100 AU. The scattering experiments are run with and wi… ▽ More The chance that a planetary system will interact with another member of its host star's nascent cluster would be greatly increased if gas giant planets form in situ on wide orbits. In this paper, we explore the outcomes of planet-planet scattering for a distribution of multiplanet systems that all have one of the planets on an initial orbit of 100 AU. The scattering experiments are run with and without stellar flybys. We convolve the outcomes with distributions for protoplanetary disk and stellar cluster sizes to generalize the results where possible. We find that the frequencies of large mutual inclinations and high eccentricities are sensitive to the number of planets in a system, but not strongly to stellar flybys. However, flybys do play a role in changing the low and moderate portions of the mutual inclination distributions, and erase dynamically cold initial conditions on average. Wide-orbit planets can be mixed throughout the planetary system, and in some cases, can potentially become hot Jupiters, which we demonstrate using scattering experiments that include a tidal damping model. If planets form on wide orbits in situ, then there will be discernible differences in the proper motion distributions of a sample of wide-orbit planets compared with a pure scattering formation mechanism. Stellar flybys can enhance the frequency of ejections in planetary systems, but auto-ionization is likely to remain the dominant source of free-floating planets. △ Less

Submitted 23 April, 2012; originally announced April 2012.

Comments: Accepted for publication by ApJ

Chondrule Formation in Bow Shocks around Eccentric Planetary Embryos

Authors: Melissa A. Morris, Aaron C. Boley, Steven J. Desch, Themis Athanassiadou

Abstract: Recent isotopic studies of Martian meteorites by Dauphas & Pourmond (2011) have established that large (~ 3000 km radius) planetary embryos existed in the solar nebula at the same time that chondrules - millimeter-sized igneous inclusions found in meteorites - were forming. We model the formation of chondrules by passage through bow shocks around such a planetary embryo on an eccentric orbit. We n… ▽ More Recent isotopic studies of Martian meteorites by Dauphas & Pourmond (2011) have established that large (~ 3000 km radius) planetary embryos existed in the solar nebula at the same time that chondrules - millimeter-sized igneous inclusions found in meteorites - were forming. We model the formation of chondrules by passage through bow shocks around such a planetary embryo on an eccentric orbit. We numerically model the hydrodynamics of the flow, and find that such large bodies retain an atmosphere, with Kelvin-Helmholtz instabilities allowing mixing of this atmosphere with the gas and particles flowing past the embryo. We calculate the trajectories of chondrules flowing past the body, and find that they are not accreted by the protoplanet, but may instead flow through volatiles outgassed from the planet's magma ocean. In contrast, chondrules are accreted onto smaller planetesimals. We calculate the thermal histories of chondrules passing through the bow shock. We find that peak temperatures and cooling rates are consistent with the formation of the dominant, porphyritic texture of most chondrules, assuming a modest enhancement above the likely solar nebula average value of chondrule densities (by a factor of 10), attributable to settling of chondrule precursors to the midplane of the disk or turbulent concentration. We calculate the rate at which a planetary embryo's eccentricity is damped and conclude that a single planetary embryo scattered into an eccentric orbit can, over ~ 10e5 years, produce ~ 10e24 g of chondrules. In principle, a small number (1-10) of eccentric planetary embryos can melt the observed mass of chondrules in a manner consistent with all known constraints. △ Less

Submitted 3 April, 2012; originally announced April 2012.

Comments: Accepted for publication in The Astrophysical Journal

Constraining the Planetary System of Fomalhaut Using High-Resolution ALMA Observations

Authors: Aaron C. Boley, Matthew J. Payne, Stuartt Corder, William Dent, Eric B. Ford, Megan Shabram

Abstract: The dynamical evolution of planetary systems leaves observable signatures in debris disks. Optical images trace micron-sized grains, which are strongly affected by stellar radiation and need not coincide with their parent body population. Observations of mm-size grains accurately trace parent bodies, but previous images lack the resolution and sensitivity needed to characterize the ring's morpholo… ▽ More The dynamical evolution of planetary systems leaves observable signatures in debris disks. Optical images trace micron-sized grains, which are strongly affected by stellar radiation and need not coincide with their parent body population. Observations of mm-size grains accurately trace parent bodies, but previous images lack the resolution and sensitivity needed to characterize the ring's morphology. Here we present ALMA 350 GHz observations of the Fomalhaut debris ring. These observations demonstrate that the parent body population is 13-19 AU wide with a sharp inner and outer boundary. We discuss three possible origins for the ring, and suggest that debris confined by shepherd planets is the most consistent with the ring's morphology. △ Less

Submitted 30 March, 2012; originally announced April 2012.

Comments: First ALMA cycle 0 science results. Accepted by ApJL

Chemistry in a gravitationally unstable protoplanetary disc

Authors: J. D. Ilee, A. C. Boley, P. Caselli, R. H. Durisen, T. W. Hartquist, J. M. C. Rawlings

Abstract: Until now, axisymmetric, alpha-disc models have been adopted for calculations of the chemical composition of protoplanetary discs. While this approach is reasonable for many discs, it is not appropriate when self-gravity is important. In this case, spiral waves and shocks cause temperature and density variations that affect the chemistry. We have adopted a dynamical model of a solar-mass star surr… ▽ More Until now, axisymmetric, alpha-disc models have been adopted for calculations of the chemical composition of protoplanetary discs. While this approach is reasonable for many discs, it is not appropriate when self-gravity is important. In this case, spiral waves and shocks cause temperature and density variations that affect the chemistry. We have adopted a dynamical model of a solar-mass star surrounded by a massive (0.39 Msun), self-gravitating disc, similar to those that may be found around Class 0 and early Class I protostars, in a study of disc chemistry. We find that for each of a number of species, e.g. H2O, adsorption and desorption dominate the changes in the gas-phase fractional abundance; because the desorption rates are very sensitive to temperature, maps of the emissions from such species should reveal the locations of shocks of varying strengths. The gas-phase fractional abundances of some other species, e.g. CS, are also affected by gas-phase reactions, particularly in warm shocked regions. We conclude that the dynamics of massive discs have a strong impact on how they appear when imaged in the emission lines of various molecular species. △ Less

Submitted 15 July, 2011; originally announced July 2011.

Comments: 10 figures and 3 tables, accepted for publication in MNRAS

Migration of Gas Giant Planets in Gravitationally Unstable Disks

Authors: Scott Michael, Richard H. Durisen, Aaron C. Boley

Abstract: Characterization of migration in gravitationally unstable disks is necessary to understand the fate of protoplanets formed by disk instability. As part of a larger study, we are using a 3D radiative hydrodynamics code to investigate how an embedded gas giant planet interacts with a gas disk that undergoes gravitational instabilities (GIs). This Letter presents results from simulations with a Jupit… ▽ More Characterization of migration in gravitationally unstable disks is necessary to understand the fate of protoplanets formed by disk instability. As part of a larger study, we are using a 3D radiative hydrodynamics code to investigate how an embedded gas giant planet interacts with a gas disk that undergoes gravitational instabilities (GIs). This Letter presents results from simulations with a Jupiter-mass planet placed in orbit at 25 AU within a 0.14 $M_{\odot}$ disk. The disk spans 5 to 40 AU around a 1 $M_{\odot}$ star and is initially marginally unstable. In one simulation, the planet is inserted prior to the eruption of GIs; in another, it is inserted only after the disk has settled into a quasi-steady GI-active state, where heating by GIs roughly balances radiative cooling. When the planet is present from the beginning, its own wake stimulates growth of a particular global mode with which it strongly interacts, and the planet plunges inward six AU in about 10$^3$ years. In both cases with embedded planets, there are times when the planet's radial motion is slow and varies in direction. At other times, when the planet appears to be interacting with strong spiral modes, migration both inward and outward can be relatively rapid, covering several AUs over hundreds of years. Migration in both cases appears to stall near the inner Lindblad resonance of a dominant low-order mode. Planet orbit eccentricities fluctuate rapidly between about 0.02 to 0.1 throughout the GI-active phases of the simulations. △ Less

Submitted 8 June, 2011; originally announced June 2011.

Comments: Submitted to ApJL

The Heavy Element Composition of Disk Instability Planets Can Range From Sub- to Super-Nebular

Authors: Aaron C. Boley, Ravit Helled, Matthew J. Payne

Abstract: Transit surveys combined with Doppler data have revealed a class of gas giant planets that are massive and highly enriched in heavy elements (e.g., HD149026b, GJ436b, and HAT-P-20b). It is tempting to consider these planets as validation of core accretion plus gas capture because it is often assumed that disk instability planets should be of nebular composition. We show in this paper, to the contr… ▽ More Transit surveys combined with Doppler data have revealed a class of gas giant planets that are massive and highly enriched in heavy elements (e.g., HD149026b, GJ436b, and HAT-P-20b). It is tempting to consider these planets as validation of core accretion plus gas capture because it is often assumed that disk instability planets should be of nebular composition. We show in this paper, to the contrary, that gas giants that form by disk instability can have a variety of heavy element compositions, ranging from sub- to super-nebular values. High levels of enrichment can be achieved through one or multiple mechanisms, including enrichment at birth, planetesimal capture, and differentiation plus tidal stripping. As a result, the metallicity of an individual gas giant cannot be used to discriminate between gas giant formation modes. △ Less

Submitted 5 April, 2011; originally announced April 2011.

Comments: Accepted by ApJ

On The Possibility of Enrichment and Differentiation in Gas Giants During Birth by Disk Instability

Authors: Aaron C. Boley, Richard H. Durisen

Abstract: We investigate the coupling between rock-size solids and gas during the formation of gas giant planets by disk fragmentation in the outer regions of massive disks. In this study, we use three-dimensional radiative hydrodynamics simulations and model solids as a spatial distribution of particles. We assume that half of the total solid fraction is in small grains and half in large solids. The former… ▽ More We investigate the coupling between rock-size solids and gas during the formation of gas giant planets by disk fragmentation in the outer regions of massive disks. In this study, we use three-dimensional radiative hydrodynamics simulations and model solids as a spatial distribution of particles. We assume that half of the total solid fraction is in small grains and half in large solids. The former are perfectly entrained with the gas and set the opacity in the disk, while the latter are allowed to respond to gas drag forces, with the back reaction on the gas taken into account. To explore the maximum effects of gas-solid interactions, we first consider 10cm-size particles. We then compare these results to a simulation with 1 km-size particles, which explores the low-drag regime. We show that (1) disk instability planets have the potential to form large cores due to aerodynamic capturing of rock-size solids in spiral arms before fragmentation; (2) that temporary clumps can concentrate tens of $M_{\oplus}$ of solids in very localized regions before clump disruption; (3) that the formation of permanent clumps, even in the outer disk, is dependent on the grain-size distribution, i.e., the opacity; (4) that nonaxisymmetric structure in the disk can create disk regions that have a solids-to-gas ratio greater than unity; (5) that the solid distribution may affect the fragmentation process; (6) that proto-gas giants and proto-brown dwarfs can start as differentiated objects prior to the H$_2$ collapse phase; (7) that spiral arms in a gravitationally unstable disk are able to stop the inward drift of rock-size solids, even redistributing them to larger radii; and, (8) that large solids can form spiral arms that are offset from the gaseous spiral arms. We conclude that planet embryo formation can be strongly affected by the growth of solids during the earliest stages of disk accretion. △ Less

Submitted 17 September, 2010; v1 submitted 14 May, 2010; originally announced May 2010.

Comments: Accepted by ApJ. 55 pages including 24 figures. In response to comments from the referee, we have included a new simulation with km-size objects and have revised some discussions and interpretations. Major conclusions remain unchanged, and new conclusions have been added in response to the new run

The Properties of Prestellar Discs in Isolated and Multiple Prestellar Systems

Authors: T. Hayfield, L. Mayer, J. Wadsley, A. C. Boley

Abstract: We present high-resolution 3D smoothed particle hydrodynamics simulations of the formation and evolution of protostellar discs in a turbulent molecular cloud. Using a piecewise polytropic equation of state, we perform two sets of simulations. In both cases we find that isolated systems undergo a fundamentally different evolution than members of binary or multiple systems. When formed, isolated sys… ▽ More We present high-resolution 3D smoothed particle hydrodynamics simulations of the formation and evolution of protostellar discs in a turbulent molecular cloud. Using a piecewise polytropic equation of state, we perform two sets of simulations. In both cases we find that isolated systems undergo a fundamentally different evolution than members of binary or multiple systems. When formed, isolated systems must accrete mass and increase their specific angular momentum, leading to the formation of massive, extended discs, which undergo strong gravitational instabilities and are susceptible to disc fragmentation. Fragments with initial masses of 5.5 M_jup, 7.4 M_jup and 12 M_jup are produced in our simulations. In binaries and small clusters, we observe that due to competition for material from the parent core, members do not accrete significant amounts of high specific angular momentum gas relative to isolated systems. We find that discs in multiple systems are strongly self-gravitating but that they are stable against fragmentation due to disc truncation and mass profile steeping by tides, accretion of high specific angular momentum gas by other members, and angular momentum being redirected into members' orbits. In general, we expect disc fragmentation to be less likely in clusters and to be more a feature of isolated systems. △ Less

Submitted 7 July, 2011; v1 submitted 12 March, 2010; originally announced March 2010.

Comments: 15 pages, 21 figures. Accepted for publication in MNRAS

Clumps in the Outer Disk by Disk Instability: Why They are Initially Gas Giants and the Legacy of Disruption

Authors: Aaron C. Boley, Tristen Hayfield, Lucio Mayer, Richard H. Durisen

Abstract: We explore the initial conditions for fragments in the extended regions ($r\gtrsim 50$ AU) of gravitationally unstable disks. We combine analytic estimates for the fragmentation of spiral arms with 3D SPH simulations to show that initial fragment masses are in the gas giant regime. These initial fragments will have substantial angular momentum, and should form disks with radii of a few AU. We sh… ▽ More We explore the initial conditions for fragments in the extended regions ($r\gtrsim 50$ AU) of gravitationally unstable disks. We combine analytic estimates for the fragmentation of spiral arms with 3D SPH simulations to show that initial fragment masses are in the gas giant regime. These initial fragments will have substantial angular momentum, and should form disks with radii of a few AU. We show that clumps will survive for multiple orbits before they undergo a second, rapid collapse due to H$_2$ dissociation and that it is possible to destroy bound clumps by transporting them into the inner disk. The consequences of disrupted clumps for planet formation, dust processing, and disk evolution are discussed. We argue that it is possible to produce Earth-mass cores in the outer disk during the earliest phases of disk evolution. △ Less

Submitted 11 January, 2010; v1 submitted 24 September, 2009; originally announced September 2009.

Comments: Accepted for publication in Icarus. The arguments have been greatly expanded (from v1) to address comments by the referees

Resolving Globular Cluster Formation within a Cosmological Context

Authors: Aaron C. Boley, George Lake, Justin Read, Romain Teyssier

Abstract: We place constraints on the formation redshifts for blue globular clusters (BGCs), independent of the details of hydrodynamics and population III star formation. The observed radial distribution of BGCs in the Milky Way Galaxy suggests that they formed in biased dark matter halos at high redshift. As a result, simulations of a ~1 Mpc box up to z~10 must resolve BGC formation in LCDM. We find tha… ▽ More We place constraints on the formation redshifts for blue globular clusters (BGCs), independent of the details of hydrodynamics and population III star formation. The observed radial distribution of BGCs in the Milky Way Galaxy suggests that they formed in biased dark matter halos at high redshift. As a result, simulations of a ~1 Mpc box up to z~10 must resolve BGC formation in LCDM. We find that most halo stars could be produced from destroyed BGCs and other low-mass clusters that formed at high redshift. We present a proof-of-concept simulation that captures the formation of globular-like star clusters. △ Less

Submitted 21 October, 2009; v1 submitted 9 August, 2009; originally announced August 2009.

Comments: Accepted for publication in ApJL

The Two Modes of Gas Giant Planet Formation

Authors: Aaron C. Boley

Abstract: I argue for two modes of gas giant planet formation and discuss the conditions under which each mode operates. Gas giant planets at disk radii $r>100$ AU are likely to form in situ by disk instability, while core accretion plus gas capture remains the dominant formation mechanism for $r<100$ AU. During the mass accretion phase, mass loading can push disks toward fragmentation conditions at large… ▽ More I argue for two modes of gas giant planet formation and discuss the conditions under which each mode operates. Gas giant planets at disk radii $r>100$ AU are likely to form in situ by disk instability, while core accretion plus gas capture remains the dominant formation mechanism for $r<100$ AU. During the mass accretion phase, mass loading can push disks toward fragmentation conditions at large $r$. Massive, extended disks can fragment into clumps of a few to tens of Jupiter masses. This is confirmed by radiation hydrodynamics simulations. The two modes of gas giant formation should lead to a bimodal distribution of gas giant semi-major axes. Because core accretion is expected to be less efficient in low-metallicity systems, the ratio of gas giants at large $r$ to planets at small $r$ should increase with decreasing metallicity. △ Less

Submitted 23 February, 2009; originally announced February 2009.

Comments: Submitted to ApJL after addressing referee's comments

Journal ref: Astrophys.J.695:L53-L57,2009

Gravitational Instabilities, Chondrule Formation, and the FU Orionis Phenomenon

Authors: Aaron C. Boley, Richard H. Durisen

Abstract: Using analytic arguments and numerical simulations, we examine whether chondrule formation and the FU Orionis phenomenon can be caused by the burst-like onset of gravitational instabilities (GIs) in dead zones. At least two scenarios for bursting dead zones can work, in principle. If the disk is on the verge of fragmention, GI activation near $r\sim4$ to 5 AU can produce chondrule-forming shocks… ▽ More Using analytic arguments and numerical simulations, we examine whether chondrule formation and the FU Orionis phenomenon can be caused by the burst-like onset of gravitational instabilities (GIs) in dead zones. At least two scenarios for bursting dead zones can work, in principle. If the disk is on the verge of fragmention, GI activation near $r\sim4$ to 5 AU can produce chondrule-forming shocks, at least under extreme conditions. Mass fluxes are also high enough during the onset of GIs to suggest that the outburst is related to an FU Orionis phenomenon. This situation is demonstrated by numerical simulations. In contrast, as supported by analytic arguments, if the burst takes place close to $r\sim1$ AU, then even low pitch angle spiral waves can create chondrule-producing shocks and outbursts. We also study the stability of the massive disks in our simulations against fragmentation and find that although disk evolution is sensitive to changes in opacity, the disks we study do not fragment, even at high resolution and even for extreme assumptions. △ Less

Submitted 10 June, 2008; originally announced June 2008.

Comments: To appear in ApJ