<p>I have to admit that power-laws (in general) used to be my <a href="https://en.m.wikipedia.org/wiki/Shtick" rel="nofollow noreferrer">shtick</a> so I am happy to shed some light on their general importance in physics which obviously also hold for astronomy.</p>
<p>The main idea of a <a href="https://en.wikipedia.org/wiki/Power_law" rel="nofollow noreferrer">power law</a> is nicely written in Wikipedia, but the essential part is the I highlighted in the following quote:</p>
<blockquote>
<p>[A Powerlaw is] a functional relationship between two quantities, where a relative change in one quantity results in a proportional relative change in the other quantity, <strong>independent of the initial size of those quantities</strong></p>
</blockquote>
<p>The (mathematically) nice part is exactly this, that over a <em>really large range</em> of the x-values, the y-values follow the same dependency. Usually, &quot;really large range&quot; means values spreading over 3 to even 10 powers of ten.</p>
<p>When fitting power laws in practice, there are effects on the sides of the scale, meaning for small and large <span class="math-container">$x$</span>-values, where usually attributed to &quot;finite size effects&quot;.</p>
<h3>Further reading</h3>
<ul>
<li>There is the research flavor <a href="https://en.wikipedia.org/wiki/Self-organized_criticality" rel="nofollow noreferrer">Self-organized criticality (SOC)</a> which <em>&quot;is considered to be one of the mechanisms by which complexity arises in nature&quot;</em>. One of the founding fathers (and mothers) was the late <a href="https://en.wikipedia.org/wiki/Per_Bak" rel="nofollow noreferrer">Per Bak</a>. His book <em><a href="https://books.google.de/books/about/How_Nature_Works.html?id=Bth4QgAACAAJ&amp;redir_esc=y" rel="nofollow noreferrer">How Nature Works</a></em> is accessible for laypersons and <a href="https://en.wiktionary.org/wiki/IMHO" rel="nofollow noreferrer">IMHO</a> a good read.</li>
</ul>

CC BY-SA 4.0
2021-06-18T07:38:18.213Z
2021-06-21T06:10:42.247Z
2021-06-21T06:10:42.247Z
B--rian
<p>The <a href="https://en.wikipedia.org/wiki/Polytrope" rel="nofollow noreferrer">polytrope model</a> of stars made early numerical calculations of stellar structure possible using anything from mechanical calculators to early electronic computers, and in certain cases, even analytically!</p>
<blockquote>
<p>In astrophysics, a polytrope refers to a solution of the Lane–Emden equation in which the pressure depends upon the density in the form</p>
<p><span class="math-container">$$P=K\rho ^{{(n+1)/n}}$$</span></p>
<p>where <span class="math-container">$P$</span> is pressure, <span class="math-container">$\rho$</span> is density and <span class="math-container">$K$</span> is a constant of proportionality. The constant <span class="math-container">$n$</span> is known as the polytropic index; note however that the polytropic index has an alternative definition as with <span class="math-container">$n$</span> as the exponent.</p>
</blockquote>
<p>See answers to:</p>
<ul>
<li>Math SE: <a href="https://math.stackexchange.com/q/4148101/284619">Could there be exact solutions to the Lane-Emden equation for real n≥0 other than 0, 1, or 5?</a></li>
<li><a href="https://astronomy.stackexchange.com/q/29743/7982">Has stellar evolution ever been modeled analytically?</a> especially <a href="https://astronomy.stackexchange.com/a/29820/7982">@ProfRob's answer</a>!</li>
<li><a href="https://astronomy.stackexchange.com/q/16633/7982">Does Eddington&#39;s variable polytropic index better fit data from the Standard Solar Model?</a></li>
<li><a href="https://astronomy.stackexchange.com/q/43913/7982">Gravity inside a star?</a></li>
</ul>

