Showing 1–50 of 342 results for author: Graham, R

Thermally activated dynamics of annealed glasses near the yielding transition under cyclic shear

Authors: Ian R Graham, Paulo E Arratia, Robert A Riggleman

Abstract: While experiments and simulations have provided a rich picture of the dynamic heterogeneity in glasses at constant temperature or under steady shear, the dynamics of glasses under oscillatory shear remain comparatively less explored. Recent work has shown that oscillatory shear protocols can embed a ``memory'' into a glass's structure, whereby the material will exhibit dynamics that are encoded by… ▽ More While experiments and simulations have provided a rich picture of the dynamic heterogeneity in glasses at constant temperature or under steady shear, the dynamics of glasses under oscillatory shear remain comparatively less explored. Recent work has shown that oscillatory shear protocols can embed a ``memory'' into a glass's structure, whereby the material will exhibit dynamics that are encoded by the oscillatory shear protocol applied. However, most of the computational work studying the memory effect has been performed in the zero temperature limit, and the effects of thermalization are poorly characterized. In this work, we use nonequilibrium molecular dynamics simulations to study the dynamics of a model two-dimensional glass former at low, non-zero temperatures under oscillatory shear. While we show that the systems' dynamics are independent of sample preparation for either small or larger strain amplitudes, the dynamics become distinct near the yield point when the deformation is applied at finite temperature. We then characterize the dynamic heterogeneity using two metrics, one derived from the vibrational modes and one that exploits machine learning to identify regions prone to rearrangement. This analysis provides evidence that the dynamics below and above yield emerge from distinct structural origins that may be important for developing improved constitutive models that can predict memory in disordered solids. △ Less

Submitted 1 October, 2024; originally announced October 2024.

Substantial extension of the lifetime of the terrestrial biosphere

Authors: R. J. Graham, Itay Halevy, Dorian Abbot

Abstract: Approximately one billion years (Gyr) in the future, as the Sun brightens, Earth's carbonate-silicate cycle is expected to drive CO$_2$ below the minimum level required by vascular land plants, eliminating most macroscopic land life. Here, we couple global-mean models of temperature- and CO$_2$-dependent plant productivity for C$_3$ and C$_4$ plants, silicate weathering, and climate to re-examine… ▽ More Approximately one billion years (Gyr) in the future, as the Sun brightens, Earth's carbonate-silicate cycle is expected to drive CO$_2$ below the minimum level required by vascular land plants, eliminating most macroscopic land life. Here, we couple global-mean models of temperature- and CO$_2$-dependent plant productivity for C$_3$ and C$_4$ plants, silicate weathering, and climate to re-examine the time remaining for terrestrial plants. If weathering is weakly temperature-dependent (as recent data suggest) and/or strongly CO$_2$-dependent, we find that the interplay between climate, productivity, and weathering causes the future luminosity-driven CO$_2$ decrease to slow and temporarily reverse, averting plant CO$_2$ starvation. This dramatically lengthens plant survival from 1 Gyr up to $\sim$1.6-1.86 Gyr, until extreme temperatures halt photosynthesis, suggesting a revised kill mechanism for land plants and potential doubling of the future lifespan of Earth's land macrobiota. An increased future lifespan for the complex biosphere may imply that Earth life had to achieve a smaller number of ``hard steps'' (unlikely evolutionary transitions) to produce intelligent life than previously estimated. These results also suggest that complex photosynthetic land life on Earth and exoplanets may be able to persist until the onset of the moist greenhouse transition. △ Less

Submitted 16 September, 2024; originally announced September 2024.

Comments: 19 pages, 6 figures. Accepted for publication in Planetary Science Journal

Resilience of Snowball Earth to Stochastic Events

Authors: Guillaume Chaverot, Andrea Zorzi, Xuesong Ding, Jonathan Itcovitz, Bowen Fan, Siddharth Bhatnagar, Aoshuang Ji, Robert J. Graham, Tushar Mittal

Abstract: Earth went through at least two periods of global glaciation (i.e., ``Snowball Earth'' states) during the Neoproterozoic, the shortest of which (the Marinoan) may not have lasted sufficiently long for its termination to be explained by the gradual volcanic build-up of greenhouse gases in the atmosphere. Large asteroid impacts and supervolcanic eruptions have been suggested as stochastic geological… ▽ More Earth went through at least two periods of global glaciation (i.e., ``Snowball Earth'' states) during the Neoproterozoic, the shortest of which (the Marinoan) may not have lasted sufficiently long for its termination to be explained by the gradual volcanic build-up of greenhouse gases in the atmosphere. Large asteroid impacts and supervolcanic eruptions have been suggested as stochastic geological events that could cause a sudden end to global glaciation via a runaway melting process. Here, we employ an energy balance climate model to simulate the evolution of Snowball Earth's surface temperature after such events. We find that even a large impactor (diameters of $d \sim 100\,\mathrm{km}$) and the supervolcanic Toba eruption ($74\,\mathrm{kyr}$ ago), are insufficient to terminate a Snowball state unless background CO$_2$ has already been driven to high levels by long-term outgassing. We suggest, according to our modeling framework, that Earth's Snowball states would have been resilient to termination by stochastic events. △ Less

Submitted 15 July, 2024; originally announced July 2024.

Comments: Accepted for publication in Geophysical Research Letters (GRL)

Carbon cycle instability for high-$\mathrm{CO_2}$ exoplanets: Implications for habitability

Authors: R. J. Graham, R. T. Pierrehumbert

Abstract: Implicit in the definition of the classical circumstellar habitable zone (HZ) is the hypothesis that the carbonate-silicate cycle can maintain clement climates on exoplanets with land and surface water across a range of instellations by adjusting atmospheric $\mathrm{CO_2}$ partial pressure ($p\mathrm{CO_2}$). This hypothesis is made by analogy to the Earth system, but it is an open question wheth… ▽ More Implicit in the definition of the classical circumstellar habitable zone (HZ) is the hypothesis that the carbonate-silicate cycle can maintain clement climates on exoplanets with land and surface water across a range of instellations by adjusting atmospheric $\mathrm{CO_2}$ partial pressure ($p\mathrm{CO_2}$). This hypothesis is made by analogy to the Earth system, but it is an open question whether silicate weathering can stabilize climate on planets in the outer reaches of the HZ, where instellations are lower than those received by even the Archean Earth and $\mathrm{CO_2}$ is thought likely to dominate atmospheres. Since weathering products are carried from land to ocean by the action of water, silicate weathering is intimately coupled to the hydrologic cycle, which intensifies with hotter temperatures under Earth-like conditions. Here, we use global climate model (GCM) simulations to demonstrate that the hydrologic cycle responds counterintuitively to changes in climate on planets with $\mathrm{CO_2}$-$\mathrm{H_2O}$ atmospheres at low instellations and high $p\mathrm{CO_2}$, with global evaporation and precipitation decreasing as $p\mathrm{CO_2}$ and temperatures increase at a given instellation. Within the MAC weathering formulation, weathering then decreases with increasing $p\mathrm{CO_2}$ for a range of instellations and $p\mathrm{CO_2}$ typical of the outer reaches of the HZ, resulting in an unstable carbon cycle that may lead to either runaway $\mathrm{CO_2}$ accumulation or depletion of $\mathrm{CO_2}$ to colder (possibly Snowball) conditions. While the behavior of the system has not been completely mapped out, the results suggest that silicate weathering could fail to maintain habitable conditions in the outer reaches of the nominal HZ. △ Less

Submitted 8 May, 2024; originally announced May 2024.

Comments: Accepted for publication in ApJ

A Gauss-Bonnet formula for the renormalized area of minimal submanifolds of Poincaré-Einstein manifolds

Authors: Jeffrey S. Case, C Robin Graham, Tzu-Mo Kuo, Aaron J. Tyrrell, Andrew Waldron

Abstract: Assuming the extrinsic $Q$-curvature admits a decomposition into the Pfaffian, a scalar conformal submanifold invariant, and a tangential divergence, we prove that the renormalized area of an even-dimensional minimal submanifold of a Poincaré-Einstein manifold can be expressed as a linear combination of its Euler characteristic and the integral of a scalar conformal submanifold invariant. We deriv… ▽ More Assuming the extrinsic $Q$-curvature admits a decomposition into the Pfaffian, a scalar conformal submanifold invariant, and a tangential divergence, we prove that the renormalized area of an even-dimensional minimal submanifold of a Poincaré-Einstein manifold can be expressed as a linear combination of its Euler characteristic and the integral of a scalar conformal submanifold invariant. We derive such a decomposition of the extrinsic $Q$-curvature in dimensions two and four, thereby recovering and generalizing results of Alexakis-Mazzeo and Tyrrell, respectively. We also conjecture such a decomposition for general natural submanifold scalars whose integral over compact submanifolds is conformally invariant, and verify our conjecture in dimensions two and four. Our results also apply to the area of a compact even-dimensional minimal submanifold of an Einstein manifold. △ Less

Submitted 25 March, 2024; originally announced March 2024.

Comments: 47 pages

Design and Validation of a Very Low-Power Phasor Measurement Unit

Authors: Zachary J. Lythgoe, Thomas F. Long, Michael J. Buchholz, Anthony R. Livernois, Kebba Kanuteh, David R. Allee, Anamitra Pal, Ian R. Graham, Zachary D. Drummond

Abstract: Phasor measurement units (PMUs) provide a high-resolution view of the power system at the locations where they are placed. As such, it is desirable to place them in bulk in low voltage distribution circuits. However, the power consumption of a PMU/micro-PMU is in the order of Watts (W) that results in them requiring an external power supply, which in turn increases the overall cost. This work deta… ▽ More Phasor measurement units (PMUs) provide a high-resolution view of the power system at the locations where they are placed. As such, it is desirable to place them in bulk in low voltage distribution circuits. However, the power consumption of a PMU/micro-PMU is in the order of Watts (W) that results in them requiring an external power supply, which in turn increases the overall cost. This work details the hardware design of a PMU capable of measuring and reporting voltage and current phasors for a single-phase system at an average power consumption of only 30.8 mW -- one to two orders of magnitude lower than existing academic and commercial PMUs. This enables the proposed PMU to run for two weeks using an 11-Wh battery or indefinitely if paired with an inexpensive solar panel. A test-bench developed in accordance with the 2018 IEC/IEEE 60255-118-1 PMU Standard confirms the accuracy of this PMU. Given its low power consumption, the proposed design is expected to accelerate adoption of PMUs in modern distribution grids. △ Less

Submitted 22 February, 2024; originally announced February 2024.

Comments: 6 pages, 4 figures, Presented at Texas Power and Energy Conference (TPEC) 2024

Sex-Specific Variances in Anatomy and Blood Flow of the Left Main Coronary Bifurcation: Implications for Coronary Artery Disease Risk

Authors: Ramtin Gharleghi, Mingzi Zhang, Dona Adikari, Lucy McGrath-Cadell, Robert M. Graham, Jolanda Wentzel, Mark Webster, Chris Ellis, Sze-Yuan Ooi, Susann Beier

Abstract: Studies have shown marked sex disparities in Coronary Artery Diseases (CAD) epidemiology, yet the underlying mechanisms remain unclear. We explored sex disparities in the coronary anatomy and the resulting haemodynamics in patients with suspected, but no significant CAD. Left Main (LM) bifurcations were reconstructed from CTCA images of 127 cases (42 males and 85 females, aged 38 to 81). Detailed… ▽ More Studies have shown marked sex disparities in Coronary Artery Diseases (CAD) epidemiology, yet the underlying mechanisms remain unclear. We explored sex disparities in the coronary anatomy and the resulting haemodynamics in patients with suspected, but no significant CAD. Left Main (LM) bifurcations were reconstructed from CTCA images of 127 cases (42 males and 85 females, aged 38 to 81). Detailed shape parameters were measured for comparison, including bifurcation angles, curvature, and diameters, before solving the haemodynamic metrics using CFD. The severity and location of the normalised vascular area exposed to physiologically adverse haemodynamics were statistically compared between sexes for all branches. We found significant differences between sexes in potentially adverse haemodynamics. Females were more likely than males to exhibit adversely low Time Averaged Endothelial Shear Stress along the inner wall of a bifurcation (16.8% vs. 10.7%). Males had a higher percentage of areas exposed to both adversely high Relative Residence Time (6.1% vs 4.2%, p=0.001) and high Oscillatory Shear Index (4.6% vs 2.3%, p<0.001). However, the OSI values were generally small and should be interpreted cautiously. Males had larger arteries (M vs F, LM: 4.0mm vs 3.3mm, LAD: 3.6mm 3.0mm, LCX:3.5mm vs 2.9mm), and females exhibited higher curvatures in all three branches (M vs F, LM: 0.40 vs 0.46, LAD: 0.45 vs 0.51, LCx: 0.47 vs 0.55, p<0.001) and larger inflow angle of the LM trunk (M: 12.9° vs F: 18.5°, p=0.025). Haemodynamic differences were found between male and female patients, which may contribute, at least in part, to differences in CAD risk. This work may facilitate a better understanding of sex differences in the clinical presentation of CAD, contributing to improved sex-specific screening, especially relevant for women with CAD who currently have worse predictive outcomes. △ Less

