<figure><img src="image_1.png" /></figure>

HAL
open science

# Existence of slowly rotating bipolytropes with prolate cores

Clément Staelen, Jean-Marc Huré, Anaïs Meunier, Pauline Noé

► To cite this version:

Clément Staelen, Jean-Marc Huré, Anaïs Meunier, Pauline Noé. Existence of slowly rotating bipolytropes with prolate cores. Astronomy & Astrophysics - A&A, 2024, 691, pp.A276. 10.1051/0004-6361/202450619. hal-04799405

HAL Id: hal-04799405

https://hal.science/hal-04799405v1

Submitted on 22 Nov 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Existence of slowly rotating bipolytropes with prolate cores

C. Staelen<fnref n="1" />,*, J.-M. Huré<fnref n="1" />, A. Meunier<fnref n="2" />, and P. Noé<fnref n="1" />,<fnref n="3" />

Received 6 May 2024 / Accepted 8 October 2024

## ABSTRACT

We report the existence of hydrostatic equilibrium states for a composite body made of two rigidly rotating, homogeneous layers bounded by spheroidal surfaces, where the core has a prolate shape. These new configurations require an oblate envelope that spins faster than the core. No solution with a prolate envelope is found. For some parameters, the prolate core can even be at rest. Numerical experiments based on the self-consistent field method support this result in the case of heterogeneous layers with polytropic equations of state. The possible cancellation of the first gravitational moment, $J_2$, is discussed.

**Key words.** gravitation – methods: analytical – planets and satellites: interiors – stars: interiors – stars: rotation

## 1. Introduction

A common property of self-gravitating systems carrying angular momentum is the natural propensity to flatten, which is indis- putably observed for planets, stars, and galaxies. The centrifugal force is actually extremely efficient in pushing matter away from the rotation axis. As calculated by Maclaurin for a rigidly rotating, homogeneous body, the flattening takes the form of a spheroid (i.e. an ellipsoid of revolution) spinning around its minor axis. Oblateness is also observed for compressible matter (Hachisu 1986). Prolate configurations corresponding to an elongation along the rotation axis have always been ruled out when only gravity, pressure forces, and rotation are considered. Prolateness can be provoked by additional mechanisms or ingre- dients, such as meridional circulation (Fujisawa & Eriguchi 2014) or magnetic fields (e.g. Oganesyan & Abramyan 1972; Lander & Jones 2009; Kawamura et al. 2011; Huré & Staelen 2024).

While prolateness is not permitted in Maclaurin's theory, we show in this article that a homogeneous prolate spheroid can be a figure of equilibrium if it is embedded in an enve- lope externally bounded by an oblate spheroidal surface. This effect, which is demonstrated in the slow-rotating and rigidly rotating limits, is purely hydrodynamical and requires asyn- chronous motion between the two layers. Based upon numerical simulations, we reach a similar conclusion when the two com- ponents are inhomogeneous. This result, which contributes to the theory of figures, provides new tools for studying the for- mation and evolution of piecewise homogeneous systems like rocky and/or icy planets, planetesimals, and asteroids in the Solar System (e.g. Rambaux et al. 2015). Furthermore, it could also help explain the presence of prolate molecular clouds in the interstellar medium (e.g. Fiege & Pudritz 2000; Cai & Taam 2010).

## 2. Rotating bodies with spheroid bounding layers

The hydrostatic equilibrium of a rigidly rotating fluid mass made of several layers bounded by spheroidal surfaces was examined in detail by Hamy (1889, 1890). Assuming the mass density was monotonic with depth, he demonstrated that all the surfaces are necessarily confocal to each other (i.e. all ellipses share the same foci) and all spheroids necessarily rotate asynchronously (see also Véronet 1912; Moulton 1916; Montalvo et al. 1983). These conditions constitute Hamy's 'no-go' theorem.

Out of confocality, the interfaces between layers cannot be exact spheroids, but numerical experiments show that the depar- ture is weak for a wide range of configurations. This led Huré (2022) to investigate the approximation of spheroidal layers for fast rotators. The nested spheroidal figures of equilibrium (NSFoE) theory reported therein assumes that hydrostatic equi- librium is satisfied at the pole and at the equator of each layer. The NSFoE theory is expected to be reliable when all spheroids are close to confocal with each other. As a consequence, this the- ory goes beyond the slowly rotating regime usually considered (Clairaut 1743; Hamy 1889) and reliably describes fast rotators. With this approach, any relative motion of layers is possible, from global rotation, in which all the rotation rates ($\Omega_i$) are the same, to a fully asynchronous regime, in which the $\Omega_i$'s are all different.

In all these theories, only oblate configurations have been obtained or reported<fnref n="1" />. Relaxing this assumption opens up to new states. In this work, we return to the slowly rotating limit, which enables us to bypass the complexity of the general formulae and to derive analytical criteria for the existence of prolate subcom- ponents.

<footnote n="1">¹ Univ. Bordeaux, CNRS, LAB, UMR 5804, F-33600 Pessac, France</footnote>

<footnote n="2">² Univ. Luxembourg, Department of Mathematics, L-4364 Esch-sur-Alzette, Luxembourg</footnote>

<footnote n="3">³ Univ. Strasbourg, Télécom Physique Strasbourg, F-67400 Illkirch-Graffenstaden, France</footnote>

<footnote>* Corresponding author; clement.staelen@u-bordeaux.fr</footnote>

<footnote n="1">¹ Note that prolate subsystems were already considered by D'Alembert (1747) in the context of a rigid core supporting a shallow ocean.</footnote>
<figure><img src="image_3.png" /><figcaption>**Fig. 1.** Heterogeneous system made of two homogeneous layers bound by spheroidal surfaces, namely the core (black) and the envelope (white). Four core-envelope configurations are analysed. O stands for oblate and P for prolate; the first letter refers to the core and the second to the envelope.</figcaption></figure>

### 3. Equation set relevant to a two-layer configuration

