Showing 1–50 of 59 results for author: Barnes, W

Modeling Time-Variable Elemental Abundances in Coronal Loop Simulations

Authors: Jeffrey W. Reep, John Unverferth, Will T. Barnes, Sherry Chhabra

Abstract: Numerous recent X-ray observations of coronal loops in both active regions (ARs) and solar flares have shown clearly that elemental abundances vary with time. Over the course of a flare, they have been found to move from coronal values towards photospheric values near the flare peak, before slowly returning to coronal values during the gradual phase. Coronal loop models typically assume that the e… ▽ More Numerous recent X-ray observations of coronal loops in both active regions (ARs) and solar flares have shown clearly that elemental abundances vary with time. Over the course of a flare, they have been found to move from coronal values towards photospheric values near the flare peak, before slowly returning to coronal values during the gradual phase. Coronal loop models typically assume that the elemental abundances are fixed, however. In this work, we introduce a time-variable abundance factor into the 0D ebtel++ code that models the changes due to chromospheric evaporation in order to understand how this affects coronal loop cooling. We find that for strong heating events ($\gtrsim$ 1 erg s$^{-1}$ cm$^{-3}$), the abundances quickly tend towards photospheric values. For smaller heating rates, the abundances fall somewhere between coronal and photospheric values, causing the loop to cool more quickly than the time-fixed photospheric cases (typical flare simulations) and more slowly than time-fixed coronal cases (typical AR simulations). This suggests heating rates in quiescent AR loops no larger than $\approx 0.1$ erg s$^{-1}$ cm$^{-3}$ to be consistent with recent measurements of abundance factors $f \gtrsim 2$. △ Less

Submitted 23 July, 2024; v1 submitted 6 June, 2024; originally announced June 2024.

Comments: Accepted to ApJL. Comments and criticisms are more than welcome! The paper is fully reproducible, with scripts to run the simulations and generate the figures available on this Github repository: https://github.com/jwreep/ebtel_abundances

CHIANTI -- an atomic database for emission lines -- Paper XVIII. Version 11, advanced ionization equilibrium models: density and charge transfer effects

Authors: R. P. Dufresne, G. Del Zanna, P. R. Young, K. P. Dere, E. Deliporanidou, W. T. Barnes, E. Landi

Abstract: Version 11 of the CHIANTI database and software package is presented. Advanced ionization equilibrium models have been added for low charge states of seven elements (C, N, O, Ne, Mg, Si and S), and represent a significant improvement especially when modelling the solar transition region. The models include the effects of higher electron density and charge transfer on ionization and recombination r… ▽ More Version 11 of the CHIANTI database and software package is presented. Advanced ionization equilibrium models have been added for low charge states of seven elements (C, N, O, Ne, Mg, Si and S), and represent a significant improvement especially when modelling the solar transition region. The models include the effects of higher electron density and charge transfer on ionization and recombination rates. As an illustration of the difference these models make, a synthetic spectrum is calculated for an electron pressure of 7$\times 10^{15}$ cm$^{-3}$ K and compared with an active region observation from HRTS. Increases are seen of factors of two to five in the predicted radiances of the strongest lines in the UV from Si IV, C IV, and N V, compared to the previous modelling using the coronal approximation. Much better agreement (within 20\%) with the observation is found for the majority of the lines. The new atomic models better equip both those who are studying the transition region and those who are interpreting emission from higher density astrophysical and laboratory plasma. In addition to the advanced models, several ion datasets have been added or updated, and data for the radiative recombination energy loss rate have been updated. △ Less

Submitted 25 March, 2024; originally announced March 2024.

Comments: 19 pages, 5 figures, 1 table, submitted to ApJ

Determining the nanoflare heating frequency of an X-ray Bright Point observed by MaGIXS

Authors: Biswajit Mondal, P. S. Athiray, Amy R. Winebarger, Sabrina L. Savage, Ken Kobayashi, Stephen Bradshaw, Will Barnes, Patrick R. Champey, Peter Cheimets, Jaroslav Dudik, Leon Golub, Helen E. Mason, David E. McKenzie, Christopher S. Moore, Chad Madsen, Katharine K. Reeves, Paola Testa, Genevieve D. Vigil, Harry P. Warren, Robert W. Walsh, Giulio Del Zanna

Abstract: Nanoflares are thought to be one of the prime candidates that can heat the solar corona to its multi-million kelvin temperature. Individual nanoflares are difficult to detect with the present generation instruments, however their presence can be inferred by comparing simulated nanoflare-heated plasma emissions with the observed emission. Using HYDRAD coronal loop simulations, we model the emission… ▽ More Nanoflares are thought to be one of the prime candidates that can heat the solar corona to its multi-million kelvin temperature. Individual nanoflares are difficult to detect with the present generation instruments, however their presence can be inferred by comparing simulated nanoflare-heated plasma emissions with the observed emission. Using HYDRAD coronal loop simulations, we model the emission from an X-ray bright point (XBP) observed by the Marshall Grazing Incidence X-ray Spectrometer (MaGIXS), along with nearest-available observations from the Atmospheric Imaging Assembly (AIA) onboard Solar Dynamics Observatory (SDO) and X-Ray Telescope (XRT) onboard Hinode observatory. The length and magnetic field strength of the coronal loops are derived from the linear-force-free extrapolation of the observed photospheric magnetogram by Helioseismic and Magnetic Imager (HMI) onboard SDO. Each loop is assumed to be heated by random nanoflares, whose magnitude and frequency are determined by the loop length and magnetic field strength. The simulation results are then compared and matched against the measured intensity from AIA, XRT, and MaGIXS. Our model results indicate the observed emissions from the XBP under study could be well matched by a distribution of nanoflares with average delay times 1500 s to 3000 s, which suggest that the heating is dominated by high-frequency events. Further, we demonstrate the high sensitivity of MaGIXS and XRT to diagnose the heating frequency using this method, while AIA passbands are found to be the least sensitive. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Comments: Accepted for publication in the Astrophysical Journal (ApJ)

Thermal Evolution of an Active Region through Quiet and Flaring Phases as Observed by NuSTAR XRT, and AIA

Authors: Jessie Duncan, Reed B. Masek, Albert Y. Shih, Lindsay Glesener, Will Barnes, Katharine K. Reeves, Yixian Zhang, Iain G. Hannah, Brian W. Grefenstette

Abstract: Solar active regions contain a broad range of temperatures, with the thermal plasma distribution often observed to peak in the few millions of kelvin. Differential emission measure (DEM) analysis can allow instruments with diverse temperature responses to be used in concert to estimate this distribution. NuSTAR HXR observations are uniquely sensitive to the highest-temperature components of the co… ▽ More Solar active regions contain a broad range of temperatures, with the thermal plasma distribution often observed to peak in the few millions of kelvin. Differential emission measure (DEM) analysis can allow instruments with diverse temperature responses to be used in concert to estimate this distribution. NuSTAR HXR observations are uniquely sensitive to the highest-temperature components of the corona, and thus extremely powerful for examining signatures of reconnection-driven heating. Here, we use NuSTAR diagnostics in combination with EUV and SXR observations (from SDO/AIA and Hinode/XRT) to construct DEMs over 170 distinct time intervals during a five-hour observation of an alternately flaring and quiet active region (NOAA designation AR 12712). This represents the first HXR study to examine the time evolution of the distribution of thermal plasma in an active region. During microflares, we find that the initial microflare-associated plasma heating is dominantly heating of material that is already relatively hot, followed later on by broader heating of initially-cooler material. During quiescent times, we show that the amount of extremely hot (>10 MK) material in this region is significantly (~3 orders of magnitude) less than that found in the quiescent active region observed in HXRs by FOXSI-2 (Ishikawa et al. 2017). This result implies there can be radically different high-temperature thermal distributions in different active regions, and strongly motivates future HXR DEM studies covering a large number of these regions. △ Less

Submitted 8 December, 2023; originally announced December 2023.

A large topographic feature on the surface of the trans-Neptunian object (307261) 2002 MS$_4$ measured from stellar occultations

Authors: F. L. Rommel, F. Braga-Ribas, J. L. Ortiz, B. Sicardy, P. Santos-Sanz, J. Desmars, J. I. B. Camargo, R. Vieira-Martins, M. Assafin, B. E. Morgado, R. C. Boufleur, G. Benedetti-Rossi, A. R. Gomes-Júnior, E. Fernández-Valenzuela, B. J. Holler, D. Souami, R. Duffard, G. Margoti, M. Vara-Lubiano, J. Lecacheux, J. L. Plouvier, N. Morales, A. Maury, J. Fabrega, P. Ceravolo , et al. (179 additional authors not shown)

Abstract: This work aims at constraining the size, shape, and geometric albedo of the dwarf planet candidate 2002 MS4 through the analysis of nine stellar occultation events. Using multichord detection, we also studied the object's topography by analyzing the obtained limb and the residuals between observed chords and the best-fitted ellipse. We predicted and organized the observational campaigns of nine st… ▽ More This work aims at constraining the size, shape, and geometric albedo of the dwarf planet candidate 2002 MS4 through the analysis of nine stellar occultation events. Using multichord detection, we also studied the object's topography by analyzing the obtained limb and the residuals between observed chords and the best-fitted ellipse. We predicted and organized the observational campaigns of nine stellar occultations by 2002 MS4 between 2019 and 2022, resulting in two single-chord events, four double-chord detections, and three events with three to up to sixty-one positive chords. Using 13 selected chords from the 8 August 2020 event, we determined the global elliptical limb of 2002 MS4. The best-fitted ellipse, combined with the object's rotational information from the literature, constrains the object's size, shape, and albedo. Additionally, we developed a new method to characterize topography features on the object's limb. The global limb has a semi-major axis of 412 $\pm$ 10 km, a semi-minor axis of 385 $\pm$ 17 km, and the position angle of the minor axis is 121 $^\circ$ $\pm$ 16$^\circ$. From this instantaneous limb, we obtained 2002 MS4's geometric albedo and the projected area-equivalent diameter. Significant deviations from the fitted ellipse in the northernmost limb are detected from multiple sites highlighting three distinct topographic features: one 11 km depth depression followed by a 25$^{+4}_{-5}$ km height elevation next to a crater-like depression with an extension of 322 $\pm$ 39 km and 45.1 $\pm$ 1.5 km deep. Our results present an object that is $\approx$138 km smaller in diameter than derived from thermal data, possibly indicating the presence of a so-far unknown satellite. However, within the error bars, the geometric albedo in the V-band agrees with the results published in the literature, even with the radiometric-derived albedo. △ Less

Submitted 23 August, 2023; v1 submitted 15 August, 2023; originally announced August 2023.

Journal ref: A&A 678, A167 (2023)

The SunPy Project: An Interoperable Ecosystem for Solar Data Analysis

Authors: The SunPy Community, Will Barnes, Steven Christe, Nabil Freij, Laura Hayes, David Stansby, Jack Ireland, Stuart Mumford, Daniel Ryan, Albert Shih

Abstract: The SunPy Project is a community of scientists and software developers creating an ecosystem of Python packages for solar physics. The project includes the sunpy core package as well as a set of affiliated packages. The sunpy core package provides general purpose tools to access data from different providers, read image and time series data, and transform between commonly used coordinate systems.… ▽ More The SunPy Project is a community of scientists and software developers creating an ecosystem of Python packages for solar physics. The project includes the sunpy core package as well as a set of affiliated packages. The sunpy core package provides general purpose tools to access data from different providers, read image and time series data, and transform between commonly used coordinate systems. Affiliated packages perform more specialized tasks that do not fall within the more general scope of the sunpy core package. In this article, we give a high-level overview of the SunPy Project, how it is broader than the sunpy core package, and how the project curates and fosters the affiliated package system. We demonstrate how components of the SunPy ecosystem, including sunpy and several affiliated packages, work together to enable multi-instrument data analysis workflows. We also describe members of the SunPy Project and how the project interacts with the wider solar physics and scientific Python communities. Finally, we discuss the future direction and priorities of the SunPy Project. △ Less

