> For the complete documentation index, see [llms.txt](https://slowdiveptg.gitbook.io/notes/llms.txt). Markdown versions of documentation pages are available by appending `.md` to page URLs; this page is available as [Markdown](https://slowdiveptg.gitbook.io/notes/stellar-evolution/chapter-8.-opacity.md). # Chapter 8. Opacity In this chapter we deal with the **opacity** $\kappa(\rho,T)$ \[cm$^2$/g], which corresponds to the efficiency of energy transport via radiation. In optically thick gas, energy is transported via photon diffusion. As is discussed in [Chapter 5](/notes/stellar-evolution/chapter-5.-energy-transport.md), we have $$ F\_\nu=-\frac{c}{3\rho\kappa\_\nu}\frac{\text de\_\text{int}}{\text dr} $$ where $e*\text{int}$ is given by the blackbody radiation luminosity $B*\nu$ $$ e\_\text{int}=\frac{4\pi}cB\_\nu(T)=\frac{8\pi h\nu^3}{c^3}\frac1{e^{h\nu/k\_BT}-1} $$ Thus $$ F\_\nu=-\frac{4\pi}{3\rho\kappa\_\nu}\frac{\partial B\_\nu}{\partial T}\frac{\partial T}{\partial r} $$ The total flux is $$ F=-\frac{4ac}{3\rho\kappa\_R}T^3\frac{\partial T}{\partial r} $$ where $\kappa\_R$ is the Rosseland mean opacity. Obviously, the efficiency of energy transported outward is inversely proportional to the opacity. We can further define optical depth $\tau$, $$ \tau\_\nu=\int\rho\kappa\_\nu\text dr\sim\rho\kappa\_\nu r $$ So when $\tau\gg1$ so that energy is transported via photon diffusion, the flux is approximately $$ F\sim\frac{acT^4}{\rho\kappa r}\sim\frac{\sigma\_{SB}T^4}{\tau} $$ ## Electron Scattering If an electromagnetic wave passes an electron, the electric field makes the electron oscillate. The oscillating electron represents a classical dipole that radiates in other directions, i.e. the electron scatters part of the energy of the incoming waves. From classical electrodynamics, elastic scattering of non-relativistic electrons, namely **Thomson scattering**, yields a cross section of $$ \sigma\_T=\frac{8\pi}{3}{r\_\text e^2} $$ where $r\_\text e$ is the classical electron radius. Thus the opacity is $$ \kappa^\text{es}*\nu=\frac{8\pi}{3}\frac{r*\text e^2}{\mu\_em\_\text p}\simeq0.35\text{ cm}^2\text{/g} $$ This opacity is $\nu$-independent, but when $T>10^8$ K, $\kappa*\nu$ slightly decreases ($\sim20\\%$ at $10^8$ K) for higher $T$. At that time the kinetic energy of electron is over $\sim10$ keV, $\gtrsim1\\%$ of $m*\text{e}c^2$. Therefore relativistic corrections (**Compton scattering**) become important. ## Absorption Due to Free-free Transitions If during its thermal motion a free electron passes an ion, the two charged particles form a system which can absorb and emit radiation. This mechanism is only effective as long as electron and ion are sufficiently close. Let us first consider the emissivity $\varepsilon\_\nu^\text{ff}$ \[erg/s/cm$^3$/Hz]. In a super, super rough estimation, it should be proportional to the densities of both electron and ions, as well as the duration of the fly-by, which inversely proportional to the magnitude of their relative velocity. $$ \varepsilon\_\nu^\text{ff}\propto \frac{n\_\text{e}n\_\text{ion}}{|v|} $$ Assuming the thermal equilibrium, electrons obey the Boltzmann distribution, $|v|\propto T^{1/2}$. As a result, $$ \varepsilon\_\nu^\text{ff}\propto \rho^2T^{-1/2} $$ Thus the opacity due to free-free transitions is $$ \kappa\_\nu^\text{ff}=\frac{\alpha\_\nu}{\rho}=\frac{\varepsilon\_\nu^\text{ff}}{\rho B\_\nu(T)} $$ The second equation is given by Kirchhoff's law of thermal radiation, which states that > For a body of any arbitrary material emitting and absorbing thermal electromagnetic radiation at every wavelength in thermodynamic equilibrium, the ratio of its emissive power to its dimensionless coefficient of absorption is equal to a universal function only of radiative wavelength and temperature. That universal function describes the perfect black-body emissive power. And thus $$ B\_\nu(T)=\frac{\varepsilon\_\nu^\text{ff}}{\alpha\_\nu} $$ where the **absorption coefficient** $\alpha*\nu\equiv\rho\kappa*\nu$ \[1/cm] denotes the attentuation of radiation per unit distance. So, $$ \kappa\_\nu^\text{ff}\propto\rho^{-1}{\varepsilon^{\text{ff}}\_\nu}\nu^{-3}\propto\rho T^{-1/2}\nu^{-3} $$ And the Rosseland mean opacity is $$ \kappa^\text{ff}=\frac{\int\frac{\partial B\_\nu(T)}{\partial T}\text d\nu}{\int\frac{1}{\kappa\_\nu}\frac{\partial B\_\nu(T)}{\partial T}\text d\nu}\propto\rho T^{-7/2} $$ This $\rho T^{-7/2}$ dependence is known as the **Kramer's law**. Precise derivation gives, $$ \kappa^\text{ff}=4\times10^{22}\left(\frac{\rho}{\text{g}\cdot\text{cm}^3}\right)\left(\frac{T}{\text K}\right)^{-7/2}\text{ cm}^2\text{/g} $$ ## Bound-free Transitions The best-known example is ionization of atoms. For a hydrogen atom at its ground state, the ionization process is $$ \ce{H + \gamma -> H+ + e-} $$ Here the energy of $\gamma$ must be higher than 13.6 eV, otherwise this reaction will not occur. In other words, for $h\nuh\nu*0$, $\alpha*\nu\propto\nu^{-3}$. Different excitation states have different $\nu*0$. The opacity $\kappa^\text{bf}*{\nu}$ for a mixture of hydrogen atoms in different states of excitation is a superposition of $\alpha*\nu$ for different stages of excitation. The resulting $\kappa^\text{bf}*{\nu}$􏰁􏰂 is a sawtooth function. ![](https://1509032923-files.gitbook.io/~/files/v0/b/gitbook-legacy-files/o/assets%2F-MPMxe8Bu9WDT3p-DA_8%2Fsync%2F0787287188f64a2a222d575540bcbf810e546305.png?generation=1608873191103259\&alt=media) The relative scale of $\kappa^\text{bf}*{\nu}$ comparing to $\kappa^\text{ff}*{\nu}$ is approximately $$ \kappa^\text{bf}*{\nu}\simeq10^3Z\kappa^\text{bf}*{\nu} $$ where $Z$ is the **metalicity**. $Z\_\odot\sim0.02$. ## Bound-bound Transitions $\kappa^\text{bb}\_\nu$ is the superposition for opacity due to numerous possiblities of line adsorption. The dependency on $\nu$ is thus really messy. ![](https://1509032923-files.gitbook.io/~/files/v0/b/gitbook-legacy-files/o/assets%2F-MPMxe8Bu9WDT3p-DA_8%2Fsync%2F580e07aeedc593393ee94d1515ffa99c862be63a.png?generation=1608873191336682\&alt=media) ## The Negative Hydrogen Ion ($\ce{H-}$) At a sufficiently low temperature, protons and electrons start to recombine. When it gets even cooler, $\ce{H-}$ begins to effectively form via $$ \ce{H + e- -> H- + \gamma} $$ $\ce {H-}$ can also contribute to the opacity through photo-dissociation $$ \ce{H- + \gamma -> H + e-} $$ With Saha equation, we can estimate the $\ce{H-}$ fraction as a function of temperature and density. Then we can obtain the Rosseland mean opacity $$ \kappa^{\ce{H-}}\simeq2.5\times10^{-31}\left(\frac{Z}{Z\_\odot}\right)\left(\frac{\rho}{\text{g}\cdot\text{cm}^3}\right)\left(\frac{T}{\text K}\right)^{9}\text{ cm}^2\text{/g} $$ when $3000\