CC BY-SA 4.0
2021-06-21T02:05:39.857Z
2021-06-21T02:05:39.857Z
uhoh
<h3>Scale invariance and self-similarity</h3>
<p>Power laws basically mean that there is no preferred scale, i.e. that a physical property is <em>scale invariant</em>. Any deviation from a power law means that the Universe somehow thinks that the scale where it breaks down has some special significance. In other words, a power law describes <em>self-similarity</em><span class="math-container">$^\dagger$</span>.</p>
<p>You can see this mathematically by considering a scaling of the independent variable <span class="math-container">$x$</span> in a power law <span class="math-container">$y(x)=x^n$</span> by a constant factor <span class="math-container">$\lambda$</span>; the result is just a scaling of the dependent variable by another factor
<span class="math-container">$$
y(\lambda x) = (\lambda x)^n = \lambda^n x^n = \lambda^n y.
$$</span></p>
<p>Of course, nothing is infinitely self-similar, but many processes are to certain degree. And if you have a process that depends on many sub-processes that are themselves described by power laws within some domain, the result will evolve toward a power law as well.</p>
<h3>Power laws in astronomy…</h3>
<p>In astronomy, you often see that power laws are a good description of some relation out to a certain scale, after which you see an &quot;exponential cut-off&quot;, i.e. where the relation transitions to the steeper exponential.</p>
<p>A notable example of this is the <em>halo mass function</em> (HMF) of structures in the Universe, where the distributions of (dark matter halo) masses are well described by power law from the lowest masses out to a &quot;characteristic&quot; mass, after which the number density of object decreases rapidly with mass (<a href="https://ui.adsabs.harvard.edu/abs/1974ApJ...187..425P" rel="nofollow noreferrer">Press &amp; Schechter 1974</a>). In fact, many relations in astronomy can be traced back to this hierarchical formalism which has its origin in the assumption that overdensities are Gaussian distributed. Scaling relations between masses, sizes, star formation rates, luminosities, feedback mechanisms, etc. then result in new power laws.</p>
<p>Assuming a constant mass-to-light ratio <span class="math-container">$M/L$</span> of galaxies, the <em>luminosity function</em> (LF) of galaxies should follow closely the MF. However, observationally we find that the LF is suppressed at both higher and lower mass scales than the characteristic mass. The reason is thought to be feedback preventing star formation by heating the gas and/or pushing it out of the galaxy; at the high-mass end, you have powerful active galactic nuclei (AGN), while at the low-mass end, where the gravitational potential is smaller, stellar activity can do the job.</p>
<p><a href="https://i.sstatic.net/i4qd4.png" rel="nofollow noreferrer"><img src="https://i.sstatic.net/i4qd4.png" alt="LF_HMF" /></a>
<sup>Predicted LF assuming a <a href="https://ui.adsabs.harvard.edu/abs/2001MNRAS.321..372J" rel="nofollow noreferrer">Jenkins et al (2001)</a> with a constant <span class="math-container">$M/L$</span> (cyan), and observed <span class="math-container">$K$</span> band LF (yellow). Figure from <a href="https://ui.adsabs.harvard.edu/abs/2003ApJ...599...38B" rel="nofollow noreferrer">Benson et al. (2003)</a> with my own annotations.</sup></p>
<p>The Milky Way, it turns out, is rather close to being such a characteristic mass and luminosity, which is why you sometimes hear that the Milky Way is a &quot;typical&quot; galaxy.</p>
<p>The LF is typically described by a <a href="https://en.wikipedia.org/wiki/Luminosity_function_(astronomy)" rel="nofollow noreferrer">Schechter function</a>. In recent years though we are seeing evidence that the exponential cut-off is a bit too drastic, and that the high-mass end may be better described by another power law, albeit with a steep index, i.e. a so-called <em>broken power law</em> (see e.g. <a href="https://ui.adsabs.harvard.edu/abs/2018ApJ...855..105O" rel="nofollow noreferrer">Oesch et al. 2018</a>).</p>
<h3>…and in nature in general</h3>
<p>Power laws are common not only in astronomy, but in all of nature, and even in &quot;artificial&quot; settings, such as linguistics (the occurrence of words in a language), economics (e.g. income distributions), and sociology (e.g. city sizes).</p>
<p>In some cases, the reason can be attributed to the fact that even non-power laws often behave like power laws in limited intervals. Other functions that occur in nature, such as <span class="math-container">$\exp(x)$</span> and <span class="math-container">$\sin(x)$</span>, can be Taylor-expanded, and near &quot;critical points&quot; the first term dominates.</p>
<p>I'm actually not sure if the reason for the ubiquity of power laws in nature is entirely understood. I think you can explain it in <em>many</em> cases, but not in general. One example that <em>can</em> be explained, mathematically, is that if stochastic processes with <em>exponential</em> growth in expectation are &quot;killed&quot; (or observed) randomly, the distribution of the killed or observed state exhibits power-law behavior in one or both tails (<a href="https://ui.adsabs.harvard.edu/abs/2002PhRvE..66f7103R" rel="nofollow noreferrer">Reed &amp; Hughes 2002</a>).</p>
<h3>But beware!</h3>
<p>That said, I think there sometimes is a tendency to assume a power law relation between variables where none is necessarily expected. A power law is a linear relation if you take the logarithm of both <span class="math-container">$x$</span> and <span class="math-container">$y$</span>, and probably the fact that deviations from some relation are suppressed in a log-log plot is sometimes a bit too appealing…</p>
<hr />
<p><span class="math-container">$\dagger$</span><sub>To be fair, scale invariance and self-similarity is not exactly the same, in that the former describes a continuous scaling, while the latter includes discrete scaling such as fractals which are identical to themselves only in certain steps and certain regions.</sub></p>

