Arxiv:1905.06367V2 [Astro-Ph.EP] 27 Aug 2019 Spheric and Surface Properties, a Comprehensive Model of Known Stellar Systems
Total Page:16
File Type:pdf, Size:1020Kb
Draft version August 29, 2019 Typeset using LATEX twocolumn style in AASTeX63 VPLanet: The Virtual Planet Simulator Rory Barnes,1, 2 Rodrigo Luger,3, 2 Russell Deitrick,4, 2 Peter Driscoll,5, 2 Thomas R. Quinn,1, 2 David P. Fleming,1, 2 Hayden Smotherman,1 Diego V. McDonald,1 Caitlyn Wilhelm,1, 2 Rodolfo Garcia,1, 2 Patrick Barth,6 Benjamin Guyer,1 Victoria S. Meadows,1, 2 Cecilia M. Bitz,7, 2 Pramod Gupta,1, 2 Shawn D. Domagal-Goldman,8, 2 and John Armstrong9, 2 1Astronomy Department, University of Washington, Box 951580, Seattle, WA 98195, USA 2NASA Virtual Planetary Laboratory Lead Team, USA 3Center for Computational Astrophysics, Flatiron Institute, New York, NY 10010, USA 4Center for Space and Habitability, University of Bern, Gesellschaftsstrasse 6, CH-3012, Bern, Switzerland 5Department of Terrestrial Magnetism, Carnegie Institution for Science, 5241 Broad Branch Rd, Washington, DC, 20015, USA 6Max Planck Institute for Astronomy, K¨onigstuhl 17, 69117 Heidelberg, Germany 7Atmospheric Sciences Department, University of Washington, Box 951640, Seattle, WA 98195, USA 8Planetary Environments Laboratory, NASA Goddard Space Flight Center, 8800 Greenbelt Road, Greenbelt, MD 20771, USA 9Department of Physics, Weber State University, Ogden, UT 84408-2508, USA ABSTRACT We describe a software package called VPLanet that simulates fundamental aspects of planetary system evolution over Gyr timescales, with a focus on investigating habitable worlds. In this initial release, eleven physics modules are included that model internal, atmospheric, rotational, orbital, stellar, and galactic processes. Many of these modules can be coupled simultaneously to simulate the evolution of terrestrial planets, gaseous planets, and stars. The code is validated by reproducing a selection of ob- servations and past results. VPLanet is written in C and designed so that the user can choose the physics modules to apply to an individual object at runtime without recompiling, i.e., a single executable can simulate the diverse phenomena that are relevant to a wide range of planetary and stellar systems. This feature is enabled by matrices and vectors of function pointers that are dynamically allocated and populated based on user input. The speed and modularity of VPLanet enables large parameter sweeps and the versatility to add/remove physical phenomena to assess their importance. VPLanet is publicly available from a repository that contains extensive documentation, numerous examples, Python scripts for plotting and data management, and infrastructure for community input and future development. 1. INTRODUCTION backs among, interior, atmospheric, stellar, orbital, and Exoplanetary systems display a diversity of morpholo- galactic processes. Below we describe a set of physical gies, including a wide range of orbital architectures and models (called modules) that simulates these phenom- planetary densities. This heterogeneity likely leads to a ena, as well as their assumptions and limitations. We wide range of evolutionary trajectories that result in dis- validate the code, including both individual modules parate planetary properties. As astronomers and astro- and a subset of module combinations, by reproducing biologists probe these worlds to determine their atmo- key observations of the Earth, Solar System bodies, and arXiv:1905.06367v2 [astro-ph.EP] 27 Aug 2019 spheric and surface properties, a comprehensive model of known stellar systems. Where observations are lacking the physical effects that sculpt a planet can help priori- or unavailable, we reproduce a selection of previously tize targets and interpret observations. Here we describe published results. The software is open source and in- a software package called VPLanet that self-consistently cludes documentation, examples, and the opportunity 1 simulates many processes that influence the evolution for community involvement. of gaseous and terrestrial planets in a range of stellar While VPLanet is designed to model an arbitrary plane- systems. Our approach allows for coupling of simple tary system, the primary motivation for creating VPLanet models by simultaneously solving ordinary and partial is to investigate the potential habitability of exoplanets differential equations (ODEs and PDEs) and explicit functions of time to track the evolution of, and feed- 1 https://github.com/VirtualPlanetaryLaboratory/vplanet 2 with a single code that can capture feedbacks across the of planetary system evolution that can be used to simu- range of physical processes that affect a planet's evolu- late the evolution of 1M⊕, 1R⊕ terrestrial planets ∼ ∼ tion. TRAPPIST-1 (Gillon et al. 2016, 2017) is a good with approximately Earth-like material properties and example of why coupled processes are needed to under- structure, with the goal of understanding their potential stand planetary evolution. In that system, the host star for habitability. dimmed during the pre-main sequence, orbital interac- Though we are interested in Earth-like planets, the tions between planets are strong, tidal heating and ro- model is not restricted to Earth-radius and Earth-mass tational braking are significant, and stellar activity can planets; rather, the model is intended to investigate drive atmospheric mass loss. As discussed below, VPLanet a wide range of planetary systems. The phenomena combines a set of theoretical models to provide a first described above are also important for uninhabitable order approach to investigate these processes and the worlds such as GJ 1132 b (Berta-Thompson et al. 2015) feedbacks among them. and 40 Eri b (Ma et al. 2018). VPLanet can be used to We define habitable to mean a planet that supports simulate many such planets (and moons) to infer their liquid surface water. Since astrochemical and planet for- histories. mation studies find that most small planets form with The individual modules of VPLanet can simulate many significant inventories of water and bioessential elements aspects of planetary evolution, but the ability to cou- (e.g., van Dishoeck et al. 2014; Morbidelli et al. 2018), ple modules together facilitates novel investigations, as the most pertinent question for potentially habitable demonstrated in several previously-published studies. planets may be \do they still have water?" The pres- Deitrick et al.(2018a,b) combined the orbital, rota- ence of water is controlled by the planet's interior, at- tional, and climate modules to show that potentially mosphere, orbit, host star, and even the galaxy. Our habitable exoplanets can become globally glaciated if approach does not presume a priori that any one pro- their orbital and rotational properties evolve rapidly and cess dominates the evolution because numerous pro- with large amplitude. Fleming et al.(2018) showed that cesses can severely impact a planet's potential to sup- the coupled stellar-tidal evolution of tight binary stars port liquid water, even in the habitable zone (HZ; Kast- can lead to orbital evolution that ejects circumbinary ing et al. 1993; Kopparapu et al. 2013). For example, planets. Lincowski et al.(2018) combined stellar evo- tidal heating may be strong enough to drive a runaway lution and atmospheric escape to track water loss and greenhouse (Barnes et al. 2013), the pre-main sequence oxygen build-up on the TRAPPIST-1 planets. Finally, evolution of low-mass stars may desiccate (remove wa- Fleming et al.(2019) simulated the coupled stellar-tidal ter from) a planet (Luger & Barnes 2015), and stable evolution of binary stars to show how unresolved bina- resonant orbital oscillations can drive extreme eccentric- ries impact gyrochronology age estimates of stellar pop- ity cycles (Barnes et al. 2015). Furthermore, all these ulations. processes can operate in a given system, and hence a VPLanet is a fast and flexible code that combines a host rigorous model of planetary system evolution, including of semi-analytic models, which are all written in C. It habitability, should include as broad a range of physics can provide quick insight into an individual planetary as possible (Meadows & Barnes 2018). system with a single simulation (such as calculating tidal Coupled evolutionary models could be valuable tools heating or stellar evolution), or can perform parameter for interpreting the environments of newly discovered sweeps and generate ensembles of evolutionary tracks Earth-sized planets with nearly Earth-like levels of in- that can be compared to observations. Alternatively, cident stellar radiation (\instellation") orbiting other VPLanet can be combined with machine learning algo- stars. These planets may be capable of supporting liq- rithms to identify key parameters (Deitrick et al. 2018b). uid water, but their habitability is currently an open These capabilities can provide direct insight or can com- question given the plethora of physical processes at play, plement research with more sophisticated tools, such as i.e., a planet's presence in the star's habitable zone does 3-D global circulation models, which are too computa- not imply the planet supports liquid water. Moreover, tionally expensive to explore vast parameter spaces. For worlds like Proxima Centauri b (Anglada-Escud´eet al. example, VPLanet can be used to isolate the most impor- 2016) and TRAPPIST-1 c{g (Gillon et al. 2017) orbit tant phenomena (see 2) and explore parameter space, x bright enough host stars that their atmospheres will be and then more complicated models can target interest- probed with future ground- and space-based facilities ing regions of