Our assumptions and notations are the same as those in Huré (2022). We considered a heterogeneous body made of two homogeneous, spheroidal components, as pictured in Fig. 1. Let $\rho_i$, $\Omega_i$, $a_i$, and $b_i$ be the mass density, the rotation rate, the equatorial radius, and the polar radius of layer $i$, respectively. Index 1 is for the deepest layer (hereafter ‘the core’) and index 2 for the upper layer (hereafter ‘the envelope’). We define the fractional radius of the core $q = a_1/a_2$, the mass-density jump at the interface $\alpha = \rho_1/\rho_2$, the squared eccentricity of layer $i$, $\varepsilon_i^2 = 1 - (b_i/a_i)^2$, and the dimensionless squared rotation rate of layer $i$, $\tilde{\Omega}_i^2 = \Omega_i^2/(2\pi G\rho_2)$, where G is the gravitational constant. Cases with $q = 0$, $q = 1$, or $\alpha = 1$ have no special interest, as they all correspond to a single-component body. Asynchronous motion is allowed: $\Omega_1$ and $\Omega_2$ are not necessarily equal. As mentioned above, the NSFoE theory is reliable when the spheroids are nearly confocal (i.e. when the confocal parameter, $c = q^2\varepsilon_1^2 - \varepsilon_2^2$, is small enough). The error is typically below 1% for $|c| \le 0.5$.

The starting points of the present analysis are Eqs. (66a) and (66b) from Huré (2022), which give the rotation rates at equilibrium of each layer in the regime of slow rotations, in the case of oblate interfaces. At the first order in $\varepsilon_i^2$, we have

$$ \tilde{\Omega}_2^2 = \frac{4}{15}\varepsilon_2^2 - \frac{2}{15}(\alpha - 1)q^3(3q^2\varepsilon_1^2 - 5\varepsilon_2^2), \quad (1a) $$

$$ \alpha \tilde{\Omega}_1^2 = \tilde{\Omega}_2^2 + \frac{2}{15}(\alpha - 1)[2(\alpha - 1)\varepsilon_1^2 + 5\varepsilon_1^2 - 3\varepsilon_2^2]. \quad (1b) $$

We note that these equations are not valid for $q = 0$ or $q = 1$, as one of the layers does not exist. The equations have been established in the case of two oblate interfaces, but they still hold for prolate interfaces (i.e. $\varepsilon_i^2 < 0$). This is due to the mathematical continuity in the expression for gravitational potential when the eccentricity becomes an imaginary number.

Thus, we worked with the same equations for a prolate-oblate mixture. There is therefore a maximum of four types of configurations, as shown in Fig. 1: an oblate core and an oblate envelope (OO), an oblate core with a prolate envelope (OP), a prolate core and an oblate envelope (PO), and a prolate core with a prolate envelope (PP).

The first case is already treated in detail in Huré (2022), where it is shown that there are two families of solutions. The first is when rotations are synchronous, or $\Omega_1 = \Omega_2$. This requires that the core-envelope interface is a surface of constant pressure (type-C solutions). The mass-density jump is not free

(see Eq. (11) below). The second is when rotations are asynchronous, or $\Omega_1 \neq \Omega_2$. The gas pressure at the surface of the core varies quadratically with the cylindrical radius, R (type-V solutions), and $\alpha$ is a free parameter.

In the next section, we examine the mathematical conditions required for the existence of OP, PO, and PP configurations. Throughout this work, we assume $\alpha > 1$ (i.e. the core is denser than its envelope), as this is the case in most problems. However, the case $\alpha < 1$ is treated in Appendix A. The results of the two cases are largely the same; only the conditions of positivity of the $\Omega_i^2$'s differ (see below).

### 4. The conditions of positivity of the $\Omega_i^2$'s

The main constraint is that the rotation rates for the two layers must both be real numbers. From Eqs. (1a) and (1b), requiring that $\Omega_i^2 \ge 0$ leads to the following set of inequalities:

$$ \{[2 + 5(\alpha - 1)q^3]\varepsilon_2^2 \ge 3(\alpha - 1)q^5\varepsilon_1^2, \quad (2a) $$

$$ \{[2 + (\alpha - 1)(5q^3 - 3)]\varepsilon_2^2 \ge (\alpha - 1)[3q^5 - 2\alpha - 3]\varepsilon_1^2, \quad (2b) $$

which, in standard conditions ($\alpha > 1$), becomes

$$ \varepsilon_2^2 \ge \frac{3(\alpha - 1)q^5}{2 + 5(\alpha - 1)q^3}\varepsilon_1^2, \quad (3a) $$

$$ \varepsilon_1^2 \ge -\frac{2 + (\alpha - 1)(5q^3 - 3)}{(\alpha - 1)(3 + 2\alpha - 3q^5)}\varepsilon_2^2, \quad (3b) $$

as $3q^5 - 2\alpha - 3 < 0$. We see from Eq. (3a) that an oblate core ($\varepsilon_1^2 > 0$) is only compatible with an oblate envelope ($\varepsilon_2^2 > 0$). OP configurations are thus ruled out. So we only consider solutions with prolate cores in the following. With the condition $\varepsilon_1^2 \le 0$, Eq. (3b) leads to

$$ [2 + (\alpha - 1)(5q^3 - 3)]\varepsilon_2^2 \ge 0. \quad (4) $$

There are therefore two remaining cases, namely the PO and PP configurations.

#### 4.1. Prolate envelope solution ruled out

With $\varepsilon_2^2 \le 0$, the above inequality reads $2 + (\alpha - 1)(5q^3 - 3) \le 0$, and so Eq. (2b) becomes

$$ \frac{3(\alpha - 1)q^5}{2 + 5(\alpha - 1)q^3}\varepsilon_1^2 \le \varepsilon_2^2 \le \frac{(\alpha - 1)[3q^5 - 2(\alpha - 1) - 5]}{2 + (\alpha - 1)(5q^3 - 3)}\varepsilon_1^2, \quad (5) $$

where we take Eq. (2a) into account. As $\varepsilon_1^2 \le 0$ is assumed, we find

$$ \frac{3q^5}{2 + 5(\alpha - 1)q^3} \ge \frac{3q^5 - 2(\alpha - 1) - 5}{2 + (\alpha - 1)(5q^3 - 3)}, \quad (6) $$

which also reads

$$ 4(\alpha - 1) + 10(\alpha - 1)^2 q^3 + 10 + 25(\alpha - 1)q^3 \le 9(\alpha - 1)q^5. \quad (7) $$

As $\alpha - 1 > 0$ and given $0 < q < 1$, we have

$$ 25(\alpha - 1)q^3 \ge 9(\alpha - 1)q^5, \quad (8) $$

which implies that Eq. (7) cannot be verified. We conclude that a prolate core surrounded by a prolate envelope cannot be at equilibrium. Thus, PP configurations are ruled out as well. We conclude that in the regime of slow rotations, a body made of two rotating, homogeneous layers bound by spheroidal surfaces cannot represent an equilibrium configuration if the envelope is prolate, whatever the shape of the core.
## 4.2. Oblate envelope solution permitted

If we now impose $\varepsilon_2^2 \ge 0$, we see that Eq. (2a) is automatically verified. Furthermore, we have $2 + (\alpha - 1)(5q^3 - 3) \ge 0$ in this case, which implies

$$q \ge \left[ \frac{3(\alpha - 1) - 2}{5(\alpha - 1)} \right]^{1/3} \equiv q_{\lim} \quad (9)$$

and

$$\varepsilon_2^2 \ge \frac{(\alpha - 1)[3q^5 - 2(\alpha - 1) - 5]}{2 + (\alpha - 1)(5q^3 - 3)} \varepsilon_1^2 \equiv \varepsilon_{2, \lim}^2. \quad (10)$$

We see that these two conditions can be simultaneously satisfied in the domain of interest. There is, however, a lower limit for the relative size of the core. We find $q_{\lim} = 0$ for $\alpha = 5/3$, and $q_{\lim} = (3/5)^{1/3} \approx 0.843$ for $\alpha \to \infty$. There is also a lower limit for the eccentricity of the envelope's surface, that is, $\varepsilon_{2, \lim}^2$. Thus, we can state that in the regime of slow rotations, a body made of two rotating layers bound by spheroidal surfaces can be in hydrostatic equilibrium if an oblate envelope surrounds a prolate core.

## 4.3. Note on the global rotations

As shown above, the selection of physically relevant configurations is based on the sign of the $\Omega_i$'s, not on their relative magnitude. As outlined in Huré (2022) in relation to OO configurations only, states of global rotation, that is, $\Omega_1 = \Omega_2$, can occur either for $\alpha = 1$, which corresponds to a single-layer object, or for a special mass-density jump, $\alpha_C$, namely

$$\alpha \equiv \alpha_C = 1 + \frac{5(\varepsilon_1^2 - \varepsilon_2^2)}{2\varepsilon_1^2 + 3q^5\varepsilon_1^2 - 5q^3\varepsilon_2^2}. \quad (11)$$

This expression is obtained by equalizing $\tilde{\Omega}_1^2$ and $\tilde{\Omega}_2^2$ in Eq. (1) and assuming $\alpha \ne 1$. In this case, the rotation rate of the whole structure is given by

$$\tilde{\Omega}^2 = \frac{4}{3} \frac{\varepsilon_1^2 - \varepsilon_2^2}{2\varepsilon_1^2 + 3q^5\varepsilon_1^2 - 5q^3\varepsilon_2^2} \varepsilon_1^2 + \frac{2}{3}\varepsilon_1^2 - \frac{2}{5}\varepsilon_2^2 \ge 0. \quad (12)$$

As proved above, only OO and PO configurations are permitted. So, assuming $\alpha \ne 1$, $\varepsilon_1^2 \le 0$, and $\varepsilon_2^2 \ge 0$, Eq. (12) implies

$$(20 + 15q^3)\varepsilon_1^4 - (16 + 25q^3 + 9q^5)\varepsilon_1^2\varepsilon_2^2 + 15q^3\varepsilon_2^4 \le 0. \quad (13)$$

The left-hand side of this inequality is a second-degree polynomial in $\varepsilon_1^2$. The two roots are

$$\varepsilon_{1\pm}^2 = \frac{16 + 25q^3 + 9q^5}{10(4 + 3q^3)} \varepsilon_2^2 \left[ 1 \pm \sqrt{1 - \frac{300q^3(4 + 3q^3)}{(16 + 25q^3 + 9q^5)^2}} \right], \quad (14)$$

and it is easy to see that both roots are positive. As the prefactor of $\varepsilon_1^4$ in the second-degree polynomial is positive, Eq. (12) is verified only between the two roots $\varepsilon_{1\pm}$ (i.e. for $\varepsilon_1^2 \ge 0$). This contradicts the starting hypothesis, which means that in the limit of slow rotations, there are no equilibria with a prolate core and an oblate envelope spinning at the same rate. In other words, under our assumptions, the presence of a rotational discontinuity is mandatory for the existence of prolate cores, that is, $\Omega_1 \ne \Omega_2$. The only exception is the case $\alpha = 1$, where the distinction between the prolate core and the oblate envelope is impossible.

<figure><img src="image_4.png" /><figcaption>Fig. 2. Rotational discontinuity, $\Delta$, versus $\alpha$ defined by Eq. (16) for $\varepsilon_2^2 = 0.01$ (top) and $\varepsilon_2^2 = 0.1$ (bottom). In both cases, $\varepsilon_1^2 = -0.01$ corresponds to a prolate core surrounded by an oblate envelope. See Tables 1 and 2 for configurations A and B.</figcaption></figure>

# 5. Results

## 5.1. Rotation rates

As we focused on slow rotations (and small deviations from sphericity), the $\tilde{\Omega}_i$'s are necessarily small values. In contrast, the rotational discontinuity at the interface, defined as

$$\Delta = \frac{\tilde{\Omega}_1}{\tilde{\Omega}_2}, \quad (15)$$

can take any positive value. From Eqs. (1a) and (1b), we have

$$\Delta^2 = \frac{1}{\alpha} \left[ 1 - (\alpha - 1) \frac{3\varepsilon_2^2 - 2(\alpha - 1)\varepsilon_1^2 - 5\varepsilon_1^2}{2\varepsilon_2^2 - (\alpha - 1)q^3(3q^2\varepsilon_1^2 - 5\varepsilon_2^2)} \right], \quad (16)$$

which is a complicated function of four parameters. However, for any PO configuration, the fraction in Eq. (16) is negative, which implies that the term inside the brackets is less than unity. As a consequence, the oblate envelope always rotates faster than the prolate core.

Figure 2 displays $\Delta$ versus $\alpha$ for two pairs ($\varepsilon_1^2, \varepsilon_2^2$). We see that all curves converge to $\Delta = 1$ for $\alpha = 1$ because this case corresponds to the single-layer Maclaurin spheroid. Roughly speaking, larger mass-density jumps are permitted as the oblateness of the envelope increases. Next, we observe that $\Omega_1 \to 0$ for some special values of the parameters. This is expected from the equation set. We see from Eq. (1) that the condition $\tilde{\Omega}_1^2 = 0$ can be obtained for a special value of mass-density contrast $\alpha_{\lim} - 1 > 0$, which is the root of a second-degree polynomial, namely

$$\varepsilon_1^2 (\alpha_{\lim} - 1)^2 + \frac{5\varepsilon_1^2 - 3\varepsilon_2^2 - q^3(3q^2\varepsilon_1^2 - 5\varepsilon_2^2)}{2} (\alpha_{\lim} - 1) + \varepsilon_2^2 = 0. \quad (17)$$

It therefore follows that a prolate core can be at rest while the envelope is rotating. In a similar way, the positivity condition
$\tilde{\Omega}_1^2 \ge 0$ corresponds to a minimum value for the core's eccentricity, namely

$$ \varepsilon_1^2 \ge - \frac{[(\alpha - 1)(5q^3 - 3) + 2]\varepsilon_2^2}{(\alpha - 1)[(5 - 3q^5) + 2(\alpha - 1)]} \equiv \varepsilon_{1,\min}^2. \quad (18) $$

## 5.2. Two examples

Two-layer configurations hosting prolate cores are easily generated as soon as Eqs. (9) and (10) are simultaneously satisfied. There are four input parameters, namely $q$, $\varepsilon_1$, $\varepsilon_2$, and $\alpha$. We had the opportunity to compare the analytical solution with the equilibrium computed from the DROP code, which numerically solves the axisymmetric problem using the self-consistent field (SCF) method (Hachisu 1986). This code generates the mass-density and pressure profiles, the equilibrium surfaces, and the interfaces, which generally deviate slightly from spheroids. The performance and reliability of the simulations have already been widely demonstrated (Huré & Hersant 2017; Huré et al. 2018), in particular the multi-layer version (Basillais & Huré 2021). To illustrate the calculations reported in this section and the previous one, we considered two examples with large cores. The input and output data are listed in Table 1 for $\alpha = 1.25$, $q = 0.7$, and $\varepsilon_2^2 = 0.01$, and in Table 2 for $\alpha = 2.5$, $q = 0.8$, and $\varepsilon_2^2 = 0.1$. Regarding the SCF code, the numerical experiments were performed with 257 nodes per direction on a cylindrical ($R, Z$) grid. To quantify the accuracy of the solutions, we used the relative virial parameter $|VP/W|$ (see e.g. Kadam et al. 2016; Basillais & Huré 2021), where $VP = W + 2T + U$, $W, T$, and $U$ being the gravitational, kinetic, and internal energies, respectively. At equilibrium, we should have $VP = 0$ (see e.g. Cox & Giulì 1968). However, due to the numerical treatment (DROP is second-order accurate in the grid spacing) and the approximation of NSFoE, VP cannot be exactly null in both approaches. The absolute value of this parameter with respect to the greatest energy (i.e. $|W|$) is expected to be of the order of magnitude of the error made with the method.

As we can see by comparing columns 2 and 4 of Tables 1 and 2, the prolate cores predicted by the calculus are also observed in the numerical experiments, and the axis ratio (which is an output for DROP) coincides with the value of the analytical solution within the numerical precision. The virial parameters are small enough that we can assert that the prolateness of the core is not an artefact of the numerical treatment. For configuration A, the values of the rotation rates are also in good agreement, and this is true for the other global quantities as well. For configuration B, the accordance is a little less satisfactory. This is due to the fact that Eqs. (1a) and (1b) are first-order accurate in the eccentricities, while $\varepsilon_2^2 = 0.1$ is marginally acceptable regarding this expansion. The third columns of Tables 1 and 2 show the results obtained via the analytical approach, but with the full expressions of the rotation rates. We see that the agreement between the full NSFoE theory and the numerical SCF method is remarkable, and this confirms the existence of PO configurations.

# 6. Discussion and perspectives

We examined the conditions of the hydrostatic equilibrium of a heterogeneous body made of two rigidly rotating, homogeneous layers bound by a spheroidal surface. Given our hypothesis, and in particular the limit of slow rotations, we find a new family of solutions that involves a prolate core. We have shown that:

- The existence of a prolate core requires an asynchronous motion between the two layers.

<table><caption>**Table 1.** Example of a configuration and the main output data obtained with a prolate core and an oblate envelope.</caption><thead><tr><th>Config. A</th><th>DROP</th><th>Full NSFoE</th><th>First-order</th></tr></thead><tbody><tr><td>ε<sub>1</sub><sup>2</sup></td><td>-0.00979</td><td>← -0.01</td><td>← -0.01</td></tr><tr><td>ε<sub>2</sub><sup>2</sup></td><td>← 0.01</td><td>← 0.01</td><td>← 0.01</td></tr><tr><td>q</td><td>0.70007</td><td>← 0.70</td><td>← 0.70</td></tr><tr><td>α</td><td>← 1.25</td><td>← 1.25</td><td>← 1.25</td></tr><tr><td>-c</td><td>1.685 × 10<sup>-2</sup></td><td>1.700 × 10<sup>-2</sup></td><td>1.700 × 10<sup>-2</sup></td></tr><tr><td>ε<sub>2,lim</sub><sup>2</sup></td><td></td><td></td><td>7.439 × 10<sup>-3</sup></td></tr><tr><td>q<sub>lim</sub></td><td></td><td></td><td>-1.00000</td></tr><tr><td>b<sub>1</sub>/a<sub>2</sub></td><td>← NSFoE</td><td>0.70349</td><td>0.70349</td></tr><tr><td>Ω̃<sub>1</sub><sup>2</sup></td><td>4.585 × 10<sup>-4</sup></td><td>4.596 × 10<sup>-4</sup></td><td>4.585 × 10<sup>-4</sup></td></tr><tr><td>Ω̃<sub>2</sub><sup>2</sup></td><td>3.411 × 10<sup>-3</sup></td><td>3.420 × 10<sup>-3</sup></td><td>3.406 × 10<sup>-3</sup></td></tr><tr><td>Δ</td><td>← NSFoE</td><td>3.666 × 10<sup>-1</sup></td><td>3.669 × 10<sup>-1</sup></td></tr><tr><td>M/(ρ<sub>2</sub>a<sub>2</sub><sup>3</sup>)</td><td>4.52913</td><td>4.52877</td><td>4.52883</td></tr><tr><td>p<sub>c</sub>/(πGρ<sub>2</sub><sup>2</sup>a<sub>2</sub><sup>2</sup>)</td><td>8.963 × 10<sup>-1</sup></td><td>8.963 × 10<sup>-1</sup></td><td>8.962 × 10<sup>-1</sup></td></tr><tr><td>|VP/W|</td><td>7 × 10<sup>-7</sup></td><td>4 × 10<sup>-6</sup></td><td>3 × 10<sup>-3</sup></td></tr></tbody></table>

**Notes.** $p_c$ is the central pressure and $M$ is the mass; ‘←’ indicates input data. The second column contains the results from the numerical solution obtained with the DROP code (Basillais & Huré 2021). The third column contains the results from the full NSFoE theory (Huré 2022) and the fourth the results obtained with the first-order approximation of $\varepsilon_2^2$ used in this work.

<table><caption>**Table 2.** Same as Table 1 but for configuration B.</caption><thead><tr><th>Config. B</th><th>DROP</th><th>Full NSFoE</th><th>This work</th></tr></thead><tbody><tr><td>ε<sub>1</sub><sup>2</sup></td><td>-0.01025</td><td>← -0.01</td><td>← -0.01</td></tr><tr><td>ε<sub>2</sub><sup>2</sup></td><td>← 0.1</td><td>← 0.1</td><td>← 0.1</td></tr><tr><td>q</td><td>0.79990</td><td>← 0.8</td><td>← 0.8</td></tr><tr><td>α</td><td>← 2.5</td><td>← 2.5</td><td>← 2.5</td></tr><tr><td>-c</td><td>1.082 × 10<sup>-2</sup></td><td>1.080 × 10<sup>-2</sup></td><td>1.080 × 10<sup>-2</sup></td></tr><tr><td>ε<sub>2,lim</sub><sup>2</sup></td><td></td><td></td><td>7.855 × 10<sup>-2</sup></td></tr><tr><td>q<sub>lim</sub></td><td></td><td></td><td>6.934 × 10<sup>-1</sup></td></tr><tr><td>b<sub>1</sub>/a<sub>2</sub></td><td>← NSFoE</td><td>8.040 × 10<sup>-1</sup></td><td>8.040 × 10<sup>-1</sup></td></tr><tr><td>Ω̃<sub>1</sub><sup>2</sup></td><td>1.925 × 10<sup>-3</sup></td><td>1.931 × 10<sup>-3</sup></td><td>1.533 × 10<sup>-3</sup></td></tr><tr><td>Ω̃<sub>2</sub><sup>2</sup></td><td>8.469 × 10<sup>-2</sup></td><td>8.491 × 10<sup>-2</sup></td><td>7.983 × 10<sup>-2</sup></td></tr><tr><td>Δ</td><td>← NSFoE</td><td>1.508 × 10<sup>-1</sup></td><td>1.385 × 10<sup>-1</sup></td></tr><tr><td>M/(ρ<sub>2</sub>a<sub>2</sub><sup>3</sup>)</td><td>7.20079</td><td>7.20687</td><td>7.21242</td></tr><tr><td>p<sub>c</sub>/(πGρ<sub>2</sub><sup>2</sup>a<sub>2</sub><sup>2</sup>)</td><td>3.10451</td><td>3.10511</td><td>3.16290</td></tr><tr><td>|VP/W|</td><td>2 × 10<sup>-5</sup></td><td>7 × 10<sup>-4</sup></td><td>6 × 10<sup>-2</sup></td></tr></tbody></table>

- The envelope is necessarily oblate.

- The envelope spins faster than the core.

- The prolate core can be at rest (for some special parameters).

- Global rotation is only compatible with an oblate core surrounded by an oblate envelope.

These results are fully confirmed by numerical experiments based on the SCF method.

## 6.1. Prolate cores in compressible systems

A natural question concerns the extension of the results reported above to the case where the mass densities are no longer uniform in each layer. This is easily checked by performing numerical simulations with appropriate equations of state. For this purpose, we used the DROP code again, and a polytropic equation of state where the gas pressure, $p$, is a power law of the mass density,
<table><caption>**Table 3.** Same as Table 1 but for fully heterogeneous layers, computed with the DROP code; the solutions obtained are shown in Figs. 3 and 4.</caption><thead><tr><th></th><th>Config. C</th><th>Config. D</th></tr></thead><tbody><tr><td>ε<sub>1</sub><sup>2</sup></td><td>-0.00996</td><td>-0.01039</td></tr><tr><td>ε<sub>2</sub><sup>2</sup></td><td>← 0.01</td><td>← 0.1</td></tr><tr><td>n<sub>1</sub></td><td>← 3.0</td><td>← 0.3</td></tr><tr><td>n<sub>2</sub></td><td>← 1.5</td><td>← 0.5</td></tr><tr><td>q</td><td>0.70001</td><td>0.79985</td></tr><tr><td>b<sub>1</sub>/a<sub>2</sub></td><td>← 0.70349</td><td>← 0.80399</td></tr><tr><td>α</td><td>← 1.25</td><td>← 2.5</td></tr><tr><td>Ω<sub>1</sub><sup>2</sup>/(2πGρ<sub>c</sub>)</td><td>3.420 × 10<sup>-5</sup></td><td>2.667 × 10<sup>-3</sup></td></tr><tr><td>Ω<sub>2</sub><sup>2</sup>/(2πGρ<sub>c</sub>)</td><td>1.710 × 10<sup>-4</sup></td><td>2.614 × 10<sup>-2</sup></td></tr><tr><td>Δ</td><td>← 0.44721</td><td>← 0.31937</td></tr><tr><td>M/(ρ<sub>c</sub>a<sub>2</sub><sup>3</sup>)</td><td>0.11012</td><td>1.86401</td></tr><tr><td>p<sub>c</sub>/(Gρ<sub>c</sub><sup>2</sup>a<sub>2</sub><sup>2</sup>)</td><td>8.352 × 10<sup>-2</sup></td><td>1.02538</td></tr><tr><td>|VP/W|</td><td>7 × 10<sup>-5</sup></td><td>6 × 10<sup>-4</sup></td></tr></tbody></table>

Notes. $\rho_c$ is the central mass density.

$\rho$, namely $p \propto \rho^{1+1/n}$, and $n$ is the polytropic index. Table 3 presents two new configurations, namely:

- Configuration C, which resembles configuration A but with $(n_1, n_2) = (3, 1.5)$. These parameters are appropriate for a radiative pressure-dominated core surrounded by a convective envelope. The mass-density jump is small, and the envelope spins faster than the core by a factor of about 2.

- Configuration D, which resembles configuration B but with $(n_1, n_2) = (0.3, 0.5)$. This case corresponds to weakly compressible layers. The mass-density jump is twice as large as for configuration C, and the envelope spins faster than the core by a factor of about 3.

Figures 3 and 4 show the mass density and the deviations of the surface bounding layers from spheres. The right panels show the deviation $\Delta Z_i(R)$ of the true surfaces from spheres, namely

$$ \Delta Z_i = Z_i(R) - \sqrt{a_i^2 - R^2}, \quad i \in \{1, 2\}, \qquad (19) $$

where $Z_i$ is the equation of the boundary of layer $i$, and $R$ is the distance to the rotation axis. Negative (positive) values correspond to oblate (prolate) surfaces. For these two configurations, we observe that the core-envelope interface is indeed prolate in shape (while isobars remain oblate).

## 6.2. Possible cancellation of the gravitational moments

A prolate core embedded in an oblate envelope leaves a mark in the gravitational potential, as it mitigates the effect of an oblate envelope. It can be shown that one of the gravitational moments, $J_{2n}$, can vanish (when $n$ is an odd number). The first coefficient, $J_2$, is of special interest because it is the main deviation from sphericity and therefore of interest for dynamical studies. For spheroidal heteroeoids, the $J_{2n}$'s are analytical (Heiskanen & Moritz 1967; Cisneros-Parra et al. 2019; Basillais & Huré 2023). We find that $J_2 = 0$ when

$$ \varepsilon_2^2 \sqrt{1 - \varepsilon_2^2} + (\alpha - 1)q^5 \varepsilon_1^2 \sqrt{1 - \varepsilon_1^2} = 0, \qquad (20) $$

which is a formally cubic equation in $\varepsilon_1^2$. In the actual conditions of slow rotation, this implies

$$ \varepsilon_1^2 \approx -\frac{\varepsilon_2^2}{(\alpha - 1)q^5} \le 0, \qquad (21) $$

<figure><img src="image_5.png" /><figcaption>**Fig. 3.** Equilibrium solution computed with the DROP code for configuration C: the normalized mass-density, $\rho/\rho_c$, map in colour code (left) and the deviation from sphericity (right). See also Table 3.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Fig. 4.** Same as Fig. 3 but for configuration D.</figcaption></figure>

which has to be considered together with Eqs. (9) and (10). This concerns a wide range of equilibria, from large cores and small mass-density jumps to small cores with very high $\alpha$'s, typically.

## 6.3. Extension to fast rotators

An interesting point concerns the extension of the results to the case of fast rotations, which might reveal the existence of very elongated prolate cores. By construction, the NSFoE theory is valid for highly flattened systems, provided the confocal parameter (see Sect. 3) is small enough. It then fully applies to structures elongated along the axis of rotation. We surveyed the $(\varepsilon_1^2, \varepsilon_2^2)$ plane, with $-1 \le \varepsilon_i^2 \le 1$, for several values of $q$ and $\alpha$. The main conclusions of this additional exploration are:

- Again, there is no solution with a prolate envelope, regardless of the shape and size of the core.

- The permitted range of prolateness for the core gets wider as $\alpha$ gets closer to 1 and/or $q$ decreases.

- The permitted range of prolateness for the core gets wider as the envelope gets more flattened.

An example of a fast-rotating envelope surrounding a prolate core, obtained from both the DROP code and NSFoE, is given in Table 4. We see that the confocal parameter is marginally acceptable in this case ($|c| \sim 0.5$), but the analytical solution is still in good agreement with the numerical experiment.
<table><caption>**Table 4.** Same as Table 1 but for a fast-rotating envelope.</caption><thead><tr><th>Config. E</th><th>DROP</th><th>Full NSFoE</th></tr></thead><tbody><tr><td>ε<sub>1</sub><sup>2</sup></td><td>-0.26394</td><td>← -0.26394</td></tr><tr><td>ε<sub>2</sub><sup>2</sup></td><td>← 0.37589</td><td>← 0.37589</td></tr><tr><td>q</td><td>0.66709</td><td>← 0.66709</td></tr><tr><td>α</td><td>← 1.25</td><td>← 1.25</td></tr><tr><td>-c</td><td>4.933 × 10<sup>-1</sup></td><td>4.933 × 10<sup>-1</sup></td></tr><tr><td>b<sub>1</sub>/a<sub>2</sub></td><td>← 0.75</td><td>0.74998</td></tr><tr><td>Ω̃<sub>1</sub><sup>2</sup></td><td>2.811 × 10<sup>-2</sup></td><td>2.934 × 10<sup>-2</sup></td></tr><tr><td>Ω̃<sub>2</sub><sup>2</sup></td><td>1.406 × 10<sup>-1</sup></td><td>1.419 × 10<sup>-1</sup></td></tr><tr><td>Δ</td><td>← 0.44721</td><td>4.547 × 10<sup>-1</sup></td></tr><tr><td>M/(ρ<sub>2</sub>a<sub>2</sub><sup>3</sup>)</td><td>3.62718</td><td>3.65865</td></tr><tr><td>p<sub>c</sub>/(πGρ<sub>2</sub><sup>2</sup>a<sub>2</sub><sup>2</sup>)</td><td>7.232 × 10<sup>-1</sup></td><td>7.236 × 10<sup>-1</sup></td></tr><tr><td>|VP/W|</td><td>4 × 10<sup>-5</sup></td><td>8 × 10<sup>-3</sup></td></tr></tbody></table>

## 6.4. Note on the problem of friction between layers

A major question concerns the time stability of the asynchronously rotating layers, even in the limit of slow rotation. The existence of a permanent difference in the velocities at the core-envelope interface is a source of friction, local heating, and energy dissipation, at the expense of the relative motion. If there is no mechanism able to maintain the difference in the rotation rates, then the shear is expected to brake the envelope and accelerate the core through a boundary layer until the point of synchronicity, that is, $\Delta \to 1$. The problem of friction between asynchronously rotating, spherical layers was studied in detail by Wavre (1925) for viscous forces at the interface of the form $\kappa|\Omega_1 - \Omega_2|^{m}R^{m}$, where $\kappa$ is a friction constant and $m$ the friction index, both positive quantities. For a Stokes-like viscosity, that is, $m = 1$, Wavre showed that the relative rotation rate $|\Omega_1 - \Omega_2|$ between the spheres undergoes an exponential decrease with a characteristic time, $\tau$, given by

$$ \frac{1}{\tau} = \frac{8}{3} \pi a_1^4 \kappa \frac{I_1 + I_2}{I_1 I_2}, \qquad (22) $$

where $I_1$ and $I_2$ are the moments of inertia of layers 1 and 2 with respect to the rotation axis, respectively. As the deviation of the core-envelope interface from a sphere is assumed to be weak, Eq. (22) corresponds to the zeroth-order term in $\epsilon_1^2$, which is enough to obtain an approximate order of magnitude of $\tau$. The main difficulty of this estimation of the characteristic time is the value of $\kappa$. For a fluid core embedded in a precessing mantle, we have $\kappa \propto I_1 \sqrt{\nu\Omega_2/a_1^5}$ (Busse 1968). However, this expression is probably not relevant for non-precessing, asynchronously rotating layers. Moreover, the viscosity is often unknown in astrophysical contexts. Even for the outer core of the Earth, the uncertainty on the kinematic viscosity, $\nu$, covers 13 orders of magnitude (Lumb & Aldridge 1991), the best estimate being $\nu \sim 10^{-6} \text{ m}^2 \cdot \text{s}^{-1}$ (Gans 1972; Poirier 1988). To evaluate the order of magnitude of $\tau$, we assumed that the expression of Busse (1968) is valid in our case and that

$\Omega_2$ is the initial value of the rotation rate of the envelope, namely

$$ \tau \approx \frac{1}{2.62} \frac{a_1}{(2\pi G \rho_2 \nu^2)^{1/4}} \frac{1}{\sqrt{\Omega_2}} \frac{1 - q^5}{1 + (\alpha - 1)q^5}. \qquad (23) $$

For the parameters of configuration A in Table 1, taking $\nu = 10^{-6} \text{ m}^2 \cdot \text{s}^{-1}$ and assuming a planetary core with $a_1 = 1000 \text{ km}$ and a mass density of the envelope $\rho_2 = 4000 \text{ kg} \cdot \text{m}^{-3}$ (typically silicates), we have $\tau \sim 10^3 \text{ yr}$. Such a small value is expected, as this situation corresponds roughly to a slightly differentiated rocky planet; the friction between the components is thus strong. The situation is probably more favourable in the interstellar medium, where the components are mostly gaseous and so the friction is weaker. Unfortunately, the viscosity is even less constrained in such contexts.

*Acknowledgements.* We are grateful to the referee, F. Chambat, for his suggestions and remarks to improve the paper, in particular for clarifications on the classical works of Clairaut, Poincaré and Hamy, and for pointing out the monograph of D'Alembert. We also thank P. Auclair-Desrotour for interesting references on the problem of friction at planetary core-mantle boundaries.

## References

* Basillais, B., & Huré, J.-M. 2021, MNRAS, 506, 3773

* Basillais, B., & Huré, J.-M. 2023, MNRAS, 520, 1504

* Busse, F. H. 1968, J. Fluid Mech., 33, 739

* Cai, M. J., & Taam, R. E. 2010, ApJ, 709, L79

* Cisneros-Parra, J. U., Martinez-Herrera, F. J., & Montalvo-Castro, J. D. 2019, ApJS, 241, 8

* Clairaut, A. C. 1743, Théorie de la figure de la Terre tirée des principes de l'hydrostatique (Paris: David Fils)

* Cox, J. P., & Giuli, R. T. 1968, Principles of stellar structure Volume I: Physical principles (New York: Gordon and Breach)

* D'Alembert, J. 1747, Réflexions sur la cause générale des vents (Paris: David l'Aîné)

* Fiege, J. D., & Pudritz, R. E. 2000, ApJ, 534, 291

* Fujisawa, K., & Eriguchi, Y. 2014, MNRAS, 438, L61

* Gans, R. F. 1972, J. Geophys. Res., 77, 360

* Hachisu, I. 1986, ApJS, 61, 479

* Hamy, M. 1889, Annales de l'Observatoire de Paris, 19, F.1.

* Hamy, M. 1890, Journal de mathématiques pures et appliquées 4e série, 6, 69

* Heiskanen, W., & Moritz, H. 1967, Physical Geodesy (San Francisco: W.H. Freeman)

* Huré, J.-M. 2022, MNRAS, 512, 4031

* Huré, J.-M., & Hersant, F. 2017, MNRAS, 464, 4761

* Huré, J.-M., & Staelen, C. 2024, Phys. Rev. D, 110

* Huré, J.-M., Hersant, F., & Nasello, G. 2018, MNRAS, 475, 63

* Kadam, K., Motl, P. M., Frank, J., Clayton, G. C., & Marcello, D. C. 2016, MNRAS, 462, 2237

* Kawamura, T., Taniguchi, K., Yoshida, S., & Eriguchi, Y. 2011, MNRAS, 416, L75

* Lander, S. K., & Jones, D. I. 2009, MNRAS, 395, 2162

* Lumb, L. I., & Aldridge, K. D. 1991, J. Geomagn. Geoelec., 43, 93

* Montalvo, D., Martinez, F. J., & Cisneros, J. 1983, Rev. Mex. Astron. Astrofis., 5, 293

* Moulton, E. J. 1916, Trans. Am. Math. Soc., 17, 100

* Oganesyan, R. S., & Abramyan, M. G. 1972, Astrophysics, 8, 352

* Poirier, J. P. 1988, Geophys. J. Int., 92, 99

* Rambaux, N., Chambat, F., & Castillo-Rogez, J. C. 2015, A&A, 584, A127

* Véronet, A. 1912, Journal de mathématiques pures et appliquées, 8, 331

* Wavre, R. 1925, Archive des sciences physiques et naturelles 5me période, 7, 133
# Appendix A: The case $\alpha < 1$

Considering $\alpha < 1$, the system of inequalities (1) becomes

$$ \varepsilon_1^2 \geq -\frac{2 - 5(1 - \alpha)q^3}{3(1 - \alpha)q^5} \varepsilon_2^2, \quad (A.1a) $$

$$ \frac{2 - (1 - \alpha)[5q^3 - 3]}{(1 - \alpha)[5 - 2(1 - \alpha) - 3q^5]} \varepsilon_2^2 \geq \varepsilon_1^2 \quad (A.1b) $$

Combining the two inequalities of the system, we readily obtain

$$ \frac{2 - (1 - \alpha)[5q^3 - 3]}{(1 - \alpha)[5 - 2(1 - \alpha) - 3q^5]} \varepsilon_2^2 \geq -\frac{2 - 5(1 - \alpha)q^3}{3(1 - \alpha)q^5} \varepsilon_2^2 \quad (A.2) $$

As before, we now have two cases to consider on the sign of $\varepsilon_2^2$.

## A.1. Prolate envelope ruled out

We first assume $\varepsilon_2^2 < 0$. We thus have

$$ -\frac{2 - 5(1 - \alpha)q^3}{3(1 - \alpha)q^5} \geq \frac{2 - (1 - \alpha)[5q^3 - 3]}{(1 - \alpha)[5 - 2(1 - \alpha) - 3q^5]} \quad (A.3) $$

which leads to

$$ 10 - [4 + 25q^3 - 9q^5](1 - \alpha) + 10(1 - \alpha)^2 \leq 0. \quad (A.4) $$

As $4 + 25q^3 - 9q^5 < 20$, the discriminant of the second-degree polynomial in $1 - \alpha$ in the left hand side is negative, except for $q = 1$. In this case, however, we see that the inequality becomes $\alpha^2 \leq 0$, which is not possible. Thus, in a slightly flattened, heterogeneous system made of two homogeneous components, there exists no equilibrium with a prolate envelope.

## A.2. Prolate core, oblate envelope permitted

In this study, we are interested in prolate subsystems. We then impose $\varepsilon_1^2 \leq 0$. Thus, Eq. (A.1b) is automatically verified and Eq. (A.1a) reads

$$ \varepsilon_2^2 \geq -\frac{3(1 - \alpha)q^5}{2 - 5(1 - \alpha)q^3} \varepsilon_1^2, \quad (A.5) $$

where we have necessarily $2 - 5(1 - \alpha)q^3 \geq 0$, meaning that we have an upper limit on $q$ (whereas we had a lower limit in the $\alpha > 1$ case), namely

$$ q \leq \left( \frac{2}{5(1 - \alpha)} \right)^{1/3} \equiv q'_{\text{lim}}. \quad (A.6) $$

Thus, in a slightly flattened, heterogeneous system made of two homogeneous components, there exists equilibria with an oblate envelop and a less dense, prolate core.