Submitted 30 November, 2023; originally announced November 2023.

Comments: 14 pages, 5 figures

Utilitarian Algorithm Configuration

Authors: Devon R. Graham, Kevin Leyton-Brown, Tim Roughgarden

Abstract: We present the first nontrivial procedure for configuring heuristic algorithms to maximize the utility provided to their end users while also offering theoretical guarantees about performance. Existing procedures seek configurations that minimize expected runtime. However, very recent theoretical work argues that expected runtime minimization fails to capture algorithm designers' preferences. Here… ▽ More We present the first nontrivial procedure for configuring heuristic algorithms to maximize the utility provided to their end users while also offering theoretical guarantees about performance. Existing procedures seek configurations that minimize expected runtime. However, very recent theoretical work argues that expected runtime minimization fails to capture algorithm designers' preferences. Here we show that the utilitarian objective also confers significant algorithmic benefits. Intuitively, this is because mean runtime is dominated by extremely long runs even when they are incredibly rare; indeed, even when an algorithm never gives rise to such long runs, configuration procedures that provably minimize mean runtime must perform a huge number of experiments to demonstrate this fact. In contrast, utility is bounded and monotonically decreasing in runtime, allowing for meaningful empirical bounds on a configuration's performance. This paper builds on this idea to describe effective and theoretically sound configuration procedures. We prove upper bounds on the runtime of these procedures that are similar to theoretical lower bounds, while also demonstrating their performance empirically. △ Less

Submitted 31 October, 2023; originally announced October 2023.

Extrinsic GJMS operators for submanifolds

Authors: Jeffrey S. Case, C Robin Graham, Tzu-Mo Kuo

Abstract: We derive extrinsic GJMS operators and $Q$-curvatures associated to a submanifold of a conformal manifold. The operators are conformally covariant scalar differential operators on the submanifold with leading part a power of the Laplacian in the induced metric. Upon realizing the conformal manifold as the conformal infinity of an asymptotically Poincaré--Einstein space and the submanifold as the b… ▽ More We derive extrinsic GJMS operators and $Q$-curvatures associated to a submanifold of a conformal manifold. The operators are conformally covariant scalar differential operators on the submanifold with leading part a power of the Laplacian in the induced metric. Upon realizing the conformal manifold as the conformal infinity of an asymptotically Poincaré--Einstein space and the submanifold as the boundary of an asymptotically minimal submanifold thereof, these operators arise as obstructions to smooth extension as eigenfunctions of the Laplacian of the induced metric on the minimal submanifold. We derive explicit formulas for the operators of orders 2 and 4. We prove factorization formulas when the original submanifold is a minimal submanifold of an Einstein manifold. We also show how to reformulate the construction in terms of the ambient metric for the conformal manifold, and use this to prove that the operators defined by the factorization formulas are conformally invariant for all orders in all dimensions. △ Less

Submitted 19 March, 2024; v1 submitted 20 June, 2023; originally announced June 2023.

Comments: 33 pages, v2 has minor editing, Remark 4.3 added, references updated

MSC Class: 53A31 (Primary) 53A55; 53A10; 58J70; 53C18 (Secondary)

Testing the Interaction Between a Substellar Companion and a Debris Disk in the HR 2562 System

Authors: Stella Yimiao Zhang, Gaspard Duchêne, Robert J. De Rosa, Megan Ansdell, Quinn Konopacky, Thomas Esposito, Eugene Chiang, Malena Rice, Brenda Matthews, Paul Kalas, Bruce Macintosh, Franck Marchis, Stan Metchev, Jenny Patience, Julien Rameau, Kimberly Ward-Duong, Schuyler Wolff, Michael P. Fitzgerald, Vanessa P. Bailey, Travis S. Barman, Joanna Bulger, Christine H. Chen, Jeffrey K. Chilcotte, Tara Cotten, René Doyon , et al. (29 additional authors not shown)

Abstract: The HR 2562 system is a rare case where a brown dwarf companion resides in a cleared inner hole of a debris disk, offering invaluable opportunities to study the dynamical interaction between a substellar companion and a dusty disk. We present the first ALMA observation of the system as well as the continued GPI monitoring of the companion's orbit with 6 new epochs from 2016 to 2018. We update the… ▽ More The HR 2562 system is a rare case where a brown dwarf companion resides in a cleared inner hole of a debris disk, offering invaluable opportunities to study the dynamical interaction between a substellar companion and a dusty disk. We present the first ALMA observation of the system as well as the continued GPI monitoring of the companion's orbit with 6 new epochs from 2016 to 2018. We update the orbital fit and, in combination with absolute astrometry from GAIA, place a 3$σ$ upper limit of 18.5 $M_J$ on the companion's mass. To interpret the ALMA observations, we used radiative transfer modeling to determine the disk properties. We find that the disk is well resolved and nearly edge on. While the misalignment angle between the disk and the orbit is weakly constrained due to the short orbital arc available, the data strongly support a (near) coplanar geometry for the system. Furthermore, we find that the models that describe the ALMA data best have an inner radius that is close to the companion's semi-major axis. Including a posteriori knowledge of the system's SED further narrows the constraints on the disk's inner radius and place it at a location that is in reasonable agreement with, possibly interior to, predictions from existing dynamical models of disk truncation by an interior substellar companion. HR\,2562 has the potential over the next few years to become a new testbed for dynamical interaction between a debris disk and a substellar companion. △ Less

Submitted 9 February, 2023; originally announced February 2023.

Comments: 22 pages, 17 figures

Exploring the Relationship Between Softness and Excess Entropy in Glass-forming Systems

Authors: Ian R. Graham, Paulo E. Arratia, Robert A. Riggleman

Abstract: We explore the relationship between a machine-learned structural quantity (softness) and excess entropy in simulations of supercooled liquids. Excess entropy is known to scale well the dynamical properties of liquids, but this quasi-universal scaling is known to breakdown in the supercooled and glassy regimes. Using numerical simulations, we test whether a local form of the excess entropy can lead… ▽ More We explore the relationship between a machine-learned structural quantity (softness) and excess entropy in simulations of supercooled liquids. Excess entropy is known to scale well the dynamical properties of liquids, but this quasi-universal scaling is known to breakdown in the supercooled and glassy regimes. Using numerical simulations, we test whether a local form of the excess entropy can lead to predictions that derive from softness, which has been shown to correlate well with the tendency for individual particles to rearrange. To that end, we explore leveraging softness to compute excess entropy in the traditional fashion over softness groupings. Our results show that by computing the excess entropy over softness-binned groupings, we can build a strong quantitative relationship between the rearrangement barriers across the explored systems. △ Less

Submitted 1 February, 2023; originally announced February 2023.

A Decade of Linear and Circular Polarimetry with the POLISH2 Polarimeter

Authors: Sloane J. Wiktorowicz, Agnieszka Słowikowska, Larissa A. Nofi, Nicole Rider, Angie Wolfgang, Ninos Hermis, Daniel Jontof-Hutter, Amanda J. Bayless, Gary M. Cole, Kirk B. Crawford, Valeri V. Tsarev, Michael C. Owens, Ernest G. Jaramillo, Geoffrey A. Maul, James R. Graham, Maxwell A. Millar-Blanchaer, Kimberly Bott, Jon C. Mauerhan

Abstract: The POLISH2 optical polarimeter has been in operation at the Lick Observatory 3-m Shane telescope since 2011, and it was commissioned at the Gemini North 8-m in 2016. This instrument primarily targets exoplanets, asteroids, and the Crab pulsar, but it has also been used for a wide variety of planetary, galactic, and supernova science. POLISH2's photoelastic modulators, employed instead of rotating… ▽ More The POLISH2 optical polarimeter has been in operation at the Lick Observatory 3-m Shane telescope since 2011, and it was commissioned at the Gemini North 8-m in 2016. This instrument primarily targets exoplanets, asteroids, and the Crab pulsar, but it has also been used for a wide variety of planetary, galactic, and supernova science. POLISH2's photoelastic modulators, employed instead of rotating waveplates or ferro-electric liquid crystal modulators, offer the unprecedented ability to achieve sensitivity and accuracy of order 1 ppm (0.0001%), which are difficult to obtain with conventional polarimeters. Additionally, POLISH2 simultaneously measures intensity (Stokes I), linear polarization (Stokes Q and U), and circular polarization (Stokes V), which fully describe the polarization state of incident light. We document our laboratory and on-sky calibration methodology, our archival on-sky database, and we demonstrate conclusive detection of circular polarization of certain objects. △ Less

Submitted 16 November, 2022; originally announced November 2022.

Comments: 38 pages, 16 figures, accepted by ApJS

CO$_2$ ocean bistability on terrestrial exoplanets

Authors: R. J. Graham, Tim Lichtenberg, Ray Pierrehumbert

Abstract: Cycling of carbon dioxide between the atmosphere and interior of rocky planets can stabilize global climate and enable planetary surface temperatures above freezing over geologic time. However, variations in global carbon budget and unstable feedback cycles between planetary sub-systems may destabilize the climate of rocky exoplanets toward regimes unknown in the Solar System. Here, we perform cle… ▽ More Cycling of carbon dioxide between the atmosphere and interior of rocky planets can stabilize global climate and enable planetary surface temperatures above freezing over geologic time. However, variations in global carbon budget and unstable feedback cycles between planetary sub-systems may destabilize the climate of rocky exoplanets toward regimes unknown in the Solar System. Here, we perform clear-sky atmospheric radiative transfer and surface weathering simulations to probe the stability of climate equilibria for rocky, ocean-bearing exoplanets at instellations relevant for planetary systems in the outer regions of the circumstellar habitable zone. Our simulations suggest that planets orbiting G- and F-type stars (but not M-type stars) may display bistability between an Earth-like climate state with efficient carbon sequestration and an alternative stable climate equilibrium where CO$_2$ condenses at the surface and forms a blanket of either clathrate hydrate or liquid CO$_2$. At increasing instellation and with ineffective weathering, the latter state oscillates between cool, surface CO$_2$-condensing and hot, non-condensing climates. CO$_2$ bistable climates may emerge early in planetary history and remain stable for billions of years. The carbon dioxide-condensing climates follow an opposite trend in $p$CO$_2$ versus instellation compared to the weathering-stabilized planet population, suggesting the possibility of observational discrimination between these distinct climate categories. △ Less

Submitted 30 September, 2022; originally announced October 2022.

Comments: Accepted for publication in JGR-Planets. 19 pages (including references and Tables), 6 figures, corresponding author email: robert.graham@physics.ox.ac.uk

Journal ref: Journal of Geophysical Research: Planets, 127, e2022JE007456

Formalizing Preferences Over Runtime Distributions

Authors: Devon R. Graham, Kevin Leyton-Brown, Tim Roughgarden

Abstract: When trying to solve a computational problem, we are often faced with a choice between algorithms that are guaranteed to return the right answer but differ in their runtime distributions (e.g., SAT solvers, sorting algorithms). This paper aims to lay theoretical foundations for such choices by formalizing preferences over runtime distributions. It might seem that we should simply prefer the algori… ▽ More When trying to solve a computational problem, we are often faced with a choice between algorithms that are guaranteed to return the right answer but differ in their runtime distributions (e.g., SAT solvers, sorting algorithms). This paper aims to lay theoretical foundations for such choices by formalizing preferences over runtime distributions. It might seem that we should simply prefer the algorithm that minimizes expected runtime. However, such preferences would be driven by exactly how slow our algorithm is on bad inputs, whereas in practice we are typically willing to cut off occasional, sufficiently long runs before they finish. We propose a principled alternative, taking a utility-theoretic approach to characterize the scoring functions that describe preferences over algorithms. These functions depend on the way our value for solving our problem decreases with time and on the distribution from which captimes are drawn. We describe examples of realistic utility functions and show how to leverage a maximum-entropy approach for modeling underspecified captime distributions. Finally, we show how to efficiently estimate an algorithm's expected utility from runtime samples. △ Less

