Open Access
Issue
A&A
Volume 691, November 2024
Article Number A276
Number of page(s) 7
Section Astrophysical processes
DOI https://doi.org/10.1051/0004-6361/202450619
Published online 22 November 2024

© The Authors 2024

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1. Introduction

A common property of self-gravitating systems carrying angular momentum is the natural propensity to flatten, which is indisputably 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 ingredients, 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 envelope 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 asynchronous motion between the two layers. Based upon numerical simulations, we reach a similar conclusion when the two components are inhomogeneous. This result, which contributes to the theory of figures, provides new tools for studying the formation 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 departure 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 equilibrium 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 theory 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 (Ωi) are the same, to a fully asynchronous regime, in which the Ωi’s are all different.

In all these theories, only oblate configurations have been obtained or reported1. 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 subcomponents.

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 ρi, Ωi, ai, and bi 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 = a1/a2, the mass-density jump at the interface α = ρ1/ρ2, the squared eccentricity of layer i, εi2 = 1 − (bi/ai)2, and the dimensionless squared rotation rate of layer i, Ω i 2 = Ω i 2 / ( 2 π G ρ 2 ) $ {\tilde{\varOmega}}_i^2 = \varOmega_i^2/(2\uppi G\rho_2) $, where G is the gravitational constant. Cases with q = 0, q = 1, or α = 1 have no special interest, as they all correspond to a single-component body. Asynchronous motion is allowed: Ω1 and Ω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 = q2ε12 − ε22, is small enough). The error is typically below 1% for |c|≲0.5.

thumbnail 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.

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 εi2, we have

