Arxiv:2104.06479V1
Total Page:16
File Type:pdf, Size:1020Kb
Draft version April 15, 2021 Preprint typeset using LATEX style emulateapj v. 12/16/11 PRIMORDIAL GIANT PLANET OBLIQUITY DRIVEN BY A CIRCUMPLANETARY DISK Rebecca G. Martin1 and Philip J. Armitage2,3 1Department of Physics and Astronomy, University of Nevada, Las Vegas, 4505 South Maryland Parkway, Las Vegas, NV 89154, USA 2Center for Computational Astrophysics, Flatiron Institute, New York, NY 10010, USA and 3Department of Physics and Astronomy, Stony Brook University, Stony Brook, NY 11794, USA Draft version April 15, 2021 ABSTRACT Detached circumplanetary disks are unstable to tilting as a result of the stellar tidal potential. We examine how a tilted circumplanetary disk affects the evolution of the spin axis of an oblate planet. The disk is evolved using time-dependent equations for linear wave-like warp evolution, including terms representing the effect of the tidal potential and planetary oblateness. For a disk withasufficiently large mass, we find that the planet spin quickly aligns to the misaligned disk. The tilt of the planetary spin axis then increases on the same timescale as the disk. This can be an efficient mechanism for generating primordial obliquity in giant planets. We suggest that directly imaged exoplanets at large orbital radii, where the disk mass criterion is more likely to be satisfied, could have significant obliquities due to the tilt instability of their circumplanetary disks. Subject headings: accretion, accretion discs – hydrodynamics – instabilities –planets and satellites: formation – planetary systems – stars: pre-main sequence 1. INTRODUCTION generate primordial obliquity, as a byproduct of planetary in- The rotation state of planets is one of the basic observables teraction with a tilted circumplanetary disk. The physical that offers information about planet formation and planetary properties of circumplanetary disks, specifically their large system evolution. For giant planets, that open a gap in the size relative to the Hill sphere, and moderately large geomet- protoplanetary disk (Lin & Papaloizou 1986), the primordial ric thickness, favor growth of the tilt instability of disks in spin state is determined by the interaction between the grow- binary systems (Lubow 1992; Lubow & Ogilvie 2000; Mar- ing planet and its circumplanetary disk, which forms because tin et al. 2020). The tilt may be able to grow large enough the planetis much smaller than the Hill radius(Artymowicz & to excite Kodai–Lidov oscillations (Kozai 1962; Lidov 1962), Lubow 1996; Lubow et al. 1999; D’Angelo et al. 2002). The at which point eccentricity is induced and accretion becomes magnitude of the planetary angular momentum that results episodic (Martin et al. 2014, 2021). Here, we show that a cir- from circumplanetary disk accretion depends upon whether cumplanetary disk, whose tilt is being excited by an external the planet-disk interaction is hydrodynamic(Donget al. 2020) tidal potential, can couple to the planetary spin via oblateness, or magnetospheric (Koenigl 1991; Batygin 2018), butonceset and generate obliquity. The planet exerts a torque on the disk is expected to be conserved in the absence of tidal effects. In (e.g Tremaine & Davis 2014; Speedie & Zanazzi 2020) and the Solar System Jupiter and Saturn rotate rapidly, but at a an equal and opposite torque is felt by the planet. This can rate that is still only 30-40% of their nominal breakup angular change the spin axis of the planet depending upon the ratio velocity. For a small sample of very massive extrasolar plan- of the disk and planet spin angular momenta. An analogous ets, Bryan et al. (2018) inferred somewhat lower but non-zero process can realign a spinning black hole with a misaligned rotation rates. disk (Scheuer & Feiler 1996; Martin et al. 2007; Perego et al. Misalignment between the spin and orbital angular momen- 2009). We show that the torque from the planet causes the tum vectors of planets—the planetary obliquity—is common planet spin and the disk to precess around the sum of their an- in the Solar System. Jupiter’s obliquity is low (3◦), but Saturn gular momenta. They align towards each other on a relatively (27◦) and the ice giants have highly significant obliquities. short timescale. On a longer timescale, both can have increas- arXiv:2104.06479v1 [astro-ph.EP] 13 Apr 2021 Measurement of extrasolar planetary obliquities is challeng- ing tilts as a result of the tidal tilt instability, if the disk mass ing, but there are some indications of non-zero obliquity for is sufficiently large. the directly imaged planetary mass companion in the 2M0122 2. CIRCUMPLANETARY DISK AND PLANET-SPIN MODEL system (Bryan et al. 2020). Unlike the magnitude of the angu- lar momentum, the obliquity is relatively fragile to late-time We model a circumplanetary disk around a planet of mass evolutionary changes. Resonances between the precession Mp that orbits a star of mass Ms at orbital separation ap. The frequencies of planetary spin axes and orbits affect both ter- angular velocity of the planet-star system is Ωp =Ωpez where restrial and giant planets in the Solar System (Laskar & Robu- 3 the angular frequency is Ωp = G(Ms + Mp)/ap. The disk tel 1993; Ward & Hamilton 2004; Ward & Canup 2006). A q related process, where the orbital precession is driven by the has surface density Σ(r) (thatis fixed in time) anda unit vector l gravitational potential of the protoplanetary disks, can gener- to describe the tilt at each radius, d(r, t) = (ld,x, ld,y, ld,z). The ate early-time obliquity (Millholland & Batygin 2019). Giant tilt is assumed to be small so that ld,z ≈ 1 and ld,x, ld,y ≪ 1. impacts (Safronov 1966; Benz et al. 1989; Morbidelli et al. In order to explore the evolution to high inclination we would 2012), and planet-planet scattering events (Li 2021; Hong need to use hydrodynamical simulations. et al. 2021), can also lead to planet spin-orbit misalignment. A circumplanetary disk is expected to be in the wave-like In this Letter, we investigate another mechanism that can disk regime (Papaloizou & Pringle 1983) since the Shakura & Sunyaev (1973) α viscosity parameter is much smaller than 2 R. G. Martin & P. J. Armitage the disk aspect ratio H/r (e.g. Martin et al. 2019), where H is differences. The method is first order explicit in time. The the disc scale height and r is the spherical distance from the grid extends from inner radius rin up to rout. In order to ob- planet. Thus, we solve the 1D wave-like warped disk equa- tain convergence, we use a non-uniform grid. There are 25 tions in the frame of the planet-star system. The disk is in grid points spaced linearly in log(r) from the inner radius up Keplerian rotation around the planet at all radii with angular to rchange = 0.02 rH. There are 75 points that are linearly dis- tributed with radius from r up to r . frequency given by Ω = GM /r3. We solve equations (12) change out p The inclinations of the disk and planet are determined via and (13) in Lubow & Ogilvieq (2000) with an additional torque on the disk from the planet. These can be written as −1 2 2 i = cos 1 − lx − ly (9) l G q 2 ∂ d 1 ∂ Σr Ω + Ωp × ld = + Ts + Tp (1) and the phase angles as ∂t ! r ∂r l and = −1 y 2 3 3 φ tan . (10) ∂G ΣH r Ω ∂ld lx ! + Ωp × G + αΩG = , (2) ∂t 4 ∂r For the disk, we calculate the angular momentum weighted where 2πG is the internal disk torque. We only need to solve vector components. We transform these back into the inertial the first two components of these vector equations since the frame in our plots by rotating about the z axis by Ωpt. disc angular momentum vector is a unit vector. We use the boundary conditions G = 0 and ∂ld/∂r = 0 at r = rin and 3. DISK EVOLUTION WITH NO PLANETARY TORQUE G r = rout. The initial conditions are chosen to be = 0 and the We model a Jupiter mass and radius planet orbiting at disk is initially flat but tilted by 10◦ with l = sin(10◦) and d,x ap = 5.2au. The circumplanetary disk has surface density l = − d,y 0. profile Σ ∝ r 3/2 distributed between the inner disk radius We work in a frame where the star is on the positive x axis r = r r = . r r T = in 2 p and the outer radius out 0 4 H, where H is the and so the stellar torque on the disk per unit area is s Hill radius. The inner disk radius is close to the peak in the T , T T = ( s,x s,y) where s,x 0 and surface density of a steady state disk that joins on to the planet 3GMs 2 at radius rp. The outer radius is chosen to be close to the Ts,y = − Σr ld,x. (3) tidal truncation radius of a circumplanetary disk (Martin & 2a3 p Lubow 2011)1. We take the disk parameters α = 0.01 and The torque per unit area from a spinning oblate planet on the H/r = 0.14. With this disk aspect ratio the disk is strongly disk is given in linear theory by unstable to tilting. As shown in Martin et al. (2020), the com- bination of the disk aspect ratio and the disk size determines − Tp,x = T(ld,y lp,y) (4) the instability. While it appears that instability requires spe- and cific parameters in the analytic model, hydrodynamic simula- tions show that that a disk with a more realistic surface density Tp,y = T(−ld,x + lp,x) (5) profile is also unstable to tilting.