Submitted 2 June, 2023; v1 submitted 25 May, 2022; originally announced May 2022.

Fermi acceleration in rotating drums

Authors: Krzysztof Burdzy, Mauricio Duarte, Carl-Erik Gauthier, C. Robin Graham, Jaime San Martin

Abstract: Consider hard balls in a bounded rotating drum. If there is no gravitation then there is no Fermi acceleration, i.e., the energy of the balls remains bounded forever. If there is gravitation, Fermi acceleration may arise. A number of explicit formulas for the system without gravitation are given. Some of these are based on an explicit realization, which we derive, of the well-known microcanonical… ▽ More Consider hard balls in a bounded rotating drum. If there is no gravitation then there is no Fermi acceleration, i.e., the energy of the balls remains bounded forever. If there is gravitation, Fermi acceleration may arise. A number of explicit formulas for the system without gravitation are given. Some of these are based on an explicit realization, which we derive, of the well-known microcanonical ensemble measure. △ Less

Submitted 17 December, 2021; originally announced December 2021.

MSC Class: 70L99

Local X-ray Transform on Asymptotically Hyperbolic Manifolds via Projective Compactification

Authors: Nikolas Eptaminitakis, C. Robin Graham

Abstract: We prove local injectivity near a boundary point for the geodesic X-ray transform for an asymptotically hyperbolic metric even mod $O(ρ^5)$ in dimensions three and higher. We prove local injectivity near a boundary point for the geodesic X-ray transform for an asymptotically hyperbolic metric even mod $O(ρ^5)$ in dimensions three and higher. △ Less

Submitted 8 January, 2022; v1 submitted 26 November, 2021; originally announced November 2021.

Comments: 30 pages. v2: Final version. Minor editorial changes. Dedication and memorial added

MSC Class: 53C65; 35R30; 53C22

Journal ref: New Zealand Journal of Mathematics 52 (2021), 733-763

Predicting the unobserved: a statistical mechanics framework for non-equilibrium material response with quantified uncertainty

Authors: Shenglin Huang, Ian R. Graham, Robert A. Riggleman, Paulo Arratia, Steve Fitzgerald, Celia Reina

Abstract: Can far-from-equilibrium material response under arbitrary loading be inferred from equilibrium data and vice versa? Can the effect of element transmutation on mechanical behavior be predicted? Remarkably, such extrapolations are possible in principle for systems governed by stochastic differential equations, thanks to a set of exact relations between probability densities for trajectories derived… ▽ More Can far-from-equilibrium material response under arbitrary loading be inferred from equilibrium data and vice versa? Can the effect of element transmutation on mechanical behavior be predicted? Remarkably, such extrapolations are possible in principle for systems governed by stochastic differential equations, thanks to a set of exact relations between probability densities for trajectories derived from the path integral formalism [1, 2, 3]. In this article, we systematically investigate inferences (in the form of ensemble-averages) drawn on system/process S based on stochastic trajectory data of system/process S~, with quantified uncertainty, to directly address the aforementioned questions. Interestingly, such inferences and their associated uncertainty do not require any simulations or experiments of S. The results are exemplified over two illustrative examples by means of numerical simulations: a one-dimensional system as a prototype for polymers and biological macromolecules, and a two-dimensional glassy system. In principle, the approach can be pushed to the extreme case where S~ is simply comprised of Brownian trajectories, i.e., equilibrium non-interacting particles, and S is a complex interacting system driven far from equilibrium. In practice, however, the "further" S~ is from S, the greater the uncertainty in the predictions, for a fixed number of realizations of system S~. △ Less

Submitted 12 October, 2021; originally announced October 2021.

A multispecies pseudoadiabat for simulating condensable-rich exoplanet atmospheres

Authors: R. J. Graham, Tim Lichtenberg, Ryan Boukrouche, Ray Pierrehumbert

Abstract: Central stages in the evolution of rocky, potentially habitable planets may play out under atmospheric conditions with a large inventory of non-dilute condensable components. Variations in condensate retention and accompanying changes in local lapse rate may substantially affect planetary climate and surface conditions, but there is currently no general theory to effectively describe such atmosphe… ▽ More Central stages in the evolution of rocky, potentially habitable planets may play out under atmospheric conditions with a large inventory of non-dilute condensable components. Variations in condensate retention and accompanying changes in local lapse rate may substantially affect planetary climate and surface conditions, but there is currently no general theory to effectively describe such atmospheres. In this article, expanding on the work by Li et al. (2018), we generalize the single-component moist pseudoadiabat derivation in Pierrehumbert (2010) to allow for multiple condensing components of arbitrary diluteness and retained condensate fraction. The introduction of a freely tunable retained condensate fraction allows for a flexible, self-consistent treatment of atmospheres with non-dilute condensable components. To test the pseudoadiabat's capabilities for simulating a diverse range of climates, we apply the formula to planetary atmospheres with compositions, surface pressures, and temperatures representing important stages with condensable-rich atmospheres in the evolution of terrestrial planets: a magma ocean planet in a runaway greenhouse state; a post-impact, late veneer-analogue planet with a complex atmospheric composition; and an Archean Earth-like planet near the outer edge of the classical circumstellar habitable zone. We find that variations in the retention of multiple non-dilute condensable species can significantly affect the lapse rate and in turn outgoing radiation and the spectral signatures of planetary atmospheres. The presented formulation allows for a more comprehensive treatment of the climate evolution of rocky exoplanets and early Earth analogues. △ Less

Submitted 29 August, 2021; originally announced August 2021.

Comments: Accepted for publication in Planetary Science Journal on 8/24/2021. For a video summary see https://bit.ly/3ztVhfF

Journal ref: The Planetary Science Journal (2021), 2, 207

The application of sub-seasonal to seasonal (S2S) predictions for hydropower forecasting

Authors: Robert M. Graham, Jethro Browell, Douglas Bertram, Christopher J. White

Abstract: Inflow forecasts play an essential role in the management of hydropower reservoirs. Forecasts help operators schedule power generation in advance to maximise economic value, mitigate downstream flood risk, and meet environmental requirements. The horizon of operational inflow forecasts is often limited in range to ~2 weeks ahead, marking the predictability barrier of deterministic weather forecast… ▽ More Inflow forecasts play an essential role in the management of hydropower reservoirs. Forecasts help operators schedule power generation in advance to maximise economic value, mitigate downstream flood risk, and meet environmental requirements. The horizon of operational inflow forecasts is often limited in range to ~2 weeks ahead, marking the predictability barrier of deterministic weather forecasts. Reliable inflow forecasts in the sub-seasonal to seasonal (S2S) range would allow operators to take proactive action to mitigate risks of adverse weather conditions, thereby improving water management and increasing revenue. This study outlines a method of deriving skilful S2S inflow forecasts using a case study reservoir in the Scottish Highlands. We generate ensemble inflow forecasts by training a linear regression model for the observed inflow onto S2S ensemble precipitation predictions from the European Centre for Medium-range Weather Forecasting (ECMWF). Subsequently, post-processing techniques from Ensemble Model Output Statistics are applied to derive calibrated S2S probabilistic inflow forecasts, without the application of a separate hydrological model. We find the S2S probabilistic inflow forecasts hold skill relative to climatological forecasts up to 6 weeks ahead. The inflow forecasts hold greater skill during winter compared with summer. The forecasts, however, struggle to predict high summer inflows, even at short lead-times. The potential for the S2S probabilistic inflow forecasts to improve water management and deliver increased economic value is confirmed using a stylised cost model. While applied to hydropower forecasting, the results and methods presented here are relevant to broader fields of water management and S2S forecasting applications. △ Less

Submitted 13 August, 2021; originally announced August 2021.

Relationships among structure, memory, and flow in sheared disordered materials

Authors: K. L. Galloway, E. G. Teich, X-g Ma, Ch. Kammer, I. R. Graham, N. C. Keim, C. Reina, D. J. Jerolmack, A. G. Yodh, P. E. Arratia

Abstract: A fundamental challenge for disordered solids is predicting macroscopic yield from the microscopic arrangements of constituent particles. Yield is accompanied by a sudden and large increase in energy dissipation due to the onset of plastic rearrangements. This suggests that one path to understanding bulk rheology is to map particle configurations to their mode of deformation. Here, we perform labo… ▽ More A fundamental challenge for disordered solids is predicting macroscopic yield from the microscopic arrangements of constituent particles. Yield is accompanied by a sudden and large increase in energy dissipation due to the onset of plastic rearrangements. This suggests that one path to understanding bulk rheology is to map particle configurations to their mode of deformation. Here, we perform laboratory experiments and numerical simulations that are designed to do just that: 2D dense colloidal systems are subjected to oscillatory shear, and particle trajectories and bulk rheology are measured. We quantify particle microstructure using excess entropy. Results reveal a direct relation between excess entropy and energy dissipation, that is insensitive to the nature of interactions among particles. We use this relation to build a physically-informed model that connects rheology to microstructure. Our findings suggest a framework for tailoring the rheological response of disordered materials by tuning microstructural properties. △ Less

Submitted 13 May, 2021; originally announced May 2021.

High pCO$_2$ reduces sensitivity to CO$_2$ perturbations on temperate, Earth-like planets throughout most of habitable zone

Authors: R. J. Graham

Abstract: The nearly logarithmic radiative impact of CO$_2$ means that planets near the outer edge of the liquid water habitable zone (HZ) require $\sim$10$^6$x more CO$_2$ to maintain temperatures conducive to standing liquid water on the planetary surface than their counterparts near the inner edge. This logarithmic radiative response also means that atmospheric CO$_2$ changes of a given mass will have sm… ▽ More The nearly logarithmic radiative impact of CO$_2$ means that planets near the outer edge of the liquid water habitable zone (HZ) require $\sim$10$^6$x more CO$_2$ to maintain temperatures conducive to standing liquid water on the planetary surface than their counterparts near the inner edge. This logarithmic radiative response also means that atmospheric CO$_2$ changes of a given mass will have smaller temperature effects on higher pCO$_2$ planets. Ocean pH is linked to atmospheric pCO$_2$ through seawater carbonate speciation and calcium carbonate dissolution/precipitation, and the response of pH to changes in pCO$_2$ also decreases at higher initial pCO$_2$. Here, we use idealized climate and ocean chemistry models to demonstrate that CO$_2$ perturbations large enough to cause catastrophic changes to surface temperature and ocean pH on low-pCO$_2$ planets in the innermost region of the HZ are likely to have much smaller effects on planets with higher pCO$_2$. Major bouts of extraterrestrial fossil fuel combustion or volcanic CO$_2$ outgassing on high-pCO$_2$ planets in the mid-to-outer HZ should have mild or negligible impacts on surface temperature and ocean pH. Owing to low pCO$_2$, Phanerozoic Earth's surface environment may be unusually volatile compared to similar planets receiving lower instellation. △ Less

Submitted 2 April, 2021; originally announced April 2021.

Comments: Submitted November 30, 2020; Reviews received March 23rd, 2021; Resubmitted April 2, 2021

Guessing about Guessing: Practical Strategies for Card Guessing with Feedback

Authors: Persi Diaconis, Ron Graham, Sam Spiro

Abstract: In simple card games, cards are dealt one at a time and the player guesses each card sequentially. We study problems where feedback (e.g. correct/incorrect) is given after each guess. For decks with repeated values (as in blackjack where suits do not matter) the optimal strategy differs from the "greedy strategy" (of guessing a most likely card each round). Further, both optimal and greedy strateg… ▽ More In simple card games, cards are dealt one at a time and the player guesses each card sequentially. We study problems where feedback (e.g. correct/incorrect) is given after each guess. For decks with repeated values (as in blackjack where suits do not matter) the optimal strategy differs from the "greedy strategy" (of guessing a most likely card each round). Further, both optimal and greedy strategies are far too complicated for real time use by human players. Our main results show that simple heuristics perform close to optimal. △ Less