Submitted 19 April, 2023; originally announced April 2023.

Comments: 15 pages, 1 figure, published in Frontiers

Journal ref: Frontiers in Astronomy and Space Sciences 10 (2023): 29

Science Platforms for Heliophysics Data Analysis

Authors: Monica G. Bobra, Will T. Barnes, Thomas Y. Chen, Mark C. M. Cheung, Laura A. Hayes, Jack Ireland, Miho Janvier, Michael S. F. Kirk, James P. Mason, Stuart J. Mumford, Paul J. Wright

Abstract: We recommend that NASA maintain and fund science platforms that enable interactive and scalable data analysis in order to maximize the scientific return of data collected from space-based instruments. We recommend that NASA maintain and fund science platforms that enable interactive and scalable data analysis in order to maximize the scientific return of data collected from space-based instruments. △ Less

Submitted 2 January, 2023; originally announced January 2023.

Comments: Heliophysics 2050 White Paper

Jupiter and Saturn as Spectral Analogs for Extrasolar Gas Giants and Brown Dwarfs

Authors: Daniel J. Coulter, Jason W. Barnes, Jonathan J. Fortney

Abstract: With the advent of direct imaging spectroscopy, the number of spectra from brown dwarfs and extrasolar gas giants is growing rapidly. Many brown dwarfs and extrasolar gas giants exhibit spectroscopic and photometric variability, which is likely the result of weather patterns. However, for the foreseeable future, point-source observations will be the only viable method to extract brown dwarf and ex… ▽ More With the advent of direct imaging spectroscopy, the number of spectra from brown dwarfs and extrasolar gas giants is growing rapidly. Many brown dwarfs and extrasolar gas giants exhibit spectroscopic and photometric variability, which is likely the result of weather patterns. However, for the foreseeable future, point-source observations will be the only viable method to extract brown dwarf and exoplanet spectra. Models have been able to reproduce the observed variability, but ground truth observations are required to verify their results. To that end, we provide visual and near-infrared spectra of Jupiter and Saturn obtained from the \emph{Cassini} VIMS instrument. We disk-integrate the VIMS spectral cubes to simulate the spectra of Jupiter and Saturn as if they were directly imaged exoplanets or brown dwarfs. We present six empirical disk-integrated spectra for both Jupiter and Saturn with phase coverage of $1.7^\circ$ to $133.5^\circ$ and $39.6^\circ$ to $110.2^\circ$, respectively. To understand the constituents of these disk-integrated spectra, we also provide end member (single feature) spectra for permutations of illumination and cloud density, as well as for Saturn's rings. In tandem, these disk-integrated and end member spectra provide the ground truth needed to analyze point source spectra from extrasolar gas giants and brown dwarfs. Lastly, we discuss the impact that icy rings, such as Saturn's, have on disk-integrated spectra and consider the feasibility of inferring the presence of rings from direct imaging spectra. △ Less

Submitted 10 August, 2022; originally announced August 2022.

Geometric Assumptions in Hydrodynamic Modeling of Coronal and Flaring Loops

Authors: Jeffrey W. Reep, Ignacio Ugarte-Urra, Harry P. Warren, Will T. Barnes

Abstract: In coronal loop modeling, it is commonly assumed that the loops are semi-circular with a uniform cross-sectional area. However, observed loops are rarely semi-circular, and extrapolations of the magnetic field show that the field strength decreases with height, implying that the cross-sectional area expands with height. We examine these two assumptions directly to understand how they affect the hy… ▽ More In coronal loop modeling, it is commonly assumed that the loops are semi-circular with a uniform cross-sectional area. However, observed loops are rarely semi-circular, and extrapolations of the magnetic field show that the field strength decreases with height, implying that the cross-sectional area expands with height. We examine these two assumptions directly to understand how they affect the hydrodynamic and radiative response of short, hot loops to strong, impulsive electron beam heating events. Both the magnitude and rate of area expansion impact the dynamics directly, and an expanding cross-section significantly lengthens the time for a loop to cool and drain, increases upflow durations, and suppresses sound waves. The standard $T \sim n^{2}$ relation for radiative cooling does not hold with expanding loops, which cool with relatively little draining. An increase in the eccentricity of loops, on the other hand, only increases the draining timescale, and is a minor effect in general. Spectral line intensities are also strongly impacted by the variation in the cross-sectional area since they depend on both the volume of the emitting region as well as the density and ionization state. With a larger expansion, the density is reduced, so the lines at all heights are relatively reduced in intensity and, because of the increase of cooling times, the hottest lines remain bright for significantly longer. Area expansion is critical to accurate modeling of the hydrodynamics and radiation, and observations are needed to constrain the magnitude, rate, and location of the expansion or lack thereof. △ Less

Submitted 25 May, 2022; v1 submitted 8 March, 2022; originally announced March 2022.

Comments: Accepted to ApJ

Static and dynamic solar coronal loops with cross-sectional area variations

Authors: P. J. Cargill, S. J. Bradshaw, J. A. Klimchuk, W. T. Barnes

Abstract: The Enthalpy Based Thermal Evolution of Loops (EBTEL) approximate model for static and dynamic coronal loops is developed to include the effect of a loop cross-sectional area which increases from the base of the transition region (TR) to the corona. The TR is defined as the part of a loop between the top of the chromosphere and the location where thermal conduction changes from an energy loss to a… ▽ More The Enthalpy Based Thermal Evolution of Loops (EBTEL) approximate model for static and dynamic coronal loops is developed to include the effect of a loop cross-sectional area which increases from the base of the transition region (TR) to the corona. The TR is defined as the part of a loop between the top of the chromosphere and the location where thermal conduction changes from an energy loss to an energy gain. There are significant differences from constant area loops due to the manner in which the reduced volume of the TR responds to conductive and enthalpy fluxes from the corona. For static loops with modest area variation the standard picture of loop energy balance is retained, with the corona and TR being primarily a balance between heating and conductive losses in the corona, and downward conduction and radiation to space in the TR. As the area at the loop apex increases, the TR becomes thicker and the density in TR and corona larger. For large apex areas, the coronal energy balance changes to one primarily between heating and radiation, with conduction playing an increasingly unimportant role, and the TR thickness becoming a significant fraction of the loop length. Approximate scaling laws are derived that give agreement with full numerical solutions for the density, but not the temperature. For non-uniform areas, dynamic loops have a higher peak temperature and are denser in the radiative cooling phase by of order 50% than the constant area case for the examples considered. They also show a final rapid cooling and draining once the temperature approaches 1 MK. Although the magnitude of the emission measure will be enhanced in the radiative phase, there is little change in the important observational diagnostic of its temperature dependence. △ Less

Submitted 17 November, 2021; originally announced November 2021.

Comments: 10 pages, 9 figures, published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, 2021;, stab3163

Science goals and new mission concepts for future exploration of Titan's atmosphere geology and habitability: Titan POlar Scout/orbitEr and In situ lake lander and DrONe explorer (POSEIDON)

Authors: Sébastien Rodriguez, Sandrine Vinatier, Daniel Cordier, Gabriel Tobie, Richard K. Achterberg, Carrie M. Anderson, Sarah V. Badman, Jason W. Barnes, Erika L. Barth, Bruno Bézard, Nathalie Carrasco, Benjamin Charnay, Roger N. Clark, Patrice Coll, Thomas Cornet, Athena Coustenis, Isabelle Couturier-Tamburelli, Michel Dobrijevic, F. Michael Flasar, Remco de Kok, Caroline Freissinet, Marina Galand, Thomas Gautier, Wolf D. Geppert, Caitlin A. Griffith , et al. (39 additional authors not shown)

Abstract: In response to ESA Voyage 2050 announcement of opportunity, we propose an ambitious L-class mission to explore one of the most exciting bodies in the Solar System, Saturn largest moon Titan. Titan, a "world with two oceans", is an organic-rich body with interior-surface-atmosphere interactions that are comparable in complexity to the Earth. Titan is also one of the few places in the Solar System w… ▽ More In response to ESA Voyage 2050 announcement of opportunity, we propose an ambitious L-class mission to explore one of the most exciting bodies in the Solar System, Saturn largest moon Titan. Titan, a "world with two oceans", is an organic-rich body with interior-surface-atmosphere interactions that are comparable in complexity to the Earth. Titan is also one of the few places in the Solar System with habitability potential. Titan remarkable nature was only partly revealed by the Cassini-Huygens mission and still holds mysteries requiring a complete exploration using a variety of vehicles and instruments. The proposed mission concept POSEIDON (Titan POlar Scout/orbitEr and In situ lake lander DrONe explorer) would perform joint orbital and in situ investigations of Titan. It is designed to build on and exceed the scope and scientific/technological accomplishments of Cassini-Huygens, exploring Titan in ways that were not previously possible, in particular through full close-up and in situ coverage over long periods of time. In the proposed mission architecture, POSEIDON consists of two major elements: a spacecraft with a large set of instruments that would orbit Titan, preferably in a low-eccentricity polar orbit, and a suite of in situ investigation components, i.e. a lake lander, a "heavy" drone (possibly amphibious) and/or a fleet of mini-drones, dedicated to the exploration of the polar regions. The ideal arrival time at Titan would be slightly before the next northern Spring equinox (2039), as equinoxes are the most active periods to monitor still largely unknown atmospheric and surface seasonal changes. The exploration of Titan northern latitudes with an orbiter and in situ element(s) would be highly complementary with the upcoming NASA New Frontiers Dragonfly mission that will provide in situ exploration of Titan equatorial regions in the mid-2030s. △ Less

Submitted 20 October, 2021; originally announced October 2021.

Comments: arXiv admin note: substantial text overlap with arXiv:1908.01374

Exploring tidal obliquity variations with SMERCURY-T

Authors: Steven M. Kreyche, Jason W. Barnes, Billy L. Quarles, John E. Chambers

Abstract: We introduce our new code, SMERCURY-T, which is based on existing codes SMERCURY (Lissauer et al. 2012) and Mercury-T (Bolmont et al. 2015). The result is a mixed-variable symplectic N-body integrator that can compute the orbital and spin evolution of a planet within a multi-planet system under the influence of tidal spin torques from its star. We validate our implementation by comparing our exper… ▽ More We introduce our new code, SMERCURY-T, which is based on existing codes SMERCURY (Lissauer et al. 2012) and Mercury-T (Bolmont et al. 2015). The result is a mixed-variable symplectic N-body integrator that can compute the orbital and spin evolution of a planet within a multi-planet system under the influence of tidal spin torques from its star. We validate our implementation by comparing our experimental results to that of a secular model. As we demonstrate in a series of experiments, SMERCURY-T allows for the study of secular spin-orbit resonance crossings and captures for planets within complex multi-planet systems. These processes can drive a planet's spin state to evolve along vastly different pathways on its road toward tidal equilibrium, as tidal spin torques dampen the planet's spin rate and evolve its obliquity. Additionally, we show the results of a scenario that exemplifies the crossing of a chaotic region that exists as the overlap of two spin-orbit resonances. The test planet experiences violent and chaotic swings in its obliquity until its eventual escape from resonance as it tidally evolves. All of these processes are and have been important over the obliquity evolution of many bodies within the Solar System and beyond, and have implications for planetary climate and habitability. SMERCURY-T is a powerful and versatile tool that allows for further study of these phenomena. △ Less

