<p>I actually just found an article that explains such a calculation in detail, using a matrix propagator method: <a href="https://www.aanda.org/articles/aa/pdf/2018/12/aa34181-18.pdf" rel="nofollow noreferrer">Padovan et al. 2018</a>
This is a calculation for gas giants or bodies in hydrostatic equilibrium. The calculation of the Love numbers for terrestrial planets that have viscoelastic behaviour requires the rheology (which includes models for the rigidity). I did not find the equations that need to be solved in that case, but a similar matrix propagator method can be used apparently.</p>
<p>I will not repeat the matrix propagation method here, as it is quite well elaborated in the paper I just mentioned, except for the equation you have to start from, which is explained in another article, but looks a bit too complicated for me at this point.</p>

CC BY-SA 4.0
2021-04-24T14:42:50.650Z
2021-04-29T13:54:12.497Z
2021-04-29T13:54:12.497Z
theWrongAlice
<p>I have a question about planetology. I'm not sure if it belongs here, but the astronomy stack exchange seemed a bit odd, too.</p>
<p>I would like to enhance my understanding of Love numbers. Love numbers quantify the deviation from equilibrium tides on a planet. For a homogeneous body, you can write, for example
<span class="math-container">$$ k_2 = \frac{3}{2}\Big(1 + \frac{19}{2}\frac{\mu}{\rho gR}\Big)^{-1}$$</span>
where <span class="math-container">$\mu$</span> is the rigidity, <span class="math-container">$\rho$</span> is the density, <span class="math-container">$R$</span> is the radius of the planet and <span class="math-container">$g$</span> is the gravity constant. I have been taught that for non-homogeneous bodies, e.g. planets with a core etc, there is no simple analytical expression and you need numerical techniques. I have been reading a <a href="https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2018JE005569" rel="nofollow noreferrer">paper</a> about exactly that, but I don't see what equations are used to obtain <span class="math-container">$k_2$</span>. In this article, the interior of the planet is described by Maxwell models and a pseudo-period Andrade model. These rheological models are used in conjunction with so called structural models for the interior. It's not clear to me what these are, either. I think they are equations for pressure and temperature. The article then says: <em>&quot;Based on these, the tidal Love number <span class="math-container">$k_2$</span> is calculated and compared against measurement.&quot;</em></p>
<p>Does anyone know what equation needs to be solved to obtain the <span class="math-container">$k_2$</span> Love number? Or what steps need to be taken to do such a calculation? Or if you know a relevant reference, anything is greatly appreciated, because I don't know where to start.</p>
<p>EDIT: I found a calculation for planets in hydrostatic equilibrium (see my own answer below). If someone knows the equations for a terrestrial planet with viscoelastic behaviour, I would still be interested in that.</p>

CC BY-SA 4.0
2021-04-24T11:51:47.130Z
This would be in fact also appropriate on the astro.se, as a part of (exo)planetary science is concerned with love numbers as well. Here is one paper that goes through the entire apparatus to calculate $k_2$ for gas giants: https://academic.oup.com/mnras/article/450/4/3952/1749966 and if you're looking for a ton of references to start, the protostars&planets 6 chapter here should be useful too: https://www2.mpia-hd.mpg.de/homes/ppvi/chapter/baraffe.pdf
AtmosphericPrisonEscape
<planetary-science><numerical-modelling><tides>
Calculating Love numbers (tidal deformation) for non-homogeneous planets
theWrongAlice
# Calculating Love numbers (tidal deformation) for non-homogeneous planets

(Asked by: username_0 on 2021-04-23)

I have a question about planetology. I'm not sure if it belongs here, but the astronomy stack exchange seemed a bit odd, too.

I would like to enhance my understanding of Love numbers. Love numbers quantify the deviation from equilibrium tides on a planet. For a homogeneous body, you can write, for example
$$ k\_2 = \frac{3}{2}\Big(1 + \frac{19}{2}\frac{\mu}{\rho gR}\Big)^{-1}$$
where $\mu$ is the rigidity, $\rho$ is the density, $R$ is the radius of the planet and $g$ is the gravity constant. I have been taught that for non-homogeneous bodies, e.g. planets with a core etc, there is no simple analytical expression and you need numerical techniques. I have been reading a [paper](https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2018JE005569) about exactly that, but I don't see what equations are used to obtain $k\_2$. In this article, the interior of the planet is described by Maxwell models and a pseudo-period Andrade model. These rheological models are used in conjunction with so called structural models for the interior. It's not clear to me what these are, either. I think they are equations for pressure and temperature. The article then says: *"Based on these, the tidal Love number $k\_2$ is calculated and compared against measurement."*

Does anyone know what equation needs to be solved to obtain the $k\_2$ Love number? Or what steps need to be taken to do such a calculation? Or if you know a relevant reference, anything is greatly appreciated, because I don't know where to start.

EDIT: I found a calculation for planets in hydrostatic equilibrium (see my own answer below). If someone knows the equations for a terrestrial planet with viscoelastic behaviour, I would still be interested in that.

**Tags:** planetary-science, numerical-modelling, tides

> **Comments:**
> **username_1**: This would be in fact also appropriate on the astro.se, as a part of (exo)planetary science is concerned with love numbers as well. Here is one paper that goes through the entire apparatus to calculate $k_2$ for gas giants: https://academic.oup.com/mnras/article/450/4/3952/1749966 and if you're looking for a ton of references to start, the protostars&planets 6 chapter here should be useful too: https://www2.mpia-hd.mpg.de/homes/ppvi/chapter/baraffe.pdf  

---

## 1 answers

### Answer by username_0 on 2021-04-24. Score: 2

I actually just found an article that explains such a calculation in detail, using a matrix propagator method: [Padovan et al. 2018](https://www.aanda.org/articles/aa/pdf/2018/12/aa34181-18.pdf)
This is a calculation for gas giants or bodies in hydrostatic equilibrium. The calculation of the Love numbers for terrestrial planets that have viscoelastic behaviour requires the rheology (which includes models for the rigidity). I did not find the equations that need to be solved in that case, but a similar matrix propagator method can be used apparently.

I will not repeat the matrix propagation method here, as it is quite well elaborated in the paper I just mentioned, except for the equation you have to start from, which is explained in another article, but looks a bit too complicated for me at this point.


---