Submitted 17 July, 2021; v1 submitted 7 December, 2020; originally announced December 2020.

Comments: 20 pages, minor typos corrected, to appear in the American Mathematical Monthly excluding Section 6

MSC Class: 60C05

Gemini Planet Imager Spectroscopy of the Dusty Substellar Companion HD 206893 B

Authors: K. Ward-Duong, J. Patience, K. Follette, R. J. De Rosa, J. Rameau, M. Marley, D. Saumon, E. L. Nielsen, A. Rajan, A. Z. Greenbaum, J. Lee, J. J. Wang, I. Czekala, G. Duchêne, B. Macintosh, S. Mark Ammons, V. P. Bailey, T. Barman, J. Bulger, C. Chen, J. Chilcote, T. Cotten, R. Doyon, T. M. Esposito, M. P. Fitzgerald , et al. (33 additional authors not shown)

Abstract: We present new near-infrared Gemini Planet Imager (GPI) spectroscopy of HD 206893 B, a substellar companion orbiting within the debris disk of its F5V star. The $J$, $H$, $K1$, and $K2$ spectra from GPI demonstrate the extraordinarily red colors of the object, confirming it as the reddest substellar object observed to date. The significant flux increase throughout the infrared presents a challengi… ▽ More We present new near-infrared Gemini Planet Imager (GPI) spectroscopy of HD 206893 B, a substellar companion orbiting within the debris disk of its F5V star. The $J$, $H$, $K1$, and $K2$ spectra from GPI demonstrate the extraordinarily red colors of the object, confirming it as the reddest substellar object observed to date. The significant flux increase throughout the infrared presents a challenging atmosphere to model with existing grids. Best-fit values vary from 1200 K to 1800 K for effective temperature and from 3.0 to 5.0 for log($g$), depending on which individual wavelength band is fit and which model suite is applied. The extreme redness of the companion can be partially reconciled by invoking a high-altitude layer of sub-micron dust particles, similar to dereddening approaches applied to the peculiar red field L-dwarf population. However, reconciling the HD 206893 B spectra with even those of the reddest low-gravity L-dwarf spectra still requires the contribution of additional atmospheric dust, potentially due to the debris disk environment in which the companion resides. Orbit fitting from four years of astrometric monitoring is consistent with a $\sim$30-year period, orbital inclination of 147$^{\circ}$, and semimajor axis of 10 au, well within the estimated disk inner radius of $\sim$50 au. As one of very few substellar companions imaged interior to a circumstellar disk, the properties of this system offer important dynamical constraints on companion-disk interaction and provide a benchmark for substellar and planetary atmospheric study. △ Less

Submitted 20 October, 2020; originally announced October 2020.

Comments: Accepted for publication in AJ (October 15, 2020). 35 pages, 20 figures, with tables and appendices presented in their entirety

Card Guessing with Partial Feedback

Authors: Persi Diaconis, Ron Graham, Xiaoyu He, Sam Spiro

Abstract: Consider the following experiment: a deck with $m$ copies of $n$ different card types is randomly shuffled, and a guesser attempts to guess the cards sequentially as they are drawn. Each time a guess is made, some amount of "feedback" is given. For example, one could tell the guesser the true identity of the card they just guessed (the complete feedback model) or they could be told nothing at all… ▽ More Consider the following experiment: a deck with $m$ copies of $n$ different card types is randomly shuffled, and a guesser attempts to guess the cards sequentially as they are drawn. Each time a guess is made, some amount of "feedback" is given. For example, one could tell the guesser the true identity of the card they just guessed (the complete feedback model) or they could be told nothing at all (the no feedback model). In this paper we explore a partial feedback model, where upon guessing a card, the guesser is only told whether or not their guess was correct. We show in this setting that, uniformly in $n$, at most $m+O(m^{3/4}\log m)$ cards can be guessed correctly in expectation. This resolves a question of Diaconis and Graham from 1981, where even the $m=2$ case was open. △ Less

Submitted 18 July, 2022; v1 submitted 10 October, 2020; originally announced October 2020.

Comments: 24 pages; minor errors corrected

MSC Class: 60C05

Journal ref: Combinator. Probab. Comp. 31 (2022) 1-20

Multiband Polarimetric Imaging of HR 4796A with the Gemini Planet Imager

Authors: Pauline Arriaga, Michael P. Fitzgerald, Gaspard Duchêne, Paul Kalas, Maxwell A. Millar-Blanchaer, Marshall D. Perrin, Christine H. Chen, Johan Mazoyer, Mark Ammons, Vanessa P. Bailey, Trafis S. Barman, Joanna Bulger, Jeffrey K. Chilcote, Tara Cotten, Robert J. De Rosa, Rene Doyon, Thomas M. Esposito, Katherine B. Follette, Benjamin L. Gerard, Stephen Goodsell, James R. Graham, Alexandra Z. Greenbaum, Pascale Hibon, Justin Hom, Li-Wei Hung , et al. (27 additional authors not shown)

Abstract: HR4796A hosts a well-studied debris disk with a long history due to its high fractional luminosity and favorable inclination lending itself well to both unresolved and resolved observations. We present new J- and K1-band images of the resolved debris disk HR4796A taken in the polarimetric mode of the Gemini Planet Imager (GPI). The polarized intensity features a strongly forward scattered brightne… ▽ More HR4796A hosts a well-studied debris disk with a long history due to its high fractional luminosity and favorable inclination lending itself well to both unresolved and resolved observations. We present new J- and K1-band images of the resolved debris disk HR4796A taken in the polarimetric mode of the Gemini Planet Imager (GPI). The polarized intensity features a strongly forward scattered brightness distribution and is undetected at the far side of the disk. The total intensity is detected at all scattering angles and also exhibits a strong forward scattering peak. We use a forward modelled geometric disk in order to extract geometric parameters, polarized fraction and total intensity scattering phase functions for these data as well as H-band data previously taken by GPI. We find the polarized phase function becomes increasingly more forward scattering as wavelength increases. We fit Mie and distribution of hollow spheres grain (DHS) models to the extracted functions. We find that while it is possible to describe generate a satisfactory model for the total intensity using a DHS model, but not with a Mie model. We find that no single grain population of DHS or Mie grains of arbitrary composition can simultaneously reproduce the polarized fraction and total intensity scattering phase functions, indicating the need for more sophisticated grain models. △ Less

Submitted 11 June, 2020; originally announced June 2020.

Comments: Accepted for publication in AJ. 3 tables, 11 figures

High-resolution observations of the solar photosphere, chromosphere and transition region. A database of coordinated IRIS and SST observations

Authors: L. H. M. Rouppe van der Voort, B. De Pontieu, M. Carlsson, J. de la Cruz Rodriguez, S. Bose, G. Chintzoglou, A. Drews, C. Froment, M. Gosic, D. R. Graham, V. H. Hansteen, V. M. J. Henriques, S. Jafarzadeh, J. Joshi, L. Kleint, P. Kohutova, T. Leifsen, J. Martinez-Sykora, D. Nobrega-Siverio, A. Ortiz, T. M. D. Pereira, A. Popovas, C. Quintero Noda, A. Sainz Dalda, G. B. Scharmer , et al. (8 additional authors not shown)

Abstract: NASA's Interface Region Imaging Spectrograph (IRIS) provides high resolution observations of the solar atmosphere through UV spectroscopy and imaging. Since the launch of IRIS in June 2013, we have conducted systematic observation campaigns in coordination with the Swedish 1-m Solar Telescope (SST) on La Palma. The SST provides complementary high-resolution observations of the photosphere and chro… ▽ More NASA's Interface Region Imaging Spectrograph (IRIS) provides high resolution observations of the solar atmosphere through UV spectroscopy and imaging. Since the launch of IRIS in June 2013, we have conducted systematic observation campaigns in coordination with the Swedish 1-m Solar Telescope (SST) on La Palma. The SST provides complementary high-resolution observations of the photosphere and chromosphere. The SST observations include spectro-polarimetric imaging in photospheric Fe I lines and spectrally-resolved imaging in the chromospheric Ca II 8542 A, H-alpha, and Ca II K lines. We present a database of co-aligned IRIS and SST datasets that is open for analysis to the scientific community. The database covers a variety of targets including active regions, sunspots, plage, quiet Sun, and coronal holes. △ Less

Submitted 29 July, 2020; v1 submitted 28 May, 2020; originally announced May 2020.

Comments: Accepted for publication in Astronomy & Astrophysics (July 21, 2020)

Journal ref: A&A 641, A146 (2020)

Thermodynamic and Energetic Limits on Continental Silicate Weathering Strongly Impact the Climate and Habitability of Wet, Rocky Worlds

Authors: R. J. Graham, R. T. Pierrehumbert

Abstract: The "liquid water habitable zone" (HZ) concept is predicated on the ability of the silicate weathering feedback to stabilize climate across a wide range of instellations. However, representations of silicate weathering used in current estimates of the effective outer edge of the HZ do not account for the thermodynamic limit on concentration of weathering products in runoff set by clay precipitatio… ▽ More The "liquid water habitable zone" (HZ) concept is predicated on the ability of the silicate weathering feedback to stabilize climate across a wide range of instellations. However, representations of silicate weathering used in current estimates of the effective outer edge of the HZ do not account for the thermodynamic limit on concentration of weathering products in runoff set by clay precipitation, nor for the energetic limit on precipitation set by planetary instellation. We find that when the thermodynamic limit is included in an idealized coupled climate/weathering model, steady-state planetary climate loses sensitivity to silicate dissolution kinetics, becoming sensitive to temperature primarily through the effect of temperature on runoff and to pCO$_2$ through an effect on solute concentration mediated by pH. This increases sensitivity to land fraction, CO$_2$ outgassing, and geological factors such as soil age and lithology, all of which are found to have a profound effect on the position of the effective outer edge of the HZ. The interplay between runoff sensitivity and the energetic limit on precipitation leads to novel warm states in the outer reaches of the HZ, owing to the decoupling of temperature and precipitation. We discuss strategies for detecting the signature of silicate weathering feedback through exoplanet observations in light of insights derived from the revised picture of weathering. △ Less

Submitted 20 May, 2020; v1 submitted 29 April, 2020; originally announced April 2020.

