THE ATMOSPHERES AND ORBITAL DYNAMICS OF HOT JUPITERS
Matthew Pryal Archbald, Pennsylvania
B.S. Astronomy & Astrophysics, Penn State University, 2015
M.S. Astronomy, University of Virginia, 2017
A Dissertation Presented to the Graduate Faculty of the University of Virginia in Candidacy for the Degree of Doctor of Philosophy
Department of Astronomy
University of Virginia December 2020
Committee Members: Phillip L. Arras Zhi-Yun Li Tim Bastian Ajay B. Limaye c Copyright by Matthew Pryal
All rights reserved
December 18, 2020 To my family and friends. Whom without, I would have finished this thesis much, much earlier. iii Abstract
The detection of the first exoplanets in the early 1990s brought along many questions about their characteristics, including on their atmospheric structure and the dynam- ics of their planetary systems. One of the most puzzling types of exoplanets were the so-called ‘hot Jupiters.’ Hot Jupiters are Jupiter sized objects, orbiting very close to their host star, that absorb a large amount of radiation, making their atmospheres very hot. Their short orbital periods and hot, extended atmospheres are some of the many characteristics that make these objects unlike any planet in our Solar System. In this thesis, I present models of the electrical conductivity of the upper atmosphere of hot Jupiters as well as the migration of hypothetical moons orbiting hot Jupiters due to tidal friction effects.
In Chapter1, I give an introduction and provide context to the research outlined in the thesis. Specifically, I give an overview of the detection of exoplanets and discuss HD 189733b, the hot Jupiter at the focus of my atmosphere model. I also introduce the interaction of a stellar wind with a planet’s magnetosphere which is central to- wards setting the current that can run through the atmosphere, which is then affected by the atmosphere’s electrical conductivity. Finally, I overview the dynamics of star- planet-moons systems and tidal friction as it relates to the hypothetical moons of hot Jupiters.
In Chapter2, I detail the star-planet interaction between the charged particles of the stellar wind with a planetary magnetic field that set the current in the atmosphere. Importantly, the electrical conductivity of the atmosphere governs the rate at which the charged particles can flow and release energy into the atmosphere through col- lisions. I also review the formulas for electrical conductivity of an atmosphere and discuss the collision rates between charged particles and other species of the atmo- sphere which is then implemented into the HD 189733b atmosphere model in Chapter 3.
In Chapter3, I present a hydrostatic model of the dayside upper atmosphere of hot Jupiter HD 189733b. With this model I compute its electrical conductivity and conductance, which determines the amount of energy that can be transferred to the atmosphere through Joule heating. I find that the conductance of hot Jupiters is much greater than the conductance of the Earth, and therefore Joule heating does not significantly impact the thermal structure of the atmosphere for Jupiter sized magnetic fields (B ∼ 10 G).
Finally, in Chapter4, I examine the orbital migration of hypothetical moons of hot iv
Jupiters. I find that, even in the case of synchronous rotation and circular orbits of the moon, gravitational perturbations from the star create a forced eccentricity in the moon that makes tidal friction never cease. This tidal friction then causes the inward migration of the moon towards the planet where it can be tidally disrupted. In the case of hot Jupiter systems, this inward migration is fast enough to explain the lack of large moons detected orbiting hot Jupiters. Additionally, I place an upper limit on the mass of moons that are able to survive this orbital migration effect within 5 Gyr, to be many orders of magnitude smaller than the mass of the Earth’s moon. v Acknowledgements
My love for Astronomy began in a ninth grade Earth and Space Science class, that I was warned by my older siblings was taught by a ‘crazy alien guy’ – who I would later find out preferred to go by Mr. Durkin. Over the course of that class, I realized that Mr. Durkin’s ‘craziness,’ was actually child-like wonder and passion for a science that can be appreciated by anyone simply by looking up at the night sky. Unknown to him, Mr. Durkin’s passion was passed to me, and over my years in Astronomy, I realized that I wanted to become that ‘crazy alien guy’ for others. It is truly unbelievable to think how much Mr. Durkin’s passion about Astronomy has impacted my path in life – taking me from a four bedroom house in tiny Archbald, PA (at one time shared by 10 people!), to Penn State University, and Rome, on cross country trips, to the University of Virginia, and Hawaii and somehow to the stars and back. I urge everyone to embrace the passions in your life and share them with as many people as you can. You never know how much of an impact your passion can have on others. At each milestone I reach in my life, I can’t help but think of all the little things that those in my life have taught me to get me to where I am. Grad school has admittedly been emotionally hard. The year 2020 has been especially stressful. But with everything I’ve been taught by family and close friends, we’ve been able to reach this finish line together. Thank you to my Mom and Dad for filling me with love and being proud of me every step of the way. For teaching me empathy, kindness, how to deal with sadness, and embrace my passions. Thanks, to my siblings Bridget, Sarah, Robbie, Lizzy, and Emma for teaching me how to be a family. To learn how to deal with anger (don’t punch walls), to be more outgoing, to love unconditionally, and handle disagreements. One of the few good things about growing up, has been becoming closer to all of you each and every year and I can’t wait for our relationships to continue to grow. Being the first person in my family to get a Ph.D., I can’t help but think of my Grammie, Grampie, Grandma, and Grandpa. Thank you for creating a better life for all those around you and laying the groundwork for your loved ones to be able to dream big today. My thanks to my family could extend forever but since acknowledgements are supposed to be shorter than the thesis, I’ll just say thank you to rest of my extended family for making my life as full as it can be. Thanks to my Brohana, my cohort, old Charlottesville friends, and the rest of those at UVA and in Charlottesville for making my time in grad school fun. Thanks for grilling and chilling, having tailgates, chatting around a fire, for going to Nelly’s or Escaf´e,going to basketball games, playing poker, for Foxfield, and Raven banquets, for Festy and for everything else in between. Grad school would have been unbearable without all of you in it. Thank you to the Awesome Astronomy Enthusiasts of 2015 and everyone else at vi
Penn State for teaching me that life’s about 8 things and 8 things only. For being a bunch of nerds who embraced our shared weirdness, for really knowing how to have fun, and for always saying hi to me in the elevator. Integrity. Thank you to my NEPA friends for being a part of some of the most fundamental parts of my life and for effectively helping me grow up – although I wouldn’t say we’re there yet, and I’m not sure if we’ll ever quite make it. Thanks to my thesis committee and thesis advisor Phil for making my research as great as it could be over the past 5.5 years. And to the rest of the UVA and Penn State Astronomy departments – thanks for your kindness and support along the way. And a special thanks to all those who I’ve had the pleasure of doing outreach with over the years. Especially to all the members of Dark Skies, Bright Kids, Astronomy on Tap - C’ville, those at the Science Museum of Virginia, and everyone else. You’ve helped make the journey worthwhile and I certainly would have been much crazier without you in my life. And last, but certainly not least (except in physical size) – thank you Molly for making one of the most stressful parts of my life, somehow the happiest. I was only able to get to this finish line because of your support these past few years. I eagerly await any stressful days ahead, as long as we’re facing them together. Also, feather.
With love, The other crazy alien guy Table of contents
List of Figures xiv
List of Tables xvi
1 Introduction1 1.1 Exoplanets...... 1 1.1.1 The Discovery of Exoplanets...... 2 1.1.2 The Detection of Exoplanets...... 3 1.1.3 Hot Jupiters...... 9 1.2 The Upper Atmosphere of HD 189733b...... 11 1.2.1 Spectroscopic Observations of Exoplanet Atmospheres.... 12 1.2.2 Chemical Species in the Atmosphere of Hot Jupiters...... 13 1.2.3 Upper Atmosphere Models of Hot Jupiters...... 15 1.3 Planetary Magnetospheres and Ionospheres...... 16 1.3.1 Overview...... 17 1.3.2 The Solar Wind and Magnetic Reconnection...... 18 1.4 The Dynamics of Star-Planet-Moon Systems...... 20 1.4.1 Introduction...... 22 1.4.2 The Lack of Moons in the Inner Solar System...... 23 1.4.3 Tidal Friction in Two Body Systems...... 23 1.4.4 The Three-Body Problem...... 27 1.4.5 Model of a Star-Planet-Moon System...... 28 1.5 Overview...... 29
2 Electric Currents and Conductivity in Planetary Atmospheres 31 2.1 Context...... 31 2.2 Historical Overview...... 32 2.3 Electric Currents Overview...... 34 2.4 Joule Heating...... 36 2.4.1 Height Integrated Currents and Conductivity...... 36 2.4.2 Parallel Currents and Alfv´enConductance...... 37 2.4.3 Reflection of Alfv´enWaves...... 39
vii viii
2.4.4 Joule Heating Calculation...... 41 2.5 Electrical Conductivity Overview...... 43 2.6 Deriving the Conductivity Formulas...... 45 2.6.1 Conductivity Across a B-field...... 45 2.6.2 Conductivity Parallel to the B-field...... 48 2.7 Height Integrated Conductivity...... 50 2.7.1 Qualitative Features of the Conductivity Profile...... 50 2.7.2 Analytic Estimate of Σ...... 51 2.8 Collision Rates...... 53 2.8.1 General Collision Parameters...... 53 2.8.2 Ion/Electron-Neutral Non-resonant Collisions...... 54 2.8.3 Resonant Charge Exchange...... 55 2.8.4 Coulumb Collisions...... 56 2.8.5 Location of Peaks in the Conductivity...... 57 2.9 Summary...... 60
3 Electrical Conductance of the Upper Atmosphere of Hot Jupiter HD 189733b 62 3.1 Summary...... 62 3.2 Introduction...... 63 3.3 The Model...... 66 3.3.1 Determining the Pressure Limits...... 67 3.3.2 Basic Atmosphere Structure...... 69 3.3.3 Differential Equations...... 76 3.3.4 Balance of Production and Loss Processes...... 77 3.3.5 Charge Neutrality...... 80 3.3.6 Thermal Balance...... 81 3.4 Results...... 90 3.4.1 Density and Temperature Profiles...... 90 3.4.2 Conductivity Profile...... 95 3.4.3 Magnetic Field Effect on Conductivity...... 98 3.4.4 Joule Heating in the Upper Atmosphere of HD 189733b.... 100 3.5 Discussion...... 102 3.5.1 Effect of Molecules on the Conductance...... 102 3.5.2 Relative Size of the Pederson Conductance...... 103 3.5.3 Significance of Joule Heating...... 105 3.5.4 Location of Joule Heating in the Atmosphere...... 108 3.6 Conclusions...... 109
4 Tidal Disruption of Moons in Star-Planet-Moon Systems 111 4.1 Summary...... 111 4.2 Introduction...... 112 ix
4.3 Tidal Evolution of Two-Body Systems...... 117 4.4 The Star-Planet-Moon System...... 119 4.5 Equations of Motion...... 120 4.6 Torque Equilibrium...... 125 4.7 Solutions at First Order...... 127 4.7.1 Ignoring Perturbations of the Star...... 127 4.7.2 Ignoring Tidal Friction in the Planet...... 128 4.7.3 First Order Gravitational Perturbations...... 129 4.7.4 First Order Tidal Friction Perturbations...... 131 4.7.5 Full First Order Perturbations...... 133 4.8 Equation for a at Second Order...... 135 4.8.1 Moon Migration Timescale...... 138 4.8.2 Comparison to Other Studies...... 139 4.9 Applications and Discussion...... 140 4.9.1 Migration of Moons around hot Jupiters...... 141 4.9.2 Migration of Moons around Mercury and Venus...... 143 4.9.3 Scaling of Migration Time with Tidal Q...... 145 4.9.4 Constraining Maximum Moon Sizes...... 147 4.10 Conclusion...... 149
5 Conclusions 150 5.1 Summary...... 150 5.2 Future Applications...... 151 List of Figures
1.1 Planetary mass (in Jupiter masses) vs semi-major axis (in astronomical units – 1 AU is the distance from the Earth to the Sun) for all of the confirmed exoplanets as of Nov 2020. Colors designate the exoplanet detection method used for each discovery as follows: red – transit; blue – radial velocity; green – microlensing; yellow – direct imaging. The images designate the location that each Solar System planet would fall on the plot in order from left to right of: Mercury, Venus, Earth, Mars, Jupiter, Saturn, Uranus, and Neptune. The plot is from exoplanets.org (Han et al. 2014)...... 4 1.2 A cartoon of the four major exoplanet detection methods (N. Bathalla – NASA Ames)...... 8 1.3 Face-on (left) and side (right) view of an exoplanet transiting its host star with an atmosphere. The face-on view shows the core of the planet in black, which blocks all the light from the star and a surrounding atmosphere which allows some transmission. The side view shows the effects that different impact parameters, b, have on the transmission of star light through the atmosphere. The more atmosphere that the light passes through, the more it gets absorbed...... 14 1.4 Observed transmission spectrum of the Na I D doublet lines in the upper atmosphere of HD 189733b from Wyttenbach et al.(2015). Top panel: out of transit spectrum from the host star normalized to the maximum flux observed in the wavelength regime noted by the x-axis. A ‘telluric’ spectrum expresses the absorption from the Earth’s atmo- sphere to identify (e.g.) water lines. Bottom panel: the transmission spectrum as a ratio of the flux observed in transit versus out of transit. The red lines show a Gaussian fit of the spectral bins, which clearly show absorption at the wavelengths expected from the Na I D doublet. 14
x xi
1.5 The left side of this plot shows the temperature of the Earth’s atmo- sphere as a function of altitude. Dashed lines designate the different regimes of the atmosphere. The right side of the plot shows where the ionosphere of the Earth is in relation to other regions of the atmosphere by plotting electron density as a function of height. The peak in the electron density designates the location of the ionosphere. Image is adapted from: https://en.wikipedia.org/wiki/Ionosphere...... 19 1.6 The shape of the magnetic field around the Earth from Dungey(1961). 21 1.7 Tides on a body. Each line represents the differential gravitational force experienced at that location. Panel (a) is for a two body system without any tidal friction. Points A denote where a bystander on the planet would experience high tide and points B denote where a by- stander would experience low tide. Subfigure (b) highlights the effects that tidal friction and tidal lag would have on the differential gravita- tional forces. The deviation from a central line connecting the center of the two bodies causes a tidal torque which drives systems towards tidal locking. Image reference: https://en.wikipedia.org/wiki/Tide.. 24 1.8 Lagrange points (left) and Roche limit (right) for a two body system. On the left, lines denote points of gravitational equipotential. Point L1 denotes the location where the gravity of the star, dominates the gravity of the planet at the Hill radius. The right subfigure exaggerates the tidal effects from a planet on a moon, as a moon migrates inwards until it is fully disrupted at the Roche limit. Image credits: (left) NASA, (right) Astronoo.com...... 26 1.9 A cartoon of the star-planet-moon model that I consider in Chapter4. Each body in the system orbits in nearly circular orbits where the only deviations are due to the mutual gravitational perturbations between all the bodies in the system...... 28
2.1 Examples of the Northern Lights (Aurora Borealis) and Southern Lights (Aurora Australis) as pictured from different locations on the Earth. Photo credit: https://en.wikipedia.org/wiki/Aurora...... 33 2.2 A schematic of two closure paths for the field-aligned current, shown in red as presented in Ganushkina et al.(2018). The left most path con- nects directly to the magnetopause while the right most path connects to the solar wind through the nightside plasma sheet...... 35 2.3 A diagram of the three vectors of current flow in an ionosphere with a given B and E field. The Pederson (σP ), Hall (σH ), and parallel (σk) conductivities are derived in Chapter 2.6...... 46 2.4 Diagram of a simple collision between a large particle and a background gas of constant volume of smaller particles (Schunk & Nagy 2004).. 54 xii
3.1 Flux density of stellar spectrum vs wavelength...... 73 3.2 Number density versus pressure. The top of this plot is the top of the atmosphere. Except when noted by the legend, solid lines represent neutral species and dashed lines represent ionized species. All other species are specifically represented in the legend with electrons denoted by ‘e’ and protons denoted by ‘p.’...... 92 3.3 (a) Temperature vs pressure as calculated by the thermal balance equa- tion until 10 µbar below which point the temperature is prescribed to be that from Huitson et al.(2012). (b) Radius vs pressure as calculated by hydrostatic balance (Equation 3.9)...... 93 3.4 The density fraction of ions, neutral atoms, and molecules as a func- tion of pressure. The density fraction is calculated by summing the densities of each particle of the given species and then dividing that value by the total number density of all species in the atmosphere. The transition from ionized to atomic to molecular becomes more clear in this representation than compared to Figure 3.2...... 93 3.5 The major heating and cooling rates in the atmosphere plotted against pressure until the prescribed temperature profile from Huitson et al. (2012) takes over at 10−5 bar. All cooling rates are expressed as dashed or dotted lines and all heating rates are expressed as solid lines. The dashed lines represent a sum of the line cooling for the specified species. Rates that are included in the model but do not significantly contribute to setting thermal balance are not plotted. A balance of the heating and cooling rates sets the temperature profile in Figure 3.3...... 94 3.6 Pederson, Hall, and parallel conductivities as a function of pressure.. 96 3.7 Individual contributions to the Pederson conductivity for electrons, protons, and ions. Each line represents the Pederson conductivity solely from the species denoted in the legend. The solid black line is a sum of all individual contributions...... 97 3.8 Individual contributions to the Hall conductivity for electrons, protons, and ions. Each line represents the Hall conductivity solely from the species denoted in the legend. Note that because the Hall conductivity is dependent on the charge of the species, the absolute value of the electron contribution is represented by the dotted line. The black line is a sum of all individual contributions...... 97 3.9 Pederson conductivity versus pressure as a function of magnetic field strength, B, given in units of Gauss...... 99 3.10 A running height integral of the Pederson conductivities in Figure 3.9 versus pressure. The final number that each line approaches at the highest pressures can be thought of as the conductance in the atmo- sphere for the given B field...... 99 xiii
3.11 Pederson conductivity vs. pressure for models with and without molecules. The solid lines show the total Pederson conductivity (black) and CII contribution (blue) to the Pederson conductivity for our standard at- mosphere model. The dashed lines show the total Pederson conduc- tivity (black) and CII contribution (blue) to the Pederson conductivity for the same atmosphere model without molecules included. Since car- bon does not get tied into CO for the dashed lines model it is able to generate an anomalous peak in the Pederson conductivity at P ≈ 10−6 bar...... 103 3.12 A comparison of the total heating in the atmosphere (shown in red) over the pressure regions where we solve the thermal balance equa- tion compared to the Joule Heating. We calculate Joule Heating from Equation 3.34 and from the Pederson conductance values at different magnetic field sizes as given in Table 3.12. The black lines represent the Joule heating rate calculated for each magnetic field tested.... 107
4.1 A cartoon of the migration and disruption of a theoretical moon in a hot Jupiter system. As show in (a) the planet and moon will be in circular orbits about the host star. In (b) tidal friction for forced eccentricity will cause inward migration of the moon. In (c) the moon will migrate within the Roche limit and become tidally disrupted. Finally in (d), the ring system caused by the tidal disruption of the moon will dissipate and no evidence of the moon will remain...... 115 4.2 Hut(1981) model for the tides raised on a primary object with mass M by an outside perturber with mass m. Hut(1981) approximates the tides raised by the companion as point masses, µ, to compute the perturbing tidal force on the companion with separation r, moving around the companion as θ˙...... 118 4.3 Geometry from the reference point of the central planet. Eccentricities are highly exaggerated to highlight the differences in longitude of peri- centers ($) between the two objects. The definitions of each variable are highlighted in Section 4.4...... 121 4.4 Migration of a hypothetical moon orbiting HD 189733b, HD 209458b, Mercury, and Venus. On the x-axis is time since the beginning of the integration in years and on the y-axis is the semi-major axis of the moon (a1) divided by the Hill radius of the individual system. Each system is integrated until the semi-major axis becomes less than the Roche limit for tidal disruption at which point a moon would theoret- ically be destroyed. If the Roche limit is within the planetary radius, then the moon will collide with the atmosphere of the planet, before being tidally disurpted...... 144 xiv
4.5 A plot of tidal destruction time (in years) of Earth-size moons as a function of star-planet semi-major axis (a2). The destruction times are calculated via Equation 4.70. Each black line represents a tidal Q as given by the legend. We assume the planet mass is equal to the mass of HD 189733b, the mass of the star is equal to the mass of HD 189733A, an apsidal motion constant of k = 0.1, ρ = ρ⊕, and that the moon is originally located at half the hill radius of the system. The dashed lines highlight the orbital location of HD 189733b and Mercury. 146 4.6 Here we highlight maximum moon mass to survive tidal disruption within 5 × 109 years versus semi-major axis of the star-planet system. Each line represents the mass of the planet from 0.1 MJup − 10 MJup as given in the legend. We calculate the lines by applying Equation 4.70 where we assume Q1 = 10, ρ1 = ρ⊕, k1 = 0.1, and that each moon is initially placed at half the hill radius of the system. Systems that fall within the colored regions above each line are additionally able to survive tidal disruption...... 148 List of Tables
2.1 Relevant neutral gas polarizabilities. These values were taken from Table 4.1 in Schunk & Nagy(2004)...... 55 2.2 Relevant collision frequencies for resonant ion-netural collisions. These values were taken from Table 4.5 in Schunk & Nagy(2004). All tem- peratures are in K and densities (n) in cm−3...... 56
3.1 Major collision frequencies. Note that SN refers to Schunk & Nagy (2004)...... 68 3.2 Adopted properties for HD 189733A and HD 189733b...... 69 3.3 Hydrogen Chemistry. In the text, these reactions are denoted as RH#. Photoionization and photodissociation rates are optically thin values in units of s−1. They are calculated with the fits for the cross sections given as the reference in the table and the ionizing flux from the MUS- CLES spectrum. Three body reaction rate coefficients are in units of cm6 s−1 and two body reaction rate coefficients are in units of cm3 s−1. Rate 13a/13b is calculated using 11a/11b, the Saha equation, and the law of mass action...... 79 3.4 Carbon and Oxygen Chemistry. In the text, these reactions are de- noted as RCO#. NIST refers to the National Institute of Standards and Technology Chemical Kinetics Database...... 81 3.5 Metal Chemistry: Photoionization/Recombination. In the text, these 4 reactions are denoted as RM#. Note: T4 ≡ (T/10 K). Note: The photoionization rates of all the metals are calculated from the cross section fits of Verner & Ferland(1996) and Verner et al.(1996).... 82 3.6 Alkali Metal Chemistry: Three Body Reactions. In the text, these reactions are denoted as RM#...... 83 3.7 Additional Molecule Chemistry. In the text, these reactions are de- noted as RM#...... 83 3.8 Relevant species and heats of formation (in kJ mol−1) as listed in Le Teuff et al.(2000)...... 85
xv xvi
3.9 Major Metal Cooling Lines. Notes: Clu refers to the collision frequency given by the source and CHIANTI refers to the Atomic Database for Spectroscopic Diagnostics of Astrophysical Plasmas which can be found at: https://www.chiantidatabase.org...... 86 4 3.10 Minor Metal Cooling Lines. Note: T4 ≡ 10 K...... 87 3.11 Stellar wind parameters near HD 189733b from Kavanagh et al.(2019). 101 3.12 Values of the Pederson conductance and the maximum Joule heating of the atmosphere of HD 189733b for different sized magnetic fields. 102
4.1 HD 189733 system with hypothetical Earth-like moon. The circulariza- tion time is calculated via Equation 4.1. We use k = 0.1 for each body for all systems and the tidal Q is taken from Goldreich & Soter(1966). Note that the initial semi-major axis (a0) for the moon is located at half the Hill radius of the star-planet system...... 141 4.2 HD 209458 system with hypothetical Earth-like moon...... 142 4.3 Mercury system with hypothetical Europa-like moon...... 142 4.4 Venus system with hypothetical Europa-like moon...... 142 4.5 The Roche limit and tidal disruption time (Roche time) for each of the four systems considered as defined above. Note that we assume 6 an initial nratio = 37.1α for each of the systems. Additionally note that the Roche limits for system A and B are within the planetary radius and therefore collision with the atmosphere will occur before tidal disruption...... 142 1
Chapter 1
Introduction
1.1 Exoplanets
Context: The study of extrasolar planets, or exoplanets, is central towards the core research in this thesis. Specifically, in Chapter3 an exoplanet atmosphere model is developed to understand the temperature, density, and conductivity profile for a type of exoplanet known as a hot Jupiter. Hot Jupiters are an exoplanet unlike anything observed in our Solar System and so in this introduction I describe the historical context of the detection of exoplanets. In doing this I outline the various exoplanet detection methods, which provide context to the place of hot Jupiters among other types of detected exoplanets and paint a picture of how the various exoplanet detection methods will inform my atmosphere model. The outline of the two major exoplanet detection methods (the transit method and radial velocity method) is intended to put in perspective how hot Jupiters were able to be detected and how we obtain information on their characteristics which are then implemented into my models. Additionally, an understanding of the types of exoplanet detection methods is needed to understand why we haven’t yet detected a moon orbiting an exoplanet, 2 which is important for the context of the research described in Chapter4. In this introductory Chapter, I also describe the specific hot Jupiter on which I base my atmosphere model – HD 189733b. Additionally, in Chapter4, I develop an analytic equation that can be applied to describe the orbital migration of hypothetical moons of hot Jupiters, or any planetary system that matches the criteria outlined within Chapter4. Therefore, background on the types of exoplanets observed and how they relate to objects in my models is needed to put into context the significance of the research described throughout this thesis.
1.1.1 The Discovery of Exoplanets
The year 1992 marked the first discovery of a planet outside our Solar System, when two exoplanets were detected orbiting a pulsar located 2300 light years away in the constellation of Virgo (Wolszczan & Frail 1992). Until then, all of our knowledge on the structure, dynamics, formation, and evolution of planetary bodies was limited entirely to the study of our own Solar System. The existence of exoplanets was entirely expected at the time, but this first detection immediately put our understanding of planets to the test. Exoplanets were originally only expected to be found around main-sequence stars (Wolszczan 2018) and not around the rapidly rotating and highly magnetized neutron stars that are pulsars, which form in the supernova explosion that marks the ‘death’ of a massive main-sequence star. This unexpected discovery foreshadowed the surprises that would come with subsequent exoplanet discoveries through the years. It wasn’t until 1995 that an exoplanet was finally observed orbiting a main- sequence star (Mayor & Queloz 1995). The discovery of exoplanet 51 Pegasi B, orbiting a Sun-like star located 50 light years away in the constellation of Pegasus, 3 opened the door to a new field of study that remains at the forefront of modern as- tronomical research. The detection and understanding of exoplanets is still rapidly changing, and as of November 2020 there have now been over 3000 exoplanets con- firmed and almost 2500 candidate exoplanets (exoplanets.org – Han et al. 2014).
1.1.2 The Detection of Exoplanets
Figure 1.1 shows the majority of detected exoplanets as of November 2020 in a plot of planetary minimum mass vs distance from the planet’s host star (semi-major axis). Images of each of the Solar System planets are over-plotted on Figure 1.1 to highlight how the detected exoplanets compare to planets in our Solar System. Additionally, the colors in Figure 1.1 designate the observational method that was used to detect each confirmed exoplanet with red, blue, green, and yellow designating the ‘transit,’ ‘radial velocity,’ ‘microlensing,’ and ‘direct imaging’ methods respectively. These methods are outlined in more detail below and a cartoon of the four major detection methods is shown in Figure 1.2.
The Transit Method
The transit method for exoplanet detection has been the most successful detection method to date, largely in part to the Kepler space telescope and more recently from the Transiting Exoplanet Survey Satellite (TESS). For a review on this detection method see Deeg & Alonso(2018). The transit method was first successfully em- ployed in 2000 to confirm the planetary status of HD 209458b (Henry et al. 2000; Charbonneau et al. 2000), which had already been detected via the radial velocity method. As a planet orbits around its host star and moves in the line of sight between its 4
Fig. 1.1.— Planetary mass (in Jupiter masses) vs semi-major axis (in astronomical units – 1 AU is the distance from the Earth to the Sun) for all of the confirmed exoplanets as of Nov 2020. Colors designate the exoplanet detection method used for each discovery as follows: red – transit; blue – radial velocity; green – microlensing; yellow – direct imaging. The images designate the location that each Solar System planet would fall on the plot in order from left to right of: Mercury, Venus, Earth, Mars, Jupiter, Saturn, Uranus, and Neptune. The plot is from exoplanets.org (Han et al. 2014). 5 host star and our telescope, light from the star is periodically blocked resulting in a decrease in observed brightness (Figure 1.2). The time interval between dimming events is the same as the time it takes the planet to revolve once around its host. The amount of light blocked from the star is dependent on the relative size of the planet to the star. For example, the transit by HD 209458b results in a decrease in visual brightness of about 2% (Henry et al. 2000), but most transits of confirmed exoplanets are much smaller than this since HD 209458b is relatively large in relation to its host star than most exoplanets are to their host. If the change in brightness of the star is large enough to be detected, the observation becomes a candidate exoplanet. Exoplanet candidates can then be confirmed via follow up observations using another method to rule out false positives (Brown 2003; Santerne et al. 2013) due to stellar eclipsing binaries, star spots, or background eclipsing stellar binaries. Exoplanets detected via this method are also typically confirmed via follow-up observations that employ the radial velocity method. By detecting an exoplanet via the transit method one can gain insights on many exoplanet characteristics such as orbital period, semi-major axis (average distance from the star), and planetary radius. Additionally, many stars can be observed si- multaneously to look for these periodic dips in stellar brightness as is done with Kepler and TESS. A drawback of the method is that planetary orbits need to be aligned with the observer’s line of sight and large enough to block a sizable amount of star light for an observable transit to occur. Additionally, rates of false positives are much higher than any other method and for some surveys can be as high as 35% (Santerne et al. 2012). For this reason, the transit method typically results in the detection of planets that are physically large in radius, very close to their host star, or a combination of the two as can be seen by the red dots in Figure 1.1. Despite 6 this, the transit method has been essential in our ability to detect large amounts of exoplanets, which has resulted in a better characterization and understanding of the types of exoplanets that exist in our Galaxy.
The Radial Velocity Method
Until the launch of Kepler in 2009, the radial velocity method for exoplanet discovery was the most productive method for detecting and confirming exoplanets. In fact, the first exoplanet detected around a Sun-like star, 51 Pegasi B, was detected with this method (Mayor & Queloz 1995). For an in depth review of the radial velocity method see Wright(2018). As an exoplanet orbits around its host star, it will gravitationally accelerate the star towards and away from the observer. This acceleration will shift the relative velocity of the star towards and away from the observer, and as a result the light emitted from the star will be Doppler-shifted to higher and lower frequencies (shorter and longer wavelengths). A common daily life example of the Doppler shift phenomenon occurs in the change of the pitch of an ambulance as it moves towards and away from a listener. See Figure 1.2 for a cartoon of this effect. The size of the shift in the stellar spectral features give insights into the planet’s minimum mass, orbital period, and semi-major axis. Only the minimum mass of a planet can be determined via this method, as the size of the Doppler shift depends on the line of sight velocity of the star, which is dependent on relative inclination of the star-planet orbit to the observer. More massive planets will gravitationally perturb their star more, resulting in a larger Doppler shift. For this reason, the radial velocity method tends to find exoplanets that are more massive than the types of exoplanets detected via the transit method. Some advantages to this method are that physical shifts of light from the Doppler 7 effect are independent of distance to the planetary system and the exoplanet orbit does not need to be precisely aligned with the observer as is required with the transit method. This results in planets being able to be detected further from their star than with the transit method. Despite the physical shift of light being distance independent, statistically significant detections do require a high signal-to-noise ratio of the light because of the many corrections that need to occur to constrain the planetary gravitational perturbation. Because of this, exoplanets detected via this method are somewhat limited in physical distance from the Earth. While exoplanets can be detected in large surveys using the radial velocity method (e.g. Troup et al. 2016 with APOGEE), comparatively larger surveys such as Kepler and TESS have detected many more exoplanets. For this reason, the radial velocity method has been used more frequently to confirm the detection of exoplanets via the transit method and to gain more insights in their planetary characteristics.
Microlensing, Direct Imaging, and Timing
Microlensing, direct imaging and timing variation round out the other major meth- ods for detecting exoplanets. Their relevance towards the exoplanets examined in this thesis is outside the scope of this work and so I will not go into more detail on them other than listing them for completion of the major detection methods. Their immediate significance can be highlighted by the fact that the first exoplanet discov- ered was detected via pulsar timing variations (Wolszczan & Frail 1992) and some of these methods are more capable of detecting planets orbiting very far from their host star as seen by the orange and green points in Figure 1.1. For a review of these exoplanet detection methods see Batista(2018) for microlensing, Pueyo(2018) for direct imaging, and Kramer(2018), Hermes(2018) and Agol & Fabrycky(2018) for 8
Fig. 1.2.— A cartoon of the four major exoplanet detection methods (N. Bathalla – NASA Ames). timing variations. Examination of Figure 1.1 highlights that the majority of exoplanets detected are unlike the planets in our Solar System. They can be closer to their host than Mercury is to our Sun, frequently exist in the mass range between the Earth and Neptune, can be more massive than Jupiter, and can have large masses and be very close to their host star. This does not indicate that our Solar System is inherently unique, but rather shows the selection biases of each of the exoplanet detection methods. These methods inherently detect very massive planets or planets very close to their host star which are unlike any of the planets in our Solar System. In order to get a full census of the types of exoplanets in our Galaxy, we would need to bridge the gaps in our observational methods to detect planets of all sizes and orbital distances. So far, the primary methods we have that can bridge this gap are through microlensing and direct 9 imaging. Further bridging this gap will be a very difficult task due to the observational limitations we currently face, but will be vital in fully characterizing planets in our Galaxy and the place that our Solar System stands in our characterization. After all, the inspiration of much Astronomical research, and especially exoplanet research, is in determining our place in the Universe. From the planets we are able to detect though, there are certainly objects unlike anything we see in our Solar System. For example, in Figure 1.1 the large amount of planets between the masses of Earth and Uranus/Neptune are typically referred to as ‘Super Earths’ or ‘Mini Neptunes.’ Most relevant to the scope of this thesis though are the planets located in the upper left quadrant of Figure 1.1. These planets are about the mass of Jupiter, but located at orbital distances much closer to their host star than Jupiter is to the Sun, making their temperatures much hotter than those observed in Jupiter. Therefore, these objects are often referred to as ‘hot Jupiters.’ These object are discussed in more detail in the following section. For an overview on all of the types of exoplanets that have been detected so far see Winn & Fabrycky (2015).
1.1.3 Hot Jupiters
The detection of Jupiter sized exoplanets orbiting very close to their host star was entirely unexpected and put into question our understanding of the origins, formation, and evolution of planets and stellar system as a whole (see Dawson & Johnson(2018) for a review on the origins of hot Jupiters). As outlined in Dawson & Johnson(2018) and references therein, a leading theory of planet formation is that planets form in the disk of gas and dust surrounding a new star. Rocky planets (like Mercury, Venus, Earth, and Mars) can form in the inner region of the disk and gas giants 10
(like Jupiter, Saturn, Uranus, and Neptune) formed in the outer region of the disk (beyond the ‘snow line’) – where the presence of more solids allowed larger cores to form which could then gravitationally attract more gases and form large atmospheres. Once the planets were assembled they would maintain their orbits about their host star where they were formed. This planetary formation model was based on the only planets we knew of at the time – those located in our own Solar System. Historically, our understanding of our place in the Universe has been upended time and time again. Earth is not the center of our Solar System. The Sun is not at the center of our Galaxy. Our Galaxy is not the only galaxy in the Universe. It turns out our understanding of exoplanets is no different – our Solar System is not the standard for observed planetary systems. For this work, and following the approach of Dawson & Johnson(2018), hot Jupiters are defined as planets with masses greater than or equal to 0.25 Jupiter masses and orbital periods of less than 10 days. The existence of these objects imme- diately brought up a fundamental question: How did they get there? There are three leading ideas: (1) they formed near the locations from the stars that they are dis- covered, (2) they formed beyond the ‘snow line’ of their stellar system and migrated inwards due to interactions with the gas disk, or (3) they migrated inwards from inter- actions with another large planet that caused large eccentricities that settled into the close orbits we observe due to tidal friction. Regardless of their formation mechanism, hot Jupiters are one of the most studied types of exoplanets and the understanding and modeling of their atmospheres is a primary focus of this thesis. 11
HD 189733b
The focus of Chapter 3 of this thesis is on modeling the upper atmosphere of the hot Jupiter HD 189733b. It was detected in 2005 via the transit method (Bouchy et al. 2005) and is the nearest transiting hot Jupiter, at a distance of about 65 light years. Its close distance to the Earth makes it one of the most observed and studied hot Jupiters to date. HD 189733b has a semi-major axis of 0.03 AU around a K2V star and an orbital period of just 53 hours – meaning that it is likely tidally locked to its host star as discussed in Chapter3. The planet has a mass about 1.2 times the mass of Jupiter and a radius about 1.1 times the radius of Jupiter.
1.2 The Upper Atmosphere of HD 189733b
Context: In Chapter3, I develop a hot Jupiter atmosphere model to understand the electrical conductivity within the upper atmosphere of HD 189733b. The species in the model, thermal processes considered, and dynamics of the atmosphere are constrained by observations of hot Jupiter atmospheres by observers. Therefore, in this section I outline how hot Jupiter atmospheres are observed, how the species in those atmospheres are able to be detected through spectroscopy, and outline the specific species that have been observed in hot Jupiters atmospheres, with an emphasis on HD 189733b – the hot Jupiter that I model. Additionally, a hot Jupiter atmosphere model cannot be fully encompassing of every physical process that occurs in the atmosphere due to computational and observational constraints and so assumptions must be made to model different effects in the atmosphere. Therefore, in this section I outline the different types of atmosphere models that have been created to model hot Jupiters. This discussion adds context to why I choose the various methods 12 outlined in Chapter3 to specifically understand the electrical conductivity in the upper atmosphere of HD 189733b.
1.2.1 Spectroscopic Observations of Exoplanet Atmospheres
Due to their immense distance from the Earth, almost all the information on exoplanet atmospheres is gathered through transmission spectroscopy (Roberge & Seager 2018) of transiting exoplanets. As an exoplanet transits its host star, part of the light interacts with particles in the planet’s atmosphere (see Figure 1.3). Stellar light can excite electrons in atoms and molecules in the planetary atmosphere at distinct wavelengths corresponding to discrete transitions of energy between two bound states. In this process, light at these discrete wavelengths in this spectrum, when compared to the out of transit stellar spectrum will appear to be ‘missing.’ By comparing the light observed in transit to the light from the host star out of transit, we can learn which particles are blocking the light and gain insights on the species in the planetary atmosphere. Figure 1.4, for example, shows a transmission spectrum for HD 189733b. This observation shows the presence of Na D doublet lines in the transmission spectrum, which requires Na atoms to be in the upper atmosphere of HD 189733b to absorb the starlight at those wavelengths. These observations also allow us to characterize the temperature and density pro- file of exoplanets, which provides us the information needed to test atmosphere models and theories of hot Jupiters (Roberge & Seager 2018). Since many exoplanets, such as hot Jupiters, orbit their stars very close to their host star, this allows us to test atmo- sphere models for many different conditions, such as for models of Na and other metals in its atmosphere. While exoplanet atmospheres have been detected via transmission spectroscopy on a handful of hot Jupiters, including the first exoplanet atmosphere 13 detection of the Na doublet line in HD 209458b (Charbonneau et al. 2002), I focus my discussion on the observations and models of HD 189733b. Additionally, HD 209458b is the most observed hot Jupiter and so I highlight the species detected in its upper atmosphere as well.
1.2.2 Chemical Species in the Atmosphere of Hot Jupiters
The atmosphere of HD 189733b was first detected from the ground by Redfield et al. (2008) from the Na D doublet resonance lines using the High Resolution Spectro- graph on the Hobby-Eberly Telescope (HET). This detection was a breakthrough at the time as it was the first detection of an exoplanet atmosphere from a ground-based observatory, which is typically difficult due to the seeing of the Earth’s atmosphere. HD 189733b is the closest hot Jupiter observed to the Earth, making this ground- based spectroscopic observation possible. Na has subsequently been detected in follow up observations by many teams including Jensen et al.(2011); Czesla et al.(2015); Cauley et al.(2016); Barnes et al.(2016), with the Space Telescope Imaging Spectro- graph (STIS) on board the Hubble Space Telescope (Huitson et al. 2012), and with the HARPS echelle spectrograph on the ESO 3.6 meter telescope (Wyttenbach et al. 2015). Additionally, the strong Hα line has also been detected with the HET by Jensen et al.(2012) and with the High Resolution Echelle Spectrometer (HiRES) on the Keck telescope in Hawaii (Cauley et al. 2015, 2016) – where they also reported observing Hβ and Hγ lines. Lyα absorption has been reported by Lecavelier Des Etangs et al. (2010) using the low-resolution spectrograph on board Hubble (STIS). Observations of metals in the atmosphere of HD 189733b have been detected as well including
OI (Ben-Jaffel & Ballester 2013), KI (Pont et al. 2013), MgI (Cauley et al. 2016), 14
Fig. 1.3.— Face-on (left) and side (right) view of an exoplanet transiting its host star with an atmosphere. The face-on view shows the core of the planet in black, which blocks all the light from the star and a surrounding atmosphere which allows some transmission. The side view shows the effects that different impact parameters, b, have on the transmission of star light through the atmosphere. The more atmosphere that the light passes through, the more it gets absorbed.
Fig. 1.4.— Observed transmission spectrum of the Na I D doublet lines in the up- per atmosphere of HD 189733b from Wyttenbach et al.(2015). Top panel: out of transit spectrum from the host star normalized to the maximum flux observed in the wavelength regime noted by the x-axis. A ‘telluric’ spectrum expresses the absorp- tion from the Earth’s atmosphere to identify (e.g.) water lines. Bottom panel: the transmission spectrum as a ratio of the flux observed in transit versus out of transit. The red lines show a Gaussian fit of the spectral bins, which clearly show absorption at the wavelengths expected from the Na I D doublet. 15 and a minimal detection of CaII (Barnes et al. 2016; Czesla et al. 2015), and HeI (Salz et al. 2018) lower in the atmosphere. In addition to atoms, molecules have also been detected including H2O(Birkby et al. 2013), CO, CO2, and a minimal detection of CH4 (Swain et al. 2008), the latter three suggesting the presence of C in the atmosphere. As for HD 209458b, there have been observations of all the atoms and molecules listed above in addition to the presence of CI and OI in its extended upper atmosphere
(Vidal-Madjar et al. 2004), more significant MgI detections (Vidal-Madjar et al. 2013) and recent FeII observations (Cubillos et al. 2020). Combined, these observations of HD 189733b and HD 209458b define the major species that must be included in atmosphere models of hot Jupiters. I use these observations to motivate the species included in my model of HD 189733b, as discussed in more detail in Chapter3.
1.2.3 Upper Atmosphere Models of Hot Jupiters
Hot Jupiter atmosphere models have typically been created to investigate their phys- ical and chemical processes and to explain the observed spectral features as discussed above. For example, models for the upper atmosphere of HD 209548b have been created to study atmospheric escape due to the intense heating from its host star (Yelle 2004; Garc´ıaMu˜noz 2007; Murray-Clay et al. 2009; Koskinen et al. 2013; Salz et al. 2016) and caused by Roche Lobe overflow (Gu et al. 2003; Li et al. 2010; Lai et al. 2010; Ehrenreich & D´esert 2011). Additionally, hydrodynamic models have been created to study the interaction of colliding planetary and stellar winds (Trem- blin & Chiang 2013; Christie et al. 2016) and magnetohydrodynamic (MHD) models have been created to further investigate this effect when considering the planetary magnetic field (Trammell et al. 2011, 2014). Comparisons to observations for HD 16
209458b of H Lyα,OI, SiIII, NaI, and KI lines have been done by (e.g.) Fortney et al. (2003); Koskinen et al.(2013); Lavvas et al.(2014) and for HD 189733b by Huang et al.(2017). Atmospheric species in these models are often calculated assuming local thermo- dynamic equilibrium (LTE) (Burrows & Sharp 1999; Seager et al. 2000; Sharp & Burrows 2007; Fortney et al. 2008; Visscher et al. 2010), slight deviations from LTE (Lodders & Fegley 2002; Fortney et al. 2003, 2006; Visscher et al. 2006; Showman et al. 2009; Visscher et al. 2010), photo-chemical models (Yelle 2004; Garc´ıaMu˜noz 2007; Koskinen et al. 2013; Lavvas et al. 2014; Huang et al. 2017), or a combination of the two (Moses et al. 2011). The details of my model for the upper atmosphere of HD 189733b are discussed in detail in Chapter3. There, I develop a model that transitions from the ionized region of the upper atmosphere, through the region dom- inated by atoms, and to deeper pressures where molecules begin to dominate the density profile. I consider photo-chemical processes to form and destroy species in the atmosphere and model thermal balance to determine the temperature profile. The resulting density profile is applied towards understanding the electrical conductivity in the upper atmosphere of HD 189733b and the implications that the conductivity profile can have on changing the thermal balance and density profile over this region. Applications of the conductivity are often dependent on the planet’s magnetic field and interactions with a stellar wind as outline below and in Chapter2.
1.3 Planetary Magnetospheres and Ionospheres
Context: In Chapter3, I apply a hot Jupiter atmosphere model to understand the electrical conductivity throughout the upper atmosphere and ionosphere of HD 189733b. Energy can be transferred to an atmosphere via collisions of charged parti- 17 cles with neutrals, which is quantified by the conductivity of the atmosphere. In that Chapter, I examine the significance of the heating from these collisions (often referred to as Joule heating) towards impacting the thermal balance of the atmosphere. Ad- ditionally, in Chapter2, I review the origin of the electrical conductivity. Therefore, understanding the flow of currents in the atmosphere is central towards understand- ing Chapter2 and Chapter3 of this thesis. The current at the top of the atmosphere is set by the interaction of charged particles with a planetary magnetic field at the magnetopause. Therefore, in this section, I introduce the concept of planetary mag- netospheres and ionospheres. I also introduce the role that the charged particles in the ‘solar wind’ plays in the interaction with the planetary magnetosphere and how through ‘magnetic reconnection’ charged particles from a star are able to enter the atmospheres of planetary systems.
1.3.1 Overview
Planetary magnetospheres and ionospheres are vital in setting the current that can flow through a planetary atmosphere and the atmospheric conductivity research out- lined in Chapter3 of this thesis. The magnetosphere of a planet is defined as the region surrounding a planet where its magnetic field pressure dominates the gas pres- sure (Kivelson & Russell 1995). The ionosphere is the relatively thin region of the atmosphere (for the Earth) of maximum charge density (Kivelson & Russell 1995). All planets that have a magnetic field, have a magnetosphere, and all planets that have a neutral-gas atmospheric envelope, have an ionosphere (Kivelson & Russell 1995). Figure 1.5 shows the location of Earth’s ionosphere in relation to other regions of the atmosphere as well as the electron density profile in the atmosphere. For more detailed reviews of magnetospheres and ionsopheres see Kivelson & Russell(1995) or 18
Schunk & Nagy(2004). The source of electrons and ions in the upper atmosphere is photoionization of neutral particles (Schunk & Nagy 2004). Ions can additionally undergo other chemical processes, such as recombination, chemical reactions with neutrals, diffusion deeper into the atmosphere, or can be transported by zonal winds or turbulence in the atmosphere to set the full profile of plasma in this region. In the ionosphere, electrons can no longer freely flow along planetary magnetic field lines due to collisions with ions and neutral species (Schunk & Nagy 2004). These collisions create a finite conductivity through this region of the atmosphere as opposed to the much larger conductivity in regions above the ionosphere. I discuss important interactions with the magnetosphere and ionosphere below and how they relate towards applications of the conductivity profile of hot Jupiter atmospheres.
1.3.2 The Solar Wind and Magnetic Reconnection
In the case of the Earth, charged particles can enter the Earth’s atmosphere along planetary magnetic field lines and interact with particles in the Earth’s ionosphere (Frahm et al. 1997; Galand & Richmond 2001; Fox et al. 2008). The source of these charged particles is the Sun and they emanate away from the Sun in what is called the solar wind. High temperatures in the outer atmosphere, or corona, of the Sun (∼ 106 K) allow charged particles to regularly form and continuously outflow radially from the Sun along its magnetic field lines (Kivelson & Russell 1995). Through a process known as ‘magnetic reconnection,’ solar magnetic field lines can connect to the Earth’s field lines, creating a source of charged particles which impact the Earth’s atmosphere (see Yamada et al. 2010 for a more recent review of this process). 19
Fig. 1.5.— The left side of this plot shows the temperature of the Earth’s atmo- sphere as a function of altitude. Dashed lines designate the different regimes of the atmosphere. The right side of the plot shows where the ionosphere of the Earth is in relation to other regions of the atmosphere by plotting electron density as a function of height. The peak in the electron density designates the location of the ionosphere. Image is adapted from: https://en.wikipedia.org/wiki/Ionosphere 20
Dungey(1961) first applied magnetic reconnection in the context of the Earth-Sun system. Without reconnection, the Earth’s magnetic field would be nearly a dipole and an entirely closed system, meaning that no charged particles would be able to freely enter into the Earth’s atmosphere. In reality, solar wind plasma interacts with the Earth’s magnetosphere and compresses the magnetic field on the day side of the Earth and produces an extended tail on the night side as shown in Figure 1.6(Dungey 1961). The boundary region where the magnetic field strength of the solar wind and the Earth is equal, is called the magnetopause. Across the magnetopause, charged particles transition from a weak magnetic field and dense/compressed region to a region of low density that is strongly magnetized in the Earth dominated regime. From this interaction of the solar wind impacting the Earth’s magnetic field at the magnetopause, the Earth’s magnetic field connects to the solar wind flow on the day side on magnetic field lines originating from very near the poles, and magnetic reconnection occurs. In this setup, plasma originating from the Sun can penetrate into the Earth’s atmosphere and interacts mainly through collisions in the ionosphere. Central towards dictating the flow of these charged particles in the ionosphere is the atmospheric conductivity. I outline the conductivity in regards to my research and detail the calculations that go into calculating atmospheric conductivity in Chapter 2.
1.4 The Dynamics of Star-Planet-Moon Systems
Context: In Chapter4, I examine the orbital migration and tidal disruption of moons in star-planet-moon systems. The research in Chapter4 was motivated by the lack of moons orbiting Mercury and Venus and the lack of detected moons around exoplanets, and so I also briefly outline previous research that has been done to explain 21
Fig. 1.6.— The shape of the magnetic field around the Earth from Dungey(1961). these observations here. Specifically, In Chapter4, I derive an analytic equation for the change in the semi-major axis of a moon in a star-planet-moon system where all bodies in the system have zero free eccentricity (meaning they orbit in nearly circular orbits). Gravitational perturbations from each body in the system, however, cause slight deviations from zero eccentricity for each body (known as the forced eccentricity) that allows tidal friction to never cease, even in the case of circularized orbits and synchronous rotation. This tidal friction effect transfers energy from the orbital motion of the bodies into heat, causing orbital migration. Therefore, in this section, I discuss the history of Solar System dynamics as it relates to my model and describe the basic physics of gravitational perturbations and tidal friction in planetary systems. I also outline the types of planetary systems to which my model will be applicable. 22 1.4.1 Introduction
The study of Solar System dynamics, or the motion of bodies due to the influence of gravitational forces, has a detailed history dating all the way back to the earliest of Astronomers (Murray & Dermott 1999). Through the pursuit of trying to understand the regular motion of objects in the night sky, Astronomers have derived laws to better explain this motion and to make order out of the movement of celestial bodies. These early breakthroughs include some of the most fundamental results in Astronomy including Kepler’s laws of planetary motion and Newton’s universal law of gravitation. A more ‘modern’ breakthrough includes the prediction of Neptune by Urbain Le Verrier in 1846 to explain the irregularities in Uranus’ orbit caused by an eternal gravitational perturber (Le Verrier 1846; for a historical overview of this prediction, see Krajnovi´c 2016). While the problems that those currently working on Solar System dynamics may not at first appear as ‘grand’ as the works of the early pioneers, significant theories and predictions are still being developed to this day as there are still many outstanding questions about the dynamics of our Solar System. For example, see the prediction of ‘Planet Nine,’ a large Neptune sized planet orbiting beyond the Kuiper Belt of our Solar System (Batygin & Brown 2016). With the advancement of computer technol- ogy, dynamical simulations of Solar System bodies have provided many answers to these questions. One of those questions that I take an approach to better understand, is on the lack of moons orbiting the rocky planets of the inner Solar System which is similar in nature to the lack of moons detected orbiting any exoplanets. 23 1.4.2 The Lack of Moons in the Inner Solar System
The gas giants of our Solar System all have many moons orbiting them compared to the planets of the inner Solar System. As of Nov 2020, Jupiter has 79 detected moons, Saturn has 82, Uranus has 27, and Neptune has 141. On the other hand, Earth only has a single very large moon, Mars has two small moons in Phobos and Deimos, and Mercury and Venus both have no moons. While this phenomenon was originally believed to be a result of the formation of our Solar System (Canup & Ward 2002), there have been theories that the lack of moons in the inner Solar System is a more recent phenomenon caused by tidal friction (Burns 1973; Ward & Reid 1973; Burns 1977). A familiar example of the effects that tidal friction can have on a moon’s orbit can be found in the Earth-Moon system. Thanks to reflectors placed on the surface of the Moon by the Apollo astronauts (Bills & Ray 1999), we know that the Moon is currently drifting away from the Earth at a rate of 3.8 cm (or 1.5 inches) per year. This migration is caused by tidal friction from tides raised on the Earth by the Moon (Bills & Ray 1999). Tidal friction is described in more physical terms below.
1.4.3 Tidal Friction in Two Body Systems
In daily life, tides are most familiar from the rise and fall of the ocean twice per day. Tides on the Earth are primarily caused by the differential acceleration felt across the surface of the Earth due to the Moon’s gravity. This differential acceleration is expressed by the vectors in Figure 1.7. In this setup, as the Moon orbits the Earth, the differential acceleration will be largest, and outward from the Earth, on both sides of the planet in the line directly between the center of the Earth and the
1https://solarsystem.nasa.gov/moons/overview/ 24
Moon (denoted by points A in Figure 1.7) – causing what we know as high tide. The differential acceleration is smallest, and inward on the Earth, on the points farthest away from the line connecting the Earth and the Moon (denoted by points B in Figure 1.7) – causing what we know as low tide. As the Moon orbits the Earth, the tidal ‘bulge’ of high tide will follow the Moon causing the rise and fall of water levels on the Earth twice a day. Everything on the Earth experiences this effect, including the solid body of the planet, but it is most easily observed in the rise and fall of water levels on the surface.
(a) (b)
Fig. 1.7.— Tides on a body. Each line represents the differential gravitational force experienced at that location. Panel (a) is for a two body system without any tidal friction. Points A denote where a bystander on the planet would experience high tide and points B denote where a bystander would experience low tide. Subfigure (b) highlights the effects that tidal friction and tidal lag would have on the differential gravitational forces. The deviation from a central line connecting the center of the two bodies causes a tidal torque which drives systems towards tidal locking. Image reference: https://en.wikipedia.org/wiki/Tide
In reality, the tides we experience on the Earth are much more complex than this basic picture. The Earth also experiences significant tides from the differential acceleration felt across its surface due to gravity from the Sun. When the tidal bulges due to the Moon and the Sun line up (during full and new moons), the sea level changes are their most extreme and referred to as ‘spring’ tides. When the tidal bulges from the Moon and Sun are exactly opposite of each other (during first or last quarter moons), the sea level changes are their smallest and referred to as ‘neap’ 25 tides. Furthermore, the time that high and low tides on the Earth occurs actually ‘lags’ behind the time that would be expected from the line connecting the Earth and the Moon in Figure 1.7. This tidal lag, in regards to the ocean, is caused by friction between water and the Earth’s surface as the water moves over the Earth’s surface during high and low tides over the course of a day. The more friction there is in the system, the greater the tidal lag. This tidal friction also inputs heat into the Earth, resulting in a loss of energy from the orbital motion and the Moon slowly migrating away from the Earth. The transfer of heat is dependent on the density of objects, so by studying this tidal lag and tidal friction we can also gain insights on the interior structure of the Earth and the Moon (Peale et al. 1979). If the Moon continues to migrate away from the Earth, it will eventually reach a lo- cation where the Sun’s gravity overcomes the Earth’s and the Moon will be ‘stripped’ from the Earth-Moon system. This is set by the ‘Hill radius’ of a two-body system (Murray & Dermott 1999) and the limit is shown by L1 in Figure 1.8. For the Earth- Moon system, the limit for tidal stripping is about four times the current distance between the Earth and the Moon. Alternatively, if the Moon was migrating inwards towards the Earth, it would eventually reach a point where the stretching of the Moon due to tides from the Earth would overcome its own self-gravitation and rip the Moon apart (see Figure 1.8). The limit where acceleration from tides overcomes self-gravitation is known as the Roche limit and in the Earth-Moon system would occur at about 2% of the current distance to the Moon. Saturn’s ring system and other ring systems of the gas giants may have been caused by destroyed moons that got too close to their planet (Canup 2010; Black & Mittal 2015; Hesselbrock & Minton 2019). Additionally, this is likely to be the fate of the Martian moon Phobos as it 26 migrates towards Mars (Yoder 1982; Black & Mittal 2015). Perhaps moons that were orbiting the inner planets during the early Solar System were stripped or destroyed due to this tidal friction effect.
Fig. 1.8.— Lagrange points (left) and Roche limit (right) for a two body system. On the left, lines denote points of gravitational equipotential. Point L1 denotes the location where the gravity of the star, dominates the gravity of the planet at the Hill radius. The right subfigure exaggerates the tidal effects from a planet on a moon, as a moon migrates inwards until it is fully disrupted at the Roche limit. Image credits: (left) NASA, (right) Astronoo.com
An important consequence of tidal lag and friction in a two-body system is that it can lead to circularization of orbits and tidal locking. Circularization is the process by which eccentric orbits evolve into circular orbits due to tidal friction. Tidal locking occurs when the gravitational forces of tides, including friction, force the rotation of an object to synchronize with its revolution about its host. The phenomenon is already seen from the surface of the Earth in that one side of the Moon always faces the Earth. If the Earth was also tidally locked to the Moon such that only one side of the Earth faced the Moon, then tides would not move over the surface of the Earth and tidal friction would stop dissipating energy. Therefore, when objects become circularized and tidally locked, one would think that tidal friction effects would completely turn 27 off. However, when a third body is introduced into the system, it can cause grav- itational perturbations that create a short-period eccentricity or non-synchronous rotation effect that can never go away and allows tidal friction to always be present. I model the effects that this tidal friction can have on three-body systems in Chapter 4 and apply that model to understand the lack of moons in the inner Solar System and lack of moons detected orbiting exoplanets.
1.4.4 The Three-Body Problem
The generalized three-body problem is a subject that has a significant amounts of literature on it (e.g. Szebehely(1967) and Marchal(1990)). I, however, seek to ana- lytically derive the migration of moons in nearly circular orbits about their host planet due to the tidal friction effect described above and employ the circular, three body problem to create my model. The circular, three-body problem is defined (e.g. Mur- ray & Dermott 1999) for a system where two of the bodies move in circular, coplanar orbits about a center of mass and the third body is too small to significantly affect the motion of the other bodies. While real bodies in the Solar System are often not in circular and coplanar orbits, the restricted three-body problem can be applicable in many unique situations, such as a star-planet-moon system. If the motions of two bodies are dominated by a central body (e.g. a star), then the secondary bodies will orbit in ellipses with only small deviations due to gravitational perturbations of the other bodies in the system. These small deviations lead to the constant tidal friction effect previously discussed and can be analyzed by the ‘disturbing function’ which is discussed in more specific terms in Chapter4. 28
Fig. 1.9.— A cartoon of the star-planet-moon model that I consider in Chapter4. Each body in the system orbits in nearly circular orbits where the only deviations are due to the mutual gravitational perturbations between all the bodies in the system.
1.4.5 Model of a Star-Planet-Moon System
My model in Chapter4 describes a three-body system with a moon in a nearly circular orbit about a planet, that itself is in a nearly circular about a star. These bodies have zero relative inclination to each other and the deviations from the circular orbit are solely due to the gravitational perturbations each body feels from each other. The system is defined such that the planet is the central mass and the moon and star are the ‘inner’ and ‘outer’ bodies respectively. A cartoon of the model is shown in Figure 1.9. 29
Other Applicable Systems
The star-planet-moon model, as discussed above, can be applied to any system that matches the criteria, not just Solar System objects. In particular, we can apply this system to exoplanet systems and determine the likelihood that they can dynamically maintain any moons. To date, there have been no strong detections of moons or- biting exoplanets (exomoons), and only a single minor detection of a Neptune sized moon orbiting a hot Jupiter (Teachey & Kipping 2018). Under current observational capabilities, it has been shown that relatively large exomoons (e.g. the size of the Galilean moons or larger), could be possible to detect, but very difficult with current exoplanet detection methods (Sartoretti & Schneider 1999; Kipping 2009a,b). There- fore, there is a chance that moons in these exoplanet systems would be dynamically unstable and they would not have any moons orbiting around them similar to the inner planets in our Solar System.
1.5 Overview
In this thesis I model two problems: (1) the upper atmosphere of hot Jupiter HD 189733b with a focus on calculating its electrical conductivity profile, as well as (2) a model for the tidal evolution of a moon in a star-planet-moon system. Specifi- cally, in Chapter2, I review the equations that are applied to calculate the electrical conductivity of my hot Jupiter atmosphere model. Additionally, I outline how the magnetosphere surrounding a planet constrains the flow of current into the atmo- sphere. I also discuss how this current leads to heating of the atmosphere through the interaction of charged particles with other species in the atmosphere. In Chapter 3, I detail the hot Jupiter atmosphere model developed and apply the electrical con- 30 ductivity equations and current flow to determine the Joule heating in the atmosphere of HD 189733b. In Chapter4, I detail the model developed to explain the tidal evolu- tion of moons in three body systems that are orbiting in coplanar and circular orbits. I also apply my derived analytical expression to theoretical moons of hot Jupiters, Mercury, and Venus to try to understand the lack of moons in these systems. In Chapter5, I summarize the work throughout the thesis and discuss applications of this work for future studies. 31
Chapter 2
Electric Currents and Conductivity in Planetary Atmospheres
2.1 Context
For Solar System planets, electric currents play an important role in the dynamics and energetics of atmospheres – most notably in the ionosphere (Paschmann et al. 2002; Schunk & Nagy 2004). One type of current system induced in the atmosphere is due to the interaction of the solar wind with the planetary magnetosphere. The size of the current which can flow through the atmosphere is dependent on the collisions of charged particles with other species in the atmosphere, which is parameterized by the electrical conductivity. A consequence of the collision of charged particles in the electric current with other species in the atmosphere is through the transfer of energy by collisions, more commonly referred to as Joule heating. Through studies of the electric currents and conductivity in Solar System planets, models can be developed for more general objects and applied, for example, to exoplanet atmosphere models. In Chapter3 I develop a hot Jupiter atmosphere model for HD 189733b with the 32 purpose of determining the electrical conductivity in its upper atmosphere. Once the conductivity profile is known, I then apply parameters for the stellar wind of its host star and magnetic field strengths to set the current at the top of the atmosphere and then determine the significance of Joule heating in HD 189733b. Therefore, in this chapter, I detail how currents in planetary atmospheres are set by the stellar wind- magnetosphere interaction, outline the calculation of Joule heating in a planetary atmosphere, and review the formulae for electrical conductivity.
2.2 Historical Overview
The most beautiful example of charged particles from a stellar wind interacting with a planetary magnetic field is in the visible auroras near the polar regions of the Earth (Figure 2.1). We now know that auroral light shows are caused by solar wind particles moving along magnetic field lines connecting the solar wind to the Earth and colliding, primarily, with oxygen in the Earth’s ionosphere. When collisions occur, electrons can be excited into a higher energy state and then release energy at very specific wavelengths (primarily corresponding to green light) as they de-excite back into the ground state. This understanding, however, eluded humans for centuries. There are written notes on auroras dating back to at least 600 B.C. (Ch. 1 – Kivelson & Russell 1995 and references therein) and even possible cave drawings as early as 10,000 B.C., but it wasn’t until the 1700s when astronomers began gaining a physical understanding of the phenomenon. Edmund Halley first suggested that auroras depended on the direction of the Earth’s magnetic field and in the 1730s, Anders Celsius published findings that magnetic needles reacted strongly to auro- ral displays. The first theory that auroras were connected to particles emanating from the Sun was by H. Becquerel in 1878. In the early 1900s Kristian Birkeland 33
Fig. 2.1.— Examples of the Northern Lights (Aurora Borealis) and Southern Lights (Aurora Australis) as pictured from different locations on the Earth. Photo credit: https://en.wikipedia.org/wiki/Aurora 34 began obtaining extensive data on auroras and concluded that large electric currents flowed along magnetic field lines during auroral storms, which was later supported by mathematician Carl Størmer. Advances in photography, spectroscopy, and the radio transmitter accelerated the understanding of auroras. In 1918, Sydney Chapman and collaborators suggested that a single charged beam of plasma from the Sun caused magnetic disturbances on the Earth which laid the foundation for our modern understanding of solar wind- magnetosphere interaction. Finally, advances in rockets and space probes in the 1950s allowed direct in situ confirmations of the solar wind, the magnetosphere, and their connection to auroras – most notably done by James Van Allen. In subsequent years, astronomers began observations of auroras on other planets in the solar system and our modern physical models of stellar wind-magnetosphere interaction were devel- oped. We are currently at the early stages of developing models for the stellar wind- magnetosphere interactions of exoplanets but the historical context of auroras is im- portant towards understanding the basic picture of how charged particles form a star can interact with the magnetic field and atmosphere of a planet. In this thesis I ap- ply the understanding of stellar wind-magnetosphere interaction to set the currents that can flow in hot Jupiter atmospheres. A general overview of electric currents is discussed in more detail below.
2.3 Electric Currents Overview
An electric current is defined as the net flow of charge from one location to another. See Ganushkina et al.(2015) for an expert overview of current systems in the Earth’s magnetosphere and Ganushkina et al.(2018) for a more elementary overview on 35 currents in planetary atmospheres. Most relevant to this thesis is the flow of electrons, protons, and ions which originate in a star and move through a planetary atmosphere while interacting with the magnetic field and other species in the atmosphere. If we assume that a planet has a strong dipolar magnetic field, stellar wind particles will collide with the magnetic field, compressing it on the side facing the star (‘the dayside’) at the magnetopause and elongating it on the side away from the star (‘the nightside’) as shown in Figure 1.6. The dayside currents are most applicable to the research outlined in Chapter3 and the discussion throughout this chapter is, therefore, focused on these currents. While there are many types of dayside and nightside currents in planetary mag- netospheres as outlined in Ganushkina et al.(2015, 2018), the type of current most applicable to the Joule heating discussed in Chapter 2.4 is the current flowing par- allel to the planetary magnetic field, appropriately called the ‘field-aligned current’ or ‘Birkeland current,’ and first suggested by Birkeland(1908). A schematic of this current with respect to the magnetopause of a planet is shown in Figure 2.2.
Fig. 2.2.— A schematic of two closure paths for the field-aligned current, shown in red as presented in Ganushkina et al.(2018). The left most path connects directly to the magnetopause while the right most path connects to the solar wind through the nightside plasma sheet. 36
This field-aligned current on the left side of Figure 2.2 connects to the dayside magnetopause and is driven by the solar wind, which acts like a generator through dayside reconnection (Cowley 2000; Iijima & Potemra 1976; Siscoe et al. 1991; Stern 1983; Xu & Kivelson 1994). This circuit then connects with one end in the ionosphere and the other in the solar wind as shown by Stern(1983). This current flows into and closes in the ionosphere as described by (e.g.) Iijima & Potemra(1976) and then returns to space (Kivelson & Russell 1995; Gombosi 2004). The stellar wind interaction with the magnetosphere can then set the current that can flow in the ionosphere that can then generate Joule heating from the conductivity profile as outlined in the next section.
2.4 Joule Heating
Joule heating is the process by which an electric current, such as the field-aligned current, generates heat. A familiar example of Joule heating is in the heating of the filament of an incandescent light bulb as current passes through it. Kivelson & Ridley(2008) discuss how to calculate the Joule heating of an atmosphere given a conductance and dependent on the stellar wind-magnetosphere interaction that sets the field-aligned current through the atmosphere. The equations introduced in Kivelson & Ridley(2008) and Cohen et al.(2018) for Joule heating, are a consequence of current flows in the ionosphere of a planet and the physical mechanism behind the currents and the heating is discussed below.
2.4.1 Height Integrated Currents and Conductivity
The region where currents in the ionosphere flow is a very thin region compared the entirety of the magnetosphere (Paschmann et al. 2002). Thus, the ionosphere can, in 37
effect, be treated as a thin conducting layer in which vertical currents or electric fields can be ignored. This formalism allows us to refer to ‘height-integrated’ currents and conductivities for the ionosphere. We save the detailed discussion on height integrated conductivities for Section 2.7 but define here the height integrated values of electrical conductivity and current to be given by:
Z Σ ≡ σdz (2.1)
Z j ≡ Jdz (2.2) where σ is the electrical conductivity in the ionosphere, Σ is the height integrated conductivity which is called the conductance, J is the current density, and j is the height integrated current density of the ionosphere over altitude z. In calculating these height-integrated values, we assume that the electric and magnetic fields do not vary much over the ionosphere. The use of the height-integrated concept of currents in the ionosphere can be used to describe the closure of parallel currents (like the field-aligned current) in this region. Parallel currents that close in the ionosphere can then transfer heat to the atmosphere through Joule heating.
2.4.2 Parallel Currents and Alfv´enConductance
Parallel currents that close in the atmosphere are created as a result of Alfv´enwaves generated by the movement of planetary magnetic field lines due to the stellar wind- magnetosphere interaction introduced in Chapter 1.3. Alfv´enwaves, which are a type of magnetohydronamic (MHD) wave generated as the response of charged particles to a restoring force from tension in a magnetic field line, travel at the Alfv´envelocity 38 given by: B v = √ (2.3) A 4πρ where B is the magnetic field strength of the field line and ρ is the mass density of the charged particles. The electric field from the current can be written in terms of the magnetic field perturbation which generates the Alfv´en wave. This perturbation can be expressed as: 1 1 k B2 δE = − δv × B = z δB × eˆ (2.4) c c ω 4πρ z
Here, the δ denotes the perturbations to the various vectors as a result of the per- turbation to the magnetic field, kz and ω are the wave numbers of the Alfv´enwave propagating along magnetic field lines, ρ is the mass density of the charged particles, and e denotes the unit vector in the z direction. We can substitute in for the Alfv´en speed above to get the perturbation to the electric field to be:
v δE = ± A δB × eˆ (2.5) c z
2 2 2 since ω = k vA for the propagating wave. The perturbation to the electric current therefore only has anx ˆ andy ˆ component. We can then determine the parallel current in the z-direction from Ampere’s law as: c 2 c ∂δBy ∂δBx c Jz = (∇ × δB)z = − = ± ∇⊥ · δE⊥ (2.6) 4π 4π ∂x ∂y 4πvA
2 We can define a quantity called the ‘Alfv´enconductance’ to be ΣA ≡ c /(4πvA) and 39 therefore the current in the z direction is:
Jz = ±ΣA∇⊥ · δE⊥ (2.7) where a positive value means the current moves along +ˆz and a negative value means the current moves along −ˆz. Note that this current is the current that is ‘incident’ on the ionosphere due to the Alfv´enwave.
2.4.3 Reflection of Alfv´enWaves
At all boundaries where two plasmas have different values of Alfv´envelocity, part of the wave is reflected and part of the wave is transmitted. This is also the case for the Alfv´enwave created from the stellar wind-magnetosphere reaction when it impacts the ionosphere which has a conductance given by ΣP . Here the ‘P’ denotes the ‘Pederson’ conductance which is discussed in more detail in the following sections. The field-aligned current density in the ionosphere can be written as:
Jz = ΣP ∇⊥ · E⊥ (2.8) when the Pederson conductance is dominant. Physically, more than just the Pederson conductivity governs the flow of charged particles through the ionosphere. The Hall conductivity can also play a role as shown by Equation 3.7 in Section 3.1.1 of Paschmann et al.(2002). However, for the simplest scenario where ionospheric conductivity is uniform over the plasma region, as we are assuming is the case for the height integrated conductivities of the ionosphere, only the Pederson current contributes to closing field-aligned currents and generating Joule heating in this region. The height-integrated approximation cannot be done to fully 40
calculate the instantaneous Joule heating over all regions of the atmosphere. In this case, full conductivity, collision, and current profiles would need to be known from the top of the atmosphere through the base of the model and an MHD model would need to be carried out to get the full effect. We, therefore, choose to follow the simplest approach outlined in Section 3.4 of Paschmann et al.(2002) and discussed in Kivelson & Ridley(2008) to calculate the Joule heating in our model, which is independent of the Hall conductivity discussed in more detail in the following sections. We can now determine the reflection coefficient of the Alfv´enwave by considering the parallel current variation generated by the difference of the Pederson and Alfv´en conductance. For the Alfv´encurrent, the incident current will have a plus sign, while the reflected (upgoing) wave will has a minus sign. The sum of the two Alfv´enwaves is thus: