This page was generated by nbsphinx from docs/notebooks/formulary/solar_plasma_beta.ipynb.
Interactive online version: Binder badge.

Plasma beta in the solar atmosphere

This notebook demonstrates plasmapy.formulary by calculating plasma \(β\) in different regions of the solar atmosphere.

Contents

  1. Introduction

  2. Solar corona

  3. Solar chromosphere

  4. Quiet solar photosphere

  5. Sunspot photosphere

Introduction

Plasma beta (\(β\)) is one of the most fundamental plasma parameters. \(β\) is the ratio of the thermal plasma pressure to the magnetic pressure:

\[β = \frac{p_{therm}}{p_{mag}}.\]

How a plasma behaves depends strongly on \(β\). When \(β ≫ 1\), the magnetic field is not strong enough to significantly alter the dynamics, so the plasma motion is more gas-like. When \(β ≪ 1\), magnetic tension and magnetic pressure dominate the dynamics.

Let’s use plasmapy.formulary.beta to compare \(β\) in different parts of the solar atmosphere.

[1]:
import astropy.units as u

from plasmapy.formulary import beta

Solar corona

Let’s start by defining some plasma parameters for an active region in the solar corona.

[2]:
T_corona = 1 * u.MK
n_corona = 1e9 * u.cm**-3
B_corona = 50 * u.G

When we use these parameters in beta, we find that \(β\) is quite small which implies that the corona is magnetically dominated.

[3]:
beta(T_corona, n_corona, B_corona)
[3]:
$0.0013879798 \; \mathrm{}$

Solar chromosphere

Next let’s calculate \(β\) for the chromosphere. Bogod et al. (2015) found that the quiet chromosphere ranges from ∼40–200 G. We can get the temperature and number density of hydrogen from model C7 of Avrett & Loeser (2007) for 1 Mm above the photosphere.

[4]:
T_chromosphere = 6225 * u.K
n_chromosphere = 2.711e13 * u.cm**-3
B_chromosphere = [40, 200] * u.G
[5]:
beta(T_chromosphere, n_chromosphere, B_chromosphere)
[5]:
$[0.36599237,~0.014639695] \; \mathrm{}$

When \(B\) is small, plasma \(β\) is not too far from 1, which means that both magnetic and plasma pressure gradient forces are important when the chromospheric magnetic field is relatively weak. When near the higher range of \(B\), \(β\) is small so that the magnetic forces are more important than plasma pressure gradient forces.

Quiet solar photosphere

Let’s specify some characteristic plasma parameters for the solar photosphere.

[6]:
T_photosphere = 5800 * u.K
n_photosphere = 1e17 * u.cm**-3
B_photosphere = 400 * u.G

When we calculate β for the photosphere, we find that it is an order of magnitude larger than 1, so plasma pressure gradient forces are more important than magnetic forces.

[7]:
beta(T_photosphere, n_photosphere, B_photosphere)
[7]:
$12.578567 \; \mathrm{}$

Sunspot photosphere

The magnetic field in the solar photosphere is strongest in sunspots, so we would expect β to be lowest there. Let’s estimate some plasma parameters for a sunspot.

[8]:
T_sunspot = 4500 * u.K
B_sunspot = 2 * u.kG

When we calculate β, we find that both pressure gradient and magnetic forces will be important.

[9]:
beta(T_sunspot, n_photosphere, B_sunspot)
[9]:
$0.39036931 \; \mathrm{}$