Submitted 7 September, 2021; originally announced September 2021.

Comments: 21 pages, 9 figures, 3 tables; accepted by The Planetary Science Journal

Understanding Heating in Active Region Cores through Machine Learning II. Classifying Observations

Authors: W. T. Barnes, S. J. Bradshaw, N. M. Viall

Abstract: Constraining the frequency of energy deposition in magnetically-closed active region cores requires sophisticated hydrodynamic simulations of the coronal plasma and detailed forward modeling of the optically-thin line-of-sight integrated emission. However, understanding which set of model inputs best matches a set of observations is complicated by the need for any proposed heating model to simulta… ▽ More Constraining the frequency of energy deposition in magnetically-closed active region cores requires sophisticated hydrodynamic simulations of the coronal plasma and detailed forward modeling of the optically-thin line-of-sight integrated emission. However, understanding which set of model inputs best matches a set of observations is complicated by the need for any proposed heating model to simultaneously satisfy multiple observable constraints. In this paper, we train a random forest classification model on a set of forward-modeled observable quantities, namely the emission measure slope, the peak temperature of the emission measure distribution, and the time lag and maximum cross-correlation between multiple pairs of AIA channels. We then use our trained model to classify the heating frequency in every pixel of active region NOAA 1158 using the observed emission measure slopes, peak temperatures, time lags, and maximum cross-correlations and are able to map the heating frequency across the entire active region. We find that high-frequency heating dominates in the inner core of the active region while intermediate frequency dominates closer to the periphery of the active region. Additionally, we assess the importance of each observed quantity in our trained classification model and find that the emission measure slope is the dominant feature in deciding with which heating frequency a given pixel is most consistent. The technique presented here offers a very promising and widely applicable method for assessing observations in terms of detailed forward models given an arbitrary number of observable constraints. △ Less

Submitted 15 July, 2021; originally announced July 2021.

Comments: 29 pages, 8 figures, accepted for publication in The Astrophysical Journal

Forecasting the Remaining Duration of an Ongoing Solar Flare

Authors: Jeffrey W. Reep, Will T. Barnes

Abstract: The solar X-ray irradiance is significantly heightened during the course of a solar flare, which can cause radio blackouts due to ionization of the atoms in the ionosphere. As the duration of a solar flare is not related to the size of that flare, it is not directly clear how long those blackouts can persist. Using a random forest regression model trained on data taken from X-ray light curves, we… ▽ More The solar X-ray irradiance is significantly heightened during the course of a solar flare, which can cause radio blackouts due to ionization of the atoms in the ionosphere. As the duration of a solar flare is not related to the size of that flare, it is not directly clear how long those blackouts can persist. Using a random forest regression model trained on data taken from X-ray light curves, we have developed a direct forecasting method that predicts how long the event will remain above background levels. We test this on a large collection of flares observed with GOES-15, and show that it generally outperforms simple linear regression. This forecast is computationally light enough to be performed in real time, allowing for the prediction to be made during the course of a flare. △ Less

Submitted 1 June, 2021; v1 submitted 5 March, 2021; originally announced March 2021.

Comments: Resubmitted to Space Weather. Comments welcome!

Titan: Earth-like on the Outside, Ocean World on the Inside

Authors: Shannon M. MacKenzie, Samuel P. D. Birch, Sarah Horst, Christophe Sotin, Erika Barth, Juan M. Lora, Melissa G. Trainer, Paul Corlies, Michael J. Malaska, Ella Sciamma-O'Brien, Alexander E. Thelen, Elizabeth P. Turtle, Jani Radebaugh, Jennifer Hanley, Anezina Solomonidou, Claire Newman, Leonardo Regoli, Sebastien Rodriguez, Benoit Seignovert, Alexander G. Hayes, Baptiste Journaux, Jordan Steckloff, Delphine Nna-Mvondo, Thomas Cornet, Maureen Palmer , et al. (8 additional authors not shown)

Abstract: Thanks to the Cassini-Huygens mission, Titan, the pale orange dot of Pioneer and Voyager encounters has been revealed to be a dynamic, hydrologically-shaped, organic-rich ocean world offering unparalleled opportunities to explore prebiotic chemistry. And while Cassini-Huygens revolutionized our understanding of each of the three layers of Titan--the atmosphere, the surface, and the interior--we ar… ▽ More Thanks to the Cassini-Huygens mission, Titan, the pale orange dot of Pioneer and Voyager encounters has been revealed to be a dynamic, hydrologically-shaped, organic-rich ocean world offering unparalleled opportunities to explore prebiotic chemistry. And while Cassini-Huygens revolutionized our understanding of each of the three layers of Titan--the atmosphere, the surface, and the interior--we are only beginning to hypothesize how these realms interact. In this paper, we summarize the current state of Titan knowledge and discuss how future exploration of Titan would address some of the next decade's most compelling planetary science questions. We also demonstrate why exploring Titan, both with and beyond the Dragonfly New Frontiers mission, is a necessary and complementary component of an Ocean Worlds Program that seeks to understand whether habitable environments exist elsewhere in our solar system. △ Less

Submitted 16 February, 2021; originally announced February 2021.

Comments: Submitted to the PSJ Focus Issue on Ocean World Exploration

Solving the Alhazen-Ptolemy Problem: Determining Specular Points on Spherical Surfaces for Radiative Transfer of Titan's Seas

Authors: William J. Miller, Jason W. Barnes, Shannon M. MacKenzie

Abstract: Given a light source, a spherical reflector, and an observer, where on the surface of the sphere will the light be directly reflected to the observer, i.e. where is the the specular point? This is known as the Alhazen-Ptolemy problem, and finding this specular point for spherical reflectors is useful in applications ranging from computer rendering to atmospheric modeling to GPS communications. Exi… ▽ More Given a light source, a spherical reflector, and an observer, where on the surface of the sphere will the light be directly reflected to the observer, i.e. where is the the specular point? This is known as the Alhazen-Ptolemy problem, and finding this specular point for spherical reflectors is useful in applications ranging from computer rendering to atmospheric modeling to GPS communications. Existing solutions rely upon finding the roots of a quartic equation and evaluating numerically which root provides the real specular point. We offer a formulation, and two solutions thereof, for which the correct root is predeterminable, thereby allowing the construction of the fully analytical solutions we present. Being faster to compute, our solutions should prove useful in cases which require repeated calculation of the specular point, such as Monte-Carlo radiative transfer, including reflections off of Titan's hydrocarbon seas. △ Less

Submitted 3 December, 2020; originally announced December 2020.

Comments: 11 pages with 5 figures. Submitted to Planetary Science Journal

Nanoflare Diagnostics from Magnetohydrodynamic Heating Profiles

Authors: Kalman J. Knizhnik, Will T. Barnes, Jeffrey W. Reep, Vadim M. Uritsky

Abstract: The nanoflare paradigm of coronal heating has proven extremely promising for explaining the presence of hot, multi-million degree loops in the solar corona. In this paradigm, localized heating events supply enough energy to heat the solar atmosphere to its observed temperatures. Rigorously modeling this process, however, has proven difficult, since it requires an accurate treatment of both the mag… ▽ More The nanoflare paradigm of coronal heating has proven extremely promising for explaining the presence of hot, multi-million degree loops in the solar corona. In this paradigm, localized heating events supply enough energy to heat the solar atmosphere to its observed temperatures. Rigorously modeling this process, however, has proven difficult, since it requires an accurate treatment of both the magnetic field dynamics and reconnection as well as the plasma's response to magnetic perturbations. In this paper, we combine fully 3D magnetohydrodynamic (MHD) simulations of coronal active region plasma driven by photospheric motions with spatially-averaged, time-dependent hydrodynamic (HD) modeling of coronal loops to obtain physically motivated observables that can be quantitatively compared with observational measurements of active region cores. We take the behavior of reconnected field lines from the MHD simulation and use them to populate the HD model to obtain the thermodynamic evolution of the plasma and subsequently the emission measure distribution. We find the that the photospheric driving of the MHD model produces only very low-frequency nanoflare heating which cannot account for the full range of active region core observations as measured by the low-temperature emission measure slope. Additionally, we calculate the spatial and temporal distributions of field lines exhibiting collective behavior, and argue that loops occur due to random energization occurring on clusters of adjacent field lines. △ Less

Submitted 31 August, 2020; originally announced September 2020.

Journal ref: Astrophysical Journal, 2020, 899, 156

The Science Case for a Titan Flagship-class Orbiter with Probes

Authors: Conor A. Nixon, James Abshire, Andrew Ashton, Jason W. Barnes, Nathalie Carrasco, Mathieu Choukroun, Athena Coustenis, Louis-Alexandre Couston, Niklas Edberg, Alexander Gagnon, Jason D. Hofgartner, Luciano Iess, Stéphane Le Mouélic, Rosaly Lopes, Juan Lora, Ralph D. Lorenz, Adrienn Luspay-Kuti, Michael Malaska, Kathleen Mandt, Marco Mastrogiuseppe, Erwan Mazarico, Marc Neveu, Taylor Perron, Jani Radebaugh, Sébastien Rodriguez , et al. (14 additional authors not shown)

Abstract: We outline a flagship-class mission concept focused on studying Titan as a global system, with particular emphasis on the polar regions. Investigating Titan from the unique standpoint of a polar orbit would enable comprehensive global maps to uncover the physics and chemistry of the atmosphere, and the topography and geophysical environment of the surface and subsurface. The mission includes two k… ▽ More We outline a flagship-class mission concept focused on studying Titan as a global system, with particular emphasis on the polar regions. Investigating Titan from the unique standpoint of a polar orbit would enable comprehensive global maps to uncover the physics and chemistry of the atmosphere, and the topography and geophysical environment of the surface and subsurface. The mission includes two key elements: (1) an orbiter spacecraft, which also acts as a data relay, and (2) one or more small probes to directly investigate Titan's seas and make the first direct measurements of their liquid composition and physical environment. The orbiter would carry a sophisticated remote sensing payload, including a novel topographic lidar, a long-wavelength surface-penetrating radar, a sub-millimeter sounder for winds and for mesospheric/thermospheric composition, and a camera and near-infrared spectrometer. An instrument suite to analyze particles and fields would include a mass spectrometer to focus on the interactions between Titan's escaping upper atmosphere and the solar wind and Saturnian magnetosphere. The orbiter would enter a stable polar orbit around 1500 to 1800 km, from which vantage point it would make global maps of the atmosphere and surface. One or more probes, released from the orbiter, would investigate Titan's seas in situ, including possible differences in composition between higher and lower latitude seas, as well as the atmosphere during the parachute descent. The number of probes, as well as the instrument complement on the orbiter and probe, remain to be finalized during a mission study that we recommend to NASA as part of the NRC Decadal Survey for Planetary Science now underway, with the goal of an overall mission cost in the "small flagship" category of ~$2 bn. International partnerships, similar to Cassini-Huygens, may also be included for consideration. △ Less

Submitted 13 August, 2020; originally announced August 2020.

Comments: 13 pages, white paper submitted to the NRC Decadal Survey for Planetary Science and Astrobiology

Tidal Currents Detected in Kraken Mare Straits from Cassini VIMS Sun Glitter Observations

Authors: Michael F. Heslar, Jason W. Barnes, Jason M. Soderblom, Benoit Seignovert, Rajani D. Dhingra, Christophe Sotin