{ Ω ~ 2 2 = 4 15 ε 2 2 2 15 ( α 1 ) q 3 ( 3 q 2 ε 1 2 5 ε 2 2 ) , α Ω ~ 1 2 = Ω ~ 2 2 + 2 15 ( α 1 ) [ 2 ( α 1 ) ε 1 2 + 5 ε 1 2 3 ε 2 2 ] . $$ \begin{aligned} \left\{ \begin{array}{ll} {\tilde{\varOmega }}_2^2 = \frac{4}{15}\varepsilon _2^2 - \frac{2}{15}(\alpha -1) q^3\left(3q^2\varepsilon _1^2-5\varepsilon _2^2\right), \\ \alpha {\tilde{\varOmega }}_1^2= {\tilde{\varOmega }}_2^2 + \frac{2}{15}(\alpha -1)\left[2(\alpha -1)\varepsilon _1^2 + 5\varepsilon _1^2-3\varepsilon _2^2\right].\end{array}\right. \end{aligned} $$(1)

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. ε i 2 < 0 $ \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 Ω 1 = Ω 2 $ {\varOmega_1 = \varOmega_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 Ω 1 Ω 2 $ \varOmega_1 \ne \varOmega_2 $. The gas pressure at the surface of the core varies quadratically with the cylindrical radius, R (type-V solutions), and α 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 α > 1 (i.e. the core is denser than its envelope), as this is the case in most problems. However, the case α < 1 is treated in Appendix A. The results of the two cases are largely the same; only the conditions of positivity of the Ω i 2 $ \varOmega_i^2 $’s differ (see below).

4. The conditions of positivity of the Ω i 2 $ \varOmega_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 Ω i 2 0 $ \varOmega_i^2\geq0 $ leads to the following set of inequalities:

{ [ 2 + 5 ( α 1 ) q 3 ] ε 2 2 3 ( α 1 ) q 5 ε 1 2 , [ 2 + ( α 1 ) ( 5 q 3 3 ) ] ε 2 2 ( α 1 ) [ 3 q 5 2 α 3 ] ε 1 2 , $$ \left \{\begin{aligned} &[2+5(\alpha -1) q^3] \varepsilon _2^2 \ge 3(\alpha -1) q^5\varepsilon _1^2, \\ &[2+(\alpha -1)(5q^3-3)] \varepsilon _2^2 \ge (\alpha -1)[3q^5-2\alpha -3]\varepsilon _1^2, \end{aligned}\right. $$(2)

which, in standard conditions (α > 1), becomes

{ ε 2 2 3 ( α 1 ) q 5 2 + 5 ( α 1 ) q 3 ε 1 2 , ε 1 2 2 + ( α 1 ) ( 5 q 3 3 ) ( α 1 ) ( 3 + 2 α 3 q 5 ) ε 2 2 , $$ \begin{aligned} \left\{ \begin{array}{ll} \varepsilon _2^2 \ge \frac{3(\alpha -1) q^5}{2+5(\alpha -1) q^3}\varepsilon _1^2, \\ \varepsilon _1^2 \ge -\frac{2+(\alpha -1) (5q^3-3)}{(\alpha -1)\left(3+2\alpha -3q^5\right)}\varepsilon _2^2, \end{array}\right. \end{aligned} $$(3)

as 3q5 − 2α − 3 < 0. We see from Eq. (3a) that an oblate core (ε12 > 0) is only compatible with an oblate envelope (ε22 > 0). OP configurations are thus ruled out. So we only consider solutions with prolate cores in the following. With the condition ε12 ≤ 0, Eq. (3b) leads to

[ 2 + ( α 1 ) ( 5 q 3 3 ) ] ε 2 2 0 . $$ \begin{aligned} \left[2 + (\alpha -1)(5q^3-3)\right]\varepsilon _2^2 \ge 0. \end{aligned} $$(4)

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

4.1. Prolate envelope solution ruled out

With ε22 ≤ 0, the above inequality reads 2 + (α − 1)(5q3 − 3)≤0, and so Eq. (2b) becomes

3 ( α 1 ) q 5 2 + 5 ( α 1 ) q 3 ε 1 2 ε 2 2 ( α 1 ) [ 3 q 5 2 ( α 1 ) 5 ] 2 + ( α 1 ) ( 5 q 3 3 ) ε 1 2 , $$ \begin{aligned} \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, \end{aligned} $$(5)

where we take Eq. (2a) into account. As ε12 ≤ 0 is assumed, we find

3 q 5 2 + 5 ( α 1 ) q 3 3 q 5 2 ( α 1 ) 5 2 + ( α 1 ) ( 5 q 3 3 ) , $$ \begin{aligned} \frac{3q^5}{2+5(\alpha -1) q^3} \ge \frac{3q^5-2(\alpha -1)-5}{2+(\alpha -1)(5q^3-3)}, \end{aligned} $$(6)

which also reads

4 ( α 1 ) + 10 ( α 1 ) 2 q 3 + 10 + 25 ( α 1 ) q 3 9 ( α 1 ) q 5 . $$ \begin{aligned} 4(\alpha -1) + 10(\alpha -1)^2q^3 + 10 + 25(\alpha -1) q^3 \le 9(\alpha -1) q^5. \end{aligned} $$(7)

As α − 1 > 0 and given 0 < q < 1, we have

25 ( α 1 ) q 3 9 ( α 1 ) q 5 , $$ \begin{aligned} {25(\alpha -1) q^3 \ge 9(\alpha -1) q^5}, \end{aligned} $$(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 ε22 ≥ 0, we see that Eq. (2a) is automatically verified. Furthermore, we have 2 + (α − 1)(5q3 − 3)≥0 in this case, which implies

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

and

ε 2 2 ( α 1 ) [ 3 q 5 2 ( α 1 ) 5 ] 2 + ( α 1 ) ( 5 q 3 3 ) ε 1 2 ε 2 , lim 2 . $$ \begin{aligned} \varepsilon _2^2 \ge \frac{(\alpha -1)\left[3q^5-2(\alpha -1)-5\right]}{2+(\alpha -1)(5q^3-3)}\varepsilon _1^2 \equiv \varepsilon _{2,\lim }^2. \end{aligned} $$(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 qlim = 0 for α = 5/3, and qlim = (3/5)1/3 ≈ 0.843 for α → ∞. There is also a lower limit for the eccentricity of the envelope’s surface, that is, ε2, lim2. 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 Ωi’s, not on their relative magnitude. As outlined in Huré (2022) in relation to OO configurations only, states of global rotation, that is, Ω 1 = Ω 2 $ \varOmega_1=\varOmega_2 $, can occur either for α = 1, which corresponds to a single-layer object, or for a special mass-density jump, αC, namely

α α C = 1 + 5 ( ε 1 2 ε 2 2 ) 2 ε 1 2 + 3 q 5 ε 1 2 5 q 3 ε 2 2 . $$ \begin{aligned} \alpha \equiv \alpha _{\rm C}= 1 + \frac{5(\varepsilon _1^2-\varepsilon _2^2)}{2\varepsilon _1^2+3q^5\varepsilon _1^2-5q^3\varepsilon _2^2}. \end{aligned} $$(11)

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

Ω ~ 2 = 4 3 ε 1 2 ε 2 2 2 ε 1 2 + 3 q 5 ε 1 2 5 q 3 ε 2 2 ε 1 2 + 2 3 ε 1 2 2 5 ε 2 2 0 . $$ \begin{aligned} {\tilde{\varOmega }}^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. \end{aligned} $$(12)

As proved above, only OO and PO configurations are permitted. So, assuming α ≠ 1, ε12 ≤ 0, and ε22 ≥ 0, Eq. (12) implies

( 20 + 15 q 3 ) ε 1 4 ( 16 + 25 q 3 + 9 q 5 ) ε 1 2 ε 2 2 + 15 q 3 ε 2 4 0 . $$ \begin{aligned} (20+15q^3)\varepsilon _1^4 - (16+25q^3+9q^5)\varepsilon _1^2\varepsilon _2^2 + 15q^3\varepsilon _2^4 \le 0. \end{aligned} $$(13)

The left-hand side of this inequality is a second-degree polynomial in ε12. The two roots are

ε 1 ± 2 = 16 + 25 q 3 + 9 q 5 10 ( 4 + 3 q 3 ) ε 2 2 [ 1 ± 1 300 q 3 ( 4 + 3 q 3 ) ( 16 + 25 q 3 + 9 q 5 ) 2 ] , $$ \begin{aligned} \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], \end{aligned} $$(14)

and it is easy to see that both roots are positive. As the prefactor of ε14 in the second-degree polynomial is positive, Eq. (12) is verified only between the two roots ε (i.e. for ε12 ≥ 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, Ω 1 Ω 2 $ \varOmega_1 \ne \varOmega_2 $. The only exception is the case α = 1, where the distinction between the prolate core and the oblate envelope is impossible.

5. Results

5.1. Rotation rates

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

Δ = Ω 1 ~ Ω 2 ~ , $$ \begin{aligned} \varDelta = \frac{\tilde{\varOmega _1}}{\tilde{\varOmega _2}}, \end{aligned} $$(15)

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

Δ 2 = 1 α [ 1 ( α 1 ) 3 ε 2 2 2 ( α 1 ) ε 1 2 5 ε 1 2 2 ε 2 2 ( α 1 ) q 3 ( 3 q 2 ε 1 2 5 ε 2 2 ) ] , $$ \begin{aligned} \varDelta ^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], \end{aligned} $$(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 Δ versus α for two pairs (ε12, ε22). We see that all curves converge to Δ = 1 $ \varDelta=1 $ for α = 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 Ω 1 0 $ \varOmega_1 \rightarrow 0 $ for some special values of the parameters. This is expected from the equation set. We see from Eq. (1) that the condition Ω 1 2 = 0 $ {\tilde{\varOmega}}_{1}^2=0 $ can be obtained for a special value of mass-density contrast αlim − 1 > 0, which is the root of a second-degree polynomial, namely

ε 1 2 ( α lim 1 ) 2 + 5 ε 1 2 3 ε 2 2 q 3 ( 3 q 2 ε 1 2 5 ε 2 2 ) 2 ( α lim 1 ) + ε 2 2 = 0 . $$ \begin{aligned} \varepsilon _{1}^2(\alpha _{\rm 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 _{\rm lim}{-}1)+\varepsilon _{2}^2=0. \end{aligned} $$(17)

thumbnail Fig. 2.

Rotational discontinuity, Δ , $ \varDelta, $ versus α defined by Eq. (16) for ε22 = 0.01 (top) and ε22 = 0.1 (bottom). In both cases, ε12 = −0.01 corresponds to a prolate core surrounded by an oblate envelope. See Tables 1 and 2 for configurations A and B.

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

ε 1 2 [ ( α 1 ) ( 5 q 3 3 ) + 2 ] ε 2 2 ( α 1 ) [ ( 5 3 q 5 ) + 2 ( α 1 ) ] ε 1 , m i n 2 . $$ \begin{aligned} \varepsilon _1^2 \ge -\frac{[(\alpha -1)(5q^3-3)+2]\varepsilon _{2}^2}{(\alpha -1)[(5-3q^5) + 2(\alpha -1)]} \equiv \varepsilon ^2_{1, \mathrm min}. \end{aligned} $$(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, ε1, ε2, and α. 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 α = 1.25, q = 0.7, and ε22 = 0.01, and in Table 2 for α = 2.5, q = 0.8, and ε22 = 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 & Giuli 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.

Table 1.

Example of a configuration and the main output data obtained with a prolate core and an oblate envelope.

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 ε22 = 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.

Table 2.

Same as Table 1 but for configuration B.

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.

  • 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, ρ, namely p ∝ ρ1 + 1/n, and n is the polytropic index. Table 3 presents two new configurations, namely:

  • Configuration C, which resembles configuration A but with (n1, n2) = (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 (n1, n2) = (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.

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.

Figures 3 and 4 show the mass density and the deviations of the surface bounding layers from spheres. The right panels shows the deviation ΔZi(R) of the true surfaces from spheres, namely

Δ Z i = Z i ( R ) a i 2 R 2 , i { 1 , 2 } , $$ \begin{aligned} \Delta Z_i = Z_i(R) - \sqrt{a_i^2-R^2}, \qquad i\in \{1,2\}, \end{aligned} $$(19)

thumbnail Fig. 3.

Equilibrium solution computed with the DROP code for configuration C: the normalized mass-density, ρ/ρc, map in colour code (left) and the deviation from sphericity (right). See also Table 3.

thumbnail Fig. 4.

Same as Fig. 3 but for configuration D.

where Zi 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, J2n, can vanish (when n is an odd number). The first coefficient, J2, is of special interest because it is the main deviation from sphericity and therefore of interest for dynamical studies. For spheroidal heteroeoids, the J2n’s are analytical (Heiskanen & Moritz 1967; Cisneros-Parra et al. 2019; Basillais & Huré 2023). We find that J2 = 0 when

ε 2 2 1 ε 2 2 + ( α 1 ) q 5 ε 1 2 1 ε 1 2 = 0 , $$ \begin{aligned} \varepsilon _2^2\sqrt{1-\varepsilon _2^2}+(\alpha -1)q^5 \varepsilon _1^2\sqrt{1-\varepsilon _1^2}=0, \end{aligned} $$(20)

which is a formally cubic equation in ε12. In the actual conditions of slow rotation, this implies

ε 1 2 ε 2 2 ( α 1 ) q 5 0 , $$ \begin{aligned} \varepsilon _1^2 \approx -\frac{\varepsilon _2^2}{(\alpha -1)q^5} \le 0, \end{aligned} $$(21)

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 α’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 (ε12, ε22) plane, with −1 ≤ εi2 ≤ 1, for several values of q and α. 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 α 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|∼0.5), but the analytical solution is still in good agreement with the numerical experiment.

Table 4.

Same as Table 1 but for a fast-rotating envelope.

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, Δ 1 $ \varDelta \rightarrow 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 κ | Ω 1 Ω 2 | m R m $ \kappa |\varOmega_1-\varOmega_2|^mR^m $, where κ 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 | Ω 1 Ω 2 | $ |\varOmega_1-\varOmega_2| $ between the spheres undergoes an exponential decrease with a characteristic time, τ, given by

1 τ = 8 3 π a 1 4 κ I 1 + I 2 I 1 I 2 , $$ \begin{aligned} \frac{1}{\tau } = \frac{8}{3}\uppi a_1^4 \kappa \frac{I_1+I_2}{I_1I_2}, \end{aligned} $$(22)

where I1 and I2 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 ε12, which is enough to obtain an approximate order of magnitude of τ. The main difficulty of this estimation of the characteristic time is the value of κ. For a fluid core embedded in a precessing mantle, we have κ I 1 ν Ω 2 / a 1 5 $ \kappa \propto I_1\sqrt{\nu \varOmega_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, ν, covers 13 orders of magnitude (Lumb & Aldridge 1991), the best estimate being ν ∼ 10−6 m2 ⋅ s−1 (Gans 1972; Poirier 1988). To evaluate the order of magnitude of τ, we assumed that the expression of Busse (1968) is valid in our case and that Ω 2 $ \varOmega_2 $ is the initial value of the rotation rate of the envelope, namely

τ 1 2.62 a 1 ( 2 π G ρ 2 ν 2 ) 1 / 4 1 Ω ~ 2 1 q 5 1 + ( α 1 ) q 5 . $$ \begin{aligned} \tau \approx \frac{1}{2.62}\frac{a_1}{(2\uppi G \rho _2 \nu ^2)^{1/4}} \frac{1}{\sqrt{{\tilde{\varOmega }}_{2}}}\frac{1-q^5}{1+(\alpha -1)q^5}. \end{aligned} $$(23)

For the parameters of configuration A in Table 1, taking ν = 10−6 m2 ⋅ s−1 and assuming a planetary core with a1 = 1000 km and a mass density of the envelope ρ2 = 4000 kg ⋅ m−3 (typically silicates), we have τ ∼ 103 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.


1

Note that prolate subsystems were already considered by D’Alembert (1747) in the context of a rigid core supporting a shallow ocean.

Acknowledgments

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

  1. Basillais, B., & Huré, J.-M. 2021, MNRAS, 506, 3773 [NASA ADS] [CrossRef] [Google Scholar]
  2. Basillais, B., & Huré, J.-M. 2023, MNRAS, 520, 1504 [NASA ADS] [CrossRef] [Google Scholar]
  3. Busse, F. H. 1968, J. Fluid Mech., 33, 739 [NASA ADS] [CrossRef] [Google Scholar]
  4. Cai, M. J., & Taam, R. E. 2010, ApJ, 709, L79 [NASA ADS] [CrossRef] [Google Scholar]
  5. Cisneros-Parra, J. U., Martinez-Herrera, F. J., & Montalvo-Castro, J. D. 2019, ApJS, 241, 8 [NASA ADS] [CrossRef] [Google Scholar]
  6. Clairaut, A. C. 1743, Théorie de la figure de la Terre tirée des principes de l’hydrostatique (Paris: David Fils) [Google Scholar]
  7. Cox, J. P., & Giuli, R. T. 1968, Principles of stellar structure Volume I: Physical principles (New York: Gordon and Breach) [Google Scholar]
  8. D’Alembert, J. 1747, Réflexions sur la cause générale des vents (Paris: David l’aîné) [Google Scholar]
  9. Fiege, J. D., & Pudritz, R. E. 2000, ApJ, 534, 291 [NASA ADS] [CrossRef] [Google Scholar]
  10. Fujisawa, K., & Eriguchi, Y. 2014, MNRAS, 438, L61 [Google Scholar]
  11. Gans, R. F. 1972, J. Geophys. Res., 77, 360 [NASA ADS] [CrossRef] [Google Scholar]
  12. Hachisu, I. 1986, ApJS, 61, 479 [NASA ADS] [CrossRef] [Google Scholar]
  13. Hamy, M. 1889, Annales de l’Observatoire de Paris, 19, F.1. [Google Scholar]
  14. Hamy, M. 1890, Journal de mathématiques pures et appliquées 4e série, 6, 69 [Google Scholar]
  15. Heiskanen, W., & Moritz, H. 1967, Physical Geodesy (San Fransisco: W.H. Freeman) [Google Scholar]
  16. Huré, J.-M. 2022, MNRAS, 512, 4031 [CrossRef] [Google Scholar]
  17. Huré, J.-M., & Hersant, F. 2017, MNRAS, 464, 4761 [CrossRef] [Google Scholar]
  18. Huré, J.-M., & Staelen, C. 2024, Phys. Rev. D, 110 [Google Scholar]
  19. Huré, J.-M., Hersant, F., & Nasello, G. 2018, MNRAS, 475, 63 [CrossRef] [Google Scholar]
  20. Kadam, K., Motl, P. M., Frank, J., Clayton, G. C., & Marcello, D. C. 2016, MNRAS, 462, 2237 [NASA ADS] [CrossRef] [Google Scholar]
  21. Kawamura, T., Taniguchi, K., Yoshida, S., & Eriguchi, Y. 2011, MNRAS, 416, L75 [CrossRef] [Google Scholar]
  22. Lander, S. K., & Jones, D. I. 2009, MNRAS, 395, 2162 [NASA ADS] [CrossRef] [Google Scholar]
  23. Lumb, L. I., & Aldridge, K. D. 1991, J. Geomagn. Geoelec., 43, 93 [NASA ADS] [CrossRef] [Google Scholar]
  24. Montalvo, D., Martinez, F. J., & Cisneros, J. 1983, Rev. Mex. Astron. Astrofis., 5, 293 [Google Scholar]
  25. Moulton, E. J. 1916, Trans. Am. Math. Soc., 17, 100 [CrossRef] [Google Scholar]
  26. Oganesyan, R. S., & Abramyan, M. G. 1972, Astrophysics, 8, 352 [Google Scholar]
  27. Poirier, J. P. 1988, Geophys. J. Int., 92, 99 [NASA ADS] [CrossRef] [Google Scholar]
  28. Rambaux, N., Chambat, F., & Castillo-Rogez, J. C. 2015, A&A, 584, A127 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Véronet, A. 1912, Journal de mathématiques pures et appliquées, 8, 331 [Google Scholar]
  30. Wavre, R. 1925, Archive des sciences physiques et naturelles 5me période, 7, 133 [Google Scholar]

Appendix A: The case α < 1

Considering α < 1, the system of inequalities (1)

{ ε 1 2 2 5 ( 1 α ) q 3 3 ( 1 α ) q 5 ε 2 2 , 2 ( 1 α ) [ 5 q 3 3 ] ( 1 α ) [ 5 2 ( 1 α ) 3 q 5 ] ε 2 2 ε 1 2 $$ \begin{aligned} \left\{ \begin{array}{ll} \varepsilon _1^2\ge -\frac{2-5(1-\alpha )q^3}{3(1-\alpha )q^5}\varepsilon _2^2,\\ \frac{2-(1-\alpha )[5q^3-3]}{(1-\alpha )\left[5-2(1-\alpha )-3q^5\right]}\varepsilon _2^2 \ge \varepsilon _1^2 \end{array}\right. \end{aligned} $$(A.1)

Combining the two inequalities of the system, we readily obtain

2 ( 1 α ) [ 5 q 3 3 ] ( 1 α ) [ 5 2 ( 1 α ) 3 q 5 ] ε 2 2 2 5 ( 1 α ) q 3 3 ( 1 α ) q 5 ε 2 2 $$ \begin{aligned} \frac{2-(1-\alpha )[5q^3-3]}{(1-\alpha )\left[5-2(1-\alpha )-3q^5\right]}\varepsilon _2^2 \ge -\frac{2-5(1-\alpha )q^3}{3(1-\alpha )q^5}\varepsilon _2^2 \end{aligned} $$(A.2)

As before, we now have two cases to consider on the sign of ε22.

A.1. Prolate envelope ruled out

We first assume ε22 < 0. We thus have

2 5 ( 1 α ) q 3 3 ( 1 α ) q 5 2 ( 1 α ) [ 5 q 3 3 ] ( 1 α ) [ 5 2 ( 1 α ) 3 q 5 ] $$ \begin{aligned} -\frac{2-5(1-\alpha )q^3}{3(1-\alpha )q^5}\ge \frac{2-(1-\alpha )[5q^3-3]}{(1-\alpha )\left[5-2(1-\alpha )-3q^5\right]} \end{aligned} $$(A.3)

which leads to

10 [ 4 + 25 q 3 9 q 5 ] ( 1 α ) + 10 ( 1 α ) 2 0 . $$ \begin{aligned} 10 - [4+25q^3-9q^5](1-\alpha ) + 10(1-\alpha )^2 \le 0. \end{aligned} $$(A.4)

As 4 + 25q3 − 9q5 < 20, the discriminant of the second-degree polynomial in 1 − α in the left hand side is negative, except for q = 1. In this case, however, we see that the inequality becomes α2 ≤ 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 ε12 ≤ 0. Thus, Eq. (A.1b) is automatically verified and Eq. (A.1a) reads

ε 2 2 3 ( 1 α ) q 5 2 5 ( 1 α ) q 3 ε 1 2 , $$ \begin{aligned} \varepsilon _2^2 \ge -\frac{3(1-\alpha )q^5}{2-5(1-\alpha )q^3}\varepsilon _1^2, \end{aligned} $$(A.5)

where we have necessarily 2 − 5(1 − α)q3 ≥ 0, meaning that we have an upper limit on q (whereas we had a lower limit in the α > 1 case), namely

q ( 2 5 ( 1 α ) ) 1 / 3 q lim . $$ \begin{aligned} q \le \left(\frac{2}{5(1-\alpha )}\right)^{1/3} \equiv q_{\lim }^{\prime }. \end{aligned} $$(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.

All Tables

Table 1.

Example of a configuration and the main output data obtained with a prolate core and an oblate envelope.

Table 2.

Same as Table 1 but for configuration B.

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.

Table 4.

Same as Table 1 but for a fast-rotating envelope.

All Figures

thumbnail 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.

In the text
thumbnail Fig. 2.

Rotational discontinuity, Δ , $ \varDelta, $ versus α defined by Eq. (16) for ε22 = 0.01 (top) and ε22 = 0.1 (bottom). In both cases, ε12 = −0.01 corresponds to a prolate core surrounded by an oblate envelope. See Tables 1 and 2 for configurations A and B.

In the text
thumbnail Fig. 3.

Equilibrium solution computed with the DROP code for configuration C: the normalized mass-density, ρ/ρc, map in colour code (left) and the deviation from sphericity (right). See also Table 3.

In the text
thumbnail Fig. 4.

Same as Fig. 3 but for configuration D.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.