Supplementary Informationvalues we use a suite of core/envelope internal structure models, with a wide range of core depths, core densities and envelope density profiles.The core is taken to be of constant density with the core density ranging between 0.8 × 10 4 < ρ core < 1.2 × 10 4 kg m −3 , and core radius extending up to 30% of the planetary radius on Neptune (0 < r core < 0.3 R, where R is the mean planetary radius), and 20% on Uranus, following the widest possible range of core densities given by [16, 17, 19]. We systematically explore this parameter space where for each case of ρ core and r core we then match an envelope represented by a 6th order polynomial with a monotonic density distribution so that ρ static = ρ core for 0 < r < r core (S1)where β = r/R is the normalized radius. The coefficients (a n ) are determined by constraints on the static density (ρ static ). The first constraint on the static density is that the density is zero at the surface, satisfied by setting the sum of the six polynomial coefficients to zero at β = 1. The second constraint is that the integrated density over the entire volume must equal to the planetary mass.The third constraint is that the density derivative at β = 1 equals the derivative of the density at the 1 bar pressure level (equal to -0.1492 and -0.2425 kg m −4 for Uranus and Neptune respectively,[31]). The first degree term in Eq. S2 is missing so that the derivative of the density goes to zero at the center for models with no core, and another constraint sets this value to zero at the core-envelope boundary for models with cores. The last constraint limits the value of J 2 to within the error estimates of the observed value of 3341.29 ± 0.72 × 10 −6 and 3408.43 ± 4.50 × 10 −6 for 1 Uranus and Neptune respectively [14, 15]. Since we are interested in determining the resulting J 4 we do not impose any constraints on J 4 , as done in most studies where the J 4 is constrained to within the observed values of J 4 (e.g., [13]). For cases with no core the density at β = 0 is not constrained, and its derivative at β = 0 is set to zero.Uranus Neptune 25