Abstract: We present Cassini VIMS observations of sun glitter -- wave-induced reflections from a liquid surface offset from a specular point -- on Kraken Mare. Sun glitter reveals rough sea surfaces around Kraken Mare, namely the coasts and narrow straits. The sun glitter observations indicate wave activity driven by the winds and tidal currents in Kraken Mare during northern summer. T104 Cassini VIMS obser… ▽ More We present Cassini VIMS observations of sun glitter -- wave-induced reflections from a liquid surface offset from a specular point -- on Kraken Mare. Sun glitter reveals rough sea surfaces around Kraken Mare, namely the coasts and narrow straits. The sun glitter observations indicate wave activity driven by the winds and tidal currents in Kraken Mare during northern summer. T104 Cassini VIMS observations show three sun glitter features in Bayta Fretum indicative of variegated wave fields. We cannot uniquely determine one source for the coastal Bayta waves, but we lean toward the interpretation of surface winds, because tidal currents should be too weak to generate capillary-gravity waves in Bayta Fretum. T105 and T110 observations reveal wave fields in the straits of Seldon Fretum, Lulworth Sinus, and Tunu Sinus that likely originate from the constriction of tidal currents. Coastlines of Bermoothes and Hufaidh Insulae adjoin rough sea surfaces, suggesting a complex interplay of wind-roughened seas and localized tidal currents. Bermoothes and Hufaidh Insulae may share characteristics of either the Torres Strait off Australia or the Aland region of Finland, summarized as an island-dense strait with shallow bathymetry that hosts complex surface circulation patterns. Hufaidh Insulae could host seafloor bedforms formed by tidal currents with an abundant sediment supply, similar to the Torres Strait. The coastlines of Hufaidh and Bermoothes Insulae likely host ria or flooded coastal inlets, suggesting the Insulae may be local peaks of primordial crust isolated by an episode of sea-level rise or tectonic uplift. △ Less

Submitted 1 July, 2020; originally announced July 2020.

Comments: Accepted to the Planetary Science Journal, 25 pages, 13 figures

KELT-9 b's Asymmetric TESS Transit Caused by Rapid Stellar Rotation and Spin-Orbit Misalignment

Authors: John P. Ahlers, Marshall C. Johnson, Keivan G Stassun, Knicole D. Colon, Jason W. Barnes, Daniel J. Stevens, Thomas Beatty, B. Scott Gaudi, Karen A. Collins, Joseph Rodriguez, George Ricker, Roland Vanderspek, David Latham, Sara Seager, Joshua Winn, Jon M. Jenkins, Douglas A. Caldwell, Robert F. Goeke, Hugh P. Osborn, Martin Paegert, Pam Rowden, Peter Tenenbaum

Abstract: KELT-9 b is an ultra hot Jupiter transiting a rapidly rotating, oblate early-A-type star in a polar orbit. We model the effect of rapid stellar rotation on KELT-9 b's transit light curve using photometry from the Transiting Exoplanet Survey Satellite (\tess) to constrain the planet's true spin-orbit angle and to explore how KELT-9 b may be influenced by stellar gravity darkening. We constrain the… ▽ More KELT-9 b is an ultra hot Jupiter transiting a rapidly rotating, oblate early-A-type star in a polar orbit. We model the effect of rapid stellar rotation on KELT-9 b's transit light curve using photometry from the Transiting Exoplanet Survey Satellite (\tess) to constrain the planet's true spin-orbit angle and to explore how KELT-9 b may be influenced by stellar gravity darkening. We constrain the host star's equatorial radius to be $1.089\pm0.017$ times as large as its polar radius and its local surface brightness to vary by $\sim38$\% between its hot poles and cooler equator. We model the stellar oblateness and surface brightness gradient and find that it causes the transit light curve to lack the usual symmetry around the time of minimum light. We take advantage of the light curve asymmetry to constrain KELT-9 b's true spin orbit angle (${87^\circ}^{+10^\circ}_{-11^\circ}$), agreeing with \citet{gaudi2017giant} that KELT-9 b is in a nearly polar orbit. We also apply a gravity darkening correction to the spectral energy distribution model from \citet{gaudi2017giant} and find that accounting for rapid rotation gives a better fit to available spectroscopy and yields a more reliable estimate for the star's polar effective temperature. △ Less

Submitted 28 April, 2020; originally announced April 2020.

Comments: Accepted for Publication in ApJ. arXiv admin note: text overlap with arXiv:1911.05025

Retrograde-rotating exoplanets experience obliquity excitations in an eccentricity-enabled resonance

Authors: Steven M. Kreyche, Jason W. Barnes, Billy L. Quarles, Jack J. Lissauer, John E. Chambers, Matthew M. Hedman

Abstract: Previous studies have shown that planets that rotate retrograde (backwards with respect to their orbital motion) generally experience less severe obliquity variations than those that rotate prograde (the same direction as their orbital motion). Here we examine retrograde-rotating planets on eccentric orbits and find a previously unknown secular spin-orbit resonance that can drive significant obliq… ▽ More Previous studies have shown that planets that rotate retrograde (backwards with respect to their orbital motion) generally experience less severe obliquity variations than those that rotate prograde (the same direction as their orbital motion). Here we examine retrograde-rotating planets on eccentric orbits and find a previously unknown secular spin-orbit resonance that can drive significant obliquity variations. This resonance occurs when the frequency of the planet's rotation axis precession becomes commensurate with an orbital eigenfrequency of the planetary system. The planet's eccentricity enables a participating orbital frequency through an interaction in which the apsidal precession of the planet's orbit causes a cyclic nutation of the planet's orbital angular momentum vector. The resulting orbital frequency follows the relationship $f = 2 \dot{\varpi} - \dotΩおめが$, where $\dot{\varpi}$ and $\dotΩおめが$ are the rates of the planet's changing longitude of periapsis and ascending node, respectively. We test this mechanism by simulating cases of a simple Earth-Jupiter system, and confirm the predicted resonance. Over the course of 100 Myr, the test Earths with rotation axis precession rates near the predicted resonant frequency experienced pronounced obliquity variations of order $10^\circ$-$30^\circ$. These variations can be significant, and suggest that while retrograde rotation is a stabilizing influence most of the time, retrograde rotators can experience large obliquity variations if they are on eccentric orbits and enter this spin-orbit resonance. △ Less

Submitted 7 April, 2020; v1 submitted 30 March, 2020; originally announced March 2020.

Comments: Accepted to The Planetary Science Journal. Figures 2 and 4b have associated video files hosted at https://github.com/SMKreyche/Kreyche-et-al.-2020-Videos

Dust Devils on Titan

Authors: Brian Jackson, Ralph D. Lorenz, Jason W. Barnes, Michelle Szurgot

Abstract: Conditions on Saturn's moon Titan suggest dust devils, which are convective, dust-laden plumes, may be active. Although the exact nature of dust on Titan is unclear, previous observations confirm an active aeolian cycle, and dust devils may play an important role in Titan's aeolian cycle, possibly contributing to regional transport of dust and even production of sand grains. The Dragonfly mission… ▽ More Conditions on Saturn's moon Titan suggest dust devils, which are convective, dust-laden plumes, may be active. Although the exact nature of dust on Titan is unclear, previous observations confirm an active aeolian cycle, and dust devils may play an important role in Titan's aeolian cycle, possibly contributing to regional transport of dust and even production of sand grains. The Dragonfly mission to Titan will document dust devil and convective vortex activity and thereby provide a new window into these features, and our analysis shows that associated winds are likely to be modest and pose no hazard to the mission. △ Less

Submitted 13 February, 2020; originally announced February 2020.

Comments: Accepted for publication in JGR: Planets; data and code available here - https://zenodo.org/record/3470280

Solar Active Region Heating Diagnostics from High Temperature Emission using the Marshall Grazing Incidence X-ray Spectrometer (MaGIXS)

Authors: P. S. Athiray, Amy R. Winebarger, Will T. Barnes, Stephen J. Bradshaw, Sabrina Savage, Harry P. Warren, Ken Kobayashi, Patrick Champey, Leon Golub, Lindsay Glesener

Abstract: The relative amount of high temperature plasma has been found to be a useful diagnostic to determine the frequency of coronal heating on sub-resolution structures. When the loops are infrequently heated, a broad emission measure (EM) over a wider range of temperatures is expected. A narrower EM is expected for high frequency heating where the loops are closer to equilibrium. The soft X-ray spectru… ▽ More The relative amount of high temperature plasma has been found to be a useful diagnostic to determine the frequency of coronal heating on sub-resolution structures. When the loops are infrequently heated, a broad emission measure (EM) over a wider range of temperatures is expected. A narrower EM is expected for high frequency heating where the loops are closer to equilibrium. The soft X-ray spectrum contains many spectral lines that provide high temperature diagnostics, including lines from Fe XVII-XIX. This region of the solar spectrum will be observed by the Marshall Grazing Incidence Spectrometer (MaGIXS) in 2020. In this paper, we derive the expected spectral lines intensity in MaGIXS to varying amounts of high temperature plasma to demonstrate that a simple line ratio of these provides a powerful diagnostic to determine the heating frequency. Similarly, we examine ratios of AIA channel intensities, filter ratios from a XRT, and energy bands from the FOXSI sounding rocket to determine their sensitivity to this parameter. We find that both FOXSI and MaGIXS provide good diagnostic capability for high-temperature plasma. We then compare the predicted line ratios to the output of a numerical model and confirm the MaGIXS ratios provide an excellent diagnostic for heating frequency. △ Less

Submitted 5 September, 2019; originally announced September 2019.

Comments: Accepted for publication in ApJ

Journal ref: The Astrophysical Journal, 2019, 884 24

Key Challenges for AAS Journals in the Next Decade

Authors: Emily M. Levesque, Lisa Prato, Christopher Sneden, Jason W. Barnes, Dawn M. Gelino, Barbara Kern, Paula Szkody, Rosemary F. G. Wyse, Leslie A. Young

Abstract: The American Astronomical Society (AAS) Journals are a vital asset of our professional society. With the push towards open access, page charges are a viable and sustainable option for continuing to effectively fund and publish the AAS Journals. However, the existing page charge model, which requires individual authors to pay page charges out of their grants or even out of pocket, is already challe… ▽ More The American Astronomical Society (AAS) Journals are a vital asset of our professional society. With the push towards open access, page charges are a viable and sustainable option for continuing to effectively fund and publish the AAS Journals. However, the existing page charge model, which requires individual authors to pay page charges out of their grants or even out of pocket, is already challenging to some researchers and could be exacerbated in the Open Access (OA) era if charges increase. A discussion of alternative models for funding page charges and publishing costs should be part of the Astro2020 decadal survey if we wish to continue supporting the sustainable and accessible publication of US research in AAS journals in the rapidly-shifting publication landscape. The AAS Publications Committee recommends that the National Academy of Sciences form a task force to develop solutions and recommendations with respect to the urgent concerns and considerations highlighted in this White Paper. △ Less

Submitted 25 July, 2019; originally announced July 2019.

Comments: 5 pages, submitted as an Astro2020 Activity, Project of State of the Profession Consideration (APC) White Paper

Understanding Heating in Active Region Cores through Machine Learning I. Numerical Modeling and Predicted Observables

Authors: W. T. Barnes, S. J. Bradshaw, N. M. Viall

