2019ApJ...878..119H

We are witness to a great and increasing interest in internal structure, composition, and evolution of exoplanets. However, direct measurements of exoplanetary mass and radius cannot be uniquely interpreted in terms of interior structure, justifying the need for an additional observable. The second degree fluid Love number, k₂, is proportional to the concentration of mass toward the body's center, hence providing valuable additional information about the internal structure. When hydrostatic equilibrium is assumed for the planetary interior, k₂ is a direct function of the planetary shape. Previous attempts were made to link the observed tidally and rotationally induced planetary oblateness in photometric light curves to k₂ using ellipsoidal shape models. Here, we construct an analytical 3D shape model function of the true planetary mean radius that properly accounts for tidal and rotational deformations. Measuring the true planetary mean radius is critical when one wishes to compare the measured k₂ to interior theoretical expectations. We illustrate the feasibility of our method and show, by applying a Differential Evolution Markov Chain to synthetic data of WASP-121b, that a precision <=65 ppm/sqrt(2)minutes is required to reliably retrieve k₂ with present understanding of stellar limb darkening (LD), therefore improving recent results based on ellipsoidal shape models. Any improvement on stellar LD would increase such performance.