Comments: 27 pages excluding references, 30 pages total, 8 figures. Submitted to ApJ, accepted 15 May 2020. (reuploaded to arxiv 5/20/2020 to include corresponding author's email address and news of acceptance)

Debris Disk Results from the Gemini Planet Imager Exoplanet Survey's Polarimetric Imaging Campaign

Authors: Thomas M. Esposito, Paul Kalas, Michael P. Fitzgerald, Maxwell A. Millar-Blanchaer, Gaspard Duchene, Jennifer Patience, Justin Hom, Marshall D. Perrin, Robert J. De Rosa, Eugene Chiang, Ian Czekala, Bruce Macintosh, James R. Graham, Megan Ansdell, Pauline Arriaga, Sebastian Bruzzone, Joanna Bulger, Christine H. Chen, Tara Cotten, Ruobing Dong, Zachary H. Draper, Katherine B. Follette, Li-Wei Hung, Ronald Lopez, Brenda C. Matthews , et al. (40 additional authors not shown)

Abstract: We report the results of a ${\sim}4$-year direct imaging survey of 104 stars to resolve and characterize circumstellar debris disks in scattered light as part of the Gemini Planet Imager Exoplanet Survey. We targeted nearby (${\lesssim}150$ pc), young (${\lesssim}500$ Myr) stars with high infrared excesses ($L_{\mathrm{IR}} / L_\star > 10^{-5}$), including 38 with previously resolved disks. Observ… ▽ More We report the results of a ${\sim}4$-year direct imaging survey of 104 stars to resolve and characterize circumstellar debris disks in scattered light as part of the Gemini Planet Imager Exoplanet Survey. We targeted nearby (${\lesssim}150$ pc), young (${\lesssim}500$ Myr) stars with high infrared excesses ($L_{\mathrm{IR}} / L_\star > 10^{-5}$), including 38 with previously resolved disks. Observations were made using the Gemini Planet Imager high-contrast integral field spectrograph in $H$-band (1.6 $μ$m) coronagraphic polarimetry mode to measure both polarized and total intensities. We resolved 26 debris disks and three protoplanetary/transitional disks. Seven debris disks were resolved in scattered light for the first time, including newly presented HD 117214 and HD 156623, and we quantified basic morphologies of five of them using radiative transfer models. All of our detected debris disks but HD 156623 have dust-poor inner holes, and their scattered-light radii are generally larger than corresponding radii measured from resolved thermal emission and those inferred from spectral energy distributions. To assess sensitivity, we report contrasts and consider causes of non-detections. Detections were strongly correlated with high IR excess and high inclination, although polarimetry outperformed total intensity angular differential imaging for detecting low inclination disks (${\lesssim} 70 °$). Based on post-survey statistics, we improved upon our pre-survey target prioritization metric predicting polarimetric disk detectability. We also examined scattered-light disks in the contexts of gas, far-IR, and millimeter detections. Comparing $H$-band and ALMA fluxes for two disks revealed tentative evidence for differing grain properties. Finally, we found no preference for debris disks to be detected in scattered light if wide-separation substellar companions were present. △ Less

Submitted 23 June, 2020; v1 submitted 28 April, 2020; originally announced April 2020.

Comments: arXiv resubmission with typographical corrections. Accepted for publication in AJ. 19 figures, 7 tables

The Gemini Planet Imager view of the HD 32297 debris disk

Authors: Gaspard Duchene, Malena Rice, Justin Hom, Joseph Zalesky, Thomas M. Esposito, Maxwell A. Millar-Blanchaer, Bin Ren, Paul Kalas, Michael Fitzgerald, Pauline Arriaga, Sebastian Bruzzone, Joanna Bulger, Christine H. Chen, Eugene Chiang, Tara Cotten, Ian Czekala, Robert J. De Rosa, Ruobing Dong, Zachary H. Draper, Katherine B. Follette, James R. Graham, Li-Wei Hung, Ronald Lopez, Bruce Macintosh, Brenda C. Matthews , et al. (38 additional authors not shown)

Abstract: We present new $H$-band scattered light images of the HD 32297 edge-on debris disk obtained with the Gemini Planet Imager (GPI). The disk is detected in total and polarized intensity down to a projected angular separation of 0.15", or 20au. On the other hand, the large scale swept-back halo remains undetected, likely a consequence of its markedly blue color relative to the parent body belt. We ana… ▽ More We present new $H$-band scattered light images of the HD 32297 edge-on debris disk obtained with the Gemini Planet Imager (GPI). The disk is detected in total and polarized intensity down to a projected angular separation of 0.15", or 20au. On the other hand, the large scale swept-back halo remains undetected, likely a consequence of its markedly blue color relative to the parent body belt. We analyze the curvature of the disk spine and estimate a radius of $\approx$100au for the parent body belt, smaller than past scattered light studies but consistent with thermal emission maps of the system. We employ three different flux-preserving post-processing methods to suppress the residual starlight and evaluate the surface brightness and polarization profile along the disk spine. Unlike past studies of the system, our high fidelity images reveal the disk to be highly symmetric and devoid of morphological and surface brightness perturbations. We find the dust scattering properties of the system to be consistent with those observed in other debris disks, with the exception of HR 4796. Finally, we find no direct evidence for the presence of a planetary-mass object in the system. △ Less

Submitted 13 April, 2020; originally announced April 2020.

Comments: Accepted for publication in the Astronomical Journal

Spectral signatures of chromospheric condensation in a major solar flare

Authors: David R. Graham, Gianna Cauzzi, Luca Zangrilli, Adam Kowalski, Paulo Simões, Joel Allred

Abstract: We study the evolution of chromospheric line and continuum emission during the impulsive phase of the X-class SOL2014-09-10T17:45 solar flare. We extend previous analyses of this flare to multiple chromospheric lines of Fe I, Fe II, Mg II, C I, and Si II, observed with IRIS, combined with radiative-hydrodynamical (RHD) modeling. For multiple flaring kernels, the lines all show a rapidly evolving d… ▽ More We study the evolution of chromospheric line and continuum emission during the impulsive phase of the X-class SOL2014-09-10T17:45 solar flare. We extend previous analyses of this flare to multiple chromospheric lines of Fe I, Fe II, Mg II, C I, and Si II, observed with IRIS, combined with radiative-hydrodynamical (RHD) modeling. For multiple flaring kernels, the lines all show a rapidly evolving double-component structure: an enhanced, emission component at rest, and a broad, highly red-shifted component of comparable intensity. The red-shifted components migrate from 25-50 km s$^{-1}$ towards the rest wavelength within $\sim$30 seconds. Using Fermi hard X-ray observations, we derive the parameters of an accelerated electron beam impacting the dense chromosphere, using them to drive a RHD simulation with the RADYN code. As in Kowalski et al. 2017a, our simulations show that the most energetic electrons penetrate into the deep chromosphere, heating it to T$\sim$10,000 K, while the bulk of the electrons dissipate their energy higher, driving an explosive evaporation, and its counterpart condensation -- a very dense (n$_e \sim 2 \times 10^{14}$ cm$^{-3}$), thin layer (30--40 km thickness), heated to 8--12,000 K, moving towards the stationary chromosphere at up to 50 km s$^{-1}$. The synthetic Fe II 2814.45A profiles closely resemble the observational data, including a continuum enhancement, and both a stationary and a highly red-shifted component, rapidly moving towards the rest wavelength. Importantly, the absolute continuum intensity, ratio of component intensities, relative time of appearance, and red-shift amplitude, are sensitive to the model input parameters, showing great potential as diagnostics. △ Less

Submitted 10 April, 2020; originally announced April 2020.

Comments: Accepted for publication by the Astrophysical Journal - 9th April 2020

HD 165054: an astrometric calibration field for high-contrast imagers in Baade's Window

Authors: Meiji M. Nguyen, Robert J. De Rosa, Jason J. Wang, Thomas M. Esposito, Paul Kalas, James R. Graham, Bruce Macintosh, Vanessa P. Bailey, Travis Barman, Joanna Bulger, Jeffrey Chilcote, Tara Cotten, Rene Doyon, Gaspard Duchêne, Michael P. Fitzgerald, Katherine B. Follette, Benjamin L. Gerard, Stephen J. Goodsell, Alexandra Z. Greenbaum, Pascale Hibon, Justin Hom, Li-Wei Hung, Patrick Ingraham, Quinn Konopacky, James E. Larkin , et al. (29 additional authors not shown)

Abstract: We present a study of the HD 165054 astrometric calibration field that has been periodically observed with the Gemini Planet Imager. HD 165054 is a bright star within Baade's Window, a region of the galactic plane with relatively low extinction from interstellar dust. HD 165054 was selected as a calibrator target due to the high number density of stars within this region ($\sim 3$ stars per square… ▽ More We present a study of the HD 165054 astrometric calibration field that has been periodically observed with the Gemini Planet Imager. HD 165054 is a bright star within Baade's Window, a region of the galactic plane with relatively low extinction from interstellar dust. HD 165054 was selected as a calibrator target due to the high number density of stars within this region ($\sim 3$ stars per square arcsecond with $H<22$), necessary because of the small field-of-view of the Gemini Planet Imager. Using nine epochs spanning over five years, we have fit a standard five-parameter astrometric model to the astrometry of seven background stars within close proximity to HD 165054 (angular separation $< 2$ arcsec). We achieved a proper motion precision of $\sim 0.3$ mas/yr, and constrained the parallax of each star to be $\lesssim 1$ mas. Our measured proper motions and parallax limits are consistent with the background stars being a part of the galactic bulge. Using these measurements we find no evidence of any systematic trend of either the plate scale or the north angle offset of GPI between 2014 and 2019. We compared our model describing the motions of the seven background stars to observations of the same field in 2014 and 2018 obtained with Keck/NIRC2, an instrument with an excellent astrometric calibration. We find that predicted position of the background sources is consistent with that measured by NIRC2, within the uncertainties of the calibration of the two instruments. In the future, we will use this field as a standard astrometric calibrator for the upgrade of GPI and potentially for other high-contrast imagers. △ Less

Submitted 6 April, 2020; originally announced April 2020.

Comments: 25 pages, 12 figures

Asymptotically Hyperbolic Manifolds with Boundary Conjugate Points but no Interior Conjugate Points

Authors: Nikolas Eptaminitakis, C. Robin Graham

Abstract: We construct non-trapping asymptotically hyperbolic manifolds with boundary conjugate points but no interior conjugate points. We construct non-trapping asymptotically hyperbolic manifolds with boundary conjugate points but no interior conjugate points. △ Less

Submitted 10 December, 2019; originally announced December 2019.

Comments: 23 pages

MSC Class: 53C20; 53C22

Imaging the 44 AU Kuiper Belt-analogue debris ring around HD 141569A with GPI polarimetry

Authors: J. S. Bruzzone, S. Metchev, G. Duchene, M. A. Millar-Blanchaer, R. Dong, J. J. Wang, J. R. Graham, J. Mazoyer, S. Wolff, S. M. Ammons, A. C. Schneider, A. Z. Greenbaum, B. C. Matthews, P. Arriaga, V. P. Bailey, T. Barman, J. Bulger, J. Chilcote, T. Cotten, R. J. De Rosa, R. Doyon, M. P. Fitzgerald, K. B. Follette, B. L. Gerard, S. J. Goodsell , et al. (31 additional authors not shown)

Abstract: We present the first polarimetric detection of the inner disk component around the pre-main sequence B9.5 star HD 141569A. Gemini Planet Imager H-band (1.65 micron) polarimetric differential imaging reveals the highest signal-to-noise ratio detection of this ring yet attained and traces structure inwards to 0.25" (28 AU at a distance of 111 pc). The radial polarized intensity image shows the east… ▽ More We present the first polarimetric detection of the inner disk component around the pre-main sequence B9.5 star HD 141569A. Gemini Planet Imager H-band (1.65 micron) polarimetric differential imaging reveals the highest signal-to-noise ratio detection of this ring yet attained and traces structure inwards to 0.25" (28 AU at a distance of 111 pc). The radial polarized intensity image shows the east side of the disk, peaking in intensity at 0.40" (44 AU) and extending out to 0.9" (100 AU). There is a spiral arm-like enhancement to the south, reminiscent of the known spiral structures on the outer rings of the disk. The location of the spiral arm is coincident with 12CO J=3-2 emission detected by ALMA, and hints at a dynamically active inner circumstellar region. Our observations also show a portion of the middle dusty ring at ~220 AU known from previous observations of this system. We fit the polarized H-band emission with a continuum radiative transfer Mie model. Our best-fit model favors an optically thin disk with a minimum dust grain size close to the blow-out size for this system: evidence of on-going dust production in the inner reaches of the disk. The thermal emission from this model accounts for virtually all of the far-infrared and millimeter flux from the entire HD 141569A disk, in agreement with the lack of ALMA continuum and CO emission beyond ~100 AU. A remaining 8-30 micron thermal excess a factor of ~2 above our model argues for a yet-unresolved warm innermost 5-15 AU component of the disk. △ Less

Submitted 26 November, 2019; originally announced November 2019.

Comments: 16 pages, 8 figures, accepted at AJ

The Gemini Planet Imager Exoplanet Survey: Dynamical Mass of the Exoplanet beta Pictoris b from Combined Direct Imaging and Astrometry

Authors: Eric L. Nielsen, Robert J. De Rosa, Jason J. Wang, Johannes Sahlmann, Paul Kalas, Gaspard Duchene, Julien Rameau, Mark S. Marley, Didier Saumon, Bruce Macintosh, Maxwell A. Millar-Blanchaer, Meiji M. Nguyen, S. Mark Ammons, Vanessa P. Bailey, Travis Barman, Joanna Bulger, Jeffrey Chilcote, Tara Cotten, Rene Doyon, Thomas M. Esposito, Michael P. Fitzgerald, Katherine B. Follette, Benjamin L. Gerard, Stephen J. Goodsell, James R. Graham , et al. (29 additional authors not shown)

Abstract: We present new observations of the planet beta Pictoris b from 2018 with GPI, the first GPI observations following conjunction. Based on these new measurements, we perform a joint orbit fit to the available relative astrometry from ground-based imaging, the Hipparcos Intermediate Astrometric Data (IAD), and the Gaia DR2 position, and demonstrate how to incorporate the IAD into direct imaging orbit… ▽ More We present new observations of the planet beta Pictoris b from 2018 with GPI, the first GPI observations following conjunction. Based on these new measurements, we perform a joint orbit fit to the available relative astrometry from ground-based imaging, the Hipparcos Intermediate Astrometric Data (IAD), and the Gaia DR2 position, and demonstrate how to incorporate the IAD into direct imaging orbit fits. We find a mass consistent with predictions of hot-start evolutionary models and previous works following similar methods, though with larger uncertainties: 12.8 [+5.3, -3.2] M_Jup. Our eccentricity determination of 0.12 [+0.04, -0.03] disfavors circular orbits. We consider orbit fits to several different imaging datasets, and find generally similar posteriors on the mass for each combination of imaging data. Our analysis underscores the importance of performing joint fits to the absolute and relative astrometry simultaneously, given the strong covariance between orbital elements. Time of conjunction is well constrained within 2.8 days of 2017 September 13, with the star behind the planet's Hill sphere between 2017 April 11 and 2018 February 16 (+/- 18 days). Following the recent radial velocity detection of a second planet in the system, beta Pic c, we perform additional two-planet fits combining relative astrometry, absolute astrometry, and stellar radial velocities. These joint fits find a significantly smaller mass for the imaged planet beta Pic b, of 8.0 +/- 2.6 M_Jup, in a somewhat more circular orbit. We expect future ground-based observations to further constrain the visual orbit and mass of the planet in advance of the release of Gaia DR4. △ Less

Submitted 25 November, 2019; originally announced November 2019.

Comments: 29 Pages, 23 Figures. Accepted to AJ

First Resolved Scattered-Light Images of Four Debris Disks in Scorpius-Centaurus with the Gemini Planet Imager

Authors: Justin Hom, Jennifer Patience, Thomas M. Esposito, Gaspard Duchêne, Kadin Worthen, Paul Kalas, Hannah Jang-Condell, Kezman Saboi, Pauline Arriaga, Johan Mazoyer, Schuyler Wolff, Maxwell A. Millar-Blanchaer, Michael P. Fitzgerald, Marshall D. Perrin, Christine H. Chen, Bruce Macintosh, Brenda C. Matthews, Jason J. Wang, James R. Graham, Franck Marchis, S. Mark Ammons, Vanessa P. Bailey, Travis Barman, Joanna Bulger, Jeffrey K. Chilcote , et al. (35 additional authors not shown)

Abstract: We present the first spatially resolved scattered-light images of four debris disks around members of the Scorpius-Centaurus (Sco-Cen) OB Association with high-contrast imaging and polarimetry using the Gemini Planet Imager (GPI). All four disks are resolved for the first time in polarized light and one disk is also detected in total intensity. The three disks imaged around HD 111161, HD 143675, a… ▽ More We present the first spatially resolved scattered-light images of four debris disks around members of the Scorpius-Centaurus (Sco-Cen) OB Association with high-contrast imaging and polarimetry using the Gemini Planet Imager (GPI). All four disks are resolved for the first time in polarized light and one disk is also detected in total intensity. The three disks imaged around HD 111161, HD 143675, and HD 145560 are symmetric in both morphology and brightness distribution. The three systems span a range of inclinations and radial extents. The disk imaged around HD 98363 shows indications of asymmetries in morphology and brightness distribution, with some structural similarities to the HD 106906 planet-disk system. Uniquely, HD 98363 has a wide co-moving stellar companion Wray 15-788 with a recently resolved disk with very different morphological properties. HD 98363 A/B is the first binary debris disk system with two spatially resolved disks. All four targets have been observed with ALMA, and their continuum fluxes range from one non-detection to one of the brightest disks in the region. With the new results, a total of 15 A/F-stars in Sco-Cen have resolved scattered light debris disks, and approximately half of these systems exhibit some form of asymmetry. Combining the GPI disk structure results with information from the literature on millimeter fluxes and imaged planets reveals a diversity of disk properties in this young population. Overall, the four newly resolved disks contribute to the census of disk structures measured around A/F-stars at this important stage in the development of planetary systems. △ Less

Submitted 21 November, 2019; originally announced November 2019.

Comments: 20 pages, 9 figures

Detection of a low-mass stellar companion to the accelerating A2IV star HR 1645

Authors: Robert J. De Rosa, Eric L. Nielsen, Julien Rameau, Gaspard Duchêne, Alexandra Z. Greenbaum, Jason J. Wang, S. Mark Ammons, Vanessa P. Bailey, Travis Barman, Joanna Bulger, Jeffrey Chilcote, Tara Cotten, Rene Doyon, Thomas M. Esposito, Michael P. Fitzgerald, Katherine B. Follette, Benjamin L. Gerard, Stephen J. Goodsell, James R. Graham, Pascale Hibon, Justin Hom, Li-Wei Hung, Patrick Ingraham, Paul Kalas Quinn Konopacky, James E. Larkin , et al. (28 additional authors not shown)

Abstract: The $\sim500$\, Myr A2IV star HR 1645 has one of the most significant low-amplitude accelerations of nearby early-type stars measured from a comparison of the {\it Hipparcos} and {\it Gaia} astrometric catalogues. This signal is consistent with either a stellar companion with a moderate mass ratio ($q\sim0.5$) on a short period ($P<1$\,yr), or a substellar companion at a separation wide enough to… ▽ More The $\sim500$\, Myr A2IV star HR 1645 has one of the most significant low-amplitude accelerations of nearby early-type stars measured from a comparison of the {\it Hipparcos} and {\it Gaia} astrometric catalogues. This signal is consistent with either a stellar companion with a moderate mass ratio ($q\sim0.5$) on a short period ($P<1$\,yr), or a substellar companion at a separation wide enough to be resolved with ground-based high contrast imaging instruments; long-period equal mass ratio stellar companions that are also consistent with the measured acceleration are excluded with previous imaging observations. The small but significant amplitude of the acceleration made HR 1645 a promising candidate for targeted searches for brown dwarf and planetary-mass companions around nearby, young stars. In this paper we explore the origin of the astrometric acceleration by modelling the signal induced by a wide-orbit M8 companion discovered with the Gemini Planet Imager, as well as the effects of an inner short-period spectroscopic companion discovered a century ago but not since followed-up. We present the first constraints on the orbit of the inner companion, and demonstrate that it is a plausible cause of the astrometric acceleration. This result demonstrates the importance of vetting of targets with measured astrometric acceleration for short-period stellar companions prior to conducting targeted direct imaging surveys for wide-orbit substellar companions. △ Less

Submitted 22 October, 2019; originally announced October 2019.

Comments: 14 pages, 9 figures. Accepted for publication in the Astronomical Journal

An updated visual orbit of the directly-imaged exoplanet 51 Eridani b and prospects for a dynamical mass measurement with Gaia

Authors: Robert J. De Rosa, Eric L. Nielsen, Jason J. Wang, S. Mark Ammons, Gaspard Duchêne, Bruce Macintosh, Meiji M. Nguyen, Julien Rameau, Vanessa P. Bailey, Travis Barman, Joanna Bulger, Jeffrey Chilcote, Tara Cotten, Rene Doyon, Thomas M. Esposito, Michael P. Fitzgerald, Katherine B. Follette, Benjamin L. Gerard, Stephen J. Goodsell, James R. Graham, Alexandra Z. Greenbaum, Pascale Hibon, Justin Hom, Li-Wei Hung, Patrick Ingraham , et al. (30 additional authors not shown)

Abstract: We present a revision to the visual orbit of the young, directly-imaged exoplanet 51 Eridani b using four years of observations with the Gemini Planet Imager. The relative astrometry is consistent with an eccentric ($e=0.53_{-0.13}^{+0.09}$) orbit at an intermediate inclination ($i=136_{-11}^{+10}$\,deg), although circular orbits cannot be excluded due to the complex shape of the multidimensional… ▽ More We present a revision to the visual orbit of the young, directly-imaged exoplanet 51 Eridani b using four years of observations with the Gemini Planet Imager. The relative astrometry is consistent with an eccentric ($e=0.53_{-0.13}^{+0.09}$) orbit at an intermediate inclination ($i=136_{-11}^{+10}$\,deg), although circular orbits cannot be excluded due to the complex shape of the multidimensional posterior distribution. We find a semi-major axis of $11.1_{-1.3}^{+4.2}$\,au and a period of $28.1_{-4.9}^{+17.2}$\,yr, assuming a mass of 1.75\,M$_{\odot}$ for the host star. We find consistent values with a recent analysis of VLT/SPHERE data covering a similar baseline. We investigated the potential of using absolute astrometry of the host star to obtain a dynamical mass constraint for the planet. The astrometric acceleration of 51~Eri derived from a comparison of the {\it Hipparcos} and {\it Gaia} catalogues was found to be inconsistent at the 2--3$σ$ level with the predicted reflex motion induced by the orbiting planet. Potential sources of this inconsistency include a combination of random and systematic errors between the two astrometric catalogs or the signature of an additional companion within the system interior to current detection limits. We also explored the potential of using {\it Gaia} astrometry alone for a dynamical mass measurement of the planet by simulating {\it Gaia} measurements of the motion of the photocenter of the system over the course of the extended eight-year mission. We find that such a measurement is only possible ($>98$\% probability) given the most optimistic predictions for the {\it Gaia} scan astrometric uncertainties for bright stars, and a high mass for the planet ($\gtrsim3.6$\,M$_{\rm Jup}$). △ Less

Submitted 22 October, 2019; originally announced October 2019.

Comments: 17 pages, 11 figures. Accepted for publication in the Astronomical Journal

Revised Astrometric Calibration of the Gemini Planet Imager

Authors: Robert J. De Rosa, Meiji M. Nguyen, Jeffrey Chilcote, Bruce Macintosh, Marshall D. Perrin, Quinn Konopacky, Jason J. Wang, Gaspard Duchêne, Eric L. Nielsen, Julien Rameau, S. Mark Ammons, Vanessa P. Bailey, Travis Barman, Joanna Bulger, Tara Cotten, Rene Doyon, Thomas M. Esposito, Michael P. Fitzgerald, Katherine B. Follette, Benjamin L. Gerard, Stephen J. Goodsell, James R. Graham, Alexandra Z. Greenbaum, Pascale Hibon, Li-Wei Hung , et al. (27 additional authors not shown)

Abstract: We present a revision to the astrometric calibration of the Gemini Planet Imager (GPI), an instrument designed to achieve the high contrast at small angular separations necessary to image substellar and planetary-mass companions around nearby, young stars. We identified several issues with the GPI Data Reduction Pipeline (DRP) that significantly affected the determination of angle of north in redu… ▽ More We present a revision to the astrometric calibration of the Gemini Planet Imager (GPI), an instrument designed to achieve the high contrast at small angular separations necessary to image substellar and planetary-mass companions around nearby, young stars. We identified several issues with the GPI Data Reduction Pipeline (DRP) that significantly affected the determination of angle of north in reduced GPI images. As well as introducing a small error in position angle measurements for targets observed at small zenith distances, this error led to a significant error in the previous astrometric calibration that has affected all subsequent astrometric measurements. We present a detailed description of these issues, and how they were corrected. We reduced GPI observations of calibration binaries taken periodically since the instrument was commissioned in 2014 using an updated version of the DRP. These measurements were compared to observations obtained with the NIRC2 instrument on Keck II, an instrument with an excellent astrometric calibration, allowing us to derive an updated plate scale and north offset angle for GPI. This revised astrometric calibration should be used to calibrate all measurements obtained with GPI for the purposes of precision astrometry. △ Less

Submitted 8 April, 2020; v1 submitted 18 October, 2019; originally announced October 2019.

Comments: 36 pages, 17 figures, accepted for publication in the Journal of Astronomical Telescopes, Instruments, and Systems. Updated version includes revisions made during the referee process

Asymmetries in adaptive optics point spread functions

Authors: Alexander Madurowicz, Bruce Macintosh, Vanessa P. Bailey, Jeffrey Chilcote, Marshall Perrin, Lisa Poyneer, Laurent Pueyo, Jean-Baptiste Ruffio, Travis Barman, Joanna Bulger, Tara Cotten, Robert J. De Rosa, Rene Doyon, Gaspard Duchêne, Thomas M. Esposito, Michael P. Fitzgerald, Katherine B. Follette, Benjamin L. Gerard, Stephen J. Goodsell, James R. Graham, Alexandra Z. Greenbaum, Pascale Hibon, Li-Wei Hung, Patrick Ingraham, Paul Kalas , et al. (23 additional authors not shown)

Abstract: An explanation for the origin of asymmetry along the preferential axis of the PSF of an AO system is developed. When phase errors from high altitude turbulence scintillate due to Fresnel propagation, wavefront amplitude errors may be spatially offset from residual phase errors. These correlated errors appear as asymmetry in the image plane under the Fraunhofer condition. In an analytic model with… ▽ More An explanation for the origin of asymmetry along the preferential axis of the PSF of an AO system is developed. When phase errors from high altitude turbulence scintillate due to Fresnel propagation, wavefront amplitude errors may be spatially offset from residual phase errors. These correlated errors appear as asymmetry in the image plane under the Fraunhofer condition. In an analytic model with an open-loop AO system, the strength of the asymmetry is calculated for a single mode of phase aberration, which generalizes to two dimensions under a Fourier decomposition of the complex illumination. Other parameters included are the spatial offset of the AO correction, which is the wind velocity in the frozen flow regime multiplied by the effective AO time delay, and propagation distance or altitude of the turbulent layer. In this model, the asymmetry is strongest when the wind is slow and nearest to the coronagraphic mask when the turbulent layer is far away, such as when the telescope is pointing low towards the horizon. A great emphasis is made about the fact that the brighter asymmetric lobe of the PSF points in the opposite direction as the wind, which is consistent analytically with the clarification that the image plane electric field distribution is actually the inverse Fourier transform of the aperture plane. Validation of this understanding is made with observations taken from the Gemini Planet Imager, as well as being reproducible in end-to-end AO simulations. △ Less

Submitted 27 September, 2019; originally announced September 2019.

Comments: 28 Pages, 13 Figures, Accepted to JATIS

The Snowball Stratosphere

Authors: R. J. Graham, Tiffany Shaw, Dorian Abbot

Abstract: According to the Snowball Earth hypothesis, Earth has experienced periods of low-latitude glaciation in its deep past. Prior studies have used general circulation models (GCMs) to examine the effects such an extreme climate state might have on the structure and dynamics of Earth's troposphere, but the behavior of the stratosphere has not been studied in detail. Understanding the snowball stratosph… ▽ More According to the Snowball Earth hypothesis, Earth has experienced periods of low-latitude glaciation in its deep past. Prior studies have used general circulation models (GCMs) to examine the effects such an extreme climate state might have on the structure and dynamics of Earth's troposphere, but the behavior of the stratosphere has not been studied in detail. Understanding the snowball stratosphere is important for developing an accurate account of the Earth's radiative and chemical properties during these episodes. Here we conduct the first analysis of the stratospheric circulation of the Snowball Earth using ECHAM6 general circulation model simulations. In order to understand the factors contributing to the stratospheric circulation, we extend the Statistical Transformed Eulerian Mean framework. We find that the stratosphere during a snowball with prescribed modern ozone levels exhibits a weaker meridional overturning circulation, reduced wave activity, stronger zonal jets, and is extremely cold relative to modern conditions. Notably, the snowball stratosphere displays no sudden stratospheric warmings. Without ozone, the stratosphere displays slightly weaker circulation, a complete lack of polar vortex, and even colder temperatures. We also explicitly quantify for the first time the cross-tropopause mass exchange rate and stratospheric mixing efficiency during the snowball and show that our values do not change the constraints on CO$_2$ inferred from geochemical proxies during the Marinoan glaciation ($\sim$635 Ma), unless the O$_2$ concentration during the snowball was orders of magnitude less than the CO$_2$ concentration. △ Less

Submitted 27 September, 2019; originally announced September 2019.

Comments: 28 pages, 11 figures, submitted to JGR-Atmospheres

An Exo-Kuiper Belt and An Extended Halo around HD 191089 in Scattered Light

Authors: Bin Ren, Élodie Choquet, Marshall D. Perrin, Gaspard Duchêne, John H. Debes, Laurent Pueyo, Malena Rice, Christine Chen, Glenn Schneider, Thomas M. Esposito, Charles A. Poteet, Jason J. Wang, S. Mark Ammons, Megan Ansdell, Pauline Arriaga, Vanessa P. Bailey, Travis Barman, Juan Sebastián Bruzzone, Joanna Bulger, Jeffrey Chilcote, Tara Cotten, Robert J. De Rosa, Rene Doyon, Michael P. Fitzgerald, Katherine B. Follette , et al. (48 additional authors not shown)

Abstract: We have obtained Hubble Space Telescope STIS and NICMOS, and Gemini/GPI scattered light images of the HD 191089 debris disk. We identify two spatial components: a ring resembling Kuiper Belt in radial extent (FWHM: ${\sim}$25 au, centered at ${\sim}$46 au), and a halo extending to ${\sim}$640 au. We find that the halo is significantly bluer than the ring, consistent with the scenario that the ring… ▽ More We have obtained Hubble Space Telescope STIS and NICMOS, and Gemini/GPI scattered light images of the HD 191089 debris disk. We identify two spatial components: a ring resembling Kuiper Belt in radial extent (FWHM: ${\sim}$25 au, centered at ${\sim}$46 au), and a halo extending to ${\sim}$640 au. We find that the halo is significantly bluer than the ring, consistent with the scenario that the ring serves as the "birth ring" for the smaller dust in the halo. We measure the scattering phase functions in the 30°-150° scattering angle range and find the halo dust is both more forward- and backward-scattering than the ring dust. We measure a surface density power law index of -0.68${\pm}$0.04 for the halo, which indicates the slow-down of the radial outward motion of the dust. Using radiative transfer modeling, we attempt to simultaneously reproduce the (visible) total and (near-infrared) polarized intensity images of the birth ring. Our modeling leads to mutually inconsistent results, indicating that more complex models, such as the inclusion of more realistic aggregate particles, are needed. △ Less

Submitted 31 July, 2019; originally announced August 2019.

Comments: 29 pages, 18 figures, ApJ accepted

Impulsive coronal heating from large-scale magnetic rearrangements: from IRIS to SDO/AIA

Authors: Fabio Reale, Paola Testa, Antonino Petralia, David R. Graham

Abstract: The Interface Region Imaging Spectrograph (IRIS) has observed bright spots at the transition region footpoints associated with heating in the overlying loops, as observed by coronal imagers. Some of these brightenings show significant blueshifts in the Si iv line at 1402.77 A (logT[K] = 4.9). Such blueshifts cannot be reproduced by coronal loop models assuming heating by thermal conduction only, b… ▽ More The Interface Region Imaging Spectrograph (IRIS) has observed bright spots at the transition region footpoints associated with heating in the overlying loops, as observed by coronal imagers. Some of these brightenings show significant blueshifts in the Si iv line at 1402.77 A (logT[K] = 4.9). Such blueshifts cannot be reproduced by coronal loop models assuming heating by thermal conduction only, but are consistent with electron beam heating, highlighting for the first time the possible importance of non-thermal electrons in the heating of non-flaring active regions. Here we report on the coronal counterparts of these brightenings observed in the hot channels of the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory. We show that the IRIS bright spots are the footpoints of very hot and transient coronal loops which clearly experience strong magnetic interactions and rearrangements, thus confirming the impulsive nature of the heating and providing important constraints for their physical interpretation. △ Less

Submitted 4 July, 2019; originally announced July 2019.

Comments: 13 pages, 9 figures, accepted by ApJ, movie available upon request to F. Reale

Decomposition of the forces on a body moving in an incompressible fluid

Authors: W R Graham

Abstract: In analysing fluid forces on a moving body, a natural approach is to seek a component due to viscosity and an `inviscid' remainder. It is also attractive to decompose the velocity field into irrotational and rotational parts, and apportion the force resultants accordingly. The `irrotational' resultants can then be identified as classical `added mass', but the remaining, `rotational', resultants ap… ▽ More In analysing fluid forces on a moving body, a natural approach is to seek a component due to viscosity and an `inviscid' remainder. It is also attractive to decompose the velocity field into irrotational and rotational parts, and apportion the force resultants accordingly. The `irrotational' resultants can then be identified as classical `added mass', but the remaining, `rotational', resultants appear not to be consistent with the physical interpretation of the rotational velocity field (as that arising from the fluid vorticity with the body stationary). The alternative presented here splits the inviscid resultants into components that are unquestionably due to independent aspects of the problem: `convective' and `accelerative'. The former are associated with the pressure field that would arise in an inviscid flow with (instantaneously) the same velocities as the real one, and with the body's velocity parameters --- angular and translational --- unchanging. The latter correspond to the pressure generated when the body accelerates from rest in quiescent fluid with its given rates of change of angular and translational velocity. They are reminiscent of the added-mass force resultants, but are simpler, and closer to the standard rigid-body inertia formulae, than the developed expressions for added-mass force and moment. Finally, the force resultants due to viscosity also include a contribution from pressure. Its presence is necessary in order to satisfy the equations governing the pressure field, and it has previously been recognised in the context of `excess' stagnation-point pressure. However, its existence does not yet seem to be widely appreciated. △ Less

Submitted 16 October, 2019; v1 submitted 7 May, 2019; originally announced May 2019.

Comments: Accepted for publication in J. Fluid Mech. Current version incorporates corrections made at proof stage

Journal ref: J. Fluid Mech. 881 (2019) 1097-1122

Performance of the Gemini Planet Imager Non-Redundant Mask and spectroscopy of two close-separation binaries HR 2690 and HD 142527

Authors: Alexandra Z. Greenbaum, Anthony Cheetham, Anand Sivaramakrishnan, Fredrik T. Rantakyrö, Gaspard Duchêne, Peter Tuthill, Robert J. De Rosa, Rebecca Oppenheimer, Bruce Macintosh, S. Mark Ammons, Vanessa P. Bailey, Travis Barman, Joanna Bulger, Andrew Cardwell, Jeffrey Chilcote, Tara Cotten, Rene Doyon, Michael P. Fitzgerald, Katherine B. Follette, Benjamin L. Gerard, Stephen J. Goodsell, James R. Graham, Pascale Hibon, Li-Wei Hung, Patrick Ingraham , et al. (29 additional authors not shown)

Abstract: The Gemini Planet Imager (GPI) contains a 10-hole non-redundant mask (NRM), enabling interferometric resolution in complement to its coronagraphic capabilities. The NRM operates both in spectroscopic (integral field spectrograph, henceforth IFS) and polarimetric configurations. NRM observations were taken between 2013 and 2016 to characterize its performance. Most observations were taken in spectr… ▽ More The Gemini Planet Imager (GPI) contains a 10-hole non-redundant mask (NRM), enabling interferometric resolution in complement to its coronagraphic capabilities. The NRM operates both in spectroscopic (integral field spectrograph, henceforth IFS) and polarimetric configurations. NRM observations were taken between 2013 and 2016 to characterize its performance. Most observations were taken in spectroscopic mode with the goal of obtaining precise astrometry and spectroscopy of faint companions to bright stars. We find a clear correlation between residual wavefront error measured by the AO system and the contrast sensitivity by comparing phase errors in observations of the same source, taken on different dates. We find a typical 5-$σ$ contrast sensitivity of $2-3~\times~10^{-3}$ at $\simλ/D$. We explore the accuracy of spectral extraction of secondary components of binary systems by recovering the signal from a simulated source injected into several datasets. We outline data reduction procedures unique to GPI's IFS and describe a newly public data pipeline used for the presented analyses. We demonstrate recovery of astrometry and spectroscopy of two known companions to HR 2690 and HD 142527. NRM+polarimetry observations achieve differential visibility precision of $σ\sim0.4\%$ in the best case. We discuss its limitations on Gemini-S/GPI for resolving inner regions of protoplanetary disks and prospects for future upgrades. We summarize lessons learned in observing with NRM in spectroscopic and polarimetric modes. △ Less

Submitted 18 April, 2019; originally announced April 2019.

Comments: Accepted to AJ, 22 pages, 14 figures

The Gemini Planet Imager Exoplanet Survey: Giant Planet and Brown Dwarf Demographics From 10-100 AU

Authors: Eric L. Nielsen, Robert J. De Rosa, Bruce Macintosh, Jason J. Wang, Jean-Baptiste Ruffio, Eugene Chiang, Mark S. Marley, Didier Saumon, Dmitry Savransky, S. Mark Ammons, Vanessa P. Bailey, Travis Barman, Celia Blain, Joanna Bulger, Jeffrey Chilcote, Tara Cotten, Ian Czekala, Rene Doyon, Gaspard Duchene, Thomas M. Esposito, Daniel Fabrycky, Michael P. Fitzgerald, Katherine B. Follette, Jonathan J. Fortney, Benjamin L. Gerard , et al. (40 additional authors not shown)

Abstract: We present a statistical analysis of the first 300 stars observed by the Gemini Planet Imager Exoplanet Survey (GPIES). This subsample includes six detected planets and three brown dwarfs; from these detections and our contrast curves we infer the underlying distributions of substellar companions with respect to their mass, semi-major axis, and host stellar mass. We uncover a strong correlation be… ▽ More We present a statistical analysis of the first 300 stars observed by the Gemini Planet Imager Exoplanet Survey (GPIES). This subsample includes six detected planets and three brown dwarfs; from these detections and our contrast curves we infer the underlying distributions of substellar companions with respect to their mass, semi-major axis, and host stellar mass. We uncover a strong correlation between planet occurrence rate and host star mass, with stars M $>$ 1.5 $M_\odot$ more likely to host planets with masses between 2-13 M$_{\rm Jup}$ and semi-major axes of 3-100 au at 99.92% confidence. We fit a double power-law model in planet mass (m) and semi-major axis (a) for planet populations around high-mass stars (M $>$ 1.5M$_\odot$) of the form $\frac{d^2 N}{dm da} \propto m^αa^β$, finding $α$ = -2.4 $\pm$ 0.8 and $β$ = -2.0 $\pm$ 0.5, and an integrated occurrence rate of $9^{+5}_{-4}$% between 5-13 M$_{\rm Jup}$ and 10-100 au. A significantly lower occurrence rate is obtained for brown dwarfs around all stars, with 0.8$^{+0.8}_{-0.5}$% of stars hosting a brown dwarf companion between 13-80 M$_{\rm Jup}$ and 10-100 au. Brown dwarfs also appear to be distributed differently in mass and semi-major axis compared to giant planets; whereas giant planets follow a bottom-heavy mass distribution and favor smaller semi-major axes, brown dwarfs exhibit just the opposite behaviors. Comparing to studies of short-period giant planets from the RV method, our results are consistent with a peak in occurrence of giant planets between ~1-10 au. We discuss how these trends, including the preference of giant planets for high-mass host stars, point to formation of giant planets by core/pebble accretion, and formation of brown dwarfs by gravitational instability. △ Less

Submitted 10 April, 2019; originally announced April 2019.

Comments: 52 pages, 18 figures. AJ in press

Slow Fibonacci Walks

Authors: Fan Chung, Ron Graham, Sam Spiro

Abstract: For a positive integer $n$, we study the number of steps to reach $n$ by a {\it Fibonacci walk} for some starting pair $a_1$ and $a_2$ satisfying the recurrence of $a_{k+2}=a_{k+1}+a_k$. The problem of slow Fibonacci walks, first suggested by Richard Stanley, is to determine the maximum number $s(n)$ of steps for such a Fibonacci walk ending at $n$. Stanley conjectured that for most $n$, there is… ▽ More For a positive integer $n$, we study the number of steps to reach $n$ by a {\it Fibonacci walk} for some starting pair $a_1$ and $a_2$ satisfying the recurrence of $a_{k+2}=a_{k+1}+a_k$. The problem of slow Fibonacci walks, first suggested by Richard Stanley, is to determine the maximum number $s(n)$ of steps for such a Fibonacci walk ending at $n$. Stanley conjectured that for most $n$, there is a slow Fibonacci walk reaching $n = a_s$ with the property that $a_{s+1}$ is the integer closest to $φn$ where $φ=(1+\sqrt{5})/2$. We prove that this is true for only a positive fraction of $n$. We give explicit formulas for the choice of the starting pairs and the determination of $s(n)$ by giving a characterization theorem. We also derive a number of density results concerning the distribution of down and up cases (that is, those $n$ with $a_{s+1}=\lfloor φn\rfloor$ or $\lceil φn \rceil$, respectively), as well as for more general `paradoxical' cases. △ Less

Submitted 19 March, 2019; originally announced March 2019.

MSC Class: 11B39

A Statistical Comparative Planetology Approach to Maximize the Scientific Return of Future Exoplanet Characterization Efforts

Authors: Jade H. Checlair, Dorian S. Abbot, Robert J. Webber, Y. Katherina Feng, Jacob L. Bean, Edward W. Schwieterman, Christopher C. Stark, Tyler D. Robinson, Eliza Kempton, Olivia D. N. Alcabes, Daniel Apai, Giada Arney, Nicolas Cowan, Shawn Domagal-Goldman, Chuanfei Dong, David P. Fleming, Yuka Fujii, R. J. Graham, Scott D. Guzewich, Yasuhiro Hasegawa, Benjamin P. C. Hayworth, Stephen R. Kane, Edwin S. Kite, Thaddeus D. Komacek, Ravi K. Kopparapu , et al. (15 additional authors not shown)

Abstract: Provided that sufficient resources are deployed, we can look forward to an extraordinary future in which we will characterize potentially habitable planets. Until now, we have had to base interpretations of observations on habitability hypotheses that have remained untested. To test these theories observationally, we propose a statistical comparative planetology approach to questions of planetary… ▽ More Provided that sufficient resources are deployed, we can look forward to an extraordinary future in which we will characterize potentially habitable planets. Until now, we have had to base interpretations of observations on habitability hypotheses that have remained untested. To test these theories observationally, we propose a statistical comparative planetology approach to questions of planetary habitability. The key objective of this approach will be to make quick and cheap measurements of critical planetary characteristics on a large sample of exoplanets, exploiting statistical marginalization to answer broad habitability questions. This relaxes the requirement of obtaining multiple types of data for a given planet, as it allows us to test a given hypothesis from only one type of measurement using the power of an ensemble. This approach contrasts with a "systems science" approach, where a few planets would be extensively studied with many types of measurements. A systems science approach is associated with a number of difficulties which may limit overall scientific return, including: the limited spectral coverage and noise of instruments, the diversity of exoplanets, and the extensive list of potential false negatives and false positives. A statistical approach could also be complementary to a systems science framework by providing context to interpret extensive measurements on planets of particular interest. We strongly recommend future missions with a focus on exoplanet characterization, and with the capability to study large numbers of planets in a homogenous way, rather than exclusively small, intense studies directed at a small sample of planets. △ Less

Submitted 12 March, 2019; originally announced March 2019.

Comments: White paper submitted in response to the solicitation of feedback for the "2020 Decadal Survey" by the National Academy of Sciences

Chern-Gauss-Bonnet formula for singular Yamabe metrics in dimension four

Authors: C. Robin Graham, Matthew J. Gursky

Abstract: We derive a formula of Chern-Gauss-Bonnet type for the Euler characteristic of a four dimensional manifold-with-boundary in terms of the geometry of the Loewner-Nirenberg singular Yamabe metric in a prescribed conformal class. The formula involves the renormalized volume and a boundary integral. It is shown that if the boundary is umbilic, then the sum of the renormalized volume and the boundary i… ▽ More We derive a formula of Chern-Gauss-Bonnet type for the Euler characteristic of a four dimensional manifold-with-boundary in terms of the geometry of the Loewner-Nirenberg singular Yamabe metric in a prescribed conformal class. The formula involves the renormalized volume and a boundary integral. It is shown that if the boundary is umbilic, then the sum of the renormalized volume and the boundary integral is a conformal invariant. Analogous results are proved for asymptotically hyperbolic metrics in dimension four for which the second elementary symmetric function of the eigenvalues of the Schouten tensor is constant. Extensions and generalizations of these results are discussed. Finally, a general result is proved identifying the infinitesimal anomaly of the renormalized volume of an asymptotically hyperbolic metric in terms of its renormalized volume coefficients, and used to outline alternate proofs of the conformal invariance of the renormalized volume plus boundary integral. △ Less

Submitted 5 February, 2019; originally announced February 2019.

Comments: 26 pages

MSC Class: 53A30 (Primary) 53C20; 53C21 (Secondary)

Dynamical Constraints on the HR 8799 Planets with GPI

Authors: Jason J. Wang, James R. Graham, Rebekah Dawson, Daniel Fabrycky, Robert J. De Rosa, Laurent Pueyo, Quinn Konopacky, Bruce Macintosh, Christian Marois, Eugene Chiang, S. Mark Ammons, Pauline Arriaga, Vanessa P. Bailey, Travis Barman, Joanna Bulger, Jeffrey Chilcote, Tara Cotten, René Doyon, Gaspard Duchêne, Thomas M. Esposito, Michael P. Fitzgerald, Katherine B. Follette, Benjamin L. Gerard, Stephen J. Goodsell, Alexandra Z. Greenbaum , et al. (30 additional authors not shown)

Abstract: The HR 8799 system uniquely harbors four young super-Jupiters whose orbits can provide insights into the system's dynamical history and constrain the masses of the planets themselves. Using the Gemini Planet Imager (GPI), we obtained down to one milliarcsecond precision on the astrometry of these planets. We assessed four-planet orbit models with different levels of constraints and found that assu… ▽ More The HR 8799 system uniquely harbors four young super-Jupiters whose orbits can provide insights into the system's dynamical history and constrain the masses of the planets themselves. Using the Gemini Planet Imager (GPI), we obtained down to one milliarcsecond precision on the astrometry of these planets. We assessed four-planet orbit models with different levels of constraints and found that assuming the planets are near 1:2:4:8 period commensurabilities, or are coplanar, does not worsen the fit. We added the prior that the planets must have been stable for the age of the system (40 Myr) by running orbit configurations from our posteriors through $N$-body simulations and varying the masses of the planets. We found that only assuming the planets are both coplanar and near 1:2:4:8 period commensurabilities produces dynamically stable orbits in large quantities. Our posterior of stable coplanar orbits tightly constrains the planets' orbits, and we discuss implications for the outermost planet b shaping the debris disk. A four-planet resonance lock is not necessary for stability up to now. However, planet pairs d and e, and c and d, are each likely locked in two-body resonances for stability if their component masses are above $6~M_{\rm{Jup}}$ and $7~M_{\rm{Jup}}$, respectively. Combining the dynamical and luminosity constraints on the masses using hot-start evolutionary models and a system age of $42 \pm 5$~Myr, we found the mass of planet b to be $5.8 \pm 0.5~M_{\rm{Jup}}$, and the masses of planets c, d, and e to be $7.2_{-0.7}^{+0.6}~M_{\rm{Jup}}$ each. △ Less

Submitted 11 September, 2018; originally announced September 2018.

Comments: 21 pages, 10 figures, accepted to AJ

The Gemini Planet Imager: Looking back over five years and forward to the future

Authors: Bruce Macintosh, Jeffrey K. Chilcote, Vanessa P. Bailey, Robert De Rosa, Eric Nielsen, Andrew Norton, Lisa Poyneer, Jason Wang, J. B. Ruffio, J. R. Graham, Christian Marois, Dmitry Savransky, Jean-Pierre Veran

Abstract: The Gemini Planet Imager (GPI), a coronagraphic adaptive optics instrument designed for spectroscopy of extrasolar planets, had first light in 2013[13]. After five years, GPI has observed more than 500 stars, producing an extensive library of science images and associated telemetry that can be analyzed to determine performance predictors. We will present a summary of on-sky performance and lessons… ▽ More The Gemini Planet Imager (GPI), a coronagraphic adaptive optics instrument designed for spectroscopy of extrasolar planets, had first light in 2013[13]. After five years, GPI has observed more than 500 stars, producing an extensive library of science images and associated telemetry that can be analyzed to determine performance predictors. We will present a summary of on-sky performance and lessons learned. The two most significant factors determining bright star contrast performance are atmospheric coherence time and the presence of dome seeing. With a possible move to Gemini North, we are planning potential upgrades including a pyramid-sensor based AO system with predictive control; we will summarize upgrade options and the science they would enable. △ Less

Submitted 18 July, 2018; originally announced July 2018.

Comments: To appear in Proceedings of the SPIE 10703 (Adaptive Optics Systems VI)