Abstract: To adequately constrain the frequency of energy deposition in active region cores in the solar corona, systematic comparisons between detailed models and observational data are needed. In this paper, we describe a pipeline for forward modeling active region emission using magnetic field extrapolations and field-aligned hydrodynamic models. We use this pipeline to predict time-dependent emission fr… ▽ More To adequately constrain the frequency of energy deposition in active region cores in the solar corona, systematic comparisons between detailed models and observational data are needed. In this paper, we describe a pipeline for forward modeling active region emission using magnetic field extrapolations and field-aligned hydrodynamic models. We use this pipeline to predict time-dependent emission from active region NOAA 1158 as observed by SDO/AIA for low-, intermediate-, and high-frequency nanoflares. In each pixel of our predicted multi-wavelength, time-dependent images, we compute two commonly-used diagnostics: the emission measure slope and the time lag. We find that signatures of the heating frequency persist in both of these diagnostics. In particular, our results show that the distribution of emission measure slopes narrows and the mean decreases with decreasing heating frequency and that the range of emission measure slopes is consistent with past observational and modeling work. Furthermore, we find that the time lag becomes increasingly spatially coherent with decreasing heating frequency while the distribution of time lags across the whole active region becomes more broad with increasing heating frequency. In a follow up paper, we train a random forest classifier on these predicted diagnostics and use this model to classify real AIA observations of NOAA 1158 in terms of the underlying heating frequency. △ Less

Submitted 7 June, 2019; originally announced June 2019.

Comments: 32 pages, 12 figures, accepted for publication in The Astrophysical Journal

Dealing With $δでるた$-Scuti Variables: Transit Light Curve Analysis of Planets Orbiting Rapidly-Rotating, Seismically Active A/F Stars

Authors: John P. Ahlers, Jason W. Barnes, Samuel A. Myers