CC BY-SA 4.0
2021-06-19T09:49:28.273Z
2021-06-21T10:02:32.333Z
2021-06-21T10:02:32.333Z
pela
<p>I often see astronomers fitting data to power laws. What about power laws makes them so useful in astronomy? Why are so many astronomical observations well-fit by power laws? I know it's a relation between two quantities, but why is a power law the go-to model?</p>

CC BY-SA 4.0
2021-06-18T03:35:56.827Z
Different but related question in History of Science and Mathematics SE: [When was the earliest use of log-log plots to demonstrate power-law behavior?](https://hsm.stackexchange.com/q/3580/3869)
uhoh
<astrophysics><data-analysis><mathematics>
Role of power laws in astronomy?
blackhole
# Role of power laws in astronomy?

(Asked by: username_0 on 2021-06-17)

I often see astronomers fitting data to power laws. What about power laws makes them so useful in astronomy? Why are so many astronomical observations well-fit by power laws? I know it's a relation between two quantities, but why is a power law the go-to model?

**Tags:** astrophysics, data-analysis, mathematics

> **Comments:**
> **username_1**: Different but related question in History of Science and Mathematics SE: [When was the earliest use of log-log plots to demonstrate power-law behavior?](https://hsm.stackexchange.com/q/3580/3869)  

---

## 3 answers

### Answer by username_3 on 2021-06-19. Score: 13 | Accepted answer

### Scale invariance and self-similarity

Power laws basically mean that there is no preferred scale, i.e. that a physical property is *scale invariant*. Any deviation from a power law means that the Universe somehow thinks that the scale where it breaks down has some special significance. In other words, a power law describes *self-similarity*$^\dagger$.

You can see this mathematically by considering a scaling of the independent variable $x$ in a power law $y(x)=x^n$ by a constant factor $\lambda$; the result is just a scaling of the dependent variable by another factor
$$
y(\lambda x) = (\lambda x)^n = \lambda^n x^n = \lambda^n y.
$$

Of course, nothing is infinitely self-similar, but many processes are to certain degree. And if you have a process that depends on many sub-processes that are themselves described by power laws within some domain, the result will evolve toward a power law as well.

### Power laws in astronomy…

In astronomy, you often see that power laws are a good description of some relation out to a certain scale, after which you see an "exponential cut-off", i.e. where the relation transitions to the steeper exponential.

A notable example of this is the *halo mass function* (HMF) of structures in the Universe, where the distributions of (dark matter halo) masses are well described by power law from the lowest masses out to a "characteristic" mass, after which the number density of object decreases rapidly with mass ([Press & Schechter 1974](https://ui.adsabs.harvard.edu/abs/1974ApJ...187..425P)). In fact, many relations in astronomy can be traced back to this hierarchical formalism which has its origin in the assumption that overdensities are Gaussian distributed. Scaling relations between masses, sizes, star formation rates, luminosities, feedback mechanisms, etc. then result in new power laws.

Assuming a constant mass-to-light ratio $M/L$ of galaxies, the *luminosity function* (LF) of galaxies should follow closely the MF. However, observationally we find that the LF is suppressed at both higher and lower mass scales than the characteristic mass. The reason is thought to be feedback preventing star formation by heating the gas and/or pushing it out of the galaxy; at the high-mass end, you have powerful active galactic nuclei (AGN), while at the low-mass end, where the gravitational potential is smaller, stellar activity can do the job.

[![LF_HMF](https://i.sstatic.net/i4qd4.png)](https://i.sstatic.net/i4qd4.png)
Predicted LF assuming a [Jenkins et al (2001)](https://ui.adsabs.harvard.edu/abs/2001MNRAS.321..372J) with a constant $M/L$ (cyan), and observed $K$ band LF (yellow). Figure from [Benson et al. (2003)](https://ui.adsabs.harvard.edu/abs/2003ApJ...599...38B) with my own annotations.

The Milky Way, it turns out, is rather close to being such a characteristic mass and luminosity, which is why you sometimes hear that the Milky Way is a "typical" galaxy.

The LF is typically described by a [Schechter function](https://en.wikipedia.org/wiki/Luminosity_function_(astronomy)). In recent years though we are seeing evidence that the exponential cut-off is a bit too drastic, and that the high-mass end may be better described by another power law, albeit with a steep index, i.e. a so-called *broken power law* (see e.g. [Oesch et al. 2018](https://ui.adsabs.harvard.edu/abs/2018ApJ...855..105O)).

### …and in nature in general

Power laws are common not only in astronomy, but in all of nature, and even in "artificial" settings, such as linguistics (the occurrence of words in a language), economics (e.g. income distributions), and sociology (e.g. city sizes).

In some cases, the reason can be attributed to the fact that even non-power laws often behave like power laws in limited intervals. Other functions that occur in nature, such as $\exp(x)$ and $\sin(x)$, can be Taylor-expanded, and near "critical points" the first term dominates.

I'm actually not sure if the reason for the ubiquity of power laws in nature is entirely understood. I think you can explain it in *many* cases, but not in general. One example that *can* be explained, mathematically, is that if stochastic processes with *exponential* growth in expectation are "killed" (or observed) randomly, the distribution of the killed or observed state exhibits power-law behavior in one or both tails ([Reed & Hughes 2002](https://ui.adsabs.harvard.edu/abs/2002PhRvE..66f7103R)).

### But beware!

That said, I think there sometimes is a tendency to assume a power law relation between variables where none is necessarily expected. A power law is a linear relation if you take the logarithm of both $x$ and $y$, and probably the fact that deviations from some relation are suppressed in a log-log plot is sometimes a bit too appealing…

---

$\dagger$To be fair, scale invariance and self-similarity is not exactly the same, in that the former describes a continuous scaling, while the latter includes discrete scaling such as fractals which are identical to themselves only in certain steps and certain regions.


---

### Answer by username_2 on 2021-06-18. Score: 8

I have to admit that power-laws (in general) used to be my [shtick](https://en.m.wikipedia.org/wiki/Shtick) so I am happy to shed some light on their general importance in physics which obviously also hold for astronomy.

The main idea of a [power law](https://en.wikipedia.org/wiki/Power_law) is nicely written in Wikipedia, but the essential part is the I highlighted in the following quote:

> [A Powerlaw is] a functional relationship between two quantities, where a relative change in one quantity results in a proportional relative change in the other quantity, **independent of the initial size of those quantities**

The (mathematically) nice part is exactly this, that over a *really large range* of the x-values, the y-values follow the same dependency. Usually, "really large range" means values spreading over 3 to even 10 powers of ten.

When fitting power laws in practice, there are effects on the sides of the scale, meaning for small and large $x$-values, where usually attributed to "finite size effects".

### Further reading

* There is the research flavor [Self-organized criticality (SOC)](https://en.wikipedia.org/wiki/Self-organized_criticality) which *"is considered to be one of the mechanisms by which complexity arises in nature"*. One of the founding fathers (and mothers) was the late [Per Bak](https://en.wikipedia.org/wiki/Per_Bak). His book *[How Nature Works](https://books.google.de/books/about/How_Nature_Works.html?id=Bth4QgAACAAJ&redir_esc=y)* is accessible for laypersons and [IMHO](https://en.wiktionary.org/wiki/IMHO) a good read.


---

### Answer by username_1 on 2021-06-21. Score: 3

The [polytrope model](https://en.wikipedia.org/wiki/Polytrope) of stars made early numerical calculations of stellar structure possible using anything from mechanical calculators to early electronic computers, and in certain cases, even analytically!

> In astrophysics, a polytrope refers to a solution of the Lane–Emden equation in which the pressure depends upon the density in the form
>
> $$P=K\rho ^{{(n+1)/n}}$$
>
> where $P$ is pressure, $\rho$ is density and $K$ is a constant of proportionality. The constant $n$ is known as the polytropic index; note however that the polytropic index has an alternative definition as with $n$ as the exponent.

See answers to:

* Math SE: [Could there be exact solutions to the Lane-Emden equation for real n≥0 other than 0, 1, or 5?](https://math.stackexchange.com/q/4148101/284619)
* [Has stellar evolution ever been modeled analytically?](https://astronomy.stackexchange.com/q/29743/7982) especially [@ProfRob's answer](https://astronomy.stackexchange.com/a/29820/7982)!
* [Does Eddington's variable polytropic index better fit data from the Standard Solar Model?](https://astronomy.stackexchange.com/q/16633/7982)
* [Gravity inside a star?](https://astronomy.stackexchange.com/q/43913/7982)


---