Abstract: We measure the bulk system parameters of the seismically active, rapidly-rotating $δでるた$-Scuti KOI-976 and constrain the orbit geometry of its transiting binary companion using a combined approach of asteroseismology and gravity-darkening light curve analysis. KOI-976 is a $1.62\pm0.2~\mathrm{M_\odot}$ star with a measured $v\sin(i)$ of $120\pm2$ km/s and seismically-induced variable signal that vari… ▽ More We measure the bulk system parameters of the seismically active, rapidly-rotating $δでるた$-Scuti KOI-976 and constrain the orbit geometry of its transiting binary companion using a combined approach of asteroseismology and gravity-darkening light curve analysis. KOI-976 is a $1.62\pm0.2~\mathrm{M_\odot}$ star with a measured $v\sin(i)$ of $120\pm2$ km/s and seismically-induced variable signal that varies by $\sim$ 0.6\% of the star's total photometric brightness. We take advantage of the star's oblate shape and seismic activity to perform three measurements of its obliquity angle relative to the plane of the sky. We first apply rotational splitting theory to the star's variable signal observed in short-cadence \emph{Kepler} photometry to constrain KOI-976's obliquity angle, and then subtract off variability from that dataset using the linear algorithm for significance reduction software {\tt LASR}. We perform gravity-darkened fits to \emph{Kepler} variability-subtracted short-cadence photometry and to \emph{Kepler's} phase-folded long-cadence photometry to obtain two more measurements of the star's obliquity. We find that the binary system transits in a grazing configuration with measured obliquity values of $36^\circ\pm17^\circ$, $46^\circ\pm16^\circ$, and $43^\circ\pm20^\circ$ respectively for the three measurements. We perform these analyses as a way to demonstrate overcoming the challenges high-mass stars can present to transit light curve fitting and to prepare for the large number of exoplanets \emph{TESS} will discover orbiting A/F stars. △ Less

Submitted 29 July, 2019; v1 submitted 3 June, 2019; originally announced June 2019.

Comments: Accepted for publication in AJ. Replacing earlier version to fix errata

The Cassini VIMS archive of Titan: from browse products to global infrared color maps

Authors: Stéphane Le Mouélic, Thomas Cornet, Sébastien Rodriguez, Christophe Sotin, Benoît Seignovert, Jason W. Barnes, Robert H. Brown, Kevin H. Baines, Bonnie J. Buratti, Roger N. Clark, Philip D. Nicholson, Jérémie Lasue, Virginia Pasek, Jason M. Soberblom

Abstract: We have analyzed the complete Visual and Infrared Mapping Spectrometer (VIMS) data archive of Titan. Our objective is to build global surface cartographic products, by combining all the data gathered during the 127 targeted flybys of Titan into synthetic global maps interpolated on a grid at 32 pixels per degree (~1.4 km/pixel at the equator), in seven infrared spectral atmospheric windows. Multis… ▽ More We have analyzed the complete Visual and Infrared Mapping Spectrometer (VIMS) data archive of Titan. Our objective is to build global surface cartographic products, by combining all the data gathered during the 127 targeted flybys of Titan into synthetic global maps interpolated on a grid at 32 pixels per degree (~1.4 km/pixel at the equator), in seven infrared spectral atmospheric windows. Multispectral summary images have been computed for each single VIMS cube in order to rapidly identify their scientific content and assess their quality. These summary images are made available to the community on a public website (vims.univ-nantes.fr). The global mapping work faced several challenges due to the strong absorbing and scattering effects of the atmosphere coupled to the changing observing conditions linked to the orbital tour of the Cassini mission. We determined a surface photometric function which accounts for variations in incidence, emergence and phase angles, and which is able to mitigate brightness variations linked to the viewing geometry of the flybys. The atmospheric contribution has been reduced using the subtraction of the methane absorption band wings, considered as proxies for atmospheric haze scattering. We present a new global three color composite map of band ratios (red: 1.59/1.27 μみゅーm; green: 2.03/1.27 μみゅーm; blue: 1.27/1.08 μみゅーm), which has also been empirically corrected from an airmass (the solar photon path length through the atmosphere) dependence. This map provides a detailed global color view of Titan's surface partially corrected from the atmosphere and gives a global insight of the spectral variability, with the equatorial dunes fields appearing in brownish tones, and several occurrences of bluish tones localized in areas such as Sinlap, Menvra and Selk craters. (...) △ Less

Submitted 18 September, 2018; originally announced September 2018.

Journal ref: Icarus 2018

Mapping polar atmospheric features on Titan with VIMS: from the dissipation of the northern cloud to the onset of a southern polar vortex

Authors: Stéphane Le Mouélic, Sébastien Rodriguez, Rozen Robidel, Baptiste Rousseau, Benoît Seignovert, Christophe Sotin, Jason W. Barnes, Robert H. Brown, Kevin H. Baines, Bonnie J. Buratti, Roger N. Clark, Philip D. Nicholson, Pascal Rannou, Thomas Cornet

Abstract: We have analyzed the complete archive of the Visual and Infrared Mapping Spectrometer (VIMS) data in order to monitor and analyze the evolution of the clouds and haze coverage at both poles of Titan during the entire Cassini mission. Our objective is to give a cartographic synopsis from a VIMS perspective, to provide a global view of the seasonal evolution of Titan's atmosphere over the poles. We… ▽ More We have analyzed the complete archive of the Visual and Infrared Mapping Spectrometer (VIMS) data in order to monitor and analyze the evolution of the clouds and haze coverage at both poles of Titan during the entire Cassini mission. Our objective is to give a cartographic synopsis from a VIMS perspective, to provide a global view of the seasonal evolution of Titan's atmosphere over the poles. We leave the detailed comparison with the Imaging Science Subsystem (ISS) and the Composite Infrared Spectrometer (CIRS) data sets to further studies. We have computed global hyperspectral mosaics for each of the 127 targeted flybys of Titan to produce synthetic color maps emphasizing the main atmospheric features. The north pole appears fully covered by a huge cloud as soon as the first observations in 2004 and up to the equinox in 2009 (Le Mouélic et al. 2012). The northern skies then became progressively clearer, after the circulation turnover in 2009, revealing the underlying lakes and seas to the optical instruments up to 2017. The reverse situation is observed over the south pole, which was mostly clear of such a high obscuring cloud during the first years of the mission, but started to develop a polar cloud in 2012. This feature grew up month after month until the end of the mission in 2017, with a poleward latitudinal extent of 75$^\circ$S in 2013 up to 58$^\circ$S in April 2017. Thanks to the spectral capabilities of VIMS, we have detected HCN spectral signatures over the north pole in almost all flybys between 2004 and 2008. These HCN signatures started then to show up over the south pole in almost all flybys between 2012 and 2017, so perfectly matching the timing and spatial extent of the northern and southern polar atmospheric features. △ Less

Submitted 2 May, 2018; originally announced May 2018.

Journal ref: Icarus, 2018

LASR-Guided Stellar Photometric Variability Subtraction: The Linear Algorithm For Significance Reduction

Authors: John P. Ahlers, Jason W. Barnes, Sarah A. Horvath, Samuel A. Myers, Matthew M. Hedman

Abstract: We develop a technique for removing stellar variability in the light curves of $δでるた$-Scuti and similar stars. Our technique, which we name the Linear Algorithm for Significance Reduction (LASR), subtracts oscillations from a time series by minimizing their statistical significance in frequency space. We demonstrate that LASR can subtract variable signals of near-arbitrary complexity and can robustly… ▽ More We develop a technique for removing stellar variability in the light curves of $δでるた$-Scuti and similar stars. Our technique, which we name the Linear Algorithm for Significance Reduction (LASR), subtracts oscillations from a time series by minimizing their statistical significance in frequency space. We demonstrate that LASR can subtract variable signals of near-arbitrary complexity and can robustly handle close frequency pairs and overtone frequencies. We demonstrate that our algorithm performs an equivalent fit as prewhitening to the straightforward variable signal of KIC 9700322. We also show that LASR provides a better fit to seismic activity than prewhitening in the case of the complex $δでるた$-Scuti KOI-976. △ Less

Submitted 10 April, 2018; originally announced April 2018.

Comments: 9 pages, 5 figures, accepted for publication in Astronomy & Astrophysics. Pseudocode and github link to code included in manuscript

Journal ref: A&A 615, A128 (2018)

Obliquity Evolution of the Potentially Habitable Exoplanet Kepler-62f

Authors: Billy Quarles, Jason W. Barnes, Jack J. Lissauer, John Chambers

Abstract: Variations in the axial tilt, or obliquity, of terrestrial planets can affect their climates and therefore their habitability. Kepler-62f is a 1.4 R$_\oplus$ planet orbiting within the habitable zone of its K2 dwarf host star (Borucki et al. 2013). We perform N-body simulations that monitor the evolution of obliquity of Kepler-62f for 10 million year timescales to explore the effects on model assu… ▽ More Variations in the axial tilt, or obliquity, of terrestrial planets can affect their climates and therefore their habitability. Kepler-62f is a 1.4 R$_\oplus$ planet orbiting within the habitable zone of its K2 dwarf host star (Borucki et al. 2013). We perform N-body simulations that monitor the evolution of obliquity of Kepler-62f for 10 million year timescales to explore the effects on model assumptions, such as the masses of the Kepler-62 planets and the possibility of outer bodies. Significant obliquity variation occurs when the rotational precession frequency overlaps with one or more of the secular orbital frequencies, but most variations are limited to $\lesssim$10$^\circ$. Moderate variations ($\sim$10$^\circ - 20^\circ$) can occur over a broader range of initial obliquities when the relative nodal longitude ($ΔでるたΩおめが$) overlaps with the frequency and phase of a given secular mode. However, we find that adding outer gas giants on long period orbits ($\gtrsim$ 1000 days) can produce large ($\sim$60$^\circ$) variations in obliquity if Kepler-62f has a very rapid (4 hr) rotation period. The possibility of giant planets on long period orbits impacts the climate and habitability of Kepler-62f through variations in the latitudinal surface flux, where the timescale for large variation can occur on million year timescales. △ Less

Submitted 23 July, 2019; v1 submitted 22 October, 2017; originally announced October 2017.

Comments: 15 pages, 17 figures, 5 tables; Revised relative to Referee Reports. Submitted to Astrobiology

Journal ref: Astrobiology, 20, 73 (2020)

Inference of Heating Properties from "Hot" Non-flaring Plasmas in Active Region Cores. II. Nanoflare Trains

Authors: W. T. Barnes, P. J. Cargill, S. J. Bradshaw

Abstract: Despite its prediction over two decades ago, the detection of faint, high-temperature ("hot") emission due to nanoflare heating in non-flaring active region cores has proved challenging. Using an efficient two-fluid hydrodynamic model, this paper investigates the properties of the emission expected from repeating nanoflares (a nanoflare train) of varying frequency as well as the separate heating o… ▽ More Despite its prediction over two decades ago, the detection of faint, high-temperature ("hot") emission due to nanoflare heating in non-flaring active region cores has proved challenging. Using an efficient two-fluid hydrodynamic model, this paper investigates the properties of the emission expected from repeating nanoflares (a nanoflare train) of varying frequency as well as the separate heating of electrons and ions. If the emission measure distribution ($\mathrm{EM}(T)$) peaks at $T = T_m$, we find that $\mathrm{EM}(T_m)$ is independent of details of the nanoflare train, and $\mathrm{EM}(T)$ above and below $T_m$ reflects different aspects of the heating. Below $T_m$ the main influence is the relationship of the waiting time between successive nanoflares to the nanoflare energy. Above $T_m$ power-law nanoflare distributions lead to an extensive plasma population not present in a monoenergetic train. Furthermore, in some cases characteristic features are present in $\mathrm{EM}(T)$. Such details may be detectable given adequate spectral resolution and a good knowledge of the relevant atomic physics. In the absence of such resolution we propose some metrics that can be used to infer the presence of "hot" plasma. △ Less

Submitted 19 October, 2016; originally announced October 2016.

Comments: 15 pages, 11 figures, accepted to ApJ

Inference of Heating Properties from "Hot" Non-flaring Plasmas in Active Region Cores I. Single Nanoflares

Authors: W. T. Barnes, P. J. Cargill, S. J. Bradshaw

Abstract: The properties expected of "hot" non-flaring plasmas due to nanoflare heating in active regions are investigated using hydrodynamic modeling tools, including a two-fluid development of the EBTEL code. Here we study a single nanoflare and show that while simple models predict an emission measure distribution extending well above 10 MK that is consistent with cooling by thermal conduction, many othe… ▽ More The properties expected of "hot" non-flaring plasmas due to nanoflare heating in active regions are investigated using hydrodynamic modeling tools, including a two-fluid development of the EBTEL code. Here we study a single nanoflare and show that while simple models predict an emission measure distribution extending well above 10 MK that is consistent with cooling by thermal conduction, many other effects are likely to limit the existence and detectability of such plasmas. These include: differential heating between electrons and ions, ionization non-equilibrium and, for short nanoflares, the time taken for the coronal density to increase. The most useful temperature range to look for this plasma, often called the "smoking gun" of nanoflare heating, lies between $10^{6.6}$ and $10^7$ K. Signatures of the actual heating may be detectable in some instances. △ Less

Submitted 16 August, 2016; originally announced August 2016.

Comments: Accepted for publication in ApJ, 16 pages, 8 figures

Obliquity Variability of a Potentially Habitable Early Venus

Authors: Jason W. Barnes, Billy Quarles, Jack J. Lissauer, John Chambers, Matthew M. Hedman

Abstract: Venus currently rotates slowly, with its spin controlled by solid-body and atmospheric thermal tides. However, conditions may have been far different 4 billion years ago, when the Sun was fainter and most of the carbon within Venus could have been in solid form, implying a low-mass atmosphere. We investigate how the obliquity would have varied for a hypothetical rapidly rotating Early Venus. The o… ▽ More Venus currently rotates slowly, with its spin controlled by solid-body and atmospheric thermal tides. However, conditions may have been far different 4 billion years ago, when the Sun was fainter and most of the carbon within Venus could have been in solid form, implying a low-mass atmosphere. We investigate how the obliquity would have varied for a hypothetical rapidly rotating Early Venus. The obliquity variation structure of an ensemble of hypothetical Early Venuses is simpler than that Earth would have if it lacked its large Moon (Lissauer et al., 2012), having just one primary chaotic regime at high prograde obliquities. We note an unexpected long-term variability of up to $\pm7^\circ$ for retrograde Venuses. Low-obliquity Venuses show very low total obliquity variability over billion-year timescales -- comparable to that of the real Moon-influenced Earth. △ Less

Submitted 29 January, 2016; originally announced February 2016.

Comments: in revision at _Astrobiology_

Structure of Titan's evaporites

Authors: D. Cordier, T. Cornet, J. W. Barnes, S. M. MacKenzie, T. Le Bahers, D. Nna-Mvondo, P. Rannou, A. G. Ferreira

Abstract: Numerous geological features that could be evaporitic in origin have been identified on the surface of Titan. Although they seem to be water-ice poor, their main properties -chemical composition, thickness, stratification- are essentially unknown. In this paper, which follows on a previous one focusing on the surface composition (Cordier et al., 2013), we provide some answers to these questions de… ▽ More Numerous geological features that could be evaporitic in origin have been identified on the surface of Titan. Although they seem to be water-ice poor, their main properties -chemical composition, thickness, stratification- are essentially unknown. In this paper, which follows on a previous one focusing on the surface composition (Cordier et al., 2013), we provide some answers to these questions derived from a new model. This model, based on the up-to-date thermodynamic theory known as "PC-SAFT", has been validated with available laboratory measurements and specifically developed for our purpose. 1-D models confirm the possibility of an acetylene and/or butane enriched central layer of evaporitic deposit. The estimated thickness of this acetylene-butane layer could explain the strong RADAR brightness of the evaporites. The 2-D computations indicate an accumulation of poorly soluble species at the deposit's margin. Among these species, HCN or aerosols similar to tholins could play a dominant role. Our model predicts the existence of chemically trimodal "bathtub rings" which is consistent with what it is observed at the south polar lake Ontario Lacus. This work also provides plausible explanations to the lack of evaporites in the south polar region and to the high radar reflectivity of dry lakebeds. △ Less

Submitted 22 December, 2015; originally announced December 2015.

Probable Spin-Orbit Aligned Super-Earth Planet Candidate KOI-2138.01

Authors: Jason W. Barnes, Johnathon P. Ahlers, Shayne A. Seubert, Howard M. Relles

Abstract: We use rotational gravity darkening in the disk of \emph{Kepler} star KOI-2138 to show that the orbit of $2.1-R_\oplus$ transiting planet candidate KOI-2138.01 has a low projected spin-orbit alignment of $λらむだ=1^\circ\pm13$. KOI-2138.01 is just the second super-Earth with a measured spin-orbit alignment after 55 Cancri e, and the first to be aligned. With a 23.55-day orbital period, KOI-2138.01 may r… ▽ More We use rotational gravity darkening in the disk of \emph{Kepler} star KOI-2138 to show that the orbit of $2.1-R_\oplus$ transiting planet candidate KOI-2138.01 has a low projected spin-orbit alignment of $λらむだ=1^\circ\pm13$. KOI-2138.01 is just the second super-Earth with a measured spin-orbit alignment after 55 Cancri e, and the first to be aligned. With a 23.55-day orbital period, KOI-2138.01 may represent the tip of a future iceberg of solar-system-like terrestrial planets having intermediate periods and low-inclination circular orbits. △ Less

Submitted 11 December, 2015; originally announced December 2015.

Comments: 5 pages, 4 figures, published in ApJL 2015 August

Journal ref: The Astrophysical Journal Letters 808 (2), L38, 2015 August 1

Meridional variation in tropospheric methane on Titan observed with AO spectroscopy at Keck and VLT

Authors: Máté Ádámkovics, Jonathan L. Mitchell, Alexander G. Hayes, Patricio M. Rojo, Paul Corlies, Jason W. Barnes, Valentin D. Ivanov, Robert H. Brown, Kevin H. Baines, Bonnie J. Buratti, Roger N. Clark, Philip D. Nicholson, Christophe Sotin

Abstract: The spatial distribution of the tropospheric methane on Titan was measured using near-infrared spectroscopy. Ground-based observations at 1.5$μみゅー{\rm m}$ (H-band) were performed during the same night using instruments with adaptive optics at both the W. M. Keck Observatory and at the Paranal Observatory on 17 July 2014 UT. The integral field observations with SINFONI on the VLT covered the entire H-… ▽ More The spatial distribution of the tropospheric methane on Titan was measured using near-infrared spectroscopy. Ground-based observations at 1.5$μみゅー{\rm m}$ (H-band) were performed during the same night using instruments with adaptive optics at both the W. M. Keck Observatory and at the Paranal Observatory on 17 July 2014 UT. The integral field observations with SINFONI on the VLT covered the entire H-band at moderate resolving power, $R=λらむだ/Δでるたλらむだ\approx1,500$, while the Keck observations were performed with NIRSPAO near 1.55254$μみゅー{\rm m}$ at higher resolution, $R\approx25,000$. The moderate resolution observations are used for flux calibration and for the determination of model parameters that can be degenerate in the interpretation of high resolution spectra. Line-by-line calculations of CH$_4$ and CH$_3$D correlated $k$ distributions from the HITRAN 2012 database were used, which incorporate revised line assignments near 1.5$μみゅー{\rm m}$. We fit the surface albedo and aerosol distributions in the VLT SINFONI observations that cover the entire H-band window and used these quantities to constrain the models of the high-resolution Keck NIRSPAO spectra when retrieving the methane abundances. Cassini VIMS images of the polar regions, acquired on 20 July 2014 UT, are used to validate the assumption that the opacity of tropospheric aerosol is relatively uniform below 10 km. We retrieved methane abundances at latitudes between 42$^{\circ}$ S and 80$^{\circ}$ N. The tropospheric methane in the Southern mid-latitudes was enhanced by a factor of $\sim$10--40% over the nominal profile that was measured using the GCMS on Huygens. The Northern hemisphere had $\sim$90% of the nominal methane abundance up to polar latitudes (80$^{\circ}$N). These measurements suggest that a source of saturated polar air is equilibrating with dryer conditions at lower latitudes. △ Less

Submitted 29 September, 2015; originally announced September 2015.

Comments: Accepted for publication in Icarus on May 22, 2015

Spin-Orbit Misalignment of Two-Planet-System KOI-89 Via Gravity Darkening

Authors: John P. Ahlers, Jason W. Barnes, Rory Barnes

Abstract: We constrain the true spin-orbit alignment of the KOI-89 system by numerically fitting the two \emph{Kepler} photometric lightcurves produced by transiting planets KOI-89.01 and KOI-89.02. The two planets have periods of 84.69 days and 207.58 days, respectively. We find that the two bodies are low-density giant planets with radii $0.45 \pm 0.03~\mathrm{R_{jup}}$ and… ▽ More We constrain the true spin-orbit alignment of the KOI-89 system by numerically fitting the two \emph{Kepler} photometric lightcurves produced by transiting planets KOI-89.01 and KOI-89.02. The two planets have periods of 84.69 days and 207.58 days, respectively. We find that the two bodies are low-density giant planets with radii $0.45 \pm 0.03~\mathrm{R_{jup}}$ and $0.43 \pm 0.05~\mathrm{R_{jup}}$ and spin-orbit misalignments $72^{\circ} \pm 3^\circ$ and $73^{\circ+11}_{-5}$, respectively. Via dynamic stability tests we demonstrate the general trend of higher system stability with the two planets close to mutual alignment and estimate their coalignment angle to $20^\circ \pm 20^\circ$ -- i.e. the planets are misaligned with the star but may be aligned with each other. From these results, we limit KOI-89's misalignment mechanisms to star-disk-binary interactions, disk warping via planet-disk interactions, planet-planet scattering, Kozai resonance, or internal gravity waves. △ Less

Submitted 3 November, 2015; v1 submitted 25 September, 2015; originally announced September 2015.

Comments: 9 pages, 7 figures, 3 tables, accepted in The Astrophysical Journal

Follow-Up Observations of PTFO 8-8695: A 3 MYr Old T-Tauri Star Hosting a Jupiter-mass Planetary Candidate

Authors: David R. Ciardi, J. C. van Eyken, J. W. Barnes, C. A. Beichman, S. J. Carey, C. J. Crockett, J. Eastman, C. M. Johns-Krull, S. B. Howell, S. R. Kane, J. N. Mclane, P. Plavchan, L. Prato, J. Stauffer, G. T. van Belle, K. von Braun

Abstract: We present Spitzer 4.5\micron\ light curve observations, Keck NIRSPEC radial velocity observations, and LCOGT optical light curve observations of PTFO~8-8695, which may host a Jupiter-sized planet in a very short orbital period (0.45 days). Previous work by \citet{vaneyken12} and \citet{barnes13} predicts that the stellar rotation axis and the planetary orbital plane should precess with a period o… ▽ More We present Spitzer 4.5\micron\ light curve observations, Keck NIRSPEC radial velocity observations, and LCOGT optical light curve observations of PTFO~8-8695, which may host a Jupiter-sized planet in a very short orbital period (0.45 days). Previous work by \citet{vaneyken12} and \citet{barnes13} predicts that the stellar rotation axis and the planetary orbital plane should precess with a period of $300 - 600$ days. As a consequence, the observed transits should change shape and depth, disappear, and reappear with the precession. Our observations indicate the long-term presence of the transit events ($>3$ years), and that the transits indeed do change depth, disappear and reappear. The Spitzer observations and the NIRSPEC radial velocity observations (with contemporaneous LCOGT optical light curve data) are consistent with the predicted transit times and depths for the $M_\star = 0.34\ M_\odot$ precession model and demonstrate the disappearance of the transits. An LCOGT optical light curve shows that the transits do reappear approximately 1 year later. The observed transits occur at the times predicted by a straight-forward propagation of the transit ephemeris. The precession model correctly predicts the depth and time of the Spitzer transit and the lack of a transit at the time of the NIRSPEC radial velocity observations. However, the precession model predicts the return of the transits approximately 1 month later than observed by LCOGT. Overall, the data are suggestive that the planetary interpretation of the observed transit events may indeed be correct, but the precession model and data are currently insufficient to confirm firmly the planetary status of PTFO~8-8695b. △ Less

Submitted 29 June, 2015; originally announced June 2015.

Comments: Accepted for publication in The Astrophysical Journal

Evidence of Titan's Climate History from Evaporite Distribution

Authors: Shannon M. MacKenzie, Jason W. Barnes, Christophe Sotin, Jason M. Soderblom, Stéphane Le Mouélic, Sebastien Rodriguez, Kevin H. Baines, Bonnie J. Buratti, Roger N. Clark, Phillip D. Nicholson, Thomas B. McCord

Abstract: Water-ice-poor, 5-$μみゅー$m-bright material on Saturn's moon Titan has previously been geomorphologically identified as evaporitic. Here we present a global distribution of the occurrences of the 5-$μみゅー$m-bright spectral unit, identified with Cassini's Visual Infrared Mapping Spectrometer (VIMS) and examined with RADAR when possible. We explore the possibility that each of these occurrences are evaporite… ▽ More Water-ice-poor, 5-$μみゅー$m-bright material on Saturn's moon Titan has previously been geomorphologically identified as evaporitic. Here we present a global distribution of the occurrences of the 5-$μみゅー$m-bright spectral unit, identified with Cassini's Visual Infrared Mapping Spectrometer (VIMS) and examined with RADAR when possible. We explore the possibility that each of these occurrences are evaporite deposits. The 5-$μみゅー$m-bright material covers 1\% of Titan's surface and is not limited to the poles (the only regions with extensive, long-lived surface liquid). We find the greatest areal concentration to be in the equatorial basins Tui Regio and Hotei Regio. Our interpretations, based on the correlation between 5-$μみゅー$m-bright material and lakebeds, imply that there was enough liquid present at some time to create the observed 5-$μみゅー$m-bright material. We address the climate implications surrounding a lack of evaporitic material at the south polar basins: if the south pole basins were filled at some point in the past, then where is the evaporite? △ Less

Submitted 12 August, 2014; originally announced August 2014.

Spin-Orbit Alignment For 110-Day-Period KOI368.01 From Gravity Darkening

Authors: John. P. Ahlers, Shayne A. Seubert, Jason W. Barnes

Abstract: We fit the \emph{Kepler} photometric light curve of the KOI-368 system using an oblate, gravity-darkened stellar model in order to constrain its spin-orbit alignment. We find that the system is relatively well-aligned with a sky-projected spin-orbit alignment of $λらむだ=10^{\circ}\pm2^{\circ}$, a stellar obliquity of $ψぷさい=3^{\circ}\pm7^{\circ}$, and a true spin-orbit alignment of… ▽ More We fit the \emph{Kepler} photometric light curve of the KOI-368 system using an oblate, gravity-darkened stellar model in order to constrain its spin-orbit alignment. We find that the system is relatively well-aligned with a sky-projected spin-orbit alignment of $λらむだ=10^{\circ}\pm2^{\circ}$, a stellar obliquity of $ψぷさい=3^{\circ}\pm7^{\circ}$, and a true spin-orbit alignment of $\varphi=11^{\circ}\pm3^{\circ}$. Although our measurement differs significantly from zero, the low value for $\varphi$ is consistent with spin-orbit alignment. We also measure various transit parameters of the KOI-368 system: $R_{\star}=2.28\pm0.02 R_{\odot}$, $R_{p}=1.83\pm0.02R_{jup}$, and $i=89.221^{\circ}\pm0.013^{\circ}$. This work shows that our gravity-darkened model can constrain long-period, well-aligned planets and M-class stars orbiting fast-rotators, allowing for measurement of a new subcategory of transiting bodies. △ Less

Submitted 3 April, 2014; originally announced April 2014.

Comments: 5 Figures, 6 Pages

Measurement of Spin-Orbit Misalignment and Nodal Precession for the Planet around Pre-Main-Sequence Star PTFO 8-8695 From Gravity Darkening

Authors: Jason W. Barnes, Julian C. van Eyken, Brian K. Jackson, David R. Ciardi, Jonathan J. Fortney

Abstract: PTFO 8-8695b represents the first transiting exoplanet candidate orbiting a pre-main-sequence star. We find that the unusual lightcurve shapes of PTFO 8-8695 can be explained by transits of a planet across an oblate, gravity-darkened stellar disk. We simultaneously and self-consistently fit two separate lightcurves observed in 2009 December and 2010 December. Our two self-consistent fits yield M_p… ▽ More PTFO 8-8695b represents the first transiting exoplanet candidate orbiting a pre-main-sequence star. We find that the unusual lightcurve shapes of PTFO 8-8695 can be explained by transits of a planet across an oblate, gravity-darkened stellar disk. We simultaneously and self-consistently fit two separate lightcurves observed in 2009 December and 2010 December. Our two self-consistent fits yield M_p = 3.0 M_Jup and M_p = 3.6 M_Jup for assumed stellar masses of M_* = 0.34 M_Sun and M_* = 0.44 M_Sun respectively. The two fits have precession periods of 293 days and 581 days. These mass determinations (consistent with previous upper limits) along with the strength of the gravity-darkened precessing model together validate PTFO 8-8695b as just the second Hot Jupiter known to orbit an M-dwarf. Our fits show a high degree of spin-orbit misalignment in the PTFO 8-8695 system: 69 +/- 2 or 73.1 +/- 0.5 degrees, in the two cases. The large misalignment is consistent with the hypothesis that planets become Hot Jupiters with random orbital plane alignments early in a system's lifetime. We predict that as a result of the highly misaligned, precessing system, the transits should disappear for months at a time over the course of the system's precession period. The precessing, gravity-darkened model also predicts other observable effects: changing orbit inclination that could be detected by radial velocity observations, changing stellar inclination that would manifest as varying v sin i, changing projected spin-orbit alignment that could be seen by the Rossiter-McLaughlin effect, changing transit shapes over the course of the precession, and differing lightcurves as a function of wavelength. Our measured planet radii of 1.64 R_Jup and 1.68 R_Jup in each case are consistent with a young, hydrogen-dominated planet that results from a hot-start formation mechanism. △ Less

Submitted 2 August, 2013; originally announced August 2013.

Comments: Accepted at ApJ

MOST Space Telescope Photometry of the 2010 January Transit of Extrasolar Planet HD80606b

Authors: Jessica E. Roberts, Jason W. Barnes, Jason F. Rowe, Jonathan F. Fortney

Abstract: We present observations of the full January 2010 transit of HD80606b from the Canadian microsatellite, Microvariability and Oscillations of Stars (MOST). By employing a space-based telescope, we monitor the entire transit thus limiting systematic errors that result from ground observations. We determine measurements for the planetary radius (R_{p}=0.987\pm0.061R_{Jup}) and inclination (i=89.283^{o… ▽ More We present observations of the full January 2010 transit of HD80606b from the Canadian microsatellite, Microvariability and Oscillations of Stars (MOST). By employing a space-based telescope, we monitor the entire transit thus limiting systematic errors that result from ground observations. We determine measurements for the planetary radius (R_{p}=0.987\pm0.061R_{Jup}) and inclination (i=89.283^{o}\pm0.024) by constraining our fits with the observed parameters of different groups. Our measured mid-transit time of 2455210.6449\pm0.0034 HJD is consistant with the 2010 Spitzer results and is 20 minutes earlier than predicted by groups who observed the June 2009 transit. △ Less

Submitted 3 December, 2012; originally announced December 2012.

Comments: 3 figures

Outcomes and Duration of Tidal Evolution in a Star-Planet-Moon System

Authors: Takashi Sasaki, Jason W. Barnes, David P. O'Brien

Abstract: We formulated tidal decay lifetimes for hypothetical moons orbiting extrasolar planets with both lunar and stellar tides. Previous work neglected the effect of lunar tides on planet rotation, and are therefore applicable only to systems in which the moon's mass is much less than that of the planet. This work, in contrast, can be applied to the relatively large moons that might be detected around n… ▽ More We formulated tidal decay lifetimes for hypothetical moons orbiting extrasolar planets with both lunar and stellar tides. Previous work neglected the effect of lunar tides on planet rotation, and are therefore applicable only to systems in which the moon's mass is much less than that of the planet. This work, in contrast, can be applied to the relatively large moons that might be detected around newly-discovered Neptune-mass and super-Earth planets. We conclude that moons are more stable when the planet/moon systems are further from the parent star, the planets are heavier, or the parent stars are lighter. Inclusion of lunar tides allows for significantly longer lifetimes for a massive moon relative to prior formulations. We expect that the semi-major axis of the planet hosting the first detected exomoon around a G-type star is 0.4-0.6 AUえーゆー and is 0.2-0.4 AUえーゆー for an M-type star. △ Less

Submitted 1 June, 2012; originally announced June 2012.

Comments: Accepted for publication in ApJ, 19 pages, 19 figures

Using Satellites to Probe Extrasolar Planet Formation

Authors: Paul Withers, Jason W. Barnes

Abstract: Planetary satellites are an integral part of the heirarchy of planetary systems. Here we make two predictions concerning their formation. First, primordial satellites, which have an array of distinguishing characteristics, form only around giant planets. If true, the size and duration of a planetary system's protostellar nebula, as well as the location of its snow line, can be constrained by knowi… ▽ More Planetary satellites are an integral part of the heirarchy of planetary systems. Here we make two predictions concerning their formation. First, primordial satellites, which have an array of distinguishing characteristics, form only around giant planets. If true, the size and duration of a planetary system's protostellar nebula, as well as the location of its snow line, can be constrained by knowing which of its planets possess primordial satellites and which do not. Second, all satellites around terrestrial planets form by impacts. If true, this greatly enhances the constraints that can be placed on the history of terrestrial planets by their satellites' compositions, sizes, and dynamics. △ Less

Submitted 4 April, 2012; originally announced April 2012.

Journal ref: Proceedings of the International Astronomical Union October 2010 6 : E1 (2 pages)

The EVIL-MC Model for Ellipsoidal Variations of Planet-Hosting Stars and Applications to the HAT-P-7 System

Authors: Brian K. Jackson, Nikole K. Lewis, Jason W. Barnes, L. Drake Deming, Adam P. Showman, Jonathan J. Fortney

Abstract: We present a new model for Ellipsoidal Variations Induced by a Low-Mass Companion, the EVIL-MC model. We employ several approximations appropriate for planetary systems to substantially increase the computational efficiency of our model relative to more general ellipsoidal variation models and improve upon the accuracy of simpler models. This new approach gives us a unique ability to rapidly and a… ▽ More We present a new model for Ellipsoidal Variations Induced by a Low-Mass Companion, the EVIL-MC model. We employ several approximations appropriate for planetary systems to substantially increase the computational efficiency of our model relative to more general ellipsoidal variation models and improve upon the accuracy of simpler models. This new approach gives us a unique ability to rapidly and accurately determine planetary system parameters. We use the EVIL-MC model to analyze Kepler Quarter 0-2 (Q0-2) observations of the HAT-P-7 system, an F-type star orbited by a nearly Jupiter-mass companion. Our analysis corroborates previous estimates of the planet-star mass ratio q = (1.10 +/- 0.06) x 10^(-3), and we have revised the planet's dayside brightness temperature to 2680 +10/-20 K. We also find a large difference between the day- and nightside planetary flux, with little nightside emission. Preliminary dynamical+radiative modeling of the atmosphere indicates this result is qualitatively consistent with high altitude absorption of stellar heating. Similar analyses of Kepler and CoRoT photometry of other planets using EVIL-MC will play a key role in providing constraints on the properties of many extrasolar systems, especially given the limited resources for follow-up and characterization of these systems. However, as we highlight, there are important degeneracies between the contributions from ellipsoidal variations and planetary emission and reflection. Consequently, for many of the hottest and brightest Kepler and CoRoT planets, accurate estimates of the planetary emission and reflection, diagnostic of atmospheric heat budgets, will require accurate modeling of the photometric contribution from the stellar ellipsoidal variation. △ Less

Submitted 19 April, 2012; v1 submitted 27 March, 2012; originally announced March 2012.

Comments: Accepted to ApJ; minor revisions to original submission; An IDL version of the EVIL-MC model is publicly available at http://www.lpl.arizona.edu/~bjackson/idl_code/index.html

Measurement of the Spin-Orbit Misalignment of KOI-13.01 from its Gravity-Darkened Kepler Transit Lightcurve

Authors: Jason W. Barnes, Ethan Linscott, Avi Shporer

Abstract: We model the asymmetry of the KOI-13.01 transit lightcurve assuming a gravity-darkened rapidly-rotating host star in order to constrain the system's spin-orbit alignment and transit parameters. We find that our model can reproduce the Kepler lightcurve for KOI-13.01 with a sky-projected alignment of λらむだ = $23 \pm 4$ degrees and with the star's north pole tilted away from the observer by $48 \pm 4$ d… ▽ More We model the asymmetry of the KOI-13.01 transit lightcurve assuming a gravity-darkened rapidly-rotating host star in order to constrain the system's spin-orbit alignment and transit parameters. We find that our model can reproduce the Kepler lightcurve for KOI-13.01 with a sky-projected alignment of λらむだ = $23 \pm 4$ degrees and with the star's north pole tilted away from the observer by $48 \pm 4$ degrees (assuming $M_* = 2.05 M_\odot$). With both these determinations, we calculate that the net misalignment between this planet's orbit normal and its star's rotational pole is $56 \pm 4$ degrees. Degeneracies in our geometric interpretation also allow a retrograde spin-orbit angle of $124 \pm 4$ degrees. This is the first spin-orbit measurement to come from gravity darkening, and is one of only a few measurements of the full (not just the sky-projected) spin-orbit misalignment of an extrasolar planet. We also measure accurate transit parameters incorporating stellar oblateness and gravity darkening: $R_* = 1.756 \pm 0.014 R_\odot$, $R_p = 1.445 \pm 0.016 R_{Jup}$, and $i = 85.9 \pm 0.4$ degrees. The new lower planetary radius falls within the planetary mass regime for plausible interior models for the transiting body. A simple initial calculation shows that KOI-13.01's circular orbit is apparently inconsistent with the Kozai mechanism having driven its spin-orbit misalignment; planet-planet scattering and stellar spin migration remain viable mechanisms. Future Kepler data will improve the precision of the KOI-13.01 transit lightcurve, allowing more precise determination of transit parameters and the opportunity to use the Photometric Rossiter-McLaughlin effect to resolve the prograde/retrograde orbit determination degeneracy. △ Less

Submitted 16 October, 2011; originally announced October 2011.

Comments: 10 pages, 5 figures, in press in ApJS Kepler Special Issue

Observations of a Stationary Mid-Latitude Cloud System on Titan

Authors: Máté Ádámkovics, Jason W. Barnes, Markus Hartung, Imke de Pater

Abstract: We report the observation of a cloud system on Titan that remained localized near 40S latitude and 60W longitude for at least 34 hours. Ground-based observations obtained with the SINFONI imaging spectrograph at the Very Large Telescope over 4 consecutive nights recorded the lifetime and altitude of the unresolved cloud system. Concomitant measurements made by Cassini/VIMS over 3 hours resolved ch… ▽ More We report the observation of a cloud system on Titan that remained localized near 40S latitude and 60W longitude for at least 34 hours. Ground-based observations obtained with the SINFONI imaging spectrograph at the Very Large Telescope over 4 consecutive nights recorded the lifetime and altitude of the unresolved cloud system. Concomitant measurements made by Cassini/VIMS over 3 hours resolved changes in the altitude and opacity of individual regions within the system during this time. Clouds are measured from 13 to 37 km altitude with optical depths per pixel ranging from $τたう$=0.13 to 7. Short timescale rise times are consistent with previous measurements of the evolution of mid-latitude clouds; however the long timescale localization of the cloud structure is unexplained. We speculate about the role of meso-scale circulation in relation to cloud formation. △ Less

Submitted 12 March, 2010; originally announced March 2010.

Comments: Accepted for publication in Icarus

Detecting the Wind-Driven Shapes of Extrasolar Giant Planets from Transit Photometry

Authors: Jason W. Barnes, Curtis S. Cooper, Adam P. Showman, William B. Hubbard

Abstract: Several processes can cause the shape of an extrasolar giant planet's shadow, as viewed in transit, to depart from circular. In addition to rotational effects, cloud formation, non-homogenous haze production and movement, and dynamical effects (winds) could also be important. When such a planet transits its host star as seen from Earth, the asphericity will introduce a deviation in the transit l… ▽ More Several processes can cause the shape of an extrasolar giant planet's shadow, as viewed in transit, to depart from circular. In addition to rotational effects, cloud formation, non-homogenous haze production and movement, and dynamical effects (winds) could also be important. When such a planet transits its host star as seen from Earth, the asphericity will introduce a deviation in the transit lightcurve relative to the transit of a perfectly spherical (or perfectly oblate) planet. We develop a theoretical framework to interpret planet shapes. We then generate predictions for transiting planet shapes based on a published theoretical dynamical model of HD189733b. Using these shape models we show that planet shapes are unlikely to introduce detectable lightcurve deviations (those >~1e-5 of the host star), but that the shapes may lead to astrophysical sources of systematic error when measuring planetary oblateness, transit time, and impact parameter. △ Less

Submitted 19 October, 2009; originally announced October 2009.

Comments: ApJ, accepted -- scheduled for publication 2009 December 1

Journal ref: Astrophys.J.706:877-884,2009

Transit Lightcurves of Extrasolar Planets Orbiting Rapidly-Rotating Stars

Authors: Jason W. Barnes

Abstract: Main-sequence stars earlier than spectral type ~F6 or so are expected to rotate rapidly due to their radiative exteriors. This rapid rotation leads to an oblate stellar figure. It also induces the photosphere to be hotter (by up to several thousand Kelvin) at the pole than at the equator as a result of a process called gravity darkening that was first predicted by von Zeipel (1924). Transits of… ▽ More Main-sequence stars earlier than spectral type ~F6 or so are expected to rotate rapidly due to their radiative exteriors. This rapid rotation leads to an oblate stellar figure. It also induces the photosphere to be hotter (by up to several thousand Kelvin) at the pole than at the equator as a result of a process called gravity darkening that was first predicted by von Zeipel (1924). Transits of extrasolar planets across such a non-uniform, oblate disk yield unusual and distinctive lightcurves that can be used to determine the relative alignment of the stellar rotation pole and the planet orbit normal. This spin-orbit alignment can be used to constrain models of planet formation and evolution. Orderly planet formation and migration within a disk that is coplanar with the stellar equator will result in spin-orbit alignment. More violent planet-planet scattering events should yield spin-orbit misaligned planets. Rossiter-McLaughlin measurements of transits of lower-mass stars show that some planets are spin-orbit aligned, and some are not. Since Rossiter-McLaughlin measurements are difficult around rapid rotators, lightcurve photometry may be the best way to determine the spin-orbit alignment of planets around massive stars. The Kepler mission will monitor ~10^4 of these stars within its sample. The lightcurves of any detected planets will allow us to probe the planet formation process around high-mass stars for the first time. △ Less

Submitted 9 September, 2009; originally announced September 2009.

Comments: 14 pages, 9 figures, accepted to ApJ

Journal ref: Astrophys.J.705:683-692,2009

Evidence for condensed-phase methane enhancement over Xanadu on Titan

Authors: M. Ádámkovics, I. de Pater, M. Hartung, J. W. Barnes

Abstract: We present evidence for condensed phase methane precipitation near Xanadu using nine nights of observations from the SINFONI integral-field spectrograph at the Very Large Telescope and imaging analysis with empirical surface subtraction. Radiative transfer models are used to support the imaging technique by simulating the spectrometer datacubes and testing for variations in both the surface refl… ▽ More We present evidence for condensed phase methane precipitation near Xanadu using nine nights of observations from the SINFONI integral-field spectrograph at the Very Large Telescope and imaging analysis with empirical surface subtraction. Radiative transfer models are used to support the imaging technique by simulating the spectrometer datacubes and testing for variations in both the surface reflectivity spectrum and atmospheric opacity. We use the models and observations together to argue against artifacts that may arise in the image analysis. High phase angle observations from Cassini/VIMS are used to test against surface scattering artifacts that may be confused with sources of atmospheric opacity. Although changes in the surface reflectivity spectrum can reproduce observations from a particular viewing geometry on a given night, multiple observations are best modeled by condensed-phase methane opacity near the surface. These observations and modeling indicate that the condensed-phase methane opacity observed with this technique occurs predominantly near Xanadu and is most likely due to precipitation. △ Less

Submitted 13 July, 2009; originally announced July 2009.

Comments: accepted for publication in Planetary and Space Science; http://dx.doi.org/10.1016/j.pss.2009.07.001