<?xml version="1.0" encoding="utf-8" standalone="yes"?><rss version="2.0" xmlns:atom="http://www.w3.org/2005/Atom"><channel><title>Physics on 孤筝の温暖小家</title><link>https://www.guzhengsvt.cn/en/categories/physics/</link><description>Recent content from 孤筝の温暖小家</description><generator>Hugo</generator><language>en</language><managingEditor>lvbowen040427@163.com (孤筝)</managingEditor><webMaster>lvbowen040427@163.com (孤筝)</webMaster><copyright>All articles on this blog are licensed under the BY-NC-SA license agreement unless otherwise stated. Please indicate the source when reprinting!</copyright><lastBuildDate>Fri, 05 Sep 2025 11:05:15 +0800</lastBuildDate><atom:link href="https://www.guzhengsvt.cn/en/categories/physics/index.xml" rel="self" type="application/rss+xml"/><item><title>Quantum Physics</title><link>https://www.guzhengsvt.cn/en/post/physics/quantum-physics/</link><pubDate>Fri, 05 Sep 2025 11:05:15 +0800</pubDate><author>lvbowen040427@163.com (孤筝)</author><guid>https://www.guzhengsvt.cn/en/post/physics/quantum-physics/</guid><description>
<![CDATA[<h1>Quantum Physics</h1><p>Author: 孤筝(lvbowen040427@163.com)</p>
        
          <h2 id="chapter-1-wave-particle-duality-and-state-description">
<a class="header-anchor" href="#chapter-1-wave-particle-duality-and-state-description"></a>
Chapter 1 Wave-Particle Duality and State Description
</h2><h3 id="11-formation-and-applications-of-quantum-mechanics">
<a class="header-anchor" href="#11-formation-and-applications-of-quantum-mechanics"></a>
1.1 Formation and Applications of Quantum Mechanics
</h3><h4 id="111-old-quantum-theory">
<a class="header-anchor" href="#111-old-quantum-theory"></a>
1.1.1 Old Quantum Theory
</h4><h5 id="photoelectric-effect-and-the-photon-hypothesis">
<a class="header-anchor" href="#photoelectric-effect-and-the-photon-hypothesis"></a>
Photoelectric Effect and the Photon Hypothesis
</h5><ul>
<li><strong>Photon energy</strong>: $E = h\nu$</li>
<li><strong>Threshold frequency</strong>: $\nu_0 = \dfrac{W_0}{h}$; no photoelectrons are emitted when $\nu < \nu_0$</li>
<li><strong>Photoelectric equation</strong>:<br>

$$
  E_k^{\text{max}} = \frac{1}{2}\mu v^2_m = h\nu - W_0
  $$</li>
<li>The photoelectric effect demonstrates the particle nature of light.</li>
</ul>
<h5 id="energy-momentum-relation-of-photons-and-wave-particle-unity">
<a class="header-anchor" href="#energy-momentum-relation-of-photons-and-wave-particle-unity"></a>
Energy-Momentum Relation of Photons and Wave-Particle Unity
</h5><ul>
<li>
<p><strong>Relativistic energy-momentum relation</strong><br>
</p>
$$
  E^2=(pc)^2+(m_0c^2)^2,\quad m_0=0\ \Rightarrow\ E=c\,\lVert\vec p\rVert
  $$</li>
<li>
<p><strong>Photon energy</strong><br>
</p>
$$
  E=h\nu=\frac{hc}{\lambda}=\hbar\omega
  $$</li>
<li>
<p><strong>Photon momentum (vector form)</strong><br>
Let $\mathbf n$ be the unit vector in the propagation direction, then
</p>
$$
  \vec p=\frac{E}{c}\,\mathbf n=\frac{h}{\lambda}\,\mathbf n=\hbar\vec k,\quad
  \vec k=\frac{2\pi}{\lambda}\,\mathbf n
  $$</li>
<li>
<p><strong>Unified correspondence of wave and particle pictures</strong><br>
</p>
$$
  E\ \longleftrightarrow\ \hbar\omega,\qquad
  \vec p\ \longleftrightarrow\ \hbar\vec k
  $$</li>
</ul>
<h5 id="bohr-structure-of-the-hydrogen-atom">
<a class="header-anchor" href="#bohr-structure-of-the-hydrogen-atom"></a>
Bohr Structure of the Hydrogen Atom
</h5><ul>
<li>Quantization of orbital angular momentum:

$$
  L = n\hbar,\quad n=1,2,3,\dots
  $$</li>
<li>Energy levels:

$$
  E_n = -\frac{13.6\ \text{eV}}{n^2}
  $$</li>
<li>This successfully explains the line spectrum of hydrogen.</li>
</ul>
<h5 id="bohrs-postulates">
<a class="header-anchor" href="#bohrs-postulates"></a>
Bohr&rsquo;s Postulates
</h5><ul>
<li>Electrons moving in stable orbits do not radiate energy.</li>
<li>Electrons absorb or emit energy when they jump between energy levels:

$$
  \Delta E = h\nu
  $$</li>
</ul>
<h5 id="compton-effect">
<a class="header-anchor" href="#compton-effect"></a>
Compton Effect
</h5><ul>
<li>The wavelength of a high-energy photon increases after scattering from an electron:

$$
  \Delta\lambda = \lambda' - \lambda = \frac{h}{m_ec}(1-\cos\theta)
  $$</li>
<li>The experiment confirms both the particle nature of light and conservation of momentum.</li>
</ul>
<h5 id="blackbody-radiation">
<a class="header-anchor" href="#blackbody-radiation"></a>
Blackbody Radiation
</h5><ul>
<li>Energy quantization assumption: the energy of the electromagnetic field takes discrete values $E=nh\nu$.</li>
<li>Planck formula:

$$
  u(\nu,T)=\frac{8\pi h\nu^3}{c^3}\frac{1}{e^{h\nu/kT}-1}
  $$</li>
<li>This successfully explains the blackbody spectrum and marks the birth of quantum theory.</li>
</ul>
<h4 id="112-wave-particle-duality-of-microscopic-particles">
<a class="header-anchor" href="#112-wave-particle-duality-of-microscopic-particles"></a>
1.1.2 Wave-Particle Duality of Microscopic Particles
</h4><h5 id="de-broglie-hypothesis">
<a class="header-anchor" href="#de-broglie-hypothesis"></a>
de Broglie Hypothesis
</h5><ul>
<li>Microscopic particles have not only particle properties but also wave properties.</li>
<li>Every particle with momentum $\vec p$ corresponds to a matter wave whose wavelength and frequency are related to its momentum and energy.</li>
</ul>
<h5 id="de-broglie-relations">
<a class="header-anchor" href="#de-broglie-relations"></a>
de Broglie Relations
</h5><ul>
<li>Wavelength:

$$
  \lambda = \frac{h}{p}
  $$</li>
<li>Vector form:

$$
  \vec p = \hbar \vec k
  $$</li>
<li>Frequency:

$$
  E = h\nu = \hbar\omega
  $$</li>
</ul>
<hr>
<h3 id="12-states-and-wave-functions">
<a class="header-anchor" href="#12-states-and-wave-functions"></a>
1.2 States and Wave Functions
</h3><h4 id="121-uncertainty-principle">
<a class="header-anchor" href="#121-uncertainty-principle"></a>
1.2.1 Uncertainty Principle
</h4><ul>
<li>The position and momentum of a microscopic particle cannot be measured simultaneously with arbitrary precision.</li>
<li>Heisenberg uncertainty relation:

$$
  \Delta x \cdot \Delta p_x \geq \frac{\hbar}{2}
  $$</li>
<li>Energy-time uncertainty relation:

$$
  \Delta E \cdot \Delta t \geq \frac{\hbar}{2}
  $$</li>
<li>Its origin lies in wave-particle duality and the non-commutativity of operators.</li>
</ul>
<h4 id="122-wave-function">
<a class="header-anchor" href="#122-wave-function"></a>
1.2.2 Wave Function
</h4><ul>
<li>To describe the state of a microscopic particle, we introduce the <strong>wave function</strong> $\psi(\vec r,t)$.</li>
<li><strong>Probabilistic interpretation</strong>: $|\psi(\vec r,t)|^2 dV$ gives the probability of finding the particle in the volume element $dV$.</li>
<li>The wave function must satisfy the superposition principle and the Schrödinger equation.</li>
<li>The total probability over all space is 1, so the probability distribution depends only on the relative magnitude of the wave function, not its absolute scale.</li>
<li>Multiplying the wave function by a constant does not change the physical state it describes.</li>
<li>Standard requirements for a wave function: single-valued, finite, and continuous.</li>
</ul>
<h4 id="123-normalization-of-the-wave-function">
<a class="header-anchor" href="#123-normalization-of-the-wave-function"></a>
1.2.3 Normalization of the Wave Function
</h4><ul>
<li><strong>Normalization condition</strong>:

$$
  \int_{-\infty}^{\infty} \psi^* (\vec r,t) \psi (\vec r,t) dV = 1
  $$</li>
<li><strong>How to normalize a wave function</strong>

$$
  \int_{-\infty}^{\infty} |\Psi(\vec r,t)|^2 dV = A^2 \int_{-\infty}^{\infty} |\psi(\vec r,t)|^2 dV = 1
  $$
where $A$ is the normalization constant.</li>
</ul>
<hr>
<h3 id="13-schrödinger-equation">
<a class="header-anchor" href="#13-schr%c3%b6dinger-equation"></a>
1.3 Schrödinger Equation
</h3><h4 id="131-wave-equation-for-a-free-particle">
<a class="header-anchor" href="#131-wave-equation-for-a-free-particle"></a>
1.3.1 Wave Equation for a Free Particle
</h4><p><strong>Concept</strong><br>
A free particle is a particle not subject to external forces. In quantum mechanics, its state is described by the wave function $\psi(\vec{r},t)$ and satisfies the Schrödinger equation.</p>
<p><strong>Schrödinger equation for a free particle</strong><br>
</p>
$$
i\hbar \frac{\partial \psi(\vec{r},t)}{\partial t} = -\frac{\hbar^2}{2m} \nabla^2 \psi(\vec{r},t)
$$<p>
where:</p>
<ul>
<li>$\hbar$: reduced Planck constant</li>
<li>$m$: mass of the particle</li>
<li>$\nabla^2$: Laplacian operator</li>
</ul>
<p><strong>Plane-wave solution</strong><br>
</p>
$$
\psi(\vec{r},t) = A e^{i(\vec{k}\cdot\vec{r} - \omega t)}
$$<p>
where:</p>
<ul>
<li>$\vec{k}$ is the wave vector, with $|\vec{k}| = k$</li>
<li>$\omega$ is the angular frequency, satisfying

$$
E = \hbar \omega = \frac{\hbar^2 k^2}{2m}
$$</li>
</ul>
<p><strong>Momentum-wave vector relation</strong>
</p>
$$
\vec{p} = \hbar \vec{k}
$$<h5 id="plane-wave-derivation-of-the-free-particle-schrödinger-equation">
<a class="header-anchor" href="#plane-wave-derivation-of-the-free-particle-schr%c3%b6dinger-equation"></a>
Plane-Wave Derivation of the Free-Particle Schrödinger Equation
</h5><p><strong>1. Assume a plane-wave form</strong>
</p>
$$
\psi(\vec{r},t) = A e^{i(\vec{k}\cdot\vec{r} - \omega t)}
$$<p>
with $\vec{k}$ the wave vector, $\omega$ the angular frequency, and $A$ the amplitude constant.</p>
<p><strong>2. Time derivative</strong>
</p>
$$
\frac{\partial \psi}{\partial t} = \frac{\partial}{\partial t} \left[ A e^{i(\vec{k}\cdot\vec{r} - \omega t)} \right] 
= -i \omega A e^{i(\vec{k}\cdot\vec{r} - \omega t)} 
= -i \omega \psi
$$<p>
Multiplying by $i\hbar$ gives
</p>
$$
i\hbar \frac{\partial \psi}{\partial t} = \hbar \omega \psi
$$<p><strong>3. Spatial Laplacian</strong>
</p>
$$
\nabla^2 \psi = \nabla^2 \left[ A e^{i \vec{k}\cdot\vec{r}} e^{-i\omega t} \right] 
= -k^2 A e^{i(\vec{k}\cdot\vec{r} - \omega t)} 
= -k^2 \psi
$$<p>
Thus,
</p>
$$
-\frac{\hbar^2}{2m} \nabla^2 \psi = \frac{\hbar^2 k^2}{2m} \psi
$$<p><strong>4. Energy relation</strong>
</p>
$$
E = \frac{\hbar^2 k^2}{2m} = \hbar \omega
$$<p><strong>5. Final equation</strong>
</p>
$$
i\hbar \frac{\partial \psi}{\partial t} = -\frac{\hbar^2}{2m} \nabla^2 \psi
$$<p><strong>Remark</strong></p>
<ul>
<li>This derivation uses only the plane-wave form and differentiation, without relying on operator definitions.</li>
<li>It corresponds to the case $V=0$.</li>
</ul>
<h4 id="133-stationary-state-schrödinger-equation-and-stationary-wave-functions">
<a class="header-anchor" href="#133-stationary-state-schr%c3%b6dinger-equation-and-stationary-wave-functions"></a>
1.3.3 Stationary-State Schrödinger Equation and Stationary Wave Functions
</h4><p><strong>Concept</strong><br>
A stationary-state wave function has separable time dependence:
</p>
$$
\psi(\vec{r},t) = \phi(\vec{r}) e^{-i E t / \hbar}
$$<p>
Let $f(t)=e^{-i E t / \hbar}$.<br>
Here $\phi(\vec{r})$ depends only on spatial coordinates, and $E$ is the total energy of the particle.</p>
<p><strong>Derivation</strong>
Starting from the time-dependent Schrödinger equation:
</p>
$$
i\hbar \frac{\partial \psi(\vec{r},t)}{\partial t} = \left[ -\frac{\hbar^2}{2m} \nabla^2 + V(\vec{r}) \right] \psi(\vec{r},t)
$$<p>
substitute $\psi(\vec{r},t) = \phi(\vec{r}) e^{-i E t / \hbar}$ to obtain
</p>
$$
\left[ -\frac{\hbar^2}{2m} \nabla^2 + V(\vec{r}) \right] \phi(\vec{r}) = E \phi(\vec{r})
$$<p><strong>Time-independent form</strong>
</p>
$$
i \hbar \frac{df}{dt}=E f , \; f= e^{-i E t / \hbar}
$$<p>
</p>
$$
-\frac{\hbar^2}{2m} \nabla^2 \phi(\vec{r}) + V(\vec{r}) \phi(\vec{r}) = E \phi(\vec{r})
$$<p><strong>Remark</strong></p>
<ul>
<li>$\phi(\vec{r})$ is called a stationary-state wave function or eigenfunction.</li>
<li>$E$ is the corresponding energy eigenvalue.</li>
</ul>
<h4 id="derivation-of-the-schrödinger-equation-from-operators">
<a class="header-anchor" href="#derivation-of-the-schr%c3%b6dinger-equation-from-operators"></a>
Derivation of the Schrödinger Equation from Operators
</h4><p><strong>1. Start from classical energy</strong>
</p>
$$
E = \frac{p^2}{2m} + V(\vec{r},t)
$$<p><strong>2. Introduce the de Broglie relations</strong>
</p>
$$
\vec{p} = \hbar \vec{k}, \quad E = \hbar \omega
$$<p>
and the plane-wave form
</p>
$$
\psi(\vec{r},t) \sim e^{i(\vec{k}\cdot\vec{r} - \omega t)}
$$<p><strong>3. Operator representation</strong>
</p>
$$
\hat{E} = i\hbar \frac{\partial}{\partial t}, \quad \hat{\vec{p}} = -i\hbar \nabla
$$<p><strong>4. Kinetic-energy and Hamiltonian operators</strong>
</p>
$$
\hat{T} = \frac{\hat{p}^2}{2m} = -\frac{\hbar^2}{2m}\nabla^2
$$<p>
</p>
$$
\hat{H} = \hat{T} + V(\vec{r},t) = -\frac{\hbar^2}{2m}\nabla^2 + V(\vec{r},t)
$$<p><strong>5. Schrödinger equation</strong>
</p>
$$
i\hbar \frac{\partial \psi(\vec{r},t)}{\partial t} = \hat{H} \psi(\vec{r},t)
$$<p>
that is,
</p>
$$
i\hbar \frac{\partial \psi(\vec{r},t)}{\partial t} = \left[ -\frac{\hbar^2}{2m}\nabla^2 + V(\vec{r},t) \right] \psi(\vec{r},t)
$$<p>For a free particle ($V=0$), this reduces to
</p>
$$
i\hbar \frac{\partial \psi}{\partial t} = -\frac{\hbar^2}{2m}\nabla^2 \psi
$$<h4 id="principle-of-superposition-of-states">
<a class="header-anchor" href="#principle-of-superposition-of-states"></a>
Principle of Superposition of States
</h4><p><strong>Concept</strong><br>
If $\psi_1$ and $\psi_2$ are two possible states of the same system, then their linear combination
</p>
$$
\psi = c_1 \psi_1 + c_2 \psi_2
$$<p>
is also a possible state, where $c_1$ and $c_2$ are complex coefficients.</p>
<p><strong>General expansion</strong>
</p>
$$
\psi(\vec{r},t) = \sum_{n} c_n \phi_n(\vec{r},t)
$$<p>
with</p>
<ul>
<li>$c_n$ the expansion coefficients, or probability amplitudes;</li>
<li>probabilities $|c_n|^2$ satisfying

$$
\sum_n |c_n|^2 = 1
$$</li>
</ul>
<p><strong>Remark</strong></p>
<ul>
<li>Superposition is one of the most fundamental principles of quantum mechanics.</li>
<li>Different eigenstates may superpose, but a measurement yields only one eigenvalue.</li>
<li>Interference in superposed states is one of the essential features that distinguishes quantum mechanics from classical mechanics.</li>
</ul>
<h2 id="chapter-2-simple-applications-of-the-schrödinger-equation">
<a class="header-anchor" href="#chapter-2-simple-applications-of-the-schr%c3%b6dinger-equation"></a>
Chapter 2 Simple Applications of the Schrödinger Equation
</h2><h3 id="21-one-dimensional-infinite-potential-well">
<a class="header-anchor" href="#21-one-dimensional-infinite-potential-well"></a>
2.1 One-Dimensional Infinite Potential Well
</h3><h4 id="211-solving-the-equation">
<a class="header-anchor" href="#211-solving-the-equation"></a>
2.1.1 Solving the Equation
</h4><p><strong>1. Potential</strong>
</p>
$$
V(x) =
\begin{cases}
0, & 0 < x < L \\
\infty, & x \leq 0 \ \text{or} \ x \geq L
\end{cases}
$$<p><strong>2. Schrödinger equation in the well</strong>
</p>
$$
-\frac{\hbar^2}{2m} \frac{d^2 \phi(x)}{dx^2} = E \phi(x)
$$<p>
which becomes
</p>
$$
\frac{d^2 \phi(x)}{dx^2} + k^2 \phi(x) = 0,\qquad
k^2 = \frac{2mE}{\hbar^2}
$$<p><strong>3. General solution</strong>
</p>
$$
\phi(x) = A \sin(kx) + B \cos(kx)
$$<p><strong>4. Boundary conditions</strong>
</p>
$$
\phi(0) = 0, \quad \phi(L) = 0
$$<p>
Hence $B=0$ and $kL = n\pi$ for $n=1,2,3,\dots$.</p>
<p><strong>5. Eigenfunctions and eigenvalues</strong>
</p>
$$
\phi_n(x) = \sqrt{\frac{2}{L}} \sin\left(\frac{n\pi x}{L}\right), \quad n=1,2,3,\dots
$$<p>
</p>
$$
E_n = \frac{n^2 \pi^2 \hbar^2}{2mL^2}, \quad n=1,2,3,\dots
$$<p><strong>Remark</strong></p>
<ul>
<li>The energy is quantized.</li>
<li>The ground-state energy is nonzero, showing the zero-point energy.</li>
</ul>
<hr>
<h3 id="22-special-functions-in-mathematical-physics">
<a class="header-anchor" href="#22-special-functions-in-mathematical-physics"></a>
2.2 Special Functions in Mathematical Physics
</h3><h4 id="221-orthogonality-and-normalization">
<a class="header-anchor" href="#221-orthogonality-and-normalization"></a>
2.2.1 Orthogonality and Normalization
</h4><p><strong>Orthogonality</strong>
</p>
$$
\int_a^b \phi_m(x)\,\phi_n(x)\,dx = 0 \quad (m \neq n)
$$<p><strong>Normalization</strong>
</p>
$$
\int_a^b |\phi_n(x)|^2 dx = 1
$$<p><strong>Orthonormality</strong>
</p>
$$
\int_a^b \phi_m(x)\,\phi_n(x)\,dx = \delta_{mn}
$$<h4 id="222-expansion-in-an-orthonormal-set">
<a class="header-anchor" href="#222-expansion-in-an-orthonormal-set"></a>
2.2.2 Expansion in an Orthonormal Set
</h4>$$
f(x) = \sum_{n=1}^{\infty} c_n \phi_n(x), \qquad
c_n = \int_a^b f(x)\,\phi_n(x)\,dx
$$<h4 id="223-fourier-series">
<a class="header-anchor" href="#223-fourier-series"></a>
2.2.3 Fourier Series
</h4>$$
f(x) = \frac{a_0}{2} + \sum_{n=1}^\infty \left[ a_n \cos\left(\frac{n\pi x}{L}\right) + b_n \sin\left(\frac{n\pi x}{L}\right) \right]
$$<p>
where
</p>
$$
a_n = \frac{1}{L} \int_{-L}^L f(x)\cos\left(\frac{n\pi x}{L}\right)\,dx, \quad
b_n = \frac{1}{L} \int_{-L}^L f(x)\sin\left(\frac{n\pi x}{L}\right)\,dx
$$<h4 id="224-constructing-orthonormal-functions">
<a class="header-anchor" href="#224-constructing-orthonormal-functions"></a>
2.2.4 Constructing Orthonormal Functions
</h4><p>The standard method is <strong>Gram-Schmidt orthogonalization</strong>:
</p>
$$
\phi_1(x) = \frac{f_1(x)}{\sqrt{\int |f_1(x)|^2 dx}}
$$<p>
</p>
$$
\phi_2(x) = \frac{f_2(x) - \int \phi_1(x) f_2(x)\,dx \,\phi_1(x)}{\sqrt{\int \left|f_2(x) - \int \phi_1(x) f_2(x)\,dx \,\phi_1(x)\right|^2 dx}}
$$<h4 id="225-legendre-polynomials-and-other-special-functions">
<a class="header-anchor" href="#225-legendre-polynomials-and-other-special-functions"></a>
2.2.5 Legendre Polynomials and Other Special Functions
</h4><p><strong>Legendre polynomials</strong>
</p>
$$
(1-x^2)\frac{d^2 y}{dx^2} - 2x \frac{dy}{dx} + l(l+1)y = 0
$$<p>
with orthogonality
</p>
$$
\int_{-1}^{1} P_l(x) P_{l'}(x)\,dx = \frac{2}{2l+1}\delta_{ll'}
$$<p><strong>Other common special functions</strong></p>
<ul>
<li><strong>Spherical harmonics</strong> $Y_l^m(\theta,\phi)$ appear in angular momentum problems.</li>
<li><strong>Bessel functions</strong> $J_n(x)$ appear in cylindrical symmetry problems.</li>
<li><strong>Hermite polynomials</strong> $H_n(x)$ appear in harmonic oscillator problems.</li>
</ul>
<p>These special functions are solutions of the Schrödinger equation under different boundary conditions and symmetries.</p>
<hr>
<h3 id="23-linear-harmonic-oscillator">
<a class="header-anchor" href="#23-linear-harmonic-oscillator"></a>
2.3 Linear Harmonic Oscillator
</h3><hr>
<h3 id="24-hydrogen-atom">
<a class="header-anchor" href="#24-hydrogen-atom"></a>
2.4 Hydrogen Atom
</h3><h4 id="241-solution-of-the-equation-separated-into--parts">
<a class="header-anchor" href="#241-solution-of-the-equation-separated-into--parts"></a>
2.4.1 Solution of the Equation (Separated into $r,\ \theta,\ \phi$ Parts)
</h4><p><strong>1. Time-independent Schrödinger equation under a Coulomb potential</strong>
</p>
$$
V(r) = -\dfrac{e^2}{4\pi\varepsilon_0 r}
$$<p>
and
</p>
$$
-\dfrac{\hbar^2}{2m}\nabla^2\Psi(r,\theta,\phi) + V(r)\Psi = E\Psi.
$$<p><strong>2. Separation of variables</strong>
</p>
$$
\Psi(r,\theta,\phi)=R(r)\,Y(\theta,\phi).
$$<p>
This leads to three equations in $r$, $\theta$, and $\phi$ after separation.</p>
<p><strong>3. Equation in $\phi$</strong>
</p>
$$
\frac{1}{\Phi(\phi)}\frac{d^2\Phi}{d\phi^2} = -m^2 \quad\Rightarrow\quad
\Phi_m(\phi)=\frac{1}{\sqrt{2\pi}} e^{i m\phi},\quad m\in\mathbb{Z}.
$$<p><strong>4. Equation in $\theta$</strong>
</p>
$$
\frac{1}{\sin\theta}\frac{d}{d\theta}\!\left(\sin\theta\frac{d\Theta}{d\theta}\right)
+\left[l(l+1)-\frac{m^2}{\sin^2\theta}\right]\Theta=0
$$<p>
with solutions proportional to the associated Legendre functions:
</p>
$$
\Theta_{l}^{m}(\theta)\propto P_l^{m}(\cos\theta).
$$<p><strong>5. Angular part: spherical harmonics</strong>
</p>
$$
Y_l^m(\theta,\phi)=N_{l}^{m}\,P_l^{m}(\cos\theta)\,e^{im\phi},
$$<p>
which satisfy
</p>
$$
\hat L^2 Y_l^m = l(l+1)\hbar^2 Y_l^m,\qquad \hat L_z Y_l^m = m\hbar Y_l^m.
$$<p><strong>6. Radial equation</strong>
Let $u(r)=rR(r)$, then
</p>
$$
-\frac{\hbar^2}{2m}\frac{d^2 u}{dr^2} + \left[ -\frac{e^2}{4\pi\varepsilon_0 r} + \frac{\hbar^2 l(l+1)}{2m r^2} \right] u = E u.
$$<p><strong>7. Energy eigenvalues</strong>
</p>
$$
E_n = -\frac{m e^4}{2(4\pi\varepsilon_0)^2 \hbar^2}\,\frac{1}{n^2} = -\frac{13.6057\ \mathrm{eV}}{n^2},\qquad n=1,2,3,\dots
$$<p>
with $l=0,1,\dots,n-1$.</p>
<p><strong>8. Wave function</strong>
</p>
$$
\Psi_{n l m}(r,\theta,\phi)=R_{n l}(r)\,Y_l^m(\theta,\phi),
$$<p>
and
</p>
$$
R_{n l}(r)=N_{n l}\left(\frac{2r}{n a_0}\right)^{l} e^{-r/(n a_0)} L_{n-l-1}^{2l+1}\!\left(\frac{2r}{n a_0}\right),
$$<p>
where $a_0=\dfrac{4\pi\varepsilon_0\hbar^2}{m e^2}$ is the Bohr radius.</p>
<h4 id="242-results-and-discussion">
<a class="header-anchor" href="#242-results-and-discussion"></a>
2.4.2 Results and Discussion
</h4><p><strong>1. Quantum numbers and their meanings</strong></p>
<ul>
<li>$n$: principal quantum number</li>
<li>$l$: orbital angular momentum quantum number</li>
<li>$m$: magnetic quantum number</li>
</ul>
<p><strong>2. Degeneracy</strong><br>
For the Coulomb potential, the energy depends only on $n$. The degeneracy of the level with principal quantum number $n$ is $n^2$.</p>
<p><strong>3. Spatial structure of the wave function</strong></p>
<ul>
<li>The angular part is given by the spherical harmonics.</li>
<li>The radial part $R_{nl}(r)$ has $n-l-1$ radial nodes.</li>
<li>The ground state $(1,0,0)$ is spherically symmetric and has no radial node.</li>
</ul>
<p><strong>4. Summary</strong><br>
The hydrogen atom is solved by separating variables in spherical coordinates. The angular equations give spherical harmonics and angular quantum numbers, while the radial equation yields the discrete energy levels and radial eigenfunctions.</p>
<h2 id="chapter-3-operator-representation-of-dynamical-variables-and-representation-theory">
<a class="header-anchor" href="#chapter-3-operator-representation-of-dynamical-variables-and-representation-theory"></a>
Chapter 3 Operator Representation of Dynamical Variables and Representation Theory
</h2><h3 id="31-relation-between-dynamical-variables-and-operators">
<a class="header-anchor" href="#31-relation-between-dynamical-variables-and-operators"></a>
3.1 Relation Between Dynamical Variables and Operators
</h3><h4 id="311-mathematical-knowledge-of-operators">
<a class="header-anchor" href="#311-mathematical-knowledge-of-operators"></a>
3.1.1 Mathematical Knowledge of Operators
</h4><ol>
<li>
<p><strong>Definition of an operator</strong><br>
An operator is a rule acting on a function space or state space. In quantum mechanics, physical quantities are represented by operators, and the wave function is the object on which they act.</p>
</li>
<li>
<p><strong>Linearity</strong><br>
If
</p>
$$
   A(c_1\psi_1 + c_2\psi_2) = c_1 A\psi_1 + c_2 A\psi_2,
   $$<p>
then $A$ is a linear operator.</p>
</li>
<li>
<p><strong>Commutation relations</strong><br>
The commutator is defined by
</p>
$$
   [A,B] = AB - BA.
   $$<p>
If $[A,B]=0$, the two operators are said to commute.</p>
</li>
<li>
<p><strong>Hermitian operators</strong><br>
If
</p>
$$
   \langle \psi | A\varphi \rangle = \langle A\psi | \varphi \rangle,
   $$<p>
then $A$ is Hermitian. Hermitian operators have real eigenvalues and represent observables.</p>
</li>
</ol>
<h4 id="312-dynamical-variables-and-operators">
<a class="header-anchor" href="#312-dynamical-variables-and-operators"></a>
3.1.2 Dynamical Variables and Operators
</h4><ol>
<li>
<p><strong>Basic idea</strong><br>
Every classical quantity $f(q,p)$ corresponds to a quantum operator $\hat f$.</p>
</li>
<li>
<p><strong>Typical operator forms in the coordinate representation</strong>
</p>
$$
   \hat{x} = x, \qquad \hat{p} = -i\hbar \frac{\partial}{\partial x}
   $$</li>
<li>
<p><strong>Fundamental commutation relation</strong>
</p>
$$
   [\hat{x}, \hat{p}] = i\hbar
   $$</li>
<li>
<p><strong>Measurement and eigenvalue equations</strong>
</p>
$$
   \hat{A}\psi_a = a\psi_a
   $$<p>
Here $a$ is a possible measurement outcome, and $\psi_a$ is the corresponding eigenfunction.</p>
</li>
</ol>
<hr>
<h3 id="32-commutation-relations-and-the-uncertainty-principle">
<a class="header-anchor" href="#32-commutation-relations-and-the-uncertainty-principle"></a>
3.2 Commutation Relations and the Uncertainty Principle
</h3><h4 id="321-commutation-relations">
<a class="header-anchor" href="#321-commutation-relations"></a>
3.2.1 Commutation Relations
</h4><ol>
<li>
<p><strong>Definition</strong>
</p>
$$
   [A,B] = AB - BA
   $$<p>
If $[A,B]=0$, the two physical quantities can have simultaneous definite values.</p>
</li>
<li>
<p><strong>Basic relation</strong>
</p>
$$
   [\hat{x}, \hat{p}_x] = i\hbar
   $$</li>
<li>
<p><strong>Three-dimensional form</strong>
</p>
$$
   [\hat{x}_i, \hat{p}_j] = i\hbar \delta_{ij}, \quad [\hat{x}_i, \hat{x}_j]=0, \quad [\hat{p}_i, \hat{p}_j]=0
   $$</li>
<li>
<p><strong>Physical meaning</strong><br>
Commutation relations determine whether two observables can be measured simultaneously with arbitrary precision.</p>
</li>
</ol>
<h4 id="322-uncertainty-principle">
<a class="header-anchor" href="#322-uncertainty-principle"></a>
3.2.2 Uncertainty Principle
</h4><ol>
<li>
<p><strong>Mathematical form</strong>
</p>
$$
   (\Delta A)^2 = \langle (A-\langle A \rangle)^2 \rangle,\qquad
   (\Delta B)^2 = \langle (B-\langle B \rangle)^2 \rangle
   $$<p>
which leads to
</p>
$$
   \Delta A \cdot \Delta B \geq \frac{1}{2}\left| \langle [A,B] \rangle \right|
   $$</li>
<li>
<p><strong>Position-momentum uncertainty</strong>
</p>
$$
   \Delta x \cdot \Delta p \geq \frac{\hbar}{2}
   $$</li>
<li>
<p><strong>Energy-time uncertainty</strong>
</p>
$$
   \Delta E \cdot \Delta t \gtrsim \hbar
   $$</li>
</ol>
<hr>
<h3 id="33-representation-theory">
<a class="header-anchor" href="#33-representation-theory"></a>
3.3 Representation Theory
</h3><h4 id="331-mathematical-basis">
<a class="header-anchor" href="#331-mathematical-basis"></a>
3.3.1 Mathematical Basis
</h4><ol>
<li>
<p><strong>Concept of representation</strong><br>
States and operators can be represented in different bases, such as the coordinate, momentum, and energy representations.</p>
</li>
<li>
<p><strong>Expansion of a state</strong>
</p>
$$
   |\psi\rangle = \sum_n c_n |\phi_n\rangle,\qquad c_n = \langle \phi_n | \psi \rangle
   $$</li>
<li>
<p><strong>Matrix elements</strong>
</p>
$$
   A_{mn} = \langle \phi_m | \hat{A} | \phi_n \rangle
   $$</li>
<li>
<p><strong>Completeness and orthogonality</strong>
</p>
$$
   \sum_n |\phi_n\rangle \langle \phi_n| = I,\qquad
   \langle \phi_m | \phi_n \rangle = \delta_{mn}
   $$</li>
</ol>
<h4 id="332-representations-of-states-and-observables">
<a class="header-anchor" href="#332-representations-of-states-and-observables"></a>
3.3.2 Representations of States and Observables
</h4><ol>
<li>
<p><strong>Coordinate representation</strong>
</p>
$$
   \psi(x) = \langle x|\psi\rangle
   $$<p>
</p>
$$
   \hat{x} \psi(x) = x \psi(x), \quad \hat{p}_x \psi(x) = -i\hbar \frac{\partial}{\partial x}\psi(x)
   $$</li>
<li>
<p><strong>Momentum representation</strong>
</p>
$$
   \phi(p) = \langle p|\psi\rangle
   $$<p>
</p>
$$
   \hat{p} \phi(p) = p \phi(p), \quad \hat{x} \phi(p) = i\hbar \frac{\partial}{\partial p}\phi(p)
   $$</li>
<li>
<p><strong>Energy representation</strong>
</p>
$$
   |\psi\rangle = \sum_n c_n |E_n\rangle, \quad c_n = \langle E_n|\psi\rangle
   $$</li>
<li>
<p><strong>Transformations between representations</strong>
</p>
$$
   \phi(p) = \frac{1}{\sqrt{2\pi\hbar}} \int_{-\infty}^{\infty} \psi(x) e^{-ipx/\hbar} dx
   $$<p>
</p>
$$
   \psi(x) = \frac{1}{\sqrt{2\pi\hbar}} \int_{-\infty}^{\infty} \phi(p) e^{ipx/\hbar} dp
   $$</li>
</ol>
<h3 id="34-orbital-angular-momentum">
<a class="header-anchor" href="#34-orbital-angular-momentum"></a>
3.4 Orbital Angular Momentum
</h3><h4 id="341-angular-momentum">
<a class="header-anchor" href="#341-angular-momentum"></a>
3.4.1 Angular Momentum
</h4><ol>
<li>
<p><strong>Definition</strong>
</p>
$$
   \vec{L} = \vec{r} \times \vec{p},\qquad
   \hat{\vec{L}} = \hat{\vec{r}} \times \hat{\vec{p}}
   $$</li>
<li>
<p><strong>Components</strong>
</p>
$$
   \hat{L}_x = y\hat{p}_z - z\hat{p}_y, \quad
   \hat{L}_y = z\hat{p}_x - x\hat{p}_z, \quad
   \hat{L}_z = x\hat{p}_y - y\hat{p}_x
   $$</li>
<li>
<p><strong>Commutation relations</strong>
</p>
$$
   [\hat{L}_x, \hat{L}_y] = i\hbar \hat{L}_z, \quad
   [\hat{L}_y, \hat{L}_z] = i\hbar \hat{L}_x, \quad
   [\hat{L}_z, \hat{L}_x] = i\hbar \hat{L}_y
   $$<p>
</p>
$$
   \hat{L}^2 = \hat{L}_x^2 + \hat{L}_y^2 + \hat{L}_z^2
   $$</li>
</ol>
<h4 id="342-conservation-of-angular-momentum">
<a class="header-anchor" href="#342-conservation-of-angular-momentum"></a>
3.4.2 Conservation of Angular Momentum
</h4><ol>
<li>
<p><strong>Condition for conservation</strong>
</p>
$$
   [\hat{H}, \hat{L}_i] = 0 \quad \Rightarrow \quad \hat{L}_i \ \text{is conserved}
   $$</li>
<li>
<p><strong>Spherically symmetric potential</strong>
</p>
$$
   [\hat{H}, \hat{L}^2] = 0, \quad [\hat{H}, \hat{L}_z] = 0
   $$</li>
</ol>
<h4 id="343-calculation-of-orbital-angular-momentum">
<a class="header-anchor" href="#343-calculation-of-orbital-angular-momentum"></a>
3.4.3 Calculation of Orbital Angular Momentum
</h4><ol>
<li>
<p><strong>Eigenvalue equations</strong>
</p>
$$
   \hat{L}^2 Y_{lm}(\theta,\varphi) = l(l+1)\hbar^2 Y_{lm}(\theta,\varphi)
   $$<p>
</p>
$$
   \hat{L}_z Y_{lm}(\theta,\varphi) = m\hbar Y_{lm}(\theta,\varphi)
   $$</li>
<li>
<p><strong>Eigenvalues</strong>
</p>
$$
   L = \sqrt{l(l+1)} \hbar,\qquad L_z = m\hbar
   $$</li>
<li>
<p><strong>Physical meaning</strong><br>
The quantum numbers $l$ and $m$ determine the magnitude of the orbital angular momentum and its $z$-component, respectively.</p>
</li>
</ol>
<h2 id="chapter-4-perturbation-theory-and-its-applications">
<a class="header-anchor" href="#chapter-4-perturbation-theory-and-its-applications"></a>
Chapter 4 Perturbation Theory and Its Applications
</h2><h3 id="41-time-independent-perturbation-theory">
<a class="header-anchor" href="#41-time-independent-perturbation-theory"></a>
4.1 Time-Independent Perturbation Theory
</h3><h4 id="411-nondegenerate-perturbation-theory">
<a class="header-anchor" href="#411-nondegenerate-perturbation-theory"></a>
4.1.1 Nondegenerate Perturbation Theory
</h4><ol>
<li>
<p><strong>Basic idea</strong>
</p>
$$
   \hat{H} = \hat{H}^{(0)} + \lambda \hat{H}'
   $$</li>
<li>
<p><strong>Energy corrections</strong>
</p>
$$
   E_n^{(1)} = \langle \psi_n^{(0)} | \hat{H}' | \psi_n^{(0)} \rangle
   $$<p>
</p>
$$
   E_n^{(2)} = \sum_{m \neq n} \frac{|\langle \psi_m^{(0)} | \hat{H}' | \psi_n^{(0)} \rangle|^2}{E_n^{(0)} - E_m^{(0)}}
   $$</li>
<li>
<p><strong>Wave-function correction</strong>
</p>
$$
   \psi_n^{(1)} = \sum_{m \neq n} \frac{\langle \psi_m^{(0)} | \hat{H}' | \psi_n^{(0)} \rangle}{E_n^{(0)} - E_m^{(0)}} \psi_m^{(0)}
   $$</li>
</ol>
<h4 id="412-degenerate-perturbation-theory">
<a class="header-anchor" href="#412-degenerate-perturbation-theory"></a>
4.1.2 Degenerate Perturbation Theory
</h4><ol>
<li>
<p><strong>Origin of the problem</strong><br>
If the zeroth-order energy corresponds to multiple orthogonal eigenstates, the state is degenerate, and the nondegenerate formulas fail.</p>
</li>
<li>
<p><strong>Method</strong>
</p>
$$
   H'_{ij} = \langle \psi_i^{(0)} | \hat{H}' | \psi_j^{(0)} \rangle
   $$<p>
Diagonalize this matrix inside the degenerate subspace.</p>
</li>
<li>
<p><strong>Result</strong><br>
The first-order energy corrections are the eigenvalues of $H'_{ij}$, and the corrected states are the corresponding linear combinations.</p>
</li>
</ol>
<hr>
<h3 id="42-time-dependent-perturbation-theory">
<a class="header-anchor" href="#42-time-dependent-perturbation-theory"></a>
4.2 Time-Dependent Perturbation Theory
</h3><ol>
<li>
<p><strong>Basic framework</strong>
</p>
$$
   \hat{H}(t) = \hat{H}^{(0)} + \hat{H}'(t)
   $$</li>
<li>
<p><strong>State expansion</strong>
</p>
$$
   |\psi(t)\rangle = \sum_n c_n(t) e^{-iE_n^{(0)}t/\hbar} |\psi_n^{(0)}\rangle
   $$</li>
<li>
<p><strong>Transition probability amplitude</strong>
</p>
$$
   c_f^{(1)}(t) = \frac{1}{i\hbar} \int_0^t \langle \psi_f^{(0)} | \hat{H}'(t') | \psi_i^{(0)} \rangle 
   e^{i\omega_{fi} t'} dt'
   $$<p>
where $\omega_{fi} = (E_f^{(0)} - E_i^{(0)})/\hbar$.</p>
</li>
<li>
<p><strong>Fermi&rsquo;s golden rule</strong>
</p>
$$
   W_{i \to f} = \frac{2\pi}{\hbar} \, |\langle f | \hat{H}' | i \rangle|^2 \, \rho(E_f)
   $$</li>
</ol>
<p><strong>Summary</strong></p>
<ul>
<li>Time-independent perturbation theory corrects energies and wave functions for static perturbations.</li>
<li>Time-dependent perturbation theory describes transitions between energy levels, such as radiation absorption and emission.</li>
</ul>
<h2 id="electron-spin">
<a class="header-anchor" href="#electron-spin"></a>
Electron Spin
</h2><h3 id="experimental-discovery-of-electron-spin">
<a class="header-anchor" href="#experimental-discovery-of-electron-spin"></a>
Experimental Discovery of Electron Spin
</h3><ol>
<li>
<p><strong>Stern-Gerlach experiment</strong><br>
Passing a beam of silver atoms through a nonuniform magnetic field produces two trajectories, revealing an intrinsic angular momentum beyond orbital angular momentum.</p>
</li>
<li>
<p><strong>Experimental conclusions</strong></p>
<ul>
<li>The spin quantum number is $s = 1/2$.</li>
<li>The two possible spin projections are $m_s = \pm 1/2$.</li>
<li>Spin contributes an additional magnetic moment:

$$
     \vec{\mu}_s = -g_s \frac{e}{2m_e} \vec{S}, \quad g_s \approx 2
     $$</li>
</ul>
</li>
</ol>
<h3 id="theory-of-electron-spin">
<a class="header-anchor" href="#theory-of-electron-spin"></a>
Theory of Electron Spin
</h3><ol>
<li>
<p><strong>Quantum description</strong></p>
<ul>
<li>Spin is intrinsic angular momentum and does not depend on spatial coordinates.</li>
<li>Its operators satisfy

$$
     [\hat{S}_i, \hat{S}_j] = i\hbar \epsilon_{ijk} \hat{S}_k
     $$</li>
</ul>
</li>
<li>
<p><strong>Physical meaning</strong></p>
<ul>
<li>Spin determines the magnetic behavior of electrons.</li>
<li>Quantized spin leads to Fermi-Dirac statistics and the Pauli exclusion principle.</li>
</ul>
</li>
</ol>
<hr>
<h3 id="spin-angular-momentum">
<a class="header-anchor" href="#spin-angular-momentum"></a>
Spin Angular Momentum
</h3><h4 id="spin-operators">
<a class="header-anchor" href="#spin-operators"></a>
Spin Operators
</h4><ol>
<li>
<p><strong>Spin components</strong>
</p>
$$
   \hat{S}_x, \hat{S}_y, \hat{S}_z
   $$<p>
satisfying
</p>
$$
   [\hat{S}_x, \hat{S}_y] = i\hbar \hat{S}_z, \quad \text{cyclic symmetry}
   $$</li>
<li>
<p><strong>Total spin operator</strong>
</p>
$$
   \hat{S}^2 = \hat{S}_x^2 + \hat{S}_y^2 + \hat{S}_z^2
   $$<p>
with
</p>
$$
   \hat{S}^2 |\chi_s\rangle = s(s+1)\hbar^2 |\chi_s\rangle
   $$</li>
</ol>
<h4 id="matrix-representation-of-eigenfunctions">
<a class="header-anchor" href="#matrix-representation-of-eigenfunctions"></a>
Matrix Representation of Eigenfunctions
</h4><ol>
<li>
<p><strong>Spin-$1/2$ particles</strong>
</p>
$$
   |\uparrow\rangle = \begin{pmatrix}1\\0\end{pmatrix}, \quad
   |\downarrow\rangle = \begin{pmatrix}0\\1\end{pmatrix}
   $$</li>
<li>
<p><strong>Pauli-matrix form of the spin operators</strong>
</p>
$$
   \hat{S}_x = \frac{\hbar}{2} \sigma_x, \quad
   \hat{S}_y = \frac{\hbar}{2} \sigma_y, \quad
   \hat{S}_z = \frac{\hbar}{2} \sigma_z
   $$<p>
where
</p>
$$
   \sigma_x = \begin{pmatrix}0 & 1\\ 1 & 0\end{pmatrix},\quad
   \sigma_y = \begin{pmatrix}0 & -i\\ i & 0\end{pmatrix},\quad
   \sigma_z = \begin{pmatrix}1 & 0\\ 0 & -1\end{pmatrix}
   $$</li>
</ol>
<h4 id="theory-of-angular-momentum-coupling">
<a class="header-anchor" href="#theory-of-angular-momentum-coupling"></a>
Theory of Angular Momentum Coupling
</h4><ol>
<li>
<p><strong>Spin-orbit coupling</strong>
</p>
$$
   \hat{H}_{\text{SO}} = \xi(r)\, \vec{L} \cdot \vec{S}
   $$<p>
which produces fine-structure splitting.</p>
</li>
<li>
<p><strong>Total angular momentum</strong>
</p>
$$
   \vec{J} = \vec{L} + \vec{S}, \quad
   \hat{J}^2 = (\hat{L}+\hat{S})^2
   $$<p>
with eigenstates $|j, m_j\rangle$ satisfying
</p>
$$
   \hat{J}^2 |j, m_j\rangle = j(j+1)\hbar^2 |j, m_j\rangle, \quad
   \hat{J}_z |j, m_j\rangle = m_j \hbar |j, m_j\rangle
   $$</li>
<li>
<p><strong>Coupling result</strong></p>
<ul>
<li>$j = l \pm s$, $m_j = -j, -j+1, ..., j$.</li>
<li>Spin-orbit coupling is an important source of the fine structure of atomic spectra.</li>
</ul>
</li>
</ol>
<h2 id="principle-of-indistinguishability">
<a class="header-anchor" href="#principle-of-indistinguishability"></a>
Principle of Indistinguishability
</h2><h3 id="systems-of-identical-particles">
<a class="header-anchor" href="#systems-of-identical-particles"></a>
Systems of Identical Particles
</h3><h4 id="concepts-and-principles">
<a class="header-anchor" href="#concepts-and-principles"></a>
Concepts and Principles
</h4><ol>
<li>
<p><strong>Definition of identical particles</strong><br>
If two particles are completely identical in physical properties such as mass, charge, and spin, and cannot be distinguished by any experiment, they are called identical particles.</p>
</li>
<li>
<p><strong>Principle of indistinguishability</strong><br>
The physical laws are invariant under exchange of identical particles. Exchanging the positions and spins of any two identical particles leaves the Hamiltonian and observables unchanged.</p>
</li>
</ol>
<h4 id="hamiltonian-of-a-system-of-identical-particles">
<a class="header-anchor" href="#hamiltonian-of-a-system-of-identical-particles"></a>
Hamiltonian of a System of Identical Particles
</h4><ol>
<li>
<p><strong>Form</strong>
</p>
$$
   \hat{H} = \sum_{i=1}^N \hat{T}_i + \sum_{i<j} V(\vec{r}_i - \vec{r}_j)
   $$</li>
<li>
<p><strong>Symmetry</strong>
</p>
$$
   [\hat{H}, \hat{P}_{ij}] = 0
   $$<p>
where $\hat{P}_{ij}$ is the exchange operator of particles $i$ and $j$.</p>
</li>
</ol>
<h4 id="wave-functions-of-identical-particle-systems">
<a class="header-anchor" href="#wave-functions-of-identical-particle-systems"></a>
Wave Functions of Identical-Particle Systems
</h4><ol>
<li>
<p><strong>Symmetry requirement</strong>
</p>
$$
   \hat{P}_{ij} \Psi(\dots, \vec{r}_i, \vec{r}_j, \dots) =
   \pm \Psi(\dots, \vec{r}_i, \vec{r}_j, \dots)
   $$<ul>
<li><strong>+</strong> for bosons, whose wave functions are symmetric</li>
<li><strong>-</strong> for fermions, whose wave functions are antisymmetric</li>
</ul>
</li>
<li>
<p><strong>Construction of many-particle wave functions</strong></p>
<ul>
<li>Bosons: symmetrized sum</li>
<li>Fermions: antisymmetrized determinant (Slater determinant)

$$
     \Psi(\vec{r}_1, \dots, \vec{r}_N) = \frac{1}{\sqrt{N!}}
     \begin{vmatrix}
     \psi_1(\vec{r}_1) & \cdots & \psi_1(\vec{r}_N) \\
     \vdots & \ddots & \vdots \\
     \psi_N(\vec{r}_1) & \cdots & \psi_N(\vec{r}_N)
     \end{vmatrix}
     $$</li>
</ul>
</li>
</ol>
<hr>
<h3 id="pauli-exclusion-principle">
<a class="header-anchor" href="#pauli-exclusion-principle"></a>
Pauli Exclusion Principle
</h3><ol>
<li>
<p><strong>Content of the principle</strong><br>
For identical fermions with half-integer spin, no two particles may occupy the same quantum state:
</p>
$$
   \Psi(\text{same quantum state}) = 0
   $$</li>
<li>
<p><strong>Physical meaning</strong><br>
It explains the arrangement of electrons in atomic orbitals and underlies atomic structure, chemical properties, and Fermi-gas behavior.</p>
</li>
<li>
<p><strong>Examples</strong></p>
<ul>
<li>In atoms, each orbital can hold at most two electrons with opposite spins.</li>
<li>In metals, electrons form a Fermi level that determines electrical and thermal properties.</li>
</ul>
</li>
</ol>

        
        <hr><p>Published on 2025-09-05 at <a href='https://www.guzhengsvt.cn/'>孤筝の温暖小家</a>, last modified on 2025-09-05</p><p>All articles on this blog are licensed under the BY-NC-SA license agreement unless otherwise stated. Please indicate the source when reprinting!</p>]]></description><category>Physics</category></item><item><title>Digital Signal Processing</title><link>https://www.guzhengsvt.cn/en/post/physics/%E6%95%B0%E5%AD%97%E4%BF%A1%E5%8F%B7%E5%A4%84%E7%90%86/</link><pubDate>Wed, 04 Sep 2024 23:44:15 +0800</pubDate><author>lvbowen040427@163.com (孤筝)</author><guid>https://www.guzhengsvt.cn/en/post/physics/%E6%95%B0%E5%AD%97%E4%BF%A1%E5%8F%B7%E5%A4%84%E7%90%86/</guid><description>
<![CDATA[<h1>Digital Signal Processing</h1><p>Author: 孤筝(lvbowen040427@163.com)</p>
        
          <h2 id="fundamental-concepts-of-digital-signal-processing">
<a class="header-anchor" href="#fundamental-concepts-of-digital-signal-processing"></a>
Fundamental Concepts of Digital Signal Processing
</h2><h3 id="signal-classification">
<a class="header-anchor" href="#signal-classification"></a>
Signal Classification
</h3><ol>
<li><strong>Continuous-time signal</strong>: Analog signal, continuous in the time domain.</li>
<li><strong>Discrete-time signal</strong>: Continuous in amplitude but discrete in time.</li>
<li><strong>Amplitude-discrete signal</strong>: Discrete in amplitude but continuous in time.</li>
<li><strong>Digital signal</strong>: Discrete in both amplitude and time.</li>
</ol>
<h4 id="differences">
<a class="header-anchor" href="#differences"></a>
Differences
</h4><p>The distinction between a discrete-time signal and a digital signal lies solely in the <strong>quantization error</strong> present in digital signals.</p>
<h3 id="implementation-methods-of-digital-signal-processing">
<a class="header-anchor" href="#implementation-methods-of-digital-signal-processing"></a>
Implementation Methods of Digital Signal Processing
</h3><p>The primary target of digital signal processing is <strong>digital signals</strong>, achieved through <strong>numerical operations</strong>.</p>
<h4 id="software-implementation">
<a class="header-anchor" href="#software-implementation"></a>
Software Implementation
</h4><p>Programs are written based on principles and algorithms and executed on general-purpose computers.</p>
<ul>
<li><strong>Advantages</strong>: Flexibility</li>
<li><strong>Disadvantages</strong>: Slow computation speed, difficult to achieve real-time processing.</li>
<li><strong>Suitable for</strong>: Algorithm research and simulation.</li>
</ul>
<h4 id="hardware-implementation">
<a class="header-anchor" href="#hardware-implementation"></a>
Hardware Implementation
</h4><p>Hardware structures are designed according to specific requirements and algorithms, using basic components such as multipliers, adders, delay units, controllers, memory, and input/output interfaces.</p>
<ul>
<li><strong>Advantages</strong>: Fast computation speed, capable of real-time processing.</li>
<li><strong>Disadvantages</strong>: Inflexibility</li>
</ul>
<p>Hardware implementation refers to selecting an appropriate <strong>DSP chip</strong>, equipped with suitable software and language for the task, to achieve specific signal processing functions.</p>
<h4 id="dedicated-chips">
<a class="header-anchor" href="#dedicated-chips"></a>
Dedicated Chips
</h4><p>Using <strong>Digital Signal Processing (DSP) chips</strong> is currently the fastest-growing and most widely applied method. DSP chips offer significant advantages over general-purpose microcontrollers:</p>
<ul>
<li><strong>Internal multipliers and accumulators</strong> tailored for digital signal processing.</li>
<li><strong>Pipeline architecture</strong>, parallel structures, and multiple buses.</li>
<li>Specialized instructions optimized for digital signal processing, enabling high-speed computation.</li>
</ul>
<p>For ultra-high-speed real-time systems where DSP chips are insufficient, <strong>Field-Programmable Gate Arrays (FPGAs)</strong> or custom ASICs (Application-Specific Integrated Circuits) should be employed.</p>
<h3 id="characteristics-of-digital-signal-processing">
<a class="header-anchor" href="#characteristics-of-digital-signal-processing"></a>
Characteristics of Digital Signal Processing
</h3><p>Compared to analog signal processing, digital signal processing offers:</p>
<ol>
<li><strong>Flexibility</strong></li>
<li><strong>High precision and stability</strong></li>
<li><strong>Ease of large-scale integration</strong></li>
<li><strong>Capability to perform functions unattainable by analog systems</strong>, such as storage, complex transformations, and operations.</li>
</ol>
<h3 id="signal-dimensionality">
<a class="header-anchor" href="#signal-dimensionality"></a>
Signal Dimensionality
</h3><p>A signal is typically a function of one or more independent variables.</p>
<ul>
<li><strong>One-dimensional signal</strong>: Only one independent variable.</li>
<li><strong>Multidimensional signal</strong>: Two or more independent variables.</li>
</ul>
<h2 id="discrete-time-signals-and-systems">
<a class="header-anchor" href="#discrete-time-signals-and-systems"></a>
Discrete-Time Signals and Systems
</h2><h3 id="discrete-time-signals">
<a class="header-anchor" href="#discrete-time-signals"></a>
Discrete-Time Signals
</h3><p>In practice, signals are generally analog. <strong>Uniform sampling</strong> (equally spaced sampling) converts them into discrete-time signals.</p>
<p>For an analog signal \( x_a(t) \), discrete time points \( t_n \).<br>
With uniform sampling interval \( T \), \( t_n = nT \):<br>
</p>
$$
x(n) = x_a(t) \big|_{t=nT} = x_a(nT), \quad -\infty  \lt  n  \lt  \infty  
$$<p><br>
Here, \( x(n) \) is called a <strong>discrete-time signal</strong>, where \( n \) is an integer, forming a sequence:<br>
</p>
$$
x(n) = \{\cdots, x_a(-2T), x_a(-T), x_a(0), x_a(T), x_a(2T), \cdots \}  
$$<p><br>
Discrete-time signals are also referred to as <strong>sequences</strong>.</p>
<h4 id="representation-methods-for-sequences">
<a class="header-anchor" href="#representation-methods-for-sequences"></a>
Representation Methods for Sequences
</h4><h5 id="set-notation">
<a class="header-anchor" href="#set-notation"></a>
Set Notation
</h5><p>A set of numbers is denoted by \( \{\cdot\} \). A discrete-time signal can be represented as an ordered set of numbers.<br>
The underlined element in the set indicates the sample value at \( n=0 \).</p>
<h5 id="formula-representation">
<a class="header-anchor" href="#formula-representation"></a>
Formula Representation
</h5><p>Example:<br>
</p>
$$
x(n) = a^{|n|}, \quad 0  \lt  a  \lt  1, \quad -\infty  \lt  n  \lt  \infty  
$$<h5 id="graphical-representation">
<a class="header-anchor" href="#graphical-representation"></a>
Graphical Representation
</h5><p>The horizontal axis represents \( n \), and the vertical axis represents the value of \( x \), with dots atop vertical lines.</p>
<h4 id="common-standard-sequences">
<a class="header-anchor" href="#common-standard-sequences"></a>
Common Standard Sequences
</h4><h5 id="unit-impulse-sequence">
<a class="header-anchor" href="#unit-impulse-sequence"></a>
Unit Impulse Sequence \( \delta(n) \)
</h5>$$
\delta(n) = 
\begin{cases} 
1 & n = 0 \\
0 & n \neq 0 
\end{cases}  
$$<p><br>
Also called the <strong>unit sample sequence</strong>, distinct from the unit impulse signal \( \delta(t) \).</p>
<h5 id="unit-step-sequence">
<a class="header-anchor" href="#unit-step-sequence"></a>
Unit Step Sequence \( u(n) \)
</h5>$$
u(n) = 
\begin{cases} 
1 & n \geq 0 \\
0 & n  \lt  0 
\end{cases}  
$$<p><br>
Relationships:<br>
</p>
$$
\delta(n) = u(n) - u(n-1)  
$$<p><br>
</p>
$$
u(n) = \sum_{k=0}^{\infty} \delta(n - k)  
$$<h5 id="rectangular-sequence">
<a class="header-anchor" href="#rectangular-sequence"></a>
Rectangular Sequence \( R_N(n) \)
</h5>$$
R_N(n) = 
\begin{cases} 
1 & 0 \leq n \leq N-1 \\
0 & \text{Otherwise} 
\end{cases}  
$$<p><br>
\( N \) is the length of the rectangular sequence. It can be expressed using the unit step sequence:<br>
</p>
$$
R_N(n) = u(n) - u(n - N)  
$$<h5 id="real-exponential-sequence">
<a class="header-anchor" href="#real-exponential-sequence"></a>
Real Exponential Sequence
</h5>$$
x(n) = a^n u(n), \quad a \text{ is a real number}  
$$<ul>
<li><strong>Convergent sequence</strong>: \( |a|  \lt  1 \)</li>
<li><strong>Divergent sequence</strong>: \( |a|  \gt  1 \)</li>
</ul>
<h5 id="sinusoidal-sequence">
<a class="header-anchor" href="#sinusoidal-sequence"></a>
Sinusoidal Sequence
</h5>$$
x(n) = \sin(\omega n)  
$$<p><br>
Here, \( \omega \) is the <strong>digital frequency</strong> (units: radians, \( rad \)), representing the rate of change (phase shift between adjacent samples).</p>
<p><strong>Analog angular frequency \( \varOmega \)</strong><br>
If the sinusoidal sequence is derived from sampling an analog signal \( x_a(t) = \sin(\varOmega t) \):<br>
</p>
$$
x(n) = x_a(t) \big|_{t=nT} = \sin(\varOmega nT) = \sin(\omega n)  
$$<p><br>
The relationship between digital and analog frequencies is:<br>
</p>
$$
\omega = \varOmega T  
$$<p><br>
Given the sampling frequency \( F_s = \frac{1}{T} \):<br>
</p>
$$
\omega = \frac{\varOmega}{F_s}  
$$<p><br>
<em>Digital frequency is the normalized analog angular frequency relative to the sampling frequency.</em></p>
<h5 id="complex-exponential-sequence">
<a class="header-anchor" href="#complex-exponential-sequence"></a>
Complex Exponential Sequence
</h5>$$
x(n) = e^{(\sigma + j \omega_0)n} = \cos(\omega_0 n) + j \sin(\omega_0 n)  
$$<p><br>
Since \( n \) is an integer, both sinusoidal and complex exponential sequences are periodic with period \( 2\pi \).</p>
<h5 id="periodic-sequence">
<a class="header-anchor" href="#periodic-sequence"></a>
Periodic Sequence
</h5><p>If for all \( n \), there exists a <strong>smallest positive integer</strong> \( N \) such that:<br>
</p>
$$
x(n) = x(n + N), \quad -\infty  \lt  n  \lt  \infty  
$$<p><br>
then the sequence \( x(n) \) is <strong>periodic</strong> with period \( N \).</p>
<h4 id="sequence-operations">
<a class="header-anchor" href="#sequence-operations"></a>
Sequence Operations
</h4><h5 id="addition-and-multiplication">
<a class="header-anchor" href="#addition-and-multiplication"></a>
Addition and Multiplication
</h5><h5 id="shifting-flipping-and-scaling">
<a class="header-anchor" href="#shifting-flipping-and-scaling"></a>
Shifting, Flipping, and Scaling
</h5><h3 id="discrete-time-systems">
<a class="header-anchor" href="#discrete-time-systems"></a>
Discrete-Time Systems
</h3><p>For a system with input \( x(n) \), output \( y(n) \), and operation \( T[\cdot] \):<br>
</p>
$$
y(n) = T[x(n)]  
$$<h4 id="linear-systems">
<a class="header-anchor" href="#linear-systems"></a>
Linear Systems
</h4><p>A system is <strong>linear</strong> if its input-output relationship satisfies the principle of superposition.</p>
<h5 id="additivity">
<a class="header-anchor" href="#additivity"></a>
Additivity
</h5>$$
y_1(n) = T[x_1(n)], \quad y_2(n) = T[x_2(n)]  
$$<p><br>
</p>
$$
T[x_1(n) + x_2(n)] = y_1(n) + y_2(n)  
$$<h5 id="homogeneity-scaling">
<a class="header-anchor" href="#homogeneity-scaling"></a>
Homogeneity (Scaling)
</h5>$$
T[a \cdot x(n)] = a \cdot y(n)  
$$<h4 id="time-invariant-systems">
<a class="header-anchor" href="#time-invariant-systems"></a>
Time-Invariant Systems
</h4><p>A system is <strong>time-invariant</strong> if its operation \( T[\cdot] \) does not change over time, i.e., the system&rsquo;s response is independent of when the input is applied.<br>
</p>
$$
y(n) = T[x(n)]  
$$<p><br>
</p>
$$
y(n - n_0) = T[x(n - n_0)]  
$$<h4 id="characteristics-of-linear-time-invariant-lti-systems">
<a class="header-anchor" href="#characteristics-of-linear-time-invariant-lti-systems"></a>
Characteristics of Linear Time-Invariant (LTI) Systems
</h4><p>Total response = Zero-input response + Zero-state response</p>
<h5 id="unit-impulse-response">
<a class="header-anchor" href="#unit-impulse-response"></a>
Unit Impulse Response
</h5><p>With zero initial state (no zero-input response):<br>
</p>
$$
h(n) = T[\delta(n)]  
$$<p><br>
For any input \( x(n) \):<br>
</p>
$$
x(n) = \sum_{m=-\infty}^{\infty} x(m) \delta(n - m)  
$$<p><br>
Thus, the output is:<br>
</p>
$$
\begin{align} 
y(n) 
&= T[x(n)] \\ 
&= T\left[\sum_{m=-\infty}^{\infty} x(m) \delta(n - m)\right] \\ 
&= \sum_{m=-\infty}^{\infty} x(m) T[\delta(n - m)] \\ 
&= \sum_{m=-\infty}^{\infty} x(m) h(n - m) \\ 
&= x(n) * h(n) 
\end{align}  
$$<p><br>
This is the <strong>convolution sum</strong>. For details, refer to <em>Signals and Systems</em>.</p>
<h4 id="causality-of-systems">
<a class="header-anchor" href="#causality-of-systems"></a>
Causality of Systems
</h4><p><strong>Definition</strong>: A system is <strong>causal</strong> if its output at time \( n \) depends only on the input at time \( n \) and prior inputs, not future inputs.</p>
<p><strong>Necessary and Sufficient Condition</strong>:<br>
The unit impulse response satisfies:<br>
</p>
$$
h(n) = 0 \quad \text{for} \quad n  \lt  0  
$$<h4 id="stability-of-systems">
<a class="header-anchor" href="#stability-of-systems"></a>
Stability of Systems
</h4><p><strong>Definition</strong>: A system is <strong>stable</strong> if every bounded input produces a bounded output (BIBO stability).</p>
<p><strong>Necessary and Sufficient Condition</strong>:<br>
The unit impulse response is absolutely summable:<br>
</p>
$$
\sum_{m=-\infty}^{\infty} |h(n)|  \lt  \infty  
$$<h4 id="linear-constant-coefficient-difference-equations">
<a class="header-anchor" href="#linear-constant-coefficient-difference-equations"></a>
Linear Constant-Coefficient Difference Equations
</h4>
        
        <hr><p>Published on 2024-09-04 at <a href='https://www.guzhengsvt.cn/'>孤筝の温暖小家</a>, last modified on 2024-09-04</p><p>All articles on this blog are licensed under the BY-NC-SA license agreement unless otherwise stated. Please indicate the source when reprinting!</p>]]></description><category>Physics</category></item><item><title>Fundamentals of Aerodynamics</title><link>https://www.guzhengsvt.cn/en/post/physics/%E7%A9%BA%E6%B0%94%E5%8A%A8%E5%8A%9B%E5%AD%A6%E5%9F%BA%E7%A1%80/</link><pubDate>Mon, 02 Sep 2024 22:13:45 +0800</pubDate><author>lvbowen040427@163.com (孤筝)</author><guid>https://www.guzhengsvt.cn/en/post/physics/%E7%A9%BA%E6%B0%94%E5%8A%A8%E5%8A%9B%E5%AD%A6%E5%9F%BA%E7%A1%80/</guid><description>
<![CDATA[<h1>Fundamentals of Aerodynamics</h1><p>Author: 孤筝(lvbowen040427@163.com)</p>
        
          <h2 id="fundamentals-of-hydrostatics">
<a class="header-anchor" href="#fundamentals-of-hydrostatics"></a>
Fundamentals of Hydrostatics
</h2><h3 id="gas-properties">
<a class="header-anchor" href="#gas-properties"></a>
Gas Properties
</h3><h3 id="compressibility">
<a class="header-anchor" href="#compressibility"></a>
Compressibility
</h3><h2 id="">
<a class="header-anchor" href="#"></a>

</h2><h3 id="bulk-modulus">
<a class="header-anchor" href="#bulk-modulus"></a>
Bulk Modulus
</h3><p>Definition: The increase in pressure required to produce a unit relative change in volume.</p>
$$
E=-\frac{dp}{dV/V}
$$<p>For a given mass of gas, volume is inversely proportional to density, i.e.,</p>
$$
\frac{d \rho}{\rho}=-\frac{dV}{V}
$$<p>Substituting back yields</p>
$$
E=\rho \frac{dp}{d \rho}
$$<p>Bulk modulus of water at room temperature: $2.1 \times 10^9N/m^2$</p>
<p>Under normal conditions, water can be considered an incompressible fluid.</p>
<h3 id="viscosity">
<a class="header-anchor" href="#viscosity"></a>
Viscosity
</h3><h2 id="">
<a class="header-anchor" href="#"></a>

</h2><h3 id="newtons-law-of-viscosity">
<a class="header-anchor" href="#newtons-law-of-viscosity"></a>
Newton&rsquo;s Law of Viscosity
</h3><p>The frictional resistance generated by fluid motion is proportional to the contact area</p>
$$
\tau =\mu \frac{du}{d \vec{n}}
$$<p>$\tau$: Frictional <strong>stress</strong>, the frictional resistance per unit area</p>
<p>$\vec{n}$: Normal direction of the contact surface</p>
<p>$\mu$: Proportionality constant, known as the <strong>coefficient of viscosity</strong> of the fluid, with units of $N \cdot s/m^2$</p>
<p>$\frac{du}{d \vec{n}}$: Velocity gradient</p>
<p>The coefficient of viscosity varies for different fluid media and changes with temperature, being largely independent of pressure.</p>
<p>The coefficient of viscosity for gases increases with rising temperature.</p>
<h3 id="sutherlands-formula">
<a class="header-anchor" href="#sutherlands-formula"></a>
Sutherland&rsquo;s Formula
</h3><p>One of the approximate formulas describing the relationship between air viscosity coefficient and temperature is Sutherland&rsquo;s Formula.</p>
$$
\frac{\mu}{\mu_0}=(\frac{T}{288.15})^{1.5}\frac{288.15+C}{T+C}
$$<p>$\mu_0$: Air viscosity coefficient at temperature $288.15K$</p>
<p>$C$: Constant, with a value of $110.4K$</p>
<h3 id="kinematic-viscosity-coefficient">
<a class="header-anchor" href="#kinematic-viscosity-coefficient"></a>
Kinematic Viscosity Coefficient
</h3>$$
\nu=\frac{\mu}{\rho}
$$<p>$\nu$: Kinematic viscosity coefficient, unit is $m^2/s$</p>
<p>$\mu$: Viscosity coefficient</p>
<p>$\rho$: Density</p>
<h3 id="thermal-conductivity">
<a class="header-anchor" href="#thermal-conductivity"></a>
Thermal Conductivity
</h3><p>Definition: When there is a temperature gradient in a certain direction within a gas, heat will transfer from the higher-temperature region to the lower-temperature region. This property is referred to as the thermal conductivity of the gas.</p>
<p>The amount of heat transferred per unit time is proportional to the heat transfer area and the temperature gradient along the direction of heat flow, expressed as:</p>
$$
q=-\lambda \frac{\partial T}{\partial \vec{n}}
$$<p>$q$: Heat per unit time through a unit area, unit: $kJ/(m^2 \cdot s)$</p>
<p>$\frac{\partial T}{\partial \vec{n}}$: Temperature gradient, unit: $K/m$</p>
<p>$\lambda$: Thermal conductivity coefficient, unit: $kJ/(m \cdot K \cdot s)$</p>
<p>The negative sign indicates that the direction of heat transfer is always opposite to the direction of the temperature gradient.</p>
<h3 id="fluid-classification">
<a class="header-anchor" href="#fluid-classification"></a>
Fluid Classification
</h3><h1 id="">
<a class="header-anchor" href="#"></a>

</h1><h3 id="continuum-hypothesis">
<a class="header-anchor" href="#continuum-hypothesis"></a>
Continuum Hypothesis
</h3><h3 id="ideal-fluid">
<a class="header-anchor" href="#ideal-fluid"></a>
Ideal Fluid
</h3><p>不考虑粘性，在这种模型中，流体微团不承受粘性力作用。常用于气体。</p>
<p>忽略粘性的气体称为理想气体。</p>
<h3 id="pressure-isotropy">
<a class="header-anchor" href="#pressure-isotropy"></a>
Pressure Isotropy
</h3><p>The pressure at a point in an ideal fluid is independent of the orientation of the pressure surface; it is merely a continuous function of spatial coordinates.</p>
<h3 id="incompressible-fluid">
<a class="header-anchor" href="#incompressible-fluid"></a>
Incompressible Fluid
</h3><p>不考虑气体压缩性或弹性，可认为体积弹性模数无穷大，或流体密度为常数。常用于液体。</p>
<p>求解不可压流体的流动规律，只需要服从力学定律，不需要考虑热力学关系。</p>
<p>对流速较低的气体，也可按不可压流体处理流动问题。</p>
<h3 id="adiabatic-fluid">
<a class="header-anchor" href="#adiabatic-fluid"></a>
Adiabatic Fluid
</h3><p>A model that does not consider the heat transfer properties of a fluid, treating the thermal conductivity coefficient of the fluid as zero. Low-speed flowing air typically has very small thermal conductivity values and can be regarded as adiabatic.</p>
<p>A gas model that disregards the heat conduction effects between gas microelements is referred to as an adiabatic gas.</p>
<h3 id="perfect-gas">
<a class="header-anchor" href="#perfect-gas"></a>
Perfect Gas
</h3><p>For any state, there exists a certain functional relationship between the pressure, density, and temperature of a gas</p>
$$
p=p(\rho,T)
$$<h3 id="equation-of-state-for-a-perfect-gas">
<a class="header-anchor" href="#equation-of-state-for-a-perfect-gas"></a>
Equation of State for a Perfect Gas
</h3>$$
p=\frac{\overline{R}}{m}\rho T
$$<p>$\overline{R}$: Universal gas constant, $8315m^2/(s^2 \cdot K)$</p>
<p>$m$: Relative molecular mass of a specific gas</p>
<p>When $R=\frac{\overline{R}}{m}$,</p>
$$
p=\rho R T
$$<p>$R$ is the gas constant, approximately $287.035m^2/(s^2 \cdot K)$ for air</p>
<h3 id="forces-on-a-fluid-particle">
<a class="header-anchor" href="#forces-on-a-fluid-particle"></a>
Forces on a Fluid Particle
</h3><h3 id="stress">
<a class="header-anchor" href="#stress"></a>
Stress
</h3><h3 id="shear-stress-friction">
<a class="header-anchor" href="#shear-stress-friction"></a>
Shear Stress (Friction)
</h3><h1 id="">
<a class="header-anchor" href="#"></a>

</h1><h3 id="volume-force">
<a class="header-anchor" href="#volume-force"></a>
Volume Force
</h3><ul>
<li>Gravity</li>
<li>Electromagnetic Force</li>
<li>Centrifugal Force</li>
</ul>
<h3 id="static-equilibrium-equation">
<a class="header-anchor" href="#static-equilibrium-equation"></a>
Static Equilibrium Equation
</h3><p>In a ==static fluid==, take a point $P$ with pressure $p$.</p>
<p>Construct a Cartesian coordinate system where the pressure at any point in the fluid is</p>
$$
p(x,y,z)
$$<p>Build a rectangular parallelepiped centered at $P$ with edges parallel to the coordinate axes and lengths $dx, dy, dz$.</p>
<p>Observing the $x$-axis direction, the forces on the two faces are respectively</p>
$$
[p(x_0,y_0,z_0)+(\frac{\partial p}{\partial x})(\frac{dx}{2})]dx dy
$$$$
[p(x_0,y_0,z_0)-(\frac{\partial p}{\partial x})(\frac{dx}{2})]dx dy
$$<p>The body force on the fluid element in the $x$-axis direction is</p>
$$
f_x \rho dx dy dz
$$<p>Here, $f_x$ is the component of the body force per unit mass in the $x$-axis direction.</p>
<p>Since the fluid is static, the fluid element is in equilibrium.</p>
<p>The force balance equation in the $x$-axis direction is</p>
$$
[p(x_0,y_0,z_0)-(\frac{\partial p}{\partial x})(\frac{dx}{2})]dx dy-[p(x_0,y_0,z_0)+(\frac{\partial p}{\partial x})(\frac{dx}{2})]dx dy+f_x \rho dx dy dz=0
$$<p>Simplifying gives</p>
$$
\frac{\partial p}{\partial x}=\rho f_x
$$$$
\frac{\partial p}{\partial y}=\rho f_y
$$$$
\frac{\partial p}{\partial z}=\rho f_z
$$<p>$\because$ The total differential of $p$ is</p>
$$
dp=\frac{\partial p}{\partial x}dx+\frac{\partial p}{\partial y}dy+\frac{\partial p}{\partial z}dz
$$<p>$\therefore$</p>
$$
dp=\rho(f_x dx+f_y dy+f_z dz)
$$<p>Define the ==body force potential function==</p>
$$
\varOmega=\varOmega(x,y,z)
$$<p>Its total differential is</p>
$$
d \varOmega=\frac{\partial \varOmega}{\partial x}dx+\frac{\partial \varOmega}{\partial y}dy+\frac{\partial \varOmega}{\partial z}dz
$$<p>Where $\frac{\partial \varOmega}{\partial x}=f_x$, $\frac{\partial \varOmega}{\partial y}=f_y$, $\frac{\partial \varOmega}{\partial z}=-f_z$.</p>
<p>From the above relations, we obtain</p>
$$
dp=-\rho d \varOmega
$$<p>Integrating both sides over $x, y, z$ gives</p>
$$
p=-\rho \varOmega+C(\text{constant})
$$$$
C=p+\rho \varOmega
$$<p>When the pressure $p_a$ at a known point A, the difference in body force potential $\varOmega_a-\varOmega$ between two points, and the density $\rho$ of the static fluid (uniform everywhere) are known, the pressure at any point can be determined from its body force potential $\varOmega$:</p>
$$
p=p_a+\rho (\varOmega_a-\varOmega)
$$<p><strong>Corollary:</strong> An isobaric surface in the fluid must also be an equipotential surface of the body force.</p>
<h2 id="atmosphere">
<a class="header-anchor" href="#atmosphere"></a>
Atmosphere
</h2><h3 id="atmospheric-layers">
<a class="header-anchor" href="#atmospheric-layers"></a>
Atmospheric Layers
</h3><h3 id="lower-atmosphere">
<a class="header-anchor" href="#lower-atmosphere"></a>
Lower Atmosphere
</h3><ul>
<li>Altitude: Sea level – 85 km</li>
<li>Characteristics: Uniform composition, nitrogen accounts for 78.1% of the total volume, oxygen accounts for 21% of the total volume</li>
</ul>
<h3 id="troposphere">
<a class="header-anchor" href="#troposphere"></a>
Troposphere
</h3><ul>
<li>Altitude
<ul>
<li>Equator: 16~18 km</li>
<li>Mid-latitudes: 10~12 km</li>
<li>Polar regions: 7~10 km</li>
</ul>
</li>
<li>Mass: Accounts for 75% of the total atmospheric mass</li>
<li>Characteristics: Features vertical air currents, storms, and thunderstorms. Temperature decreases rapidly with increasing altitude.</li>
</ul>
<h3 id="convective-boundary-layer">
<a class="header-anchor" href="#convective-boundary-layer"></a>
Convective Boundary Layer
</h3><p>Transition layer, with a thickness ranging from several hundred meters to one or two kilometers.</p>
<h3 id="stratosphere">
<a class="header-anchor" href="#stratosphere"></a>
Stratosphere
</h3><ul>
<li>Altitude: Troposphere ~32 km</li>
<li>Mass: Accounts for about one-quarter of the atmospheric mass</li>
<li>Characteristics: No weather phenomena, horizontal air flow, constant temperature (average about 216.65 K)</li>
</ul>
<h3 id="mesosphere">
<a class="header-anchor" href="#mesosphere"></a>
Mesosphere
</h3><ul>
<li>Altitude: 32~85 km</li>
<li>Mass: 1/3000</li>
<li>Temperature: First rises then falls, can drop below 106 K at 85 km.</li>
</ul>
<h3 id="upper-atmosphere">
<a class="header-anchor" href="#upper-atmosphere"></a>
Upper Atmosphere
</h3><ul>
<li>Altitude: Above 85 km</li>
<li>Characteristics: Non-uniform composition, directly absorbs solar radiation</li>
</ul>
<h3 id="thermosphere">
<a class="header-anchor" href="#thermosphere"></a>
Thermosphere
</h3><ul>
<li>Altitude: 85~500 km</li>
<li>Temperature: Increases with altitude, reaching up to 1370 K during the day at 500 km.</li>
<li>Characteristics: Directly exposed to solar shortwave radiation</li>
</ul>
<h3 id="outer-atmosphere">
<a class="header-anchor" href="#outer-atmosphere"></a>
Outer Atmosphere
</h3><ul>
<li>Altitude: 500+ km, where the atmosphere gradually merges with interstellar space</li>
<li>Mass: $1/10^{11}$</li>
<li>Characteristics: The atmosphere is too thin to be defined by temperature. Air molecules can escape into space.</li>
</ul>
<h3 id="upper-atmosphere-and-ionosphere">
<a class="header-anchor" href="#upper-atmosphere-and-ionosphere"></a>
Upper Atmosphere and Ionosphere
</h3><ul>
<li>The upper atmosphere is dissociated into electrons and ions by solar shortwave radiation, forming the ionosphere.</li>
<li>Above 100 km altitude, the air becomes a good conductor.</li>
<li>Above 150 km, the air is too thin to transmit sound.</li>
</ul>
<h3 id="d-layer">
<a class="header-anchor" href="#d-layer"></a>
D Layer
</h3><ul>
<li>高度：60~80 km</li>
</ul>
<h3 id="e-layer">
<a class="header-anchor" href="#e-layer"></a>
E Layer
</h3><ul>
<li>高度：100~120 km</li>
</ul>
<h3 id="layer">
<a class="header-anchor" href="#layer"></a>
$F_1$ Layer
</h3><ul>
<li>Height: 180~220 km</li>
</ul>
<h3 id="layer-1">
<a class="header-anchor" href="#layer-1"></a>
$F_2$ Layer
</h3><ul>
<li>Altitude: 300~350 km</li>
</ul>
<h3 id="flow-field">
<a class="header-anchor" href="#flow-field"></a>
Flow Field
</h3><p>Flow field: The space filled with moving fluid</p>
<p>Flow parameters: Physical quantities used to characterize fluid motion, such as velocity, density, pressure, etc.</p>
<p>Fluid mechanics methods: Lagrangian method, Eulerian method</p>
<h3 id="lagrange-method">
<a class="header-anchor" href="#lagrange-method"></a>
Lagrange Method
</h3><p>Focuses on <strong>particles</strong> (motion)</p>
<ul>
<li>Studies the variation of motion parameters and trajectories of individual particles in the flow field over time.</li>
<li>Synthesizes the changes in motion parameters of all fluid particles to derive the overall flow field motion characteristics.</li>
</ul>
<h3 id="eulerian-method">
<a class="header-anchor" href="#eulerian-method"></a>
Eulerian Method
</h3><p>Focuses on <strong>spatial points</strong> (fixed)</p>
<ul>
<li>Studies how the motion parameters of fluid particles change over time as they pass through fixed spatial points.</li>
<li>By combining the variation of motion parameters at all spatial points in the flow field, the overall flow field motion can be understood.</li>
</ul>
<p>In the Eulerian method, the motion parameters of the flow field are generally functions of spatial coordinates and time.</p>
<p>Taking velocity as an example:</p>
$$
v = v(x, y, z, t)  
$$<p>The four variables are independent.</p>
<p>In a general three-dimensional space, a Cartesian coordinate system is established, and scalar parameters are decomposed into the $x$, $y$, and $z$ directions for separate analysis.</p>
$$
v_x = v_x(x, y, z, t)  
$$<p>Differentiating to obtain the acceleration components:</p>
$$
a_x = \frac{d v_x}{dt} = \frac{\partial v_x}{\partial t} + \frac{\partial v_x}{\partial x}\frac{dx}{dt} + \frac{\partial v_x}{\partial y}\frac{dy}{dt} + \frac{\partial v_x}{\partial z}\frac{dz}{dt}  
$$<p>$\because$ $\frac{dx}{dt} = v_x$, $\frac{dy}{dt} = v_y$, $\frac{dz}{dt} = v_z$</p>
<p>$\therefore$</p>
$$
a_x = \frac{\partial v_x}{\partial t} + v_x \frac{\partial v_x}{\partial x} + v_y \frac{\partial v_x}{\partial y} + v_z \frac{\partial v_x}{\partial z}  
$$<p>From this, it can be seen that acceleration is a function of time and position:</p>
$$
a_x = a_x(t, x, y, z)  
$$<h3 id="local-acceleration">
<a class="header-anchor" href="#local-acceleration"></a>
Local Acceleration
</h3><p>The first term on the right side of the equation represents the rate of change of velocity of a fluid particle at a fixed spatial point over time, known as <strong>local acceleration</strong>. (Relationship between velocity and time)</p>
<p>Local acceleration is caused by the temporal variability of velocity in the flow field.</p>
<h3 id="convective-acceleration">
<a class="header-anchor" href="#convective-acceleration"></a>
Convective Acceleration
</h3><p>The latter three terms reflect the rate of change in velocity as a fluid particle moves along the velocity vector direction from one point in space to an adjacent point at the same instant, known as <strong>convective acceleration</strong>. (Relationship between velocity and displacement)</p>
<p>Convective acceleration is caused by the non-uniformity of the flow field.</p>
<h3 id="unsteady-flow-field">
<a class="header-anchor" href="#unsteady-flow-field"></a>
Unsteady Flow Field
</h3><p>In the flow field, the physical quantities at least at one <strong>spatial point</strong> vary with time.</p>
<h3 id="steady-flow-field">
<a class="header-anchor" href="#steady-flow-field"></a>
Steady Flow Field
</h3><p>In a <strong>steady flow field</strong>, the physical quantities at <strong>any spatial point</strong> do not change with time.</p>
<h1 id="">
<a class="header-anchor" href="#"></a>

</h1><h3 id="pathline">
<a class="header-anchor" href="#pathline"></a>
Pathline
</h3><p>The collection of all spatial points traversed by a marked moving fluid particle in a flow field over a period of time is called the pathline of that fluid particle.</p>
<h3 id="streamlines">
<a class="header-anchor" href="#streamlines"></a>
Streamlines
</h3><p>In a flow field, curves that are <strong>tangent to the velocity vector</strong> at every spatial point are called streamlines.</p>
<p>Streamlines are curves formed by different fluid particles at the same instant, indicating the direction of velocity for those particles at that moment.</p>
<p><strong>Characteristics:</strong></p>
<ul>
<li>
<p>In a steady flow field, the streamlines do not change over time.</p>
</li>
<li>
<p>Unsteady flow fields have streamlines that vary with time.</p>
</li>
<li>
<p>In a steady flow field, the streamline passing through a given spatial point coincides with the <strong>trajectories of all fluid particles</strong> passing through that point.</p>
</li>
<li>
<p>Generally, streamlines do not intersect (at the same instant and spatial point, there cannot be two velocity directions).</p>
</li>
<li>
<p>At points where the velocity is zero, streamlines can intersect. Such points are typically called <strong>stagnation points</strong>.</p>
</li>
<li>
<p>At points where the velocity is infinite, streamlines can intersect. Such points are typically called <strong>singularities</strong>.</p>
</li>
<li>
<p>Streamlines can be tangent to each other, and beyond the point of tangency, the two lines coincide.</p>
</li>
<li>
<p>Every point in the flow field has a streamline passing through it. The collection of all streamlines is called the <strong>streamline pattern</strong> or simply the <strong>flow pattern</strong>.</p>
</li>
</ul>
<h3 id="streamline-differential-equation">
<a class="header-anchor" href="#streamline-differential-equation"></a>
Streamline Differential Equation
</h3><p>Let the velocity at a point $M(x,y,z)$ on a streamline be $\vec{v}$, and the infinitesimal segment length of the streamline at point $M$ be $ds$, decomposed into $v_x,v_y,v_z$ and $dx,dy,dz$ in the Cartesian coordinate system.</p>
<p>At any point on the streamline, the direction of velocity is the same as the tangent direction of the streamline, so</p>
$$
\cos(\vec{v},\vec{i})=\frac{v_x}{v}=\frac{dx}{ds}
$$<p>Here, $\vec{i}$ is the unit normal vector in the $x$-axis direction, and similarly for the $y$ and $z$ axes.</p>
$$
\frac{dx}{v_x}=\frac{dy}{v_y}=\frac{dz}{v_z}
$$<p>The above equation is the differential equation of the streamline.</p>
<p>When the velocity distribution is known, the shape of the streamline passing through any point in the flow field can be determined.</p>
<h3 id="stream-tube">
<a class="header-anchor" href="#stream-tube"></a>
Stream Tube
</h3><p>In a flow field, a <strong>closed</strong> curve C that is <strong>not a streamline</strong> is considered. Drawing streamlines through every point on C, the collection of these streamlines forms a tubular surface known as a stream tube.</p>
<h3 id="analysis-of-fluid-particle-motion">
<a class="header-anchor" href="#analysis-of-fluid-particle-motion"></a>
Analysis of Fluid Particle Motion
</h3><h3 id="motion-forms">
<a class="header-anchor" href="#motion-forms"></a>
Motion Forms
</h3><h3 id="rigid-body-motion">
<a class="header-anchor" href="#rigid-body-motion"></a>
Rigid Body Motion
</h3><ul>
<li>Translational Motion</li>
<li>Rotational Motion About an Axis</li>
</ul>
<h3 id="fluid-motion">
<a class="header-anchor" href="#fluid-motion"></a>
Fluid Motion
</h3><ul>
<li>Translational motion</li>
<li>Rotational motion about an axis</li>
<li>Deformation motion
<ul>
<li>Linear deformation</li>
<li>Shear deformation</li>
</ul>
</li>
</ul>
<h3 id="two-dimensional-analysis">
<a class="header-anchor" href="#two-dimensional-analysis"></a>
Two-Dimensional Analysis
</h3><p>![Pasted image 20240902212258.png][1]</p>
<p>Consider an arbitrary rectangular fluid element ABCD in the flow field, with side lengths of $\delta_x,\delta_y$ respectively, both being infinitesimal quantities.</p>
<p>Let $v_x,v_y$ be the component velocities of the fluid element at point A, and assume the component velocities are continuous functions of spatial coordinates. The velocities at points B and D can be expressed using Taylor series expansions around point A.</p>
<p>$\because$ The side lengths of the fluid element are sufficiently small</p>
<p>$\therefore$ Higher-order infinitesimals can be neglected</p>
$$
v_{Bx}=v_x+\frac{\partial v_x}{\partial x}\delta_x
$$$$
v_{By}=v_y+\frac{\partial v_y}{x}\delta_x
$$<p>When the fluid element moves, in addition to the overall motion, point B also moves relative to point A.</p>
<p>The relative velocity in the $x$-direction is $v_{Bx}-v_x=\frac{\partial v_x}{\partial x}\delta_x$, and similarly for the $y$-direction $\frac{\partial v_x}{\partial x}\delta_x$.</p>
<p>The relative velocity of D with respect to A is $v_{Dx}-v_x=\frac{\partial v_x}{\partial y}\delta_y,v_{Dy}-v_y=\frac{\partial v_y}{\partial y}\delta_y$.</p>
<h2 id="">
<a class="header-anchor" href="#"></a>

</h2><h3 id="linear-deformation-motion">
<a class="header-anchor" href="#linear-deformation-motion"></a>
Linear Deformation Motion
</h3><p>![Pasted image 20240902212840.png][2]</p>
<p>The relative velocities $\frac{\partial v_x}{\partial x}\delta_x$ and $\frac{\partial v_y}{\partial y}\delta_y$ represent the linear deformation rates of rectangle ABCD&rsquo;s edges. During time interval $dt$:</p>
$$
AB'=AB+\frac{\partial v_x}{\partial x}\delta_x dt
$$$$
AD'=AD+\frac{\partial v_y}{\partial y}\delta_y dt
$$<p>The relative rate of area change is:</p>
$$
\frac{d(\delta S)}{\delta S \cdot dt}=\frac{AB' \cdot CD'-AB \cdot CD}{AB \cdot CD \cdot dt}
$$<p>After neglecting higher-order terms:</p>
$$
\frac{d(\delta S)}{\delta S \cdot dt}=\frac{\partial v_x}{\partial x}+\frac{\partial v_y}{\partial y}
$$<p>Extending to three-dimensional space similarly yields:</p>
$$
\frac{d(\delta V)}{\delta V \cdot dt}=\frac{\partial v_x}{\partial x}+\frac{\partial v_y}{\partial y}+\frac{\partial v_z}{\partial z}
$$<h3 id="angular-deformation-motion">
<a class="header-anchor" href="#angular-deformation-motion"></a>
Angular Deformation Motion
</h3><p>![Pasted image 20240902214219.png][3]</p>
<p>The relative velocities $\frac{\partial v_y}{\partial x}\delta x,\frac{\partial v_x}{\partial y}\delta y$ represent the rotation of edges AB and AD around point A.</p>
<p>Defining counterclockwise rotation as positive,</p>
<p>The angular velocity of edge AB is</p>
$$
\frac{d\alpha_1}{dt}=\frac{\partial v_y}{\partial x}\delta_x / \delta_x=\frac{\partial v_y}{\partial x}
$$<p>Similarly, the angular velocity of edge AD is</p>
$$
\frac{d \alpha_2}{dt}=-\frac{\partial v_x}{\partial y}
$$<h3 id="">
<a class="header-anchor" href="#"></a>

</h3><h3 id="microelement-rotation-angular-velocity-about-the--axis">
<a class="header-anchor" href="#microelement-rotation-angular-velocity-about-the--axis"></a>
Microelement Rotation Angular Velocity about the $z$-axis
</h3><p><strong>Definition</strong>: The <strong>average value</strong> of the angular velocities of two <strong>mutually perpendicular lines</strong> in the projection of a microelement on the $xOy$ plane rotating about the $z$-axis. (Half the sum of angular velocities)</p>
$$
\epsilon_z=\frac{1}{2}(\frac{\partial v_y}{\partial x}-\frac{\partial v_x}{\partial y})
$$<h3 id="">
<a class="header-anchor" href="#"></a>

</h3><h3 id="angular-deformation-rate">
<a class="header-anchor" href="#angular-deformation-rate"></a>
Angular Deformation Rate
</h3><p><strong>Definition</strong>: Half of the <strong>change in angle per unit time</strong> between two mutually perpendicular lines in the projection of a fluid element on the $xOy$ plane. (Half of the angular velocity difference)</p>
$$
\gamma_z=\frac{1}{2}\left(\frac{\partial v_y}{\partial x}+\frac{\partial v_x}{\partial y}\right)
$$<p>When extended to three-dimensional space, the angular velocities and angular deformation rates of the fluid element along the three axes can be similarly derived.</p>
<p>Omitted.</p>
<h3 id="divergence">
<a class="header-anchor" href="#divergence"></a>
Divergence
</h3><p>Definition: The sum of the directional derivatives of each velocity component in its respective direction is called the three-degree of the velocity vector.</p>
$$
div \vec{v}=\frac{\partial v_x}{\partial x}+\frac{\partial v_y}{\partial y}+\frac{\partial v_z}{\partial z}
$$<p>Physical meaning: It quantifies the <strong>relative volume change rate</strong> of a fluid element during motion.</p>
<p>==Assumption==: The density of the fluid remains unchanged (the fluid motion is considered incompressible).</p>
<p>The volume flow rate emitted from a point is defined as</p>
$$
\lim_{\delta V \to 0}\frac{Volume outflow - Volume inflow}{\delta V \cdot dt}=\frac{\partial v_x}{\partial x}+\frac{\partial v_y}{\partial y}+\frac{\partial v_z}{\partial z}
$$<p>This equals the <strong>net volume outflow</strong> per <strong>unit time</strong> from a <strong>unit volume</strong> control volume at a certain point in space, which is also equal to the relative volume change rate of the fluid element during motion.</p>
<h3 id="curl">
<a class="header-anchor" href="#curl"></a>
Curl
</h3><p>Definition: Twice the angular velocity of rotation.</p>
$$
\vec{\omega}=curl \vec{v}=(\frac{\partial v_z}{\partial y}-\frac{\partial v_y}{\partial z})\vec{i}+(\frac{\partial v_x}{\partial z}-\frac{\partial v_z}{\partial x})\vec{j}+(\frac{\partial v_y}{\partial x}-\frac{\partial v_x}{\partial y})\vec{k}
$$<h1 id="">
<a class="header-anchor" href="#"></a>

</h1><h3 id="velocity-potential">
<a class="header-anchor" href="#velocity-potential"></a>
Velocity Potential
</h3><p>In fluid mechanics, fluid motion can be classified based on whether fluid particles have rotational motion:</p>
<ul>
<li>Rotational flow</li>
<li>Irrotational flow</li>
</ul>
<p>When the flow is considered irrotational, $\omega=0$,</p>
$$
\begin{cases}
\frac{\partial v_z}{\partial y}=\frac{\partial v_y}{\partial z} \\
\frac{\partial v_x}{\partial z}=\frac{\partial v_z}{\partial x} \\
\frac{\partial v_y}{\partial x}=\frac{\partial v_x}{\partial y} \\
\end{cases}
$$<p>The above system of equations is the necessary and sufficient condition for $v_xdx+v_ydy+v_zdz$ to constitute the total differential of some function $\phi(x,y,z)$. That is,</p>
$$
d \phi=v_xdx+v_ydy+v_zdz=\frac{\partial \phi}{\partial x}dx+\frac{\partial \phi}{\partial y}dy+\frac{\partial \phi}{\partial z}dz
$$<p>$\phi$ is called the velocity potential or velocity potential function.</p>
$$
\begin{cases}
v_x=\frac{\partial \phi}{\partial x}\\
v_y=\frac{\partial \phi}{\partial y}\\
v_z=\frac{\partial \phi}{\partial z}\\
\end{cases}
$$<p>When using <strong>cylindrical coordinates</strong>,</p>
$$
\phi=\phi(r,\theta,z)
$$$$
\begin{cases}
v_r=\frac{\partial \phi}{\partial r}\\
v_\theta=\frac{\partial \phi}{\partial \theta}\\
v_z=\frac{\partial \phi}{\partial z}\\
\end{cases}
$$<hr>
<h3 id="scalars">
<a class="header-anchor" href="#scalars"></a>
Scalars
</h3><h3 id="pressure">
<a class="header-anchor" href="#pressure"></a>
Pressure
</h3><h3 id="density">
<a class="header-anchor" href="#density"></a>
Density
</h3><h3 id="temperature">
<a class="header-anchor" href="#temperature"></a>
Temperature
</h3><h3 id="viscosity-coefficient">
<a class="header-anchor" href="#viscosity-coefficient"></a>
Viscosity Coefficient
</h3><h3 id="vector">
<a class="header-anchor" href="#vector"></a>
Vector
</h3><h3 id="flow-velocity">
<a class="header-anchor" href="#flow-velocity"></a>
Flow Velocity
</h3><h1 id="">
<a class="header-anchor" href="#"></a>

</h1><h3 id="shear-stress">
<a class="header-anchor" href="#shear-stress"></a>
Shear Stress
</h3><h3 id="ideal-gas-equation-of-state">
<a class="header-anchor" href="#ideal-gas-equation-of-state"></a>
Ideal Gas Equation of State
</h3>$$
pV=nRT
$$<p>$R=8.314J \cdot mol^{-1} \cdot K^{-1}$ is the universal molar gas constant.</p>
$$
p=\frac{n \cdot M}{V}\frac{R}{M} T
$$$$
n(\text{amount of substance}) \cdot M(\text{molar mass})=m(\text{mass})
$$$$
p=\rho R' T
$$$$
R'=\frac{R}{M}
$$<p>$R'$ is the ==specific gas constant==.</p>
<p>For ideal air $R'=287J/(kg \cdot K)$</p>
<h3 id="aerodynamics-and-moments">
<a class="header-anchor" href="#aerodynamics-and-moments"></a>
Aerodynamics and Moments
</h3><h3 id="aerodynamic-force--resultant">
<a class="header-anchor" href="#aerodynamic-force--resultant"></a>
Aerodynamic Force $R$: Resultant
</h3><p><strong>The force exerted by air on an object</strong></p>
<ul>
<li>Pressure $p$: Pressure</li>
<li>Shear stress $\tau$: Shear stress</li>
</ul>
<p>The resultant of pressure and shear stress is the force exerted by air on the object, known as the aerodynamic force.</p>
<p><strong>Wind-axis system</strong></p>
<ul>
<li>Lift $L$: Lift, the vertical component</li>
<li>Drag $D$: Drag, the horizontal component</li>
</ul>
<h3 id="freestream">
<a class="header-anchor" href="#freestream"></a>
Freestream
</h3>$$
V_{\infty}
$$<p>Freestream refers to the undisturbed incoming flow in front of the aircraft, i.e., the natural flow of air without interference from the aircraft or other objects.</p>
<p>The directions of lift and drag are determined by the freestream direction.</p>
<h3 id="angle-of-attack-aoa">
<a class="header-anchor" href="#angle-of-attack-aoa"></a>
Angle of Attack (AOA)
</h3>$$
\alpha
$$<p><strong>Angle of attack</strong> (abbreviated as AOA, commonly represented by the <a href="https://zh.wikipedia.org/wiki/%E5%B8%8C%E8%87%98%E5%AD%97%E6%AF%8D" title="Greek letter">Greek letter</a> α) is a term in <a href="https://zh.wikipedia.org/wiki/%E7%A9%BA%E6%B0%A3%E5%8B%95%E5%8A%9B%E5%AD%B8" title="aerodynamics">aerodynamics</a>. It refers to the angle between the <a href="https://zh.wikipedia.org/wiki/%E7%BF%BC%E5%BC%A6" title="chord line">chord line</a> of an <a href="https://zh.wikipedia.org/wiki/%E6%A9%9F%E7%BF%BC" title="airfoil">airfoil</a> and the direction of the <strong>free stream</strong> (or <strong>relative wind flow</strong>). For an aircraft, the angle of attack is defined as the angle between the aircraft&rsquo;s longitudinal axis and the direction of the relative wind. A <strong>positive</strong> angle of attack occurs when the airfoil is tilted <strong>upward</strong>, while a <strong>negative</strong> angle of attack occurs when it is tilted <strong>downward</strong>.</p>
<p><strong>Body Axis System</strong></p>
<ul>
<li>Normal force $N$: Perpendicular to the airfoil direction</li>
<li>Axial force $A$: Parallel to the airfoil direction</li>
</ul>
<h3 id="力矩-moment">
<a class="header-anchor" href="#%e5%8a%9b%e7%9f%a9-moment"></a>
力矩 $M$：Moment
</h3><p>The moment that makes the aircraft pitch up is positive, and the moment that makes the aircraft pitch down is negative.</p>
<h3 id="dynamic-pressure">
<a class="header-anchor" href="#dynamic-pressure"></a>
Dynamic Pressure $q$
</h3><p>The dynamic pressure generated by free stream $V_{\infty},\rho_{\infty}$</p>
$$
q_{\infty}=\frac{1}{2}\rho_{\infty}V_{\infty}^2
$$<p>Unit is $Pa$, same as pressure</p>
<h3 id="characteristic-geometric-dimension">
<a class="header-anchor" href="#characteristic-geometric-dimension"></a>
Characteristic Geometric Dimension $S$
</h3><p>For three-dimensional objects, it represents the area; for two-dimensional objects, it represents the perimeter.</p>
<h3 id="dimensionless-parameters">
<a class="header-anchor" href="#dimensionless-parameters"></a>
Dimensionless Parameters
</h3><p>Three-dimensional objects commonly use uppercase $C$, while two-dimensional objects commonly use lowercase $c$.</p>
<h2 id="">
<a class="header-anchor" href="#"></a>

</h2><h3 id="lift-coefficient">
<a class="header-anchor" href="#lift-coefficient"></a>
Lift Coefficient
</h3>$$
C_L=\frac{L}{q_{\infty}S}
$$<h2 id="">
<a class="header-anchor" href="#"></a>

</h2><h3 id="drag-coefficient">
<a class="header-anchor" href="#drag-coefficient"></a>
Drag Coefficient
</h3>$$
C_D=\frac{D}{q_{\infty}S}
$$<h2 id="">
<a class="header-anchor" href="#"></a>

</h2><h3 id="normal-force-coefficient">
<a class="header-anchor" href="#normal-force-coefficient"></a>
Normal Force Coefficient
</h3>$$
C_N=\frac{N}{q_{\infty}S}
$$<h2 id="">
<a class="header-anchor" href="#"></a>

</h2><h3 id="axial-force-coefficient">
<a class="header-anchor" href="#axial-force-coefficient"></a>
Axial Force Coefficient
</h3>$$
C_A=\frac{A}{q_{\infty}S}
$$<h2 id="">
<a class="header-anchor" href="#"></a>

</h2><h3 id="aerodynamic-coefficients">
<a class="header-anchor" href="#aerodynamic-coefficients"></a>
Aerodynamic Coefficients
</h3>$$
C_R=\frac{R}{q_{\infty}S}
$$<h2 id="">
<a class="header-anchor" href="#"></a>

</h2><h3 id="moment-coefficient">
<a class="header-anchor" href="#moment-coefficient"></a>
Moment Coefficient
</h3>$$
C_M=\frac{\vec{M}}{\vec{r} \times \vec{q_{\infty}}S}
$$<h2 id="">
<a class="header-anchor" href="#"></a>

</h2><h3 id="pressure-coefficient">
<a class="header-anchor" href="#pressure-coefficient"></a>
Pressure Coefficient
</h3><p>$p$：某点静压 (Static pressure at a point)</p>
<p>$p_{\infty}$：自由来流静压 (Freestream static pressure)</p>
$$
C_p=\frac{p-p_{\infty}}{q_{\infty}}
$$<h3 id="friction-coefficient">
<a class="header-anchor" href="#friction-coefficient"></a>
Friction Coefficient
</h3><p>$\tau$: Shear stress at a point, which is the derivative of shear stress with respect to area. Dimensions are the same as pressure.</p>
$$
C_f=\frac{\tau}{q_{\infty}}
$$<h2 id="two-centers">
<a class="header-anchor" href="#two-centers"></a>
Two Centers
</h2><h3 id="center-of-pressure-cop">
<a class="header-anchor" href="#center-of-pressure-cop"></a>
Center of Pressure (COP)
</h3><p>Center of Pressure (COP): The point of intersection between the line of action of the resultant force of fluid pressure on a plane or curved surface within the fluid and that plane or surface. The aerodynamic force $R$ produces zero moment $\vec{0}$ about this point.</p>
<h3 id="aerodynamic-center-ac">
<a class="header-anchor" href="#aerodynamic-center-ac"></a>
Aerodynamic Center (AC)
</h3><p><strong>Aerodynamic Center</strong> (English: aerodynamic center, abbreviated as AC) in <a href="https://zh.wikipedia.org/wiki/%E7%A9%BA%E6%B0%94%E5%8A%A8%E5%8A%9B%E5%AD%A6" title="Aerodynamics">aerodynamics</a> refers to a fixed point on an <a href="https://zh.wikipedia.org/wiki/%E7%BF%BC%E5%9E%8B" title="Airfoil">airfoil</a> where the pitching moment does not vary with changes in <a href="https://zh.wikipedia.org/wiki/%E8%BF%8E%E8%A7%92" title="Angle of Attack">angle of attack</a>, i.e.,</p>
$$
\frac{d C_M}{d \alpha}=0
$$<h3 id="the-difference-between-aerodynamic-center-and-center-of-pressure">
<a class="header-anchor" href="#the-difference-between-aerodynamic-center-and-center-of-pressure"></a>
The Difference Between Aerodynamic Center and Center of Pressure
</h3><p>The center of pressure is a special point where the force system is synthesized, resulting in <strong>zero net moment</strong> at that point, <strong>with the center of pressure located behind the aerodynamic center</strong>. On the other hand, the aerodynamic center is the point where <strong>the net moment remains unchanged</strong>.</p>
<p>The position of the center of pressure changes with variations in the angle of attack. As the angle of attack increases, the lift increases, causing the center of pressure to move forward. This simultaneously reduces the distance between the center of pressure and the aerodynamic center. The increased lift multiplied by the shortened moment arm precisely equals the unchanged moment, which is exactly what the definition of the aerodynamic center requires.</p>

        
        <hr><p>Published on 2024-09-02 at <a href='https://www.guzhengsvt.cn/'>孤筝の温暖小家</a>, last modified on 2024-09-02</p><p>All articles on this blog are licensed under the BY-NC-SA license agreement unless otherwise stated. Please indicate the source when reprinting!</p>]]></description><category>Physics</category></item><item><title>Thermodynamics and Statistical Physics</title><link>https://www.guzhengsvt.cn/en/post/physics/%E7%83%AD%E5%8A%9B%E5%AD%A6%E4%B8%8E%E7%BB%9F%E8%AE%A1%E7%89%A9%E7%90%86/</link><pubDate>Thu, 06 Jun 2024 00:00:59 +0800</pubDate><author>lvbowen040427@163.com (孤筝)</author><guid>https://www.guzhengsvt.cn/en/post/physics/%E7%83%AD%E5%8A%9B%E5%AD%A6%E4%B8%8E%E7%BB%9F%E8%AE%A1%E7%89%A9%E7%90%86/</guid><description>
<![CDATA[<h1>Thermodynamics and Statistical Physics</h1><p>Author: 孤筝(lvbowen040427@163.com)</p>
        
          <h2 id="physical-quantities">
<a class="header-anchor" href="#physical-quantities"></a>
Physical Quantities
</h2><h3 id="internal-energy-u">
<a class="header-anchor" href="#internal-energy-u"></a>
Internal Energy U
</h3>$$
\Delta U=Q+W
$$<p>
Q is the heat absorbed by the system, W is the work done by the system</p>
<h3 id="entropy-s">
<a class="header-anchor" href="#entropy-s"></a>
Entropy S
</h3><p>Entropy is a measure of the total amount of energy that <strong>cannot perform work</strong> in terms of dynamics. That is, when the total entropy increases, the ability to perform work decreases, and the measure of entropy is an indicator of energy degradation.<br>
In a ==reversible process==<br>
</p>
$$
\Delta S=\frac {Q}{T}
$$<p><br>
Q: In a reversible process, the system absorbs heat at a constant temperature.</p>
<h3 id="焓-h">
<a class="header-anchor" href="#%e7%84%93-h"></a>
焓 H
</h3>$$
H=U+pV
$$<h3 id="free-energy-f">
<a class="header-anchor" href="#free-energy-f"></a>
Free Energy F
</h3>$$
F=U-TS
$$<p>
Internal energy of the system - energy that cannot do work = energy that can do work (free energy)</p>
<h3 id="gibbs-function-free-enthalpy-g">
<a class="header-anchor" href="#gibbs-function-free-enthalpy-g"></a>
Gibbs Function (Free Enthalpy) G
</h3>$$
G=F+pV=U-TS+pV=H-TS
$$<h3 id="chemical-potential">
<a class="header-anchor" href="#chemical-potential"></a>
Chemical Potential $\mu$
</h3>$$
\mu=(\frac{\partial G}{\partial n})_{T,p}=G_m
$$<p>
Where:</p>
<ul>
<li>n: Amount of substance ($mol$)</li>
<li>$G_m$: Molar Gibbs function</li>
</ul>
<h2 id="equation">
<a class="header-anchor" href="#equation"></a>
Equation
</h2><h3 id="state-function-total-differential-reversible-process">
<a class="header-anchor" href="#state-function-total-differential-reversible-process"></a>
State Function Total Differential (Reversible Process)
</h3><h1 id="">
<a class="header-anchor" href="#"></a>

</h1><h3 id="total-differential-of-internal-energy-fundamental-equation-of-thermodynamics">
<a class="header-anchor" href="#total-differential-of-internal-energy-fundamental-equation-of-thermodynamics"></a>
Total Differential of Internal Energy (Fundamental Equation of Thermodynamics)
</h3>$$
dU=dQ+dW=TdS-pdV
$$<p>
i.e., $U=U(S,V)$
</p>
$$
(\frac{\partial U}{\partial S})_V=T,(\frac{\partial U}{\partial V})_S=-p
$$<p>
$\because$ The second-order partial derivatives of U are equal<br>
$\therefore$ We obtain the <strong>first Maxwell relation</strong>:
</p>
$$
(\frac{\partial T}{\partial V})_S=-(\frac{\partial p}{\partial S})_V
$$<h3 id="enthalpy-total-differential">
<a class="header-anchor" href="#enthalpy-total-differential"></a>
Enthalpy Total Differential
</h3>$$
dH=dU+d(pV)=TdS-pdV+pdV+Vdp=TdS+Vdp
$$<h3 id="differential-of-free-energy">
<a class="header-anchor" href="#differential-of-free-energy"></a>
Differential of Free Energy
</h3>$$
dF=-SdT-pdV
$$<h3 id="gibbs-function-total-differential">
<a class="header-anchor" href="#gibbs-function-total-differential"></a>
Gibbs Function Total Differential
</h3>$$
dG=-SdT+Vdp
$$<h3 id="total-differential-of-entropy">
<a class="header-anchor" href="#total-differential-of-entropy"></a>
Total Differential of Entropy
</h3>$$
dS=\frac{dU+pdV}{T}
$$<h3 id="maxwell-relations">
<a class="header-anchor" href="#maxwell-relations"></a>
Maxwell Relations
</h3><p>The Maxwell Relations are important equations in thermodynamics that describe the relationships between state variables.<br>
Here are the four Maxwell Relations:</p>
<ol>
<li>
<p><strong>First Maxwell Relation</strong>:
</p>
$$
\left ( \frac{\partial T}{\partial V} \right)_S = - \left ( \frac{\partial P}{\partial S} \right)_V
$$</li>
<li>
<p><strong>Second Maxwell Relation</strong>:
</p>
$$
\left ( \frac{\partial T}{\partial P} \right)_S = \left ( \frac{\partial V}{\partial S} \right)_P
$$</li>
<li>
<p><strong>Third Maxwell Relation</strong>:
</p>
$$
\left ( \frac{\partial S}{\partial V} \right)_T = \left ( \frac{\partial P}{\partial T} \right)_V
$$</li>
<li>
<p><strong>Fourth Maxwell Relation</strong>:
</p>
$$
\left ( \frac{\partial S}{\partial P} \right)_T = -\left ( \frac{\partial V}{\partial T} \right)_P
$$</li>
</ol>
<p><strong>Corollary</strong>:
</p>
$$
C_p-C_V=T(\frac{\partial p}{\partial T})_V(\frac{\partial V}{\partial T})_p
$$<p>For an ideal gas $pV=nRT$, substituting yields:
</p>
$$
C_p-C_V=nR, \text{which is a constant.}
$$<h3 id="equation-of-state">
<a class="header-anchor" href="#equation-of-state"></a>
Equation of State
</h3><p>Definition: The functional relationship equation between temperature and state parameters.
</p>
$$
f(p,V,T)=0
$$<h3 id="ideal-gas-equation-of-state">
<a class="header-anchor" href="#ideal-gas-equation-of-state"></a>
Ideal Gas Equation of State
</h3>$$
pV=nRT=NkT
$$<ul>
<li>n: Amount of substance</li>
<li>R: Gas constant</li>
<li>N: Total number of particles in the system</li>
<li>k: Boltzmann constant</li>
</ul>
<h3 id="van-der-waals-equation">
<a class="header-anchor" href="#van-der-waals-equation"></a>
Van der Waals Equation
</h3>$$
(p+\frac{an^2}{V^2})(V-nb)=nRT
$$<ul>
<li>a、b are constants</li>
</ul>
<h3 id="common-coefficients">
<a class="header-anchor" href="#common-coefficients"></a>
Common Coefficients
</h3><h3 id="coefficient-of-volume-expansion">
<a class="header-anchor" href="#coefficient-of-volume-expansion"></a>
Coefficient of Volume Expansion $\alpha$
</h3>$$
\alpha=\frac{1}{V}(\frac{\partial V}{\partial T})_p=\frac{1}{T}, \text{ideal gas}
$$<h1 id="">
<a class="header-anchor" href="#"></a>

</h1><h3 id="pressure-coefficient">
<a class="header-anchor" href="#pressure-coefficient"></a>
Pressure Coefficient $\beta$
</h3>$$
\beta=\frac{1}{V}(\frac{\partial p}{\partial T})_V=\frac{1}{T}, \text{ideal gas}
$$<h3 id="isothermal-compressibility">
<a class="header-anchor" href="#isothermal-compressibility"></a>
Isothermal Compressibility $\kappa_T$
</h3>$$
\kappa_T=-\frac{1}{V}(\frac{\partial V}{\partial P})_T=\frac{1}{p}, \text{ideal gas}
$$<h3 id="heat-capacity-c">
<a class="header-anchor" href="#heat-capacity-c"></a>
Heat Capacity C
</h3>$$
C=\lim_{\Delta T \to 0}\frac{\Delta Q}{\Delta T}
$$<h3 id="constant-volume-heat-capacity">
<a class="header-anchor" href="#constant-volume-heat-capacity"></a>
Constant Volume Heat Capacity $C_V$
</h3>$$
C_V=(\frac{\partial U}{\partial T})_V=T(\frac{\partial S}{\partial T})_V
$$<h3 id="heat-capacity-at-constant-pressure">
<a class="header-anchor" href="#heat-capacity-at-constant-pressure"></a>
Heat Capacity at Constant Pressure $C_p$
</h3>$$
C_p=(\frac{\partial H}{\partial T})_p=T(\frac{\partial S}{\partial T})_p
$$<h3 id="multipartite-heat-capacity">
<a class="header-anchor" href="#multipartite-heat-capacity"></a>
Multipartite Heat Capacity $C_n$
</h3>$$
C_n=T(\frac{\partial S}{\partial T})_n
$$<h3 id="adiabatic-equation-for-ideal-gas">
<a class="header-anchor" href="#adiabatic-equation-for-ideal-gas"></a>
Adiabatic Equation for Ideal Gas
</h3><p>$\gamma$：[[#Adiabatic Index]]
C: Constant
</p>
$$
pV^\gamma=C_1
$$<p>
</p>
$$
TV^{\gamma-1}=C_2
$$<p>
</p>
$$
P^{1-\gamma}T^\gamma=C_3
$$<h3 id="the-second-law-of-thermodynamics">
<a class="header-anchor" href="#the-second-law-of-thermodynamics"></a>
==The Second Law of Thermodynamics==
</h3>$$
dS\ge \frac{dQ}{T}
$$<ol>
<li>Two adiabatic lines cannot intersect.</li>
<li>T: External temperature</li>
<li>The external temperature equals the system temperature (T) if and only if the process is reversible.</li>
</ol>
<h3 id="efficiency-of-reversible-heat-engine">
<a class="header-anchor" href="#efficiency-of-reversible-heat-engine"></a>
Efficiency of Reversible Heat Engine $\eta$
</h3><p>[[#Carnot&rsquo;s Theorem]]
</p>
$$
\eta=1-\frac{T_1}{T_2}
$$<ul>
<li>$T_2$: Temperature of the high-temperature heat reservoir</li>
<li>$T_1$: Temperature of the low-temperature heat reservoir</li>
<li>$0 \lt \eta \lt 1$</li>
</ul>
<h3 id="energy-equation-in-throttling-process">
<a class="header-anchor" href="#energy-equation-in-throttling-process"></a>
Energy Equation in Throttling Process
</h3><p>For the throttling process, the following energy balance equation can be written:
</p>
$$
H_1 = H_2
$$<p>where, $H_1$ and $H_2$ are the enthalpies before and after throttling, respectively.
<strong>Proof</strong>:
The system is adiabatic, $\Delta U=W$
$U_2-U_1=p_1V_1-p_2V_2$
$U_2+p_2V_2=U_1+p_1V_1$
$H_2=H_1$</p>
<h3 id="joule-thomson-coefficient">
<a class="header-anchor" href="#joule-thomson-coefficient"></a>
Joule-Thomson Coefficient
</h3><p>[[#Throttling Process]]
</p>
$$
\mu=(\frac{\partial T}{\partial p})_H=\frac{V}{C_p}(T\alpha-1)
$$<h3 id="fundamental-equations-of-thermodynamics-for-open-systems-">
<a class="header-anchor" href="#fundamental-equations-of-thermodynamics-for-open-systems-"></a>
Fundamental Equations of Thermodynamics for Open Systems ($+\mu dn$)
</h3><h1 id="">
<a class="header-anchor" href="#"></a>

</h1><h3 id="total-differential-of-internal-energy">
<a class="header-anchor" href="#total-differential-of-internal-energy"></a>
Total Differential of Internal Energy
</h3>$$
dU=dQ+dW+\mu dn=TdS-pdV+\mu dn
$$<h3 id="enthalpy-total-differential-1">
<a class="header-anchor" href="#enthalpy-total-differential-1"></a>
Enthalpy Total Differential
</h3>$$
dH=dU+d(pV)+\mu dn=TdS+Vdp+\mu dn
$$<h3 id="differential-of-free-energy-1">
<a class="header-anchor" href="#differential-of-free-energy-1"></a>
Differential of Free Energy
</h3>$$
dF=-SdT-pdV+\mu dn
$$<h3 id="gibbs-function-total-differential-1">
<a class="header-anchor" href="#gibbs-function-total-differential-1"></a>
Gibbs Function Total Differential
</h3>$$
dG=-SdT+Vdp+\mu dn
$$<h3 id="grand-potential">
<a class="header-anchor" href="#grand-potential"></a>
Grand Potential
</h3>$$
J=-pV
$$<h3 id="the-clausius-clapeyron-equation">
<a class="header-anchor" href="#the-clausius-clapeyron-equation"></a>
The Clausius-Clapeyron Equation
</h3>$$
\frac{dp}{dT}=\frac{S_m^\beta-S_m^\alpha}{V_m^\beta-V_m^\alpha}=\frac{L}{T(V_m^\beta-V_m^\alpha)}=\frac{s^\beta-s^\alpha}{v^\beta-v^\alpha}
$$<p>
Where:</p>
<ul>
<li>$L=T(S_m^\beta-S_m^\alpha)$ is the latent heat of phase transition.</li>
<li>$S_m$ is the molar entropy</li>
<li>$s$ is the specific entropy (entropy per unit mass)</li>
<li>$v$ is the specific volume (volume per unit mass)</li>
</ul>
<h3 id="ehrenfest-equations">
<a class="header-anchor" href="#ehrenfest-equations"></a>
Ehrenfest Equations
</h3><p>In [[#second-order phase transitions]], the Clapeyron equation becomes an indeterminate form of $\frac{0}{0}$. By applying L&rsquo;Hôpital&rsquo;s rule and taking partial derivatives of the numerator and denominator with respect to $T$, we obtain
</p>
$$
\frac{d𝑝}{d𝑇}=\frac{\Delta c_𝑝}{𝑇𝑣\Delta\alpha}
$$<p>
Alternatively, taking partial derivatives with respect to $𝑝$ yields
</p>
$$
\frac{d𝑝}{d𝑇}=\frac{\Delta\alpha}{\Delta\kappa_T}
$$<p>
These two equations are known as the <strong>Ehrenfest equations</strong>.</p>
<h3 id="first-order-phase-transition-equation">
<a class="header-anchor" href="#first-order-phase-transition-equation"></a>
First-Order Phase Transition Equation
</h3><h3 id="formula">
<a class="header-anchor" href="#formula"></a>
Formula
</h3><p>For first-order phase transitions, at the transition temperature $T_c$, the Gibbs free energy $G$ itself is continuous, but its first derivatives with respect to temperature and pressure exhibit discontinuities:</p>
<ul>
<li>$\left ( \frac{\partial G}{\partial T} \right)_P = -S$</li>
<li>$\left ( \frac{\partial G}{\partial P} \right)_T = V$</li>
</ul>
<p>At the transition point $T_c$, these first derivatives undergo jumps:
</p>
$$
\Delta S = S_2 - S_1 \neq 0
$$<p>
</p>
$$
\Delta V = V_2 - V_1 \neq 0
$$<p>
Here, $S$ represents entropy, $V$ represents volume, and the subscripts 1 and 2 denote the two phases before and after the phase transition.</p>
<h3 id="proof">
<a class="header-anchor" href="#proof"></a>
Proof
</h3><p>The characteristic of a first-order phase transition lies in the <strong>discontinuity</strong> of <strong>certain first-order derivatives of the thermodynamic potential function</strong> at the transition point. Taking entropy as an example:</p>
<ol>
<li><strong>Continuity of Gibbs Free Energy</strong>:
At the transition temperature $T_c$, the Gibbs free energy $G$ is continuous:

$$
   G_1 (T_c, P) = G_2 (T_c, P)
   $$</li>
<li><strong>Discontinuity of Entropy</strong>:
Taking the partial derivative of $G$ with respect to temperature yields the expression for entropy:

$$
   S = -\left ( \frac{\partial G}{\partial T} \right)_P
   $$
At $T_c$, the jump in entropy is:

$$
   \Delta S = S_2 - S_1 = -\left ( \frac{\partial G_2}{\partial T} \right)_P + \left ( \frac{\partial G_1}{\partial T} \right)_P \neq 0
   $$</li>
</ol>
<p>This demonstrates the discontinuity of entropy at the phase transition point.</p>
<h3 id="second-order-phase-transition-equations">
<a class="header-anchor" href="#second-order-phase-transition-equations"></a>
Second-Order Phase Transition Equations
</h3><h2 id="">
<a class="header-anchor" href="#"></a>

</h2><h3 id="formula-1">
<a class="header-anchor" href="#formula-1"></a>
Formula
</h3><p>For second-order phase transitions, at the critical temperature $T_c$, the Gibbs free energy $G$ and its first derivatives (such as entropy and volume) are continuous, but its second derivatives (such as heat capacity, compressibility, and thermal expansion coefficient) are discontinuous:</p>
<ul>
<li>$\left ( \frac{\partial^2 G}{\partial T^2} \right)_P = \frac{\partial S}{\partial T} = \frac{C_P}{T}$</li>
<li>$\left ( \frac{\partial^2 G}{\partial P^2} \right)_T = \frac{\partial V}{\partial P} = -\kappa_T V$</li>
</ul>
<p>At the phase transition point $T_c$, these second derivatives exhibit jumps:
</p>
$$
\Delta C_P = C_{P 2} - C_{P 1} \neq 0
$$<p>
</p>
$$
\Delta \kappa_T = \kappa_{T 2} - \kappa_{T 1} \neq 0
$$<p>
Here, $C_P$ is the constant-pressure heat capacity, and $\kappa_T$ is the isothermal compressibility.</p>
<h3 id="proof-1">
<a class="header-anchor" href="#proof-1"></a>
Proof
</h3><p>The characteristics of second-order phase transitions lie in <strong>the discontinuity of certain second derivatives of the thermodynamic potential function at the phase transition point</strong>. Taking heat capacity as an example:</p>
<ol>
<li><strong>Continuity of Entropy</strong>:
At the phase transition temperature $T_c$, the entropy $S$ is continuous:

$$
   S_1 (T_c, P) = S_2 (T_c, P)
   $$</li>
<li><strong>Discontinuity of Heat Capacity</strong>:
Taking the partial derivative of entropy with respect to temperature yields the expression for heat capacity:

$$
   C_P = T \left ( \frac{\partial S}{\partial T} \right)_P
   $$
At $T_c$, the jump in heat capacity is:

$$
   \Delta C_P = C_{P 2} - C_{P 1} = T \left ( \left ( \frac{\partial S_2}{\partial T} \right)_P - \left ( \frac{\partial S_1}{\partial T} \right)_P \right) \neq 0
   $$</li>
</ol>
<p>This indicates that heat capacity exhibits discontinuity at the phase transition point.</p>
<h3 id="differential-of-the-gibbs-function-for-multicomponent-systems">
<a class="header-anchor" href="#differential-of-the-gibbs-function-for-multicomponent-systems"></a>
Differential of the Gibbs Function for Multicomponent Systems
</h3>$$
dG=-SdT+Vdp+\sum_i\mu_i dn_i
$$<h3 id="gibbs-phase-rule">
<a class="header-anchor" href="#gibbs-phase-rule"></a>
Gibbs Phase Rule
</h3>$$
f=k+2-\varphi
$$<p>
Where:</p>
<ul>
<li>$f$: Degrees of freedom in a multicomponent multiphase system, i.e., the number of independently variable intensive properties</li>
<li>$k$: Number of components</li>
<li>$\varphi$: Number of phases</li>
</ul>
<h3 id="de-broglie-relations">
<a class="header-anchor" href="#de-broglie-relations"></a>
de Broglie Relations
</h3><p>[[#de Broglie Wave]]
</p>
$$
\begin{cases}
\varepsilon=\hbar\omega\\
\vec{p}=\hbar\vec{k}
\end{cases}
$$<h3 id="planck-constant">
<a class="header-anchor" href="#planck-constant"></a>
Planck Constant
</h3>$$
\begin{align}
h & =6.626069934(89)×10^{-34}J\cdot s \\
& =4.135667662(25)×10^{-15}eV\cdot s
\end{align}
$$<p>
Reduced Planck Constant (Dirac Constant):
</p>
$$
\hbar \equiv \frac{h}{2\pi}=1.054571800(13)\times10^{-34}J\cdot s
$$<h3 id="uncertainty-relation">
<a class="header-anchor" href="#uncertainty-relation"></a>
Uncertainty Relation
</h3>$$
\Delta q\Delta p\ge h
$$<p>
<strong>Or</strong>
</p>
$$
\Delta q\Delta p\ge \frac{\hbar}{2}
$$<p>
==Generally==, take
</p>
$$
\Delta q\Delta p\approx h
$$<h3 id="three-dimensional-free-particle-quantum-state-count">
<a class="header-anchor" href="#three-dimensional-free-particle-quantum-state-count"></a>
Three-Dimensional Free Particle Quantum State Count
</h3><p>==Temporarily disregarding quantum spin==</p>
<blockquote>
<p>When particles move within a macroscopically sized container, their momentum and energy values are <strong>quasi-continuous</strong>.</p>
</blockquote>
<h3 id="cartesian-momentum-representation">
<a class="header-anchor" href="#cartesian-momentum-representation"></a>
Cartesian Momentum Representation
</h3><p>Find the <strong>number of quantum states for free particles</strong> within a volume $V=L^3$, in the momentum range from $p_x$ to $p_x+dp_x$, $p_y$ to $p_y+dp_y$, and $p_z$ to $p_z+dp_z$.<br>
Solution 1:<br>
From the relationship between momentum and quantum numbers for [[#Three-Dimensional Free Particles]], the possible number of $p_x$ states is:<br>
</p>
$$
dn_x=\frac{L}{2\pi\hbar}dp_x
$$<p><br>
This means $dn_x$ is characterized by $dp_x$.<br>
The same applies to the other two directions.<br>
</p>
$$
dn_xdn_ydn_z=\frac{V}{h^3}dp_xdp_ydp_z
$$<p><br>
Solution 2:<br>
From the definition of [[#Phase Cell]], the phase cell size for a three-dimensional free particle is $\Delta q_1\cdots \Delta q_r\Delta p_1\cdots \Delta p_r\approx h^3$. In this problem, the $\mu$-space volume is $d\Omega=Vdp_xdp_ydp_z$.<br>
The number of phase cells that can fit into this volume is the number of quantum states for free particles.<br>
</p>
$$
dn_xdn_ydn_z=\frac{d\Omega}{h^3}=\frac{Vdp_xdp_ydp_z}{h^3}
$$<h3 id="momentum-representation-in-spherical-coordinates">
<a class="header-anchor" href="#momentum-representation-in-spherical-coordinates"></a>
Momentum Representation in Spherical Coordinates
</h3>$$
\displaylines{p_x=p\sin \theta \cos \varphi 
\\p_y=p\sin \theta \sin \varphi 
\\p_z=p\cos \theta}
$$<p>
</p>
$$
dp_xdp_ydp_z=p^2\sin \theta dp d\theta d\varphi
$$<p>
Within volume $V$, for momentum magnitude ranging from $p$ to $p+dp$, direction from $\theta$ to $d\theta$, and $\varphi$ to $\varphi+d\varphi$, the number of possible states for a free particle is
</p>
$$
\frac{Vp^2\sin \theta dp d\theta d\varphi}{h^3}
$$<p>
<em>There are three quantum numbers $p,\theta,\varphi$ in total, corresponding to three degrees of freedom.</em></p>
<h3 id="energy-representation">
<a class="header-anchor" href="#energy-representation"></a>
Energy Representation
</h3><p>Integrating over $\theta,\varphi$ ($0 \lt \theta \lt \pi,0 \lt \varphi \lt 2\pi$) gives:
</p>
$$
\frac{4\pi Vp^2}{h^3}dp
$$<p>
Here the momentum can be in any direction, with 1 degree of freedom.
$\because$ $\varepsilon=\frac{p^2}{2m}$
$\therefore$ Within volume $V$, the number of possible quantum states for free particles in the energy range from $\varepsilon$ to $\varepsilon+d\varepsilon$ is
</p>
$$
D(\varepsilon)d\varepsilon=\frac{2\pi V}{h^3}(2m)^{3/2}\varepsilon^{1/2}d\varepsilon
$$<p>
where $D(\varepsilon)$ represents the number of possible states per unit energy interval, known as the ==density of states==.</p>
<h3 id="considering-quantum-spin">
<a class="header-anchor" href="#considering-quantum-spin"></a>
Considering Quantum Spin
</h3><p>[[#Spin Angular Momentum]]
Number of quantum states = Number of momentum quantum states obtained above * Possible spin states of the particle
eg. If the particle&rsquo;s spin quantum number is $\frac{1}{2}$, $S_z=2s+1=1$, the above result should be multiplied by 2.</p>
<h2 id="definitions-and-theorems">
<a class="header-anchor" href="#definitions-and-theorems"></a>
Definitions and Theorems
</h2><h3 id="basic-definitions">
<a class="header-anchor" href="#basic-definitions"></a>
Basic Definitions
</h3><h1 id="">
<a class="header-anchor" href="#"></a>

</h1><h3 id="equilibrium-state">
<a class="header-anchor" href="#equilibrium-state"></a>
Equilibrium State
</h3><p>The definition of equilibrium state focuses on two points:</p>
<ol>
<li>There is <strong>no macroscopic</strong> exchange of energy or matter between the system and its surroundings (distinguished from <strong>steady state</strong>). Microscopic exchanges are allowed.</li>
<li>The macroscopic properties of all parts of the system remain unchanged over a long period of time. Equilibrium state is <strong>thermodynamic equilibrium</strong>, with <strong>fluctuations</strong> present.</li>
</ol>
<p><strong>Relaxation time</strong>: The time required for a system to return to equilibrium from a non-equilibrium state.</p>
<h3 id="methods-for-describing-equilibrium-properties-macroscopic-description">
<a class="header-anchor" href="#methods-for-describing-equilibrium-properties-macroscopic-description"></a>
Methods for Describing Equilibrium Properties: Macroscopic Description
</h3><ul>
<li>State parameters: Macroscopic variables that determine equilibrium properties</li>
<li>State functions: Macroscopic variables determined by state parameters</li>
</ul>
<h3 id="extensive-and-intensive-quantities">
<a class="header-anchor" href="#extensive-and-intensive-quantities"></a>
Extensive and Intensive Quantities
</h3><p>In thermodynamics, the <strong>macroscopic variables</strong> of a homogeneous system can be divided into two categories: extensive and intensive quantities.</p>
<h3 id="广延量extensive-variables">
<a class="header-anchor" href="#%e5%b9%bf%e5%bb%b6%e9%87%8fextensive-variables"></a>
广延量（Extensive Variables）
</h3><blockquote>
<p>广延量是指<strong>依赖于系统规模或大小的热力学变量</strong>。这些变量随着系统的大小成比例地增加或减少。</p>
</blockquote>
<ul>
<li><strong>质量（Mass）</strong>：整个系统的质量。</li>
<li><strong>体积（Volume, V）</strong>：系统占据的空间体积。</li>
<li><strong>能量（Energy, E）</strong>：包括内部能量、动能和势能。</li>
<li><strong>熵（Entropy, S）</strong>：表示系统微观状态的不确定性和混乱度。</li>
<li><strong>粒子数（Number of particles, N）</strong>：系统中粒子的总数。</li>
<li><strong>热量（Heat, Q）</strong>：系统所包含或传递的热量。</li>
<li><strong>电荷（Charge, Q）</strong>：系统中的电荷总量。</li>
</ul>
<h3 id="强度量intensive-variables">
<a class="header-anchor" href="#%e5%bc%ba%e5%ba%a6%e9%87%8fintensive-variables"></a>
强度量（Intensive Variables）
</h3><blockquote>
<p>强度量是指<strong>与系统的大小无关</strong>的热力学变量，这些变量<strong>在系统内的任一点都保持一致</strong>，不随系统规模而变化。</p>
</blockquote>
<ul>
<li><strong>温度（Temperature, T）</strong>：表示系统的热状态。</li>
<li><strong>压强（Pressure, P）</strong>：系统单位面积上所受的力。</li>
<li><strong>化学势（Chemical potential, μ）</strong>：系统中每增加一个粒子所需要的自由能。</li>
<li><strong>浓度（Concentration）</strong>：单位体积内粒子的数量。</li>
<li><strong>电场（Electric field, E）</strong>：单位电荷所受的力。</li>
<li><strong>磁场（Magnetic field, H）</strong>：单位磁偶极子所受的力矩。</li>
</ul>
<h3 id="the-relationship-between-extensive-and-intensive-quantities">
<a class="header-anchor" href="#the-relationship-between-extensive-and-intensive-quantities"></a>
The Relationship Between Extensive and Intensive Quantities
</h3><p>Extensive and intensive quantities are closely related. Typically, <strong>extensive quantities can be calculated through the product or integral of intensive quantities</strong>. For example:</p>
<ul>
<li><strong>Internal energy</strong> (U): Can be determined by the system&rsquo;s temperature, volume, and particle count.</li>
<li><strong>Volume</strong>: Can be calculated by integrating the density of each unit volume within the system.</li>
<li><strong>Heat and work</strong>: Are the products of intensive quantities (such as temperature and pressure) and extensive quantities (such as entropy and volume change).</li>
</ul>
<h3 id="temperature">
<a class="header-anchor" href="#temperature"></a>
Temperature
</h3><blockquote>
<p><strong>Law of Thermal Equilibrium</strong>: If object A is in thermal equilibrium with objects B and C respectively, then B and C must also be in thermal equilibrium when they are in thermal contact.</p>
</blockquote>
<ul>
<li>Definition of temperature: Objects in mutual thermal equilibrium must possess a physical quantity intrinsic to the object itself, which is defined as temperature.</li>
<li>Temperature scale: A numerical representation method for temperature.</li>
<li>Three essential elements: Thermometric property of the thermometric substance, fixed points, and scale division.</li>
</ul>
<h3 id="work-in-thermodynamic-processes">
<a class="header-anchor" href="#work-in-thermodynamic-processes"></a>
Work in Thermodynamic Processes
</h3><p><strong>Quasi-static process</strong>: $W=\int_{V_1}^{V_2}pdV$, where $V_1$ is the initial volume of the system and $V_2$ is the final volume.</p>
<p><strong>Reversible process: A quasi-static process without dissipation</strong></p>
<blockquote>
<p>In a reversible process, the work done by the surroundings on the system can be expressed using the system&rsquo;s state parameters.</p>
</blockquote>
<p><strong>State space: A space constructed with independent state parameters as coordinate axes.</strong></p>
<blockquote>
<p>In state space, a point represents an equilibrium state of the system, and a curve represents a quasi-static process.</p>
</blockquote>
<p><strong>Process equation: The functional relationship between independent state parameters during a quasi-static process.</strong></p>
<p>Expressions for work (considering only reversible processes):</p>
<ul>
<li>Fluid volume change process: −𝑝d𝑉</li>
<li>Surface film area change process: 𝜎d𝐴</li>
<li>Elastic wire length change: 𝐹d𝐿</li>
<li>Polarization work (ignoring dielectric volume change):</li>
<li>Magnetization work (ignoring magnetic medium volume change):</li>
</ul>
<p><strong>Work in special non-static processes</strong>:</p>
<ul>
<li>Isochoric (constant volume) $𝑊=0$</li>
<li>Isobaric (constant pressure) $𝑊=−𝑝 \Delta 𝑉$</li>
</ul>
<h3 id="carnots-theorem">
<a class="header-anchor" href="#carnots-theorem"></a>
Carnot&rsquo;s Theorem
</h3><p><strong>Among all heat engines operating between two fixed temperatures, reversible heat engines have the highest efficiency.</strong>
Corollary: If two reversible heat engines operate at the same temperatures, their efficiencies are equal.
[[#Reversible Heat Engine Efficiency $ eta$]]</p>
<h3 id="principle-of-entropy-increase">
<a class="header-anchor" href="#principle-of-entropy-increase"></a>
Principle of Entropy Increase
</h3><blockquote>
<p>Under <strong>adiabatic conditions</strong>, processes that decrease entropy are impossible (entropy never decreases).</p>
</blockquote>
<p><strong>Proof</strong>: From the second law of thermodynamics:
</p>
$$
dS\ge \frac{dQ}{T}
$$<p>
For an adiabatic system, $dQ$=0, thus $dS\ge 0$.</p>
<p><strong>Corollary</strong>: In an isolated (adiabatic) system, <strong>irreversible processes</strong> always proceed in the direction of increasing entropy, i.e., $dS \gt 0$.</p>
<h3 id="throttling-process">
<a class="header-anchor" href="#throttling-process"></a>
Throttling Process
</h3><p>In thermodynamics, the throttling process (or Joule-Thomson process) is a typical adiabatic irreversible process characterized by a significant pressure drop as the fluid passes through a throttling device (such as an orifice, valve, or porous plug), while the enthalpy remains constant. Below is a detailed analysis of this process:</p>
<h3 id="characteristics-of-throttling-process">
<a class="header-anchor" href="#characteristics-of-throttling-process"></a>
Characteristics of Throttling Process
</h3><ol>
<li><strong>Adiabatic Nature</strong>: The throttling process is typically considered adiabatic because there is insufficient time for heat exchange during rapid flow.</li>
<li><strong>Isenthalpic Property</strong>: The most important characteristic is that the throttling process is isenthalpic, meaning the enthalpy ($H$) remains unchanged before and after throttling. This property plays a key role in the throttling process. [[#Energy Equation in Throttling Process]]</li>
<li><strong>Irreversibility</strong>: The throttling process is an <strong>irreversible process</strong>, and the <strong>entropy of the system usually increases</strong>.</li>
<li><strong>Pressure Drop</strong>: During the throttling process, the pressure of the fluid drops significantly.</li>
</ol>
<h3 id="relationship-between-enthalpy-temperature-and-pressure">
<a class="header-anchor" href="#relationship-between-enthalpy-temperature-and-pressure"></a>
Relationship Between Enthalpy, Temperature, and Pressure
</h3><p>For an ideal gas, enthalpy is solely a function of temperature, so the temperature remains constant during a throttling process (a special case of the Joule-Thomson effect). However, for real gases, the situation is more complex—the temperature may increase or decrease depending on the gas&rsquo;s ==Joule-Thomson coefficient== ([[#Joule-Thomson Coefficient]]), defined as:
</p>
$$
\mu_{JT} = \left ( \frac{\partial T}{\partial P} \right)_H
$$<ul>
<li>If $\mu_{JT}  \gt  0$, the gas temperature decreases during throttling.</li>
<li>If $\mu_{JT}  \lt  0$, the gas temperature increases during throttling.</li>
</ul>
<h3 id="applications-of-throttling-process">
<a class="header-anchor" href="#applications-of-throttling-process"></a>
Applications of Throttling Process
</h3><p>The throttling process is crucial in many industrial applications, such as:</p>
<ul>
<li><strong>Refrigeration Systems</strong>: Utilizing throttling valves (e.g., expansion valves) to reduce the temperature and pressure of refrigerants, thereby achieving cooling effects.</li>
<li><strong>Gas Separation</strong>: Cooling certain gases to their liquefaction temperature through throttling to separate different gas components.</li>
</ul>
<h3 id="characteristic-functions">
<a class="header-anchor" href="#characteristic-functions"></a>
Characteristic Functions
</h3><p>Characteristic functions (also known as thermodynamic potential functions or thermodynamic potentials) are used to describe the state of a thermodynamic system. Various thermodynamic properties and relationships can be derived from these functions.<br>
There are four main characteristic functions: internal energy, Helmholtz free energy, enthalpy, and Gibbs free energy.</p>
<ul>
<li>Internal energy $U(S,V)$: For adiabatic isochoric systems, U remains constant</li>
<li>Enthalpy $H(S,p)$: For adiabatic isobaric systems, H remains constant</li>
<li>Free energy $F(T,V)$: For isothermal isochoric systems, F remains constant</li>
<li>Gibbs function $G(T,p)$: For isothermal isobaric systems, G remains constant</li>
</ul>
<h3 id="radiation">
<a class="header-anchor" href="#radiation"></a>
Radiation
</h3><h3 id="equilibrium-radiation">
<a class="header-anchor" href="#equilibrium-radiation"></a>
Equilibrium Radiation
</h3><p>Equilibrium Radiation refers to the electromagnetic radiation characteristics of a system in a state of <strong>thermodynamic equilibrium</strong>. In this state, the processes of radiation and absorption within the system reach a dynamic balance, and the spectrum and intensity of electromagnetic radiation are solely determined by the system&rsquo;s temperature.</p>
<h3 id="characteristics-of-equilibrium-radiation">
<a class="header-anchor" href="#characteristics-of-equilibrium-radiation"></a>
Characteristics of Equilibrium Radiation
</h3><ol>
<li><strong>Thermal Equilibrium</strong>:
<ul>
<li>The system is in a state of thermal equilibrium, meaning the temperature is uniform across all parts, with no net heat flow.</li>
</ul>
</li>
<li><strong>Blackbody Radiation</strong>:
<ul>
<li>Under thermal equilibrium, the radiation exhibits blackbody characteristics. Blackbody radiation is an idealized form of radiation that perfectly absorbs and re-emits all frequencies, with its spectrum depending solely on temperature.</li>
</ul>
</li>
</ol>
<h3 id="radiative-flux-density">
<a class="header-anchor" href="#radiative-flux-density"></a>
Radiative Flux Density
</h3><p>辐射通量密度（Radiative Flux Density），也称为辐射出射度或辐射强度，是指单位时间内通过单位面积的辐射能量。它是描述辐射场的重要物理量，在热力学、气象学、天文学等领域有广泛应用。</p>
<h3 id="definition">
<a class="header-anchor" href="#definition"></a>
Definition
</h3><p>The radiant flux density E is expressed as:
</p>
$$
E = \frac{d\Phi}{dA}
$$<p>
where:</p>
<ul>
<li>E is the radiant flux density, measured in watts per square meter ($W/m^2$).</li>
<li>$d\Phi$ is the radiant energy flux passing through area $dA$, measured in watts ($W$).</li>
<li>$dA$ is the area, measured in square meters ($m^2$).</li>
</ul>
<p>==Or==
</p>
$$
J_u=\frac{c}{4}\frac{U}{V}=\frac{1}{4}cu
$$<p>
where:</p>
<ul>
<li>$c$: speed of light</li>
<li>$u$: radiant energy density, $u=\frac{U}{V}$, the equilibrium radiant energy per unit volume</li>
</ul>
<h3 id="applications-of-the-stefan-boltzmann-law">
<a class="header-anchor" href="#applications-of-the-stefan-boltzmann-law"></a>
Applications of the Stefan-Boltzmann Law
</h3><p>For an ideal black body, its radiant flux density is proportional to the fourth power of its temperature, as described by the Stefan-Boltzmann law:
</p>
$$
E = \sigma T^4
$$<p>
where:</p>
<ul>
<li>$E$ is the radiant flux density of the black body.</li>
<li>$\sigma$ is the Stefan-Boltzmann constant, with a value of $5.67 \times 10^{-8} \, W \, m^{-2} \, K^{-4}$.</li>
<li>$T$ is the absolute temperature of the black body, in Kelvin ($K$).</li>
</ul>
<h3 id="directionality-of-radiant-flux-density">
<a class="header-anchor" href="#directionality-of-radiant-flux-density"></a>
Directionality of Radiant Flux Density
</h3><p>辐射通量密度是一个矢量量，考虑其方向性时称为辐射强度 $I$，表示为单位立体角上的辐射通量密度：
</p>
$$
I = \frac{d\Phi}{dA \cos \theta \, d\Omega}
$$<p>
其中：</p>
<ul>
<li>$I$ 是辐射强度，单位是瓦特每平方米每立体角（$W/m^2/sr$）。</li>
<li>$d\Omega$ 是立体角，单位是球面度（sr）。</li>
<li>$\theta$ 是辐射方向与法线之间的夹角。</li>
</ul>
<h3 id="blackbody-radiation-cavity-radiation">
<a class="header-anchor" href="#blackbody-radiation-cavity-radiation"></a>
Blackbody Radiation (Cavity Radiation)
</h3><p>黑体辐射（Blackbody Radiation）是指一个理想化的物体（称为黑体）在热平衡状态下发出的电磁辐射。黑体是一个理想化的概念，它具有完全吸收和完全辐射所有频率电磁波的能力。黑体辐射的特性只取决于黑体的温度，而与其材料或表面性质无关。</p>
<h2 id="">
<a class="header-anchor" href="#"></a>

</h2><h3 id="characteristics-of-blackbody-radiation">
<a class="header-anchor" href="#characteristics-of-blackbody-radiation"></a>
Characteristics of Blackbody Radiation
</h3><ol>
<li><strong>Perfect Absorption</strong>:
<ul>
<li>A blackbody can completely absorb all incident electromagnetic radiation, regardless of wavelength or direction. Thus, a blackbody exhibits no reflection or transmission at any wavelength.</li>
</ul>
</li>
<li><strong>Perfect Emission</strong>:
<ul>
<li>A blackbody emits electromagnetic radiation at any temperature, radiating across all wavelengths. This radiation is solely determined by the blackbody&rsquo;s temperature.</li>
</ul>
</li>
</ol>
<h3 id="practical-applications-of-blackbody-radiation">
<a class="header-anchor" href="#practical-applications-of-blackbody-radiation"></a>
Practical Applications of Blackbody Radiation
</h3><ol>
<li><strong>Cosmic Background Radiation</strong>:
<ul>
<li>The cosmic microwave background radiation closely approximates blackbody radiation at a temperature of 2.725 K, serving as crucial evidence for the Big Bang theory.</li>
</ul>
</li>
<li><strong>Stellar Spectra</strong>:
<ul>
<li>The radiation from stars approximates blackbody radiation, allowing their surface temperatures to be deduced by analyzing their spectra.</li>
</ul>
</li>
<li><strong>Infrared Thermometry</strong>:
<ul>
<li>Based on the principle of blackbody radiation, the surface temperature of an object can be inferred by measuring the intensity of the infrared radiation it emits.</li>
</ul>
</li>
<li><strong>Thermal Imaging</strong>:
<ul>
<li>Thermal imaging devices utilize the principle of blackbody radiation, detecting infrared radiation emitted from an object&rsquo;s surface to create temperature-based images.</li>
</ul>
</li>
</ol>
<h3 id="criteria-for-thermodynamic-equilibrium">
<a class="header-anchor" href="#criteria-for-thermodynamic-equilibrium"></a>
Criteria for Thermodynamic Equilibrium
</h3><h3 id="isolated-system">
<a class="header-anchor" href="#isolated-system"></a>
Isolated System
</h3><p>For an isolated (adiabatic) system, S increases or remains constant, reaching $S_\max$ at equilibrium.
If a small disturbance occurs, $\Delta S \lt 0$, expand it to a second-order Taylor series $\Delta S=\delta S+\frac{1}{2}\delta^2 S$
At equilibrium, the conditions $\delta S=0,\delta^2 S \lt 0$ must be satisfied.</p>
<h3 id="isothermal-isochoric-system">
<a class="header-anchor" href="#isothermal-isochoric-system"></a>
Isothermal-Isochoric System
</h3><p>等温等容系统 F decreases or remains constant, reaching $F_\min$ at equilibrium.
If a small disturbance occurs, $\Delta F \gt 0$, which can be expanded as a second-order Taylor series $\Delta F=\delta F+\frac{1}{2}\delta^2 F$.
At equilibrium, the conditions $\delta F=0,\delta^2 F \gt 0$ are satisfied.</p>
<h3 id="isothermal-isobaric-system">
<a class="header-anchor" href="#isothermal-isobaric-system"></a>
Isothermal-Isobaric System
</h3><p>In an isothermal-isobaric system, G decreases or remains constant, reaching $G_\min$ at equilibrium.
If a slight perturbation occurs, $\Delta G \gt 0$, which can be expanded as a second-order Taylor series: $\Delta G=\delta G+\frac{1}{2}\delta^2 G$.
At equilibrium, the conditions $\delta G=0$ and $\delta^2 G>0$ must be satisfied.</p>
<h3 id="equilibrium-conditions-for-two-phase-systems-in-a-unit">
<a class="header-anchor" href="#equilibrium-conditions-for-two-phase-systems-in-a-unit"></a>
Equilibrium Conditions for Two-Phase Systems in a Unit
</h3>$$
\begin{cases}  
T^\alpha=T^\beta \\  
p^\alpha=p^\beta \\  
\mu^\alpha=\mu^\beta  
\end{cases}  
$$<h3 id="classification-of-phase-transitions">
<a class="header-anchor" href="#classification-of-phase-transitions"></a>
Classification of Phase Transitions
</h3><p>[[#Ehrenfest Equations]]
Ehrenfest (Paul Ehrenfest) classified phase transitions <strong>based on the continuity of thermodynamic potential functions and their derivatives at the phase transition point</strong>.
Ehrenfest categorized phase transitions into first-order and second-order phase transitions.</p>
<h3 id="first-order-phase-transition">
<a class="header-anchor" href="#first-order-phase-transition"></a>
First-Order Phase Transition
</h3><blockquote>
<p>A first-order phase transition refers to a phase transition where the first derivatives of the thermodynamic potential (such as Gibbs free energy), like entropy and volume, are discontinuous at the transition point.</p>
</blockquote>
<p>Typical first-order phase transitions include melting, vaporization, and sublimation.
[[#First-Order Phase Transition Equations]]</p>
<h3 id="second-order-phase-transition">
<a class="header-anchor" href="#second-order-phase-transition"></a>
Second-order Phase Transition
</h3><blockquote>
<p>A second-order phase transition refers to a phase transition process where the first derivative of the thermodynamic potential function is continuous at the transition point, but its second derivative exhibits discontinuity.</p>
</blockquote>
<p>Typical second-order phase transitions include the transition of superconductors and liquid crystal phase transitions.<br>
[[#Second-order Phase Transition Equations]]</p>
<h3 id="multicomponent-system-phase-equilibrium-conditions">
<a class="header-anchor" href="#multicomponent-system-phase-equilibrium-conditions"></a>
Multicomponent System Phase Equilibrium Conditions
</h3><h3 id="equilibrium-conditions-for-isothermal-isobaric-systems">
<a class="header-anchor" href="#equilibrium-conditions-for-isothermal-isobaric-systems"></a>
Equilibrium Conditions for Isothermal-Isobaric Systems
</h3>$$
\begin{cases}
T^\alpha=T^\beta \\
p^\alpha=p^\beta \\
\mu_i^\alpha=\mu_i^\beta
\end{cases}
$$<h3 id="the-same-applies-to-other-systems">
<a class="header-anchor" href="#the-same-applies-to-other-systems"></a>
The same applies to other systems
</h3><h3 id="gibbs-paradox">
<a class="header-anchor" href="#gibbs-paradox"></a>
Gibbs Paradox
</h3><p><strong>Description</strong>: The entropy increases abruptly from $2nR\ln 2$ to 0 when transitioning from two arbitrarily similar gases to the same gas.<br>
<strong>Reason</strong>: Indistinguishability of identical particles.</p>
<h3 id="the-third-law-of-thermodynamics">
<a class="header-anchor" href="#the-third-law-of-thermodynamics"></a>
The Third Law of Thermodynamics
</h3><h3 id="nernsts-theorem">
<a class="header-anchor" href="#nernsts-theorem"></a>
Nernst&rsquo;s Theorem
</h3><p>The change in entropy of a condensed system during an isothermal process approaches zero as the thermodynamic temperature tends to zero.
</p>
$$
\lim_{T \to 0}(\Delta S)_T=0
$$<h3 id="principle-of-the-unattainability-of-absolute-zero">
<a class="header-anchor" href="#principle-of-the-unattainability-of-absolute-zero"></a>
Principle of the Unattainability of Absolute Zero
</h3><p>It is impossible to cool an object to thermodynamic zero through <strong>a finite number of steps</strong>.</p>
<h3 id="corollaries">
<a class="header-anchor" href="#corollaries"></a>
Corollaries
</h3><ol>
<li>$\lim_{T \to 0}C_n=0$ This corollary is proven by Einstein&rsquo;s quantum statistics in the [[#solid]] section.</li>
<li>$\lim_{T \to 0}\alpha=0,\lim_{T \to 0}\beta=0$</li>
<li>At absolute zero, the entropy value can be 0</li>
</ol>
<h3 id="classical-description-of-particle-motion-states">
<a class="header-anchor" href="#classical-description-of-particle-motion-states"></a>
Classical Description of Particle Motion States
</h3><h3 id="space">
<a class="header-anchor" href="#space"></a>
$\mu$ Space
</h3><blockquote>
<p>Describes the mechanical motion state of particles using position and momentum coordinates. The $2r$ variables $q_1,q_2,...,q_r;p_1,p_2,...,p_r$ form a $2r$-dimensional space with Cartesian coordinates, known as $\mu$ space.</p>
</blockquote>
<p><strong>Degrees of Freedom of a Particle</strong>: The minimum number of coordinates required to determine the spatial position <strong>or</strong> configuration position of a particle.</p>
<h3 id="phase-cell">
<a class="header-anchor" href="#phase-cell"></a>
Phase Cell
</h3><blockquote>
<p>The volume of a particle&rsquo;s motion state in $\mu$ space</p>
</blockquote>
<p>For a particle with $r$ degrees of freedom, the size of the phase cell is $\Delta q_1\cdots \Delta q_r\Delta p_1\cdots \Delta p_r$
When it satisfies the [[#Uncertainty Relation]] $\Delta q_i\Delta p_i\approx h$, the size of the phase cell $\approx h^r$</p>
<h3 id="free-particle">
<a class="header-anchor" href="#free-particle"></a>
Free Particle
</h3><blockquote>
<p>A free particle is a particle that moves freely without the influence of any force.
In the absence of an external field, molecules of an ideal gas or free electrons in a metal can be approximately regarded as free particles.</p>
</blockquote>
<p><strong>Momentum of a free particle</strong>: $p_x=m\dot{x}$, the momentum of the particle in the x-axis direction. Here, m is the mass of the particle, and $\dot{x}$ is the velocity of the particle in the x-axis direction (the first derivative of position with respect to time).
<strong>Energy of a free particle</strong>:
</p>
$$
\varepsilon=\frac{1}{2m}\sum_ip_i^2
$$<h3 id="linear-harmonic-oscillator">
<a class="header-anchor" href="#linear-harmonic-oscillator"></a>
Linear Harmonic Oscillator
</h3><blockquote>
<p>Classical mechanics tells us that a particle with mass $m$ under the action of an elastic force $F = -Ax$ will perform simple harmonic motion along the $x$-axis near the origin, known as a linear harmonic oscillator.</p>
</blockquote>
<ul>
<li>Angular frequency of oscillation:
$$\omega=\sqrt{\frac{A}{m}}$$</li>
<li>$A$: Elastic coefficient</li>
</ul>
<h3 id="energy-of-a-linear-harmonic-oscillator-with-one-degree-of-freedom">
<a class="header-anchor" href="#energy-of-a-linear-harmonic-oscillator-with-one-degree-of-freedom"></a>
Energy of a Linear Harmonic Oscillator with One Degree of Freedom
</h3><p>The sum of kinetic and potential energy
</p>
$$
\varepsilon=\frac{p^2}{2m}+\frac{A}{2}x^2=\frac{p^2}{2m}+\frac{1}{2}m\omega^2x^2
$$<p>
Written in the standard form of an ellipse equation
</p>
$$
\frac{p^2}{2m\varepsilon}+\frac{x^2}{\frac{2\varepsilon}{m\omega^2}}=1
$$<h3 id="rotor-degree-of-freedom-2">
<a class="header-anchor" href="#rotor-degree-of-freedom-2"></a>
Rotor (Degree of Freedom: 2)
</h3><blockquote>
<p>Consider the motion of a particle \( P \) with mass \( m \) attached to the origin \( O \) by a light rod of fixed length.</p>
</blockquote>
<ul>
<li>In Cartesian coordinates, the position of the particle is determined by coordinates \( x, y, z \).</li>
</ul>
<p>Particle energy (kinetic energy):
</p>
$$
\varepsilon=\frac{1}{2}m(\dot{x}^2+\dot{y}^2+\dot{z}^2)
$$<ul>
<li>Using spherical coordinates \( r, \theta, \varphi \) to describe the particle&rsquo;s position: \( x = r\sin \theta \cos \varphi, y = r\sin \theta \sin \varphi, z = r\cos \theta \)</li>
</ul>
<p>Particle energy:
</p>
$$
\begin{align}
\varepsilon &=\frac{1}{2}m(\dot{r}^2+r^2\dot{\theta}^2+r^2\sin^2\theta\dot{\varphi}^2)\\
&=\frac{1}{2}m(r^2\dot{\theta}^2+r^2\sin^2\dot{\varphi}^2)
\end{align}
$$<ul>
<li>Introducing conjugate momenta \( p_\theta = mr^2\dot{\theta}, p_\varphi = mr^2\sin^2\theta \dot{\varphi} \), and the moment of inertia \( I \) of the particle relative to \( O \).</li>
<li>\( 0  \lt  \theta  \lt  \pi, 0  \lt  \varphi  \lt  2\pi \), two degrees of freedom.</li>
</ul>
$$
\varepsilon=\frac{1}{2I}(p_\theta^2+\frac{1}{\sin^2\theta}p_\varphi^2)
$$<h3 id="diatomic-molecules">
<a class="header-anchor" href="#diatomic-molecules"></a>
Diatomic Molecules
</h3><blockquote>
<p>The two-body problem can be reduced to a single-body problem. In statistical physics, the rotation of a diatomic molecule about its center of mass is treated as a rotor.</p>
</blockquote>
<p>According to classical mechanics, in the absence of external forces, the total angular momentum of the rotor $\vec L=\vec r\times \vec p$ is a conserved quantity, meaning neither its magnitude nor direction changes over time. Since $\vec r$ is perpendicular to $\vec L$, the motion of the particle lies in a plane perpendicular to $\vec L$. If we choose the $z$-axis to be parallel to $\vec L$, the particle&rsquo;s motion must lie within the $xy$-plane. This is equivalent to fixing $\theta=\frac{\pi}{2}, p_\theta=0$.</p>
<p>The particle&rsquo;s energy in this case is:
</p>
$$
\varepsilon=\frac{p_\varphi^2}{2I}=\frac{L^2}{2I}
$$<p>
where: $L^2=\vec{L}\cdot \vec{L}$<br>
<strong>Here, the particle has 1 degree of freedom</strong>, namely $\varphi$.</p>
<h3 id="quantum-description-of-particle-motion-states">
<a class="header-anchor" href="#quantum-description-of-particle-motion-states"></a>
Quantum Description of Particle Motion States
</h3><h3 id="de-broglie-wave">
<a class="header-anchor" href="#de-broglie-wave"></a>
de Broglie Wave
</h3><blockquote>
<p>Microscopic particles (photons, electrons, protons, neutrons, and even atoms, molecules, etc.) universally exhibit wave-particle duality.
A free particle with energy $\varepsilon$ and momentum $\vec{p}$ is associated with a plane wave of circular frequency $\omega$ and wave vector $\vec{k}$, known as the de Broglie wave.</p>
</blockquote>
<p>[[#de Broglie Relations]]</p>
<h3 id="quantum-state">
<a class="header-anchor" href="#quantum-state"></a>
Quantum State
</h3><blockquote>
<p>In quantum mechanics, the motion state of microscopic particles is called the quantum state.
The quantum state is characterized by a set of quantum numbers, and the <strong>number of quantum numbers</strong> is equal to the degrees of freedom of the particle.</p>
</blockquote>
<h1 id="">
<a class="header-anchor" href="#"></a>

</h1><h3 id="linear-harmonic-oscillator-energy">
<a class="header-anchor" href="#linear-harmonic-oscillator-energy"></a>
==Linear Harmonic Oscillator Energy==
</h3><p>[[#Energy of a Linear Harmonic Oscillator with 1 Degree of Freedom]]
In ==quantum physics==, the possible energy values of a linear harmonic oscillator with angular frequency $\omega$ are:
</p>
$$
\varepsilon_n=\hbar\omega(n+\frac{1}{2}),n=0,1,2,\cdots
$$<p>
Here, n is the quantum number representing the oscillator&rsquo;s motion state and energy, with only one such number—meaning the linear harmonic oscillator has 1 degree of freedom.
Clearly, the values of $\varepsilon_n$ are <strong>discrete</strong>.</p>
<h3 id="energy-levels">
<a class="header-anchor" href="#energy-levels"></a>
Energy Levels
</h3><blockquote>
<p>Discrete energies are called energy levels.</p>
</blockquote>
<p>In $\mu$ space, surfaces with equal energy are called ==equipotential surfaces==. Since these surfaces are discontinuous, each equipotential surface is referred to as an <strong>energy level</strong>.
[[#$ mu$ space]]</p>
<h3 id="rotor-energy">
<a class="header-anchor" href="#rotor-energy"></a>
Rotor Energy
</h3><p>[[#Diatomic Molecule]]
In quantum physics, the value of $L^2$ can only take discrete values:
</p>
$$
L^2=l(l_1)\hbar^2,l=0,1,2,\cdots
$$<p>
For a given angular momentum $L$, its projection $L_z$ along its eigen-direction (taken as the $z$-axis) can only take discrete values:
</p>
$$
L_z=m_l\hbar,m_l=-l,-l+1,\cdots,l
$$<p>
Clearly, there are $2l+1$ possible values for $m_l$.</p>
<blockquote>
<p>A rotor&rsquo;s motion state with <strong>2 degrees of freedom</strong> is characterized by <strong>two quantum numbers</strong> $l,m_l$.
In classical theory, the spatial orientation of the motion plane is arbitrary, whereas in quantum theory, $m_l$ can only take the above discrete values, known as <strong>spatial quantization</strong>.</p>
</blockquote>
<p>Thus, the rotor&rsquo;s energy is also quantized:
</p>
$$
\varepsilon_l=\frac{l(l+1)\hbar^2}{2I},l=0,1,2,\cdots
$$<h3 id="degeneracy">
<a class="header-anchor" href="#degeneracy"></a>
Degeneracy $\omega_l$
</h3><p>$\because$ It is observed that: the energy depends only on the quantum number $l$, and the quantum state depends on the two quantum numbers $m_l, l$.
$\therefore$ There are $2l+1$ quantum states for the energy level $\varepsilon_l$.</p>
<blockquote>
<p>Generally speaking, if a given energy level has more than one quantum state, it is called degenerate.
If an energy level has only one quantum state, it is called non-degenerate.
If an energy level has k quantum states, its degeneracy is k.</p>
</blockquote>
<p>Between two adjacent [[#Energy Levels]], there are several [[#Phase Cells]], and the number of phase cells is the <strong>degeneracy</strong> of the energy level.</p>
<h3 id="spin-angular-momentum">
<a class="header-anchor" href="#spin-angular-momentum"></a>
Spin Angular Momentum
</h3><blockquote>
<p>Certain fundamental particles possess intrinsic angular momentum, known as spin angular momentum $S$.</p>
</blockquote>
$$
S^2=s(s+1)\hbar^2
$$<p>
$s$: ==Spin quantum number==. Can be integer or half-integer, an inherent property of particles that can be summed. [[#Bosons and Fermions]]
The state of spin angular momentum is determined by its magnitude (spin quantum number $s$) and its projection along its <strong>eigenaxis</strong>. Using $z$ to denote the eigenaxis, the possible values of $S_z$ are
</p>
$$
S_z=m_s\hbar,m_s=s,s-1,\cdots,-s
$$<p>
totaling $2s+1$ values.</p>
<h3 id="energy-of-free-particles">
<a class="header-anchor" href="#energy-of-free-particles"></a>
Energy of Free Particles
</h3><h2 id="">
<a class="header-anchor" href="#"></a>

</h2><h3 id="boundary-conditions">
<a class="header-anchor" href="#boundary-conditions"></a>
Boundary Conditions
</h3><blockquote>
<p>To determine the possible states of motion of a particle, it is necessary to know the boundary conditions of the de Broglie wave at the container walls.
Typically, standing wave conditions or periodic boundary conditions are employed.</p>
</blockquote>
<p>For a <strong>one-dimensional</strong> free particle in a container of length $L$, the periodic boundary condition requires that the possible states of motion must satisfy that an integer multiple of the de Broglie wavelength $\lambda$ equals the container length $L$. That is,
</p>
$$
L=\left\vert n_x \right\vert\lambda,\left\vert n_x \right\vert=0,1,2,\cdots
$$<ul>
<li>$n_x$: The projection of the quantum number $n$ in the $x$-direction. For a one-dimensional free particle, the quantum number is simply $n_x$.</li>
</ul>
<h3 id="one-dimensional-free-particle">
<a class="header-anchor" href="#one-dimensional-free-particle"></a>
One-Dimensional Free Particle
</h3><ul>
<li>Wave vector $k_x$:

$$
k_x=\frac{2\pi}{L}n_x,n_x=0,\pm 1,\pm 2,\cdots
$$</li>
<li>Momentum $p_x$:

$$
p_x=\frac{2\pi\hbar}{L}n_x,n_x=0,\pm 1,\pm 2,\cdots
$$</li>
<li>Energy $\varepsilon_{n_x}$

$$
\varepsilon_{n_x}=\frac{p_x^2}{2m}=\frac{2\pi^2\hbar^2}{m}\cdot\frac{n_x^2}{L^2},n_x=0,\pm 1,\pm 2,\cdots
$$</li>
</ul>
<h3 id="three-dimensional-free-particle">
<a class="header-anchor" href="#three-dimensional-free-particle"></a>
Three-Dimensional Free Particle
</h3><ul>
<li>Momentum:

$$
\begin{cases}
p_x=\frac{2\pi\hbar}{L}n_x,n_x=0,\pm 1,\pm 2,\cdots\\
p_y=\frac{2\pi\hbar}{L}n_y,n_y=0,\pm 1,\pm 2,\cdots\\
p_z=\frac{2\pi\hbar}{L}n_z,n_z=0,\pm 1,\pm 2,\cdots
\end{cases}
$$</li>
<li>Energy:

$$
\varepsilon=\frac{1}{2m}(p_x^2+p_y^2+p_z^2)=\frac{2\pi^2\hbar^2}{m}\frac{n_x^2+n_y^2+n_z^2}{L^2}
$$</li>
<li>Quantum numbers: $n_x,n_y,n_z$ (three in total, corresponding to three degrees of freedom)</li>
<li>Energy levels: Determined by $n_x^2+n_y^2+n_z^2$</li>
</ul>
<h3 id="description-of-the-microscopic-motion-states-of-a-system">
<a class="header-anchor" href="#description-of-the-microscopic-motion-states-of-a-system"></a>
Description of the Microscopic Motion States of a System
</h3><h3 id="microstates">
<a class="header-anchor" href="#microstates"></a>
Microstates
</h3><blockquote>
<p>The microscopic state of a system is its mechanical motion state.<br>
It manifests as different occupation patterns where particles occupy different energy levels and phase cells.</p>
</blockquote>
<h3 id="system">
<a class="header-anchor" href="#system"></a>
System
</h3><p>==Limited to systems composed of identical and nearly independent particles==</p>
<ul>
<li><strong>System composed of identical particles</strong>: A system consisting of the same type of particles with identical intrinsic properties (same mass, charge, spin, etc.).
Examples include a free electron gas composed of free electrons, helium gas composed of $^4He$ atoms, etc. ^f1b50b</li>
<li><strong>System composed of nearly independent particles</strong>: The interaction between particles is very weak, and the average energy of interaction is much smaller than the average energy of a single particle, so the interaction between particles can be ignored, and the total energy of the system can be expressed as the sum of the energies of individual particles.
Examples include systems composed of ideal gases. ^24ebad

$$
E=\sum^{N}_{i=1}\varepsilon_i
$$
Where:</li>
<li>$\varepsilon_i$: Energy of the i-th particle</li>
<li>$N$: Total number of particles</li>
</ul>
<h3 id="classical-mechanics-description">
<a class="header-anchor" href="#classical-mechanics-description"></a>
Classical Mechanics Description
</h3><blockquote>
<p>In classical physics, identical particles [[#^f1b50b]] are distinguishable.
Determining the microscopic state of motion of a system requires determining the mechanical state of motion (individual quantum state) of each particle.</p>
</blockquote>
<p>Number of variables: $2Nr$<br>
$r$: degrees of freedom of a single particle</p>
<h3 id="bosons-and-fermions">
<a class="header-anchor" href="#bosons-and-fermions"></a>
Bosons and Fermions
</h3><blockquote>
<p>In nature, microscopic particles can be divided into two categories, known as bosons and fermions.</p>
</blockquote>
<p>[[#Spin Angular Momentum]]</p>
<h2 id="">
<a class="header-anchor" href="#"></a>

</h2><h3 id="fermion">
<a class="header-anchor" href="#fermion"></a>
Fermion
</h3><blockquote>
<p>Spin quantum number is a half-integer</p>
</blockquote>
<p>Examples: electron, proton, neutron, $\mu$ meson (all with $\frac{1}{2}$ spin); $^2H$ atom, $^3H$ nucleus (with $\frac{3}{2}$ spin), etc.</p>
<h3 id="boson">
<a class="header-anchor" href="#boson"></a>
Boson
</h3><blockquote>
<p>Spin quantum number is an integer</p>
</blockquote>
<p>Examples: photon (1), $\pi$ meson (0), $^1H$ atom (1), $^4He$ atom (4), etc.</p>
<h3 id="pauli-exclusion-principle">
<a class="header-anchor" href="#pauli-exclusion-principle"></a>
Pauli Exclusion Principle
</h3><blockquote>
<p>In a system containing multiple identical and nearly independent fermions, an individual quantum state can accommodate at most one fermion.</p>
</blockquote>
<h3 id="distribution">
<a class="header-anchor" href="#distribution"></a>
Distribution
</h3><blockquote>
<p>A set composed of the number of particles at each <strong>energy level</strong>, this set is called a distribution.</p>
</blockquote>
$$
\{a_l\}=a_1,a_2,\cdots,a_l,\cdots
$$<p>
$a_l$: The number of particles at energy level $\varepsilon_l$
Satisfies:
</p>
$$
\sum_la_l=N,\sum_la_l\varepsilon_l=E
$$<p>
A distribution encompasses a large number of possible microstates, and a microstate corresponds to a distribution.</p>
<h3 id="three-systems">
<a class="header-anchor" href="#three-systems"></a>
Three Systems
</h3><h2 id="">
<a class="header-anchor" href="#"></a>

</h2><h3 id="boltzmann-system">
<a class="header-anchor" href="#boltzmann-system"></a>
Boltzmann System
</h3><p>A system composed of <strong>identical</strong> [[#^f1b50b]]<strong>nearly independent particles</strong>[[#^24ebad]] that are ==distinguishable==, with <strong>no restriction on the number of particles</strong> occupying a single quantum state, is called a Boltzmann system.</p>
<h3 id="bose-system">
<a class="header-anchor" href="#bose-system"></a>
Bose System
</h3><p>A system composed of bosons, not constrained by the Pauli exclusion principle.
In a Bose system consisting of multiple <strong>identical and nearly independent</strong> bosons, <strong>the number of bosons occupying the same individual quantum state</strong> is unrestricted.</p>
<h3 id="fermi-system">
<a class="header-anchor" href="#fermi-system"></a>
Fermi System
</h3><p>A system composed of fermions, obeying the Pauli exclusion principle.</p>
<h3 id="three-types-of-statistics">
<a class="header-anchor" href="#three-types-of-statistics"></a>
Three Types of Statistics
</h3><h2 id="">
<a class="header-anchor" href="#"></a>

</h2><h3 id="boltzmann-statistics-classical-statistics">
<a class="header-anchor" href="#boltzmann-statistics-classical-statistics"></a>
Boltzmann Statistics (Classical Statistics)
</h3><p>Distinguishable identical particles, with no restriction on the number of particles per phase cell.</p>
<h3 id="bose-statistics">
<a class="header-anchor" href="#bose-statistics"></a>
Bose Statistics
</h3><p>Indistinguishable identical particles, with no restriction on the number of particles per phase cell.</p>
<h3 id="fermi-statistics">
<a class="header-anchor" href="#fermi-statistics"></a>
Fermi Statistics
</h3><p>Indistinguishable identical particles, with a maximum of one particle per phase cell.</p>
<h3 id="principle-of-equal-probability">
<a class="header-anchor" href="#principle-of-equal-probability"></a>
Principle of Equal Probability
</h3><blockquote>
<p>For an <strong>isolated</strong> system in <strong>equilibrium</strong>, all possible microstates of the system are equally probable.</p>
</blockquote>
<h3 id="most-probable-distribution">
<a class="header-anchor" href="#most-probable-distribution"></a>
Most Probable Distribution
</h3><blockquote>
<p>According to the principle of equal a priori probabilities, for an isolated system in equilibrium, the distribution with the greatest number of microscopic states has the highest probability of occurrence, known as the most probable distribution.</p>
</blockquote>
<h3 id="boltzmann-distribution">
<a class="header-anchor" href="#boltzmann-distribution"></a>
Boltzmann Distribution
</h3><blockquote>
<p>The most probable distribution of particles in a Boltzmann system is called the <strong>Maxwell-Boltzmann distribution</strong>, or simply the Boltzmann distribution.</p>
</blockquote>
<ul>
<li>Approximate equation: $\ln m!=m (\ln m-1),m\gg 1$</li>
<li>Denote $\Omega_{m.B.}$ as $\Omega$<br>
Under Boltzmann statistics, this distribution contains the most microstates, with the maximum number of microstates $\Omega$.<br>

$$
a_l=\omega_le^{-\alpha-\beta\varepsilon_l}=e^{-\alpha-\beta\varepsilon_l}\frac{\Delta\omega_l}{h_0^r}
$$<br>
The latter is the ==classical expression of the Boltzmann distribution==.<br>

$$
N=\sum_la_l=\sum_l\omega_le^{-\alpha-\beta\varepsilon_l}=\sum_se^{-\alpha-\beta\varepsilon_s}
$$<br>

$$
E=\sum_l\varepsilon_la_l=\sum_l\omega_l\varepsilon_le^{-\alpha-\beta\varepsilon_l}=\sum_s\varepsilon_se^{-\alpha-\beta\varepsilon_s}
$$<br>

$$
f_s=\frac{a_l}{\omega_l}=e^{-\alpha-\beta\varepsilon_l}
$$</li>
<li>$a_l$: The most probable number of particles (average number of particles) at energy level $\varepsilon_l$, which should satisfy the condition $a_l\gg 1$.</li>
<li>$\alpha$: Lagrange multiplier.</li>
<li>$\beta$: Lagrange multiplier, generally determined by experimental conditions.</li>
<li>$\omega_l$: The number of quantum states at energy level $\varepsilon_l$, i.e., the degeneracy. Should satisfy the condition $\omega_l\gg 1$.</li>
<li>$s$: A quantum state at energy level $\varepsilon_l$.</li>
<li>$\sum_l$: Sum over quantum number $l$.</li>
<li>$\sum_s$: Sum over all quantum states.</li>
<li>$f_s$: The average number of particles in quantum state $s$.</li>
</ul>
<h3 id="bose-distribution">
<a class="header-anchor" href="#bose-distribution"></a>
Bose Distribution
</h3><blockquote>
<p>The most probable distribution of particles in a Bose system is called the Bose-Einstein distribution, abbreviated as the Bose distribution.</p>
</blockquote>
$$
a_l=\frac{\omega_l}{e^{\alpha+\beta\varepsilon_l}-1}
$$<p>
</p>
$$
N=\sum_la_l=\sum_l\frac{\omega_l}{e^{\alpha+\beta\varepsilon_l}-1}
$$<p>
</p>
$$
E=\sum_l\varepsilon_la_l=\sum_l\frac{\varepsilon_l\omega_l}{e^{\alpha+\beta\varepsilon_l}-1}
$$<p>
</p>
$$
F_s=\frac{a_l}{\omega_l}=\frac{1}{e^{\alpha+\beta\varepsilon_l}-1}
$$<h3 id="fermi-distribution">
<a class="header-anchor" href="#fermi-distribution"></a>
Fermi Distribution
</h3><blockquote>
<p>The most probable distribution of a Fermi system.</p>
</blockquote>
$$
a_l=\frac{\omega_l}{e^{\alpha+\beta\varepsilon_l}+1}
$$<p>
</p>
$$
N=\sum_la_l=\sum_l\frac{\omega_l}{e^{\alpha+\beta\varepsilon_l}+1}
$$<p>
</p>
$$
E=\sum_l\varepsilon_la_l=\sum_l\frac{\varepsilon_l\omega_l}{e^{\alpha+\beta\varepsilon_l}+1}
$$<p>
</p>
$$
F_s=\frac{a_l}{\omega_l}=\frac{1}{e^{\alpha+\beta\varepsilon_l}+1}
$$<h3 id="relationship-among-the-three-distributions">
<a class="header-anchor" href="#relationship-among-the-three-distributions"></a>
Relationship Among the Three Distributions
</h3><p>If $\alpha$ satisfies the condition ([[#Classical Limit Condition]], non-degeneracy condition)
</p>
$$
e^\alpha \gg 1
$$<p>
then the $\pm$ in the Bose and Fermi distributions can be neglected, and both distributions transition to the Boltzmann distribution.
In this case, for all $l$, we have
</p>
$$
\frac{a_l}{\omega_l}\ll 1
$$<p>
</p>
$$
\Omega_{B.E.}\approx\frac{\Omega_{M.B.}}{N!}\approx\Omega_{F.D.}
$$<h3 id="localization">
<a class="header-anchor" href="#localization"></a>
Localization
</h3><ul>
<li>In nature, some systems can be considered as composed of localized particles.
For example, atoms or ions in a crystal oscillate slightly around their equilibrium positions.</li>
<li>Although these particles are <strong>indistinguishable by their quantum nature</strong>, they can be distinguished based on their positions, so ==localized particles== can be treated as distinguishable particles.</li>
<li>Systems composed of localized particles are called <em>localized systems</em>.
Examples include paramagnetic solids and nuclear spin systems.</li>
</ul>
<blockquote>
<p>Localized systems and Bose (Fermi) systems that satisfy classical limit conditions both obey the Boltzmann distribution.
However, their microscopic state numbers are different.</p>
</blockquote>
<h3 id="statistical-expressions-of-thermodynamic-quantities">
<a class="header-anchor" href="#statistical-expressions-of-thermodynamic-quantities"></a>
Statistical Expressions of Thermodynamic Quantities
</h3><h1 id="">
<a class="header-anchor" href="#"></a>

</h1><h3 id="internal-energy">
<a class="header-anchor" href="#internal-energy"></a>
Internal Energy
</h3>$$
U=\sum_la_l\varepsilon_l=\sum_l\varepsilon_l\omega_le^{-\alpha-\beta\varepsilon_l}
$$<h3 id="particle-partition-function">
<a class="header-anchor" href="#particle-partition-function"></a>
Particle Partition Function
</h3><p>Abbreviated as partition function
</p>
$$
Z_1=\sum_l\omega_le^{-\beta\varepsilon_l}=\sum_se^{-\beta\varepsilon_s}
$$<p>
<strong>Applicable scope: Nearly independent particle system, particles are distinguishable, and the number of particles each phase cell (quantum state) can accommodate is not limited.</strong></p>
$$
N=e^{-\alpha}\sum_l\omega_le^{-\beta\varepsilon_l}=e^{-\alpha}Z_1
$$<p>By eliminating $\alpha$ from the above three equations, we obtain the ==statistical expression of internal energy==:
</p>
$$
U=-N\frac{\partial}{\partial\beta}\ln Z_1
$$<h3 id="total-differential-of-internal-energy-1">
<a class="header-anchor" href="#total-differential-of-internal-energy-1"></a>
Total Differential of Internal Energy
</h3>$$
dU=\sum_la_ld\varepsilon_l+\sum_l\varepsilon_lda_l
$$<ul>
<li>$\sum_la_ld\varepsilon_l$: The change in internal energy caused by the alteration of energy levels due to changes in external parameters while the particle distribution remains unchanged, which represents <strong>the work done by the surroundings on the system</strong> during the process.</li>
<li>$\sum_l\varepsilon_lda_l$: The change in internal energy caused by the redistribution of particles while the energy levels remain unchanged, which represents <strong>the heat absorbed by the system from the surroundings</strong> during the process.</li>
<li>In an infinitesimal quasi-static process, the heat absorbed by the system from the surroundings equals the increase in internal energy due to <strong>the redistribution of particles among energy levels</strong>.
Clearly, heat has no microscopic counterpart and is a unique <strong>macroscopic quantity</strong> in thermal phenomena.
$dQ$ is not a total differential but merely an <strong>infinitesimal quantity</strong>.</li>
</ul>
<h3 id="entropy">
<a class="header-anchor" href="#entropy"></a>
Entropy
</h3>$$
dS=\frac{1}{T}dQ=\frac{1}{T}(dU-Ydy)
$$<p>
</p>
$$
\beta(dU-Ydy)=Nd(\ln Z_1-\beta\frac{\partial}{\partial\beta}\ln Z_1)
$$<p>
$\beta$ is a function of temperature, let
</p>
$$
\beta=\frac{1}{kT}
$$<p>
</p>
$$
k=R/N_A
$$<ul>
<li>$k$: Boltzmann constant, $k=1.381\times10^{-23}J\cdot K^{-1}$</li>
<li>$N_A$: Avogadro&rsquo;s number, $N_A=6.022\times10^{23}mol^{-1}$</li>
<li>$R$: Molar gas constant, $R=8.314J\cdot K^{-1}\cdot mol^{-1}$</li>
</ul>
$$
dS=Nkd(\ln Z_1-\beta\frac{\partial}{\partial \beta}\ln Z_1)
$$<p>
Integrating gives:
</p>
$$
S=Nk(\ln Z_1-\beta\frac{\partial}{\partial \beta}\ln Z_1)
$$<p>
The integration constant is chosen to be 0.</p>
$$
\ln Z_1=\ln N+\alpha
$$<p>
From the Boltzmann distribution
</p>
$$
a_l=\omega_le^{-\alpha-\beta\varepsilon_l}
$$<p>
we obtain
</p>
$$
\alpha+\beta\varepsilon_l=\ln\frac{\omega_l}{\alpha_l}
$$<p>
Simplifying the entropy gives
</p>
$$
\begin{align}
S&=k(N\ln N+\sum_l(\alpha+\beta\varepsilon_l)a_l)\\
&=k(N\ln N+\sum_la_l\ln \omega_l-\sum_la_l\ln a_l)\\
&=k\ln \Omega
\end{align}
$$<h3 id="boltzmann-relation">
<a class="header-anchor" href="#boltzmann-relation"></a>
Boltzmann Relation
</h3>$$
S=k\ln \Omega
$$<blockquote>
<p>The more microstates a macrostate corresponds to, the greater its disorder and the larger its entropy.</p>
</blockquote>
<ul>
<li>For distinguishable particle systems (localized systems), $\Omega=\Omega_{M.B.}$</li>
<li>For Bose (Fermi) systems satisfying the classical limit condition, $\Omega=\frac{\Omega_{M.B.}}{N!}$</li>
</ul>
<h3 id="free-energy">
<a class="header-anchor" href="#free-energy"></a>
Free Energy
</h3>$$
F=
\begin{cases}
-NkT\ln Z_1, \text{localized system}\\
-NkT\ln Z_1+kT\ln N!, \text{Bose/Fermi systems satisfying classical limit conditions}
\end{cases}
$$<h3 id="equation-of-state-for-ideal-gas">
<a class="header-anchor" href="#equation-of-state-for-ideal-gas"></a>
Equation of State for Ideal Gas
</h3>$$
p=\frac{N}{\beta}\frac{\partial}{\partial V}\ln Z_1=\frac{NkT}{V}
$$<h3 id="classical-limit-conditions">
<a class="header-anchor" href="#classical-limit-conditions"></a>
Classical Limit Conditions
</h3><p>[[#Three Distribution Relations]]
</p>
$$
e^\alpha=\frac{Z_1}{N}=\frac{V}{N}(\frac{2\pi mkT}{h^2})^{\frac{3}{2}}\gg 1
$$<ol>
<li>The smaller $N/V$, the more rarefied the gas;</li>
<li>The higher the temperature $T$;</li>
<li>The larger the molecular mass $m$
The easier it is to satisfy the classical limit conditions.
Another expression:

$$
n\lambda^3=e^{-\alpha}\ll 1
$$</li>
</ol>
<h3 id="maxwells-velocity-distribution-law">
<a class="header-anchor" href="#maxwells-velocity-distribution-law"></a>
Maxwell&rsquo;s Velocity Distribution Law
</h3><p>Under general conditions, <strong>gases satisfy the classical limit condition</strong>.</p>
<blockquote>
<p>Based on the <strong>Boltzmann distribution</strong>, the translational motion of the center of mass of gas molecules is studied, and the velocity distribution law of gas molecules is derived.</p>
</blockquote>
<h3 id="maxwell-velocity-distribution-law">
<a class="header-anchor" href="#maxwell-velocity-distribution-law"></a>
Maxwell Velocity Distribution Law
</h3>$$
f(v_x,x_y,x_z)dv_xdv_ydv_z=n(\frac{m}{2\pi kT})^{3/2}e^{-\frac{m}{2kT}(v_x^2+v_y^2+v_z^2)}dv_xdv_ydv_z
$$<p>
$f(v_x,x_y,x_z)$ should satisfy
</p>
$$
f(v_x,x_y,x_z)=\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}f(v_x,x_y,x_z)dv_xdv_ydv_z=n
$$<p>
$n=\frac{N}{V}$: number of molecules per unit volume</p>
<h3 id="maxwells-speed-distribution-law">
<a class="header-anchor" href="#maxwells-speed-distribution-law"></a>
Maxwell&rsquo;s Speed Distribution Law
</h3><p>Using the <strong>spherical polar coordinates</strong> volume element $v^2\sin \theta dv d \theta d \varphi$ to replace the Cartesian volume element $dv_xdv_ydv_z$, integrating over $\theta,\varphi$ gives the number of molecules per unit volume with speeds in the range $dv$ as:
</p>
$$
f(v)dv=4\pi n(\frac{m}{2\pi kT})^{3/2}e^{-\frac{m}{2kT}v^2}v^2dv
$$<p>
This equation is called the molecular speed distribution of gases.
</p>
$$
f(v)=4\pi n(\frac{m}{2\pi kT})^{3/2}\int_{0}^{\infty}e^{-\frac{m}{2kT}v^2}v^2dv=n
$$<h3 id="most-probable-speed">
<a class="header-anchor" href="#most-probable-speed"></a>
Most Probable Speed $v_p$
</h3><blockquote>
<p>The speed distribution function has a maximum value. The speed at which the speed distribution function reaches its maximum is called the most probable speed.
If the speed is divided into equal intervals, the interval containing $v_p$ has the highest number of molecules.</p>
</blockquote>
$$
\frac{d}{dv}(e^{-\frac{m}{2kT}v^2}v^2)=0
$$<p>
</p>
$$
v_p=\sqrt{\frac{2kT}{m}}
$$<p>
Average speed $\overline{v}$: the average value of speed $v$
</p>
$$
\overline{v}=4\pi(\frac{m}{2\pi kT})^{3/2}\int_{0}^{\infty}ve^{-\frac{m}{2kT}v^2}v^2dv=\sqrt{\frac{8kT}{\pi m}}
$$<p>
Root-mean-square speed $v_s$: the square root of the average value of $v^2$
</p>
$$
v_s^2=\overline{v_2}=4\pi(\frac{m}{2\pi kT})^{3/2}\int_{0}^{\infty}v^2e^{-\frac{m}{2kT}v^2}v^2dv=\frac{3kT}{m}
$$<p>
</p>
$$
v_s=\sqrt{\frac{3kT}{m}}=\sqrt{\frac{3RT}{M}}
$$<p>
$M$: molar mass</p>
<h3 id="wall-collision-number-effusion">
<a class="header-anchor" href="#wall-collision-number-effusion"></a>
Wall Collision Number (Effusion)
</h3><blockquote>
<p>The number of molecules colliding with a unit area of the container wall per unit time is called the wall collision number.
The number of particles escaping per unit time through a unit area hole in the container wall is called effusion.
Clearly: Wall collision number = Effusion</p>
</blockquote>
<p>$dA$ is an infinitesimal area element on the container wall, with its normal direction along the $x$-axis. Let $d\varGamma dA dt$ represent the number of molecules that collide with area $dA$ within time $dt$, with velocities in the range $dv_xdv_ydv_z$.
</p>
$$
d\varGamma dA dt=fv_xdv_xdv_ydv_zdAdt
$$<p>
</p>
$$
d\varGamma=fv_xdv_xdv_ydv_z
$$<p>
Integrating over velocities,
</p>
$$
\begin{cases}
v_x,0\to \infty\\
v_y,-\infty \to +\infty\\
v_z,-\infty \to +\infty
\end{cases}
$$<p>
we obtain
</p>
$$
\varGamma=n\sqrt{\frac{kT}{2\pi m}}=
\frac{1}{4}n\overline{v}
$$<h3 id="equipartition-theorem-classical-statistics">
<a class="header-anchor" href="#equipartition-theorem-classical-statistics"></a>
Equipartition Theorem (Classical Statistics)
</h3><blockquote>
<p>For a <strong>classical system</strong> in equilibrium at temperature $T$, the <strong>average value</strong> of each independent quadratic term in the particle energy equals $\frac{1}{2}kT$.</p>
</blockquote>
<h3 id="monatomic-ideal-gas">
<a class="header-anchor" href="#monatomic-ideal-gas"></a>
Monatomic Ideal Gas
</h3>$$
\varepsilon=\frac{1}{2m}(p_x^2+p_y^2+p_z^2)
$$<p>
The above equation has three squared terms. Applying the equipartition theorem gives:
</p>
$$
\overline{\varepsilon}=\frac{3}{2}kT
$$<p>
Thus, the internal energy of a monatomic <strong>ideal gas</strong> is:
</p>
$$
U=\frac{3}{2}NkT
$$<p>
</p>
$$
C_V=\frac{3}{2}Nk
$$<p>
Using the thermodynamic relation $C_p-C_V=Nk$, the isobaric heat capacity is obtained as:
</p>
$$
C_p=\frac{5}{2}Nk
$$<h3 id="adiabatic-index">
<a class="header-anchor" href="#adiabatic-index"></a>
Adiabatic Index
</h3>$$
\gamma=\frac{C_p}{C_V}=\frac{5}{3}=1.667
$$<p>
[[#Ideal Gas Adiabatic Equation]]</p>
<h3 id="diatomic-molecular-gas">
<a class="header-anchor" href="#diatomic-molecular-gas"></a>
Diatomic Molecular Gas
</h3><blockquote>
<p>The energy of a diatomic molecule has five quadratic terms (momentum in three directions of the molecule, the motion between the two atoms can be reduced to a rotor, two degrees of freedom, i.e., two quadratic terms).</p>
</blockquote>
$$
\overline{\varepsilon}=\frac{5}{2}kT
$$<p>
</p>
$$
U=\frac{5}{2}NkT
$$<p>
</p>
$$
C_V=\frac{5}{2}Nk
$$<p>
</p>
$$
C_p=\frac{7}{2}Nk
$$<p>
</p>
$$
\gamma=\frac{C_p}{C_V}=\frac{7}{5}=1.4
$$<h3 id="solids">
<a class="header-anchor" href="#solids"></a>
Solids
</h3><blockquote>
<p>Atoms in a solid can undergo small vibrations around their equilibrium positions.
Assuming each atom&rsquo;s vibration is an independent simple harmonic motion,
the atomic model of a solid simplifies to a <strong>three-dimensional linear oscillator</strong>.</p>
</blockquote>
<p>The energy of an atom in <strong>one degree of freedom</strong> is:
</p>
$$
\varepsilon=\frac{1}{2m}p^2+\frac{1}{2}m\omega^2q^2
$$<p>
with two quadratic terms.<br>
The average energy per atom is:
</p>
$$
\overline{\varepsilon}=3kT
$$<p>
The internal energy of the solid is:
</p>
$$
U=3NkT
$$<p>
The heat capacity at constant volume is:
</p>
$$
C_V=3Nk
$$<p>
In experiments, $C_p$ is usually measured. Using the thermodynamic relation:
</p>
$$
C_p-C_V=\frac{TV\alpha^2}{\kappa T}
$$<p>
$C_V$ can be derived.<br>
At room temperature and high temperatures, $C_V$ matches predictions, but <strong>at low temperatures, the heat capacity drops rapidly with decreasing temperature—a phenomenon unexplained by classical theory.</strong><br>
[[#Principle of Unattainability of Absolute Zero]]</p>
<h3 id="metal">
<a class="header-anchor" href="#metal"></a>
Metal
</h3><p>There are free electrons in metals. If the equipartition theorem is applied to the electrons, the heat capacity of the free electrons would be of the same order of magnitude as the heat capacity of the ionic vibrations.
Experimental results show that above $3K$, the heat capacity of free electrons is negligible compared to the heat capacity of ionic vibrations.</p>
<h3 id="quantum-statistics-of-localizednearly-independent-particle-systems-satisfying-classical-limit-conditions">
<a class="header-anchor" href="#quantum-statistics-of-localizednearly-independent-particle-systems-satisfying-classical-limit-conditions"></a>
Quantum Statistics of Localized/Nearly Independent Particle Systems Satisfying Classical Limit Conditions
</h3><h3 id="diatomic-ideal-gas">
<a class="header-anchor" href="#diatomic-ideal-gas"></a>
Diatomic Ideal Gas
</h3><p>An ideal gas is a ==non-localized system==. Since it satisfies the [[#Classical Limit Condition]], it can be discussed using the Boltzmann distribution.</p>
<blockquote>
<p>Under certain approximations, the motion of diatomic molecules can be divided into <strong>translation $t$, vibration $v$, and rotation $r$</strong>.</p>
</blockquote>
$$
\varepsilon=\varepsilon^t+\varepsilon^v+\varepsilon^r
$$<p>
</p>
$$
\omega=\omega_t+\omega_v+\omega r
$$<p>
</p>
$$
Z_1=Z_1^t\cdot Z_1^v \cdot Z_1^r
$$<p>
</p>
$$
U=-N\frac{\partial}{\partial \beta}\ln Z_1=-N\frac{\partial}{\partial \beta}(\ln Z_1^t+\ln Z_1^v+\ln Z_1^r)=U^t+U^v+U^r
$$<p>
</p>
$$
C_V=C_V^t+C_V^v+C_V^r
$$<h3 id="translational-energy">
<a class="header-anchor" href="#translational-energy"></a>
Translational Energy
</h3>$$
\begin{cases}
Z_1^t=V(\frac{2\pi m}{h^2\beta})^{3/2}\\
U^t=\frac{3N}{2\beta}=\frac{3}{2}NkT\\
C_V^t=\frac{3}{2}Nk
\end{cases}
$$<p>
Consistent with the results obtained from the equipartition theorem of classical statistics.</p>
<h3 id="vibrational-energy-and-characteristic-temperature">
<a class="header-anchor" href="#vibrational-energy-and-characteristic-temperature"></a>
Vibrational Energy and Characteristic Temperature
</h3>$$
\begin{cases}
Z_1^v=\frac{e^{-\frac{\beta\hbar\omega}{2}}}{1-e^{-\beta\hbar\omega}}\\
U^v=\frac{N\hbar\omega}{2}+\frac{N\hbar\omega}{e^{\beta\hbar\omega}-1}\\
C_V^v=Nk(\frac{\hbar \omega}{kT})^2\cdot \frac{e^{\hbar \omega /kT}}{(e^{\hbar \omega/kT}-1)^2}
\end{cases}
$$<p>$U^v$:</p>
<ul>
<li>The first term is temperature-independent and represents the ==zero-point energy== of $N$ oscillators.</li>
<li>The second term is the ==thermal excitation energy== of $N$ oscillators at temperature $T$.</li>
</ul>
<p>Introducing the ==vibrational characteristic temperature== $\theta_v$, which satisfies:
</p>
$$
k\theta_v=\hbar\omega
$$<p>
Substituting into the above system of equations yields:
</p>
$$
\begin{cases}
U_v=\frac{Nk \theta_v}{e^{\frac{\theta_v}{T}}-1}\\
C_V^v=Nk(\frac{\theta_v}{T})^2\frac{e^{\frac{\theta_v}{T}}}{(e^{\frac{\theta_v}{T}}-1)^2}
\end{cases}
$$<p>
$\because$ The ==characteristic temperature of diatomic molecules== is on the order of $10^3K$,
$\therefore$ In the room temperature range, $T\ll \theta_v$. Within this range, the contribution of vibrational degrees of freedom to heat capacity is nearly zero.<br>
At room temperature, the energy level spacing $\hbar \omega$ of diatomic molecules is much larger than $kT$. Due to the discrete energy levels, oscillators must acquire energy $\hbar \omega$ to transition to an excited state.</p>
<ul>
<li>In the case of $T \ll \theta_v$, the probability of an oscillator acquiring thermal energy $\hbar \omega$ and transitioning to an excited state is extremely low. Thus, on average, <strong>almost all oscillators remain frozen in the ground state</strong>.<br>
When the gas temperature rises, they hardly absorb any energy (unable to transition to higher energy levels). This explains why ==vibrational degrees of freedom do not participate in energy equipartition at room temperature==.</li>
<li>In the case of $T \gg \theta_v$, $C_V^v \to Nk$ (can be proven using L&rsquo;Hôpital&rsquo;s rule).</li>
</ul>
<h3 id="rotational-energy">
<a class="header-anchor" href="#rotational-energy"></a>
Rotational Energy
</h3><p>Rotational energy levels
</p>
$$
\varepsilon^r=\frac{l(l+1)\hbar^2}{2I},l=0,1,2,\cdots
$$<p>
Degeneracy
</p>
$$
\omega^r=2l+1
$$<p>
Partition function
</p>
$$
Z_1^r=\sum_{i=0}^\infty(2l+1)e^{-\beta\frac{l(l+1)\hbar^2}{2I}}
$$<p>
Characteristic rotational temperature
</p>
$$
\theta_r=\frac{\hbar^2}{2Ik} \lt 100K
$$<p>
</p>
$$
Z_1^r=\sum_{i=0}^{\infty}(2l+1)e^{-l(l+1)\frac{\theta_r}{T}}
$$<p>
==At room temperature== $\frac{\theta_r}{T}\ll 1$, integration yields
</p>
$$
Z_1^r=\frac{2I}{\beta \hbar^2}
$$<p>
From this we obtain
</p>
$$
\begin{cases}
U^r=NkT\\
C_V^r=Nk
\end{cases}
$$<p>Consistent with classical statistics.</p>
<h3 id="heteronuclear-diatomic-molecules">
<a class="header-anchor" href="#heteronuclear-diatomic-molecules"></a>
Heteronuclear Diatomic Molecules
</h3>$$
\begin{cases}
Z_1^t=V(\frac{2\pi m}{h_0^2\beta})^{3/2}\\
Z_1^v=\frac{2\pi}{h_0\beta \omega}\\
Z_1^r=\frac{8\pi^2I}{h_0^2\beta}
\end{cases}
$$<p>
</p>
$$
\begin{cases}
C_V^t=\frac{3}{2}Nk\\
C_V^v=Nk\\
C_V^r=Nk
\end{cases}
$$<h3 id="monatomic-ideal-gas-1">
<a class="header-anchor" href="#monatomic-ideal-gas-1"></a>
Monatomic Ideal Gas
</h3><h2 id="">
<a class="header-anchor" href="#"></a>

</h2><h3 id="entropy-1">
<a class="header-anchor" href="#entropy-1"></a>
Entropy
</h3><p>==Monatomic Ideal Gas==
</p>
$$
S=Nk(\ln Z_1-\beta\frac{\partial}{\partial \beta})-k \ln N!
$$<p>
Using the approximation $\ln N!=N(\ln N-1)$, it simplifies to
</p>
$$
S=\frac{3}{2}Nk \ln T+Nk \ln\frac{V}{N}+\frac{3}{2}Nk[\frac{5}{3}+\ln(\frac{2\pi mk}{h^2})]
$$<p>
The entropy in the above equation is an extensive quantity and represents absolute entropy.</p>
<h3 id="vapor-pressure-equation">
<a class="header-anchor" href="#vapor-pressure-equation"></a>
Vapor Pressure Equation
</h3><p>Sackur-Tetrode equation
</p>
$$
\ln p=-\frac{L}{RT}+\frac{5}{2}\ln T+\frac{5}{2}+\ln[k^5/2(\frac{2\pi m}{h^2})^3/2]
$$<h3 id="chemical-potential-1">
<a class="header-anchor" href="#chemical-potential-1"></a>
Chemical Potential
</h3>$$
\begin{align}
\mu &=(\frac{\partial F}{\partial N})_{T,V}\\
&=-kT \ln \frac{Z_1}{N}\\
&=kT \ln[\frac{N}{V}(\frac{h^2}{2\pi mkT})^3/2]
\end{align}
$$<p>
$\because$ For an ideal gas, $\frac{N}{V}(\frac{h^2}{2\pi mkT})^3/2\ll 1$<br>
$\therefore$ ==The chemical potential of an ideal gas is negative==.</p>
<h3 id="solids-1">
<a class="header-anchor" href="#solids-1"></a>
Solids
</h3><p>The thermal motion of atoms in a solid can be regarded as the vibration of $3N$ oscillators. <em>Einstein assumed that these $3N$ oscillators all have the same frequency.</em> Let $\omega$ denote the angular frequency of the oscillators, the energy levels of the oscillators are
</p>
$$
\varepsilon_n=\hbar \omega(n+\frac{1}{2}),n=0,1,2,\cdots
$$<p>
The partition function
</p>
$$
Z_1=\frac{e^{-\frac{\beta\hbar\omega}{2}}}{1-e^{-\beta\hbar\omega}}
$$<p>
Internal energy
</p>
$$
U=-3N \frac{\partial}{\partial \beta}\ln Z_1=3N \frac{\hbar \omega}{2}+3N\frac{\hbar \omega}{e^{\beta \hbar \omega}-1}
$$<p>
Heat capacity
</p>
$$
C_V=(\frac{\partial U}{\partial T})_V=3Nk(\frac{\hbar \omega}{kT})^2\frac{e^{\frac{\hbar \omega}{kT}}}{(e^{\frac{\hbar \omega}{kT}}-1)^2}
$$<p>
Introducing the ==Einstein characteristic temperature== $\theta_E=\frac{\hbar \omega}{k}$, the heat capacity can be rewritten as
</p>
$$
C_V=3Nk(\frac{\theta_E}{T})^2\frac{e^{\frac{\theta_E}{T}}}{(e^{\frac{\theta_E}{T}}-1)^2}
$$<ul>
<li>When $T \gg \theta_E$, the approximation $e^{\theta_E/T}-1\approx \theta_E/T$ can be made, yielding

$$
    C_V=3Nk
    $$
which is consistent with the result from the equipartition theorem.
At this point, the energy level spacing is much smaller than $kT$, and the quantization effect can be ignored.</li>
<li>When $T \ll \theta_E$, the approximation $e^{\frac{\theta_E}{T}}-1 \approx e^{\frac{\theta_E}{T}}$ holds,

$$
    C_V=3Nk(\frac{\theta_E}{T})^2e^{-\frac{\theta_E}{T}}
    $$
As $T \to 0$, $C_V \to 0$.
This proves the third law of thermodynamics[[#推论]].</li>
</ul>
<h3 id="bosefermi-systems">
<a class="header-anchor" href="#bosefermi-systems"></a>
Bose/Fermi Systems
</h3><p>[[#Bose Systems]], [[#Fermi Systems]].</p>
<h3 id="grand-partition-function-for-bosonic-systems">
<a class="header-anchor" href="#grand-partition-function-for-bosonic-systems"></a>
Grand Partition Function for Bosonic Systems
</h3>$$
\varXi=\prod_{l}\varXi_l=\prod_{l}(1-e^{-\alpha-\beta \varepsilon_l})^{-\omega_l}
$$<p>
</p>
$$
\ln \varXi=-\sum_i \omega_l \ln(1-e^{-\alpha-\beta \varepsilon_l})
$$<h3 id="average-particle-number">
<a class="header-anchor" href="#average-particle-number"></a>
Average Particle Number
</h3>$$
\overline{N}=\sum_l \frac{\omega_l}{e^{\alpha+\beta \varepsilon_l}-1}=-\frac{\partial}{\partial \alpha}\ln \varXi
$$<h3 id="internal-energy-1">
<a class="header-anchor" href="#internal-energy-1"></a>
Internal Energy
</h3>$$
U=\sum_l \frac{\varepsilon_l \omega_l}{e^{\alpha+\beta \varepsilon_l}-1}=-\frac{\partial}{\partial \beta}\ln \varXi
$$<h3 id="grand-partition-function-for-fermi-systems">
<a class="header-anchor" href="#grand-partition-function-for-fermi-systems"></a>
Grand Partition Function for Fermi Systems
</h3>$$
\varXi=\prod_{l}\varXi_l=\prod_{l}(1+e^{-\alpha-\beta \varepsilon_l})^{\omega_l}
$$<p>
</p>
$$
\ln \varXi=\sum_i \omega_l \ln(1+e^{-\alpha-\beta \varepsilon_l})
$$<h2 id="">
<a class="header-anchor" href="#"></a>

</h2><h3 id="average-particle-number-1">
<a class="header-anchor" href="#average-particle-number-1"></a>
Average Particle Number
</h3>$$
\overline{N}=\sum_l \frac{\omega_l}{e^{\alpha+\beta \varepsilon_l}+1}=-\frac{\partial}{\partial \alpha}\ln \varXi
$$<h3 id="internal-energy-2">
<a class="header-anchor" href="#internal-energy-2"></a>
Internal Energy
</h3>$$
U=\sum_l \frac{\varepsilon_l \omega_l}{e^{\alpha+\beta \varepsilon_l}+1}=-\frac{\partial}{\partial \beta}\ln \varXi
$$<h3 id="generalized-force">
<a class="header-anchor" href="#generalized-force"></a>
Generalized Force
</h3>$$
Y=-\frac{1}{\beta}\frac{\partial}{\partial y}\ln \varXi
$$<h3 id="equation-of-state-1">
<a class="header-anchor" href="#equation-of-state-1"></a>
Equation of State
</h3>$$
p=\frac{1}{\beta}\frac{\partial}{\partial V}\ln \varXi
$$<h3 id="entropy-2">
<a class="header-anchor" href="#entropy-2"></a>
Entropy
</h3>$$
\begin{align}
S &=k(\ln \varXi-\alpha \frac{\partial}{\partial \alpha}\ln \varXi-\beta \frac{\partial}{\partial \beta}\ln \varXi)\\
&=k(\ln \varXi+\alpha \overline{N}+\beta U)\\
&=k \ln \Omega
\end{align}
$$<h3 id="grand-potential-1">
<a class="header-anchor" href="#grand-potential-1"></a>
Grand Potential
</h3><p>[[#Grand Potential]]
</p>
$$
J=U-TS-\overline{N}\mu=-kT \ln \varXi
$$<h1 id="">
<a class="header-anchor" href="#"></a>

</h1><h3 id="lagrange-multiplier">
<a class="header-anchor" href="#lagrange-multiplier"></a>
Lagrange Multiplier
</h3>$$
\beta=\frac{1}{kT},\alpha=-\frac{\mu}{kT}
$$<hr>
<h2 id="tags">
<a class="header-anchor" href="#tags"></a>
title: Bose-Einstein Condensation
description:
keywords:
tags:
</h2><h3 id="system-properties">
<a class="header-anchor" href="#system-properties"></a>
System Properties
</h3><blockquote>
<p>Consider a system composed of $N$ identical, nearly independent bosons, with temperature $T$ and volume $V$. For clarity, assume the particles have zero spin.</p>
</blockquote>
<p>According to the Bose distribution, the number of particles in each energy level is
</p>
$$
a_l=\frac{\omega_l}{e^{\alpha+\beta \varepsilon_l}-1}=\frac{\omega_l}{e^{\frac{\varepsilon_l-\mu}{kT}}-1}
$$<p>
$\because$ The number of particles in an energy level cannot be negative,
$\therefore$ Let $\varepsilon_0$ denote the lowest energy level, $\varepsilon_0 \gt \mu$.
If the lowest energy level is 0, then $\mu \lt 0$.</p>
<h3 id="chemical-potential-2">
<a class="header-anchor" href="#chemical-potential-2"></a>
Chemical Potential
</h3><p>The chemical potential of a Bose system is less than 0. The chemical potential increases as the temperature decreases. When the temperature drops to a certain critical temperature $T_C$, $\mu \to 0$.
$T_C$ is determined by the following equation:
</p>
$$
\frac{2\pi}{h^3}(2m)^{3/2}\int_0^{\infty}\frac{\varepsilon^{1/2d \varepsilon}}{e^{\frac{\varepsilon}{kT_C}}-1}=n
$$<p>
</p>
$$
T_C=\frac{2\pi}{2.612^{2/3}}\frac{\hbar^2}{mk}n^{3/2}
$$<h3 id="condensation-phenomenon">
<a class="header-anchor" href="#condensation-phenomenon"></a>
Condensation Phenomenon
</h3><p>When $T \le T_C$, the ==particle number density== accumulated at the lowest energy level $\varepsilon=0$ is given by
</p>
$$
n_0(T)=n[1-(\frac{T}{T_C})^{3/2}]
$$<p>
At <strong>absolute zero</strong>, particles will occupy the lowest energy states as much as possible. For <strong>bosons</strong>, the number of particles that can occupy a single quantum state is unrestricted. Therefore, at absolute zero, all bosonic particles will be in the lowest energy level $\varepsilon=0$.</p>
<p>The above equation shows that:</p>
<blockquote>
<p>When $T  \lt  T_C$, a <strong>macroscopic number</strong> of particles condense into the lowest energy level. This phenomenon is called ==Bose-Einstein condensation==.</p>
</blockquote>
$$
U=0.770NkT(\frac{T}{T_C})^{3/2}
$$<p>
</p>
$$
C_V=(\frac{\partial U}{\partial T})_V=\frac{5}{2}\frac{U}{T}=1.925Nk(\frac{T}{T_C})^{3/2}
$$<p>
At $T=T_C$, $C_V$ reaches its maximum value of $1.925Nk$.<br>
When the <strong>temperature is sufficiently high</strong>,
</p>
$$
C_V=\frac{3}{2}Nk
$$<h3 id="critical-condition-for-bose-einstein-condensation-in-an-ideal-bose-gas">
<a class="header-anchor" href="#critical-condition-for-bose-einstein-condensation-in-an-ideal-bose-gas"></a>
Critical Condition for Bose-Einstein Condensation in an Ideal Bose Gas
</h3>$$
n(\frac{h}{\sqrt{2\pi mkT_C}})^3=n \lambda^3 \ge 2.612
$$<p>
$\lambda$: thermal de Broglie wavelength of atoms</p>
<h3 id="photon-gas-bose-gas">
<a class="header-anchor" href="#photon-gas-bose-gas"></a>
Photon Gas (Bose Gas)
</h3><blockquote>
<p>From the particle perspective, the radiation field inside an empty cavity can be regarded as a <strong>photon gas</strong>.
The radiation field inside the empty cavity can be decomposed into the superposition of infinitely many <strong>monochromatic plane waves</strong>.</p>
</blockquote>
<p>Based on the [[#de Broglie relations]], and the photon $\omega=ck$, we obtain the ==energy-momentum relation for photons==
</p>
$$
\varepsilon=cp
$$<ul>
<li>Photons are bosons with a spin quantum number of 1, and they obey Bose-Einstein statistics when in equilibrium.</li>
<li>Since the cavity walls continuously emit and absorb photons, the number of photons in the photon gas is not conserved. When deriving the Bose distribution, only the condition that $E$ is constant exists, and not the condition that $N$ is constant. Therefore, only one Lagrange multiplier $\beta$ should be introduced, with $\alpha=0$.</li>
</ul>
<h3 id="photon-gas-statistical-distribution">
<a class="header-anchor" href="#photon-gas-statistical-distribution"></a>
Photon Gas Statistical Distribution
</h3>$$
a_l=\frac{\omega_l}{e^{\beta \omega_l}-1}
$$<h3 id="thermodynamic-potential">
<a class="header-anchor" href="#thermodynamic-potential"></a>
Thermodynamic Potential
</h3><p>From
</p>
$$
\alpha=-\frac{\mu}{kT}=0
$$<p>
we obtain the thermodynamic potential of the photon gas as $\mu=0$.</p>
<h3 id="number-of-quantum-states">
<a class="header-anchor" href="#number-of-quantum-states"></a>
Number of Quantum States
</h3><p>The spin quantum number of a photon is 1, and its spin projection along the momentum direction can take two possible values, $\pm \hbar$, corresponding to left and right circular polarization. Considering that the photon spin has 2 projections, as shown in the [[#Energy Representation]] section of [[#Number of Quantum States for a Three-Dimensional Free Particle]],</p>
<blockquote>
<p>Within a cavity of volume $V$, in the momentum range from $p$ to $p+dp$,</p>
</blockquote>
<p>the number of quantum states for photons is
</p>
$$
\Omega=\frac{8\pi V}{h^3}p^2dp
$$<blockquote>
<p>Within a cavity of volume $V$, in the angular frequency range from $\omega$ to $\omega+d \omega$,</p>
</blockquote>
<p>the number of quantum states for photons is
</p>
$$
\Omega=\frac{V}{\pi^2c^3}\omega^2d \omega
$$<h3 id="energy-distribution-function">
<a class="header-anchor" href="#energy-distribution-function"></a>
Energy Distribution Function
</h3><blockquote>
<p>In a cavity of volume $V$, within the frequency range from $\omega$ to $\omega+d \omega$,</p>
</blockquote>
<p>The average number of photons
</p>
$$
\overline{N}=\frac{V}{\pi^2c^3}\frac{\omega^2d \omega}{e^{\hbar \omega/kT}-1}
$$<p>
The energy of the radiation field
</p>
$$
U(\omega,T)d \omega=\overline{N}\hbar \omega=\frac{V}{\pi^2c^3}\frac{\hbar \omega^3}{e^{\hbar \omega/kT}-1}d \omega
$$<h2 id="">
<a class="header-anchor" href="#"></a>

</h2><h3 id="rayleigh-jeans-law">
<a class="header-anchor" href="#rayleigh-jeans-law"></a>
Rayleigh-Jeans Law
</h3><p>In the low-frequency range where $\frac{\hbar \omega}{kT}\ll 1$, $e^{\frac{\hbar \omega}{kT}}\approx 1+\frac{\hbar \omega}{kT}$, the above equation can be approximated as
</p>
$$
U(\omega,T)d \omega=\frac{V}{\pi^2c^3}\omega^2kTd \omega
$$<p>
This equation is known as the ==Rayleigh-Jeans Law==, which <strong>fits well at low frequencies</strong> but leads to the ultraviolet catastrophe at high frequencies.</p>
<h3 id="wiens-law">
<a class="header-anchor" href="#wiens-law"></a>
Wien&rsquo;s Law
</h3><p>In the high-frequency range where $\frac{\hbar \omega}{kT}\gg 1$, $e^{\frac{\hbar \omega}{kT}}-1\approx e^{\frac{\hbar \omega}{kT}}$, the above equation can be approximated as
</p>
$$
U(\omega,T)d \omega=\frac{V}{\pi^2c^3}\hbar \omega^3e^{-\frac{\hbar \omega}{kT}}d \omega
$$<p>
This formula matches Wien&rsquo;s Law, providing good agreement at high frequencies but poor agreement at low frequencies.</p>
<blockquote>
<p>Inside a cavity of volume $V$</p>
</blockquote>
<p>The internal energy of equilibrium radiation is
</p>
$$
U=\frac{\pi^2k^4}{15c^3\hbar^3}VT^4
$$<h3 id="grand-partition-function">
<a class="header-anchor" href="#grand-partition-function"></a>
Grand Partition Function
</h3>$$
\ln \varXi=\frac{pi^2V}{45c^3}\frac{1}{(\beta \hbar)^3}
$$<h3 id="photon-gas-internal-energy">
<a class="header-anchor" href="#photon-gas-internal-energy"></a>
Photon Gas Internal Energy
</h3>$$
U=-\frac{\partial}{\partial \beta}\ln \varXi=\frac{\pi^2k^4V}{15c^3\hbar^3}T^4
$$<p>
Consistent with [[#Internal Energy Distribution Function]]</p>
<h3 id="pressure">
<a class="header-anchor" href="#pressure"></a>
Pressure
</h3>$$
p=-\frac{1}{\beta}\frac{\partial}{\partial V}\ln \varXi=\frac{\pi^2k^4}{45c^3\hbar^3}T^4
$$<p>
Comparing the above two equations gives
</p>
$$
p=\frac{1}{3}\frac{U}{V}
$$<h3 id="entropy-3">
<a class="header-anchor" href="#entropy-3"></a>
Entropy
</h3>$$
S=k(\ln \varXi+\beta U)=\frac{4}{45}\frac{\pi^2k^4}{c^3\hbar^3}T^3V
$$<p>
When $T \to 0$, $S \to 0$, consistent with [[#Third Law of Thermodynamics]]</p>
<h3 id="radiation-flux-density">
<a class="header-anchor" href="#radiation-flux-density"></a>
Radiation Flux Density
</h3><p>The formula for [[#Radiation Flux Density]] derived from [[#Radiation]] is
</p>
$$
J_u=\frac{c}{4}\frac{U}{V}=\frac{1}{4}cu
$$<p>
Based on the concept of [[#Wall Collision Rate (Effusion)]], the radiation flux density of a photon gas is obtained as
</p>
$$
J_u=\frac{\pi^2k^4}{60c^2\hbar^3}T^4
$$<h3 id="free-electron-gas-in-metals-fermi-gas">
<a class="header-anchor" href="#free-electron-gas-in-metals-fermi-gas"></a>
Free Electron Gas in Metals (Fermi Gas)
</h3><blockquote>
<p>The free electrons in metals form a highly degenerate Fermi gas.</p>
</blockquote>
<h3 id="average-electron-number">
<a class="header-anchor" href="#average-electron-number"></a>
Average Electron Number
</h3><p>According to the Fermi distribution, at temperature $T$, the average number of electrons occupying a quantum state with energy $\varepsilon$ is:
</p>
$$
f=\frac{1}{e^{\frac{\varepsilon-\mu}{kT}}+1}
$$<p>
The electron spin has 2 possible projections along its momentum direction. Based on the [[#Energy Representation]] in [[#Three-Dimensional Free Particle Quantum State Number]],</p>
<blockquote>
<p>Within volume $V$, the number of electron quantum states in the energy range from $\varepsilon$ to $\varepsilon+d \varepsilon$ is:</p>
</blockquote>
$$
D(\varepsilon)d\varepsilon=\frac{4\pi V}{h^3}(2m)^{3/2}\varepsilon^{1/2}d\varepsilon
$$<p>
Thus, the average number of electrons in this range is:
</p>
$$
dN=f \cdot D(\varepsilon)d\varepsilon=\frac{4\pi V}{h^3}(2m)^{3/2}\frac{\varepsilon^{1/2}d\varepsilon}{e^{\frac{\varepsilon-\mu}{kT}}+1}
$$<h3 id="total-number-of-electrons">
<a class="header-anchor" href="#total-number-of-electrons"></a>
Total Number of Electrons
</h3><p>Integrating over $d \varepsilon$ gives
</p>
$$
\frac{4\pi V}{h^3}(2m)^{3/2}\int_0^{\infty}\frac{\varepsilon^{1/2}d\varepsilon}{e^{\frac{\varepsilon-\mu}{kT}}+1}=N
$$<p>
Given the number of electrons $N$, temperature $T$, and volume $V$, the chemical potential $\mu$ can be determined from the above equation.</p>
<h3 id="system-properties-at">
<a class="header-anchor" href="#system-properties-at"></a>
System Properties at $T=0$
</h3><h3 id="average-electron-number-per-quantum-state">
<a class="header-anchor" href="#average-electron-number-per-quantum-state"></a>
Average Electron Number per Quantum State
</h3>$$
\begin{cases}
f=1,\varepsilon \lt \mu(0)\\
f=0,\varepsilon \gt \mu(0)
\end{cases}
$$<h3 id="fermi-level">
<a class="header-anchor" href="#fermi-level"></a>
Fermi Level
</h3><blockquote>
<p>The chemical potential at zero temperature is also known as the Fermi level, $\mu(0)$, which is the chemical potential when $T=0$.</p>
</blockquote>
<p>From
</p>
$$
\frac{4\pi V}{h^3}(2m)^{3/2}\int_0^{\mu(0)}\varepsilon^{1/2}d \varepsilon=N
$$<p>
we obtain
</p>
$$
\mu(0)=\frac{\hbar^2}{2m}(3\pi^2\frac{N}{V})^{2/3}
$$<h3 id="fermi-momentum">
<a class="header-anchor" href="#fermi-momentum"></a>
Fermi Momentum
</h3><p>Let
</p>
$$
\mu(0)=\frac{p_F^2}{2m}
$$<p>
we obtain
</p>
$$
p_F=(3\pi^2n)^{1/3}\hbar
$$<p>
This momentum is called ==Fermi momentum==.</p>
<h3 id="fermi-velocity">
<a class="header-anchor" href="#fermi-velocity"></a>
Fermi Velocity
</h3>$$
v_F=\frac{p_F}{m}
$$<p>
This velocity is called the ==Fermi velocity==.</p>
<h3 id="fermi-temperature">
<a class="header-anchor" href="#fermi-temperature"></a>
Fermi Temperature
</h3>$$
kT_F=\mu(0)
$$<p>
</p>
$$
T_F=\frac{\mu(0)}{k}
$$<p>
This temperature is called the ==Fermi temperature==.</p>
<h3 id="energy-of-electron-gas">
<a class="header-anchor" href="#energy-of-electron-gas"></a>
Energy of Electron Gas
</h3>$$
U(0)=\frac{4\pi V}{h^3}(2m)^{3/2}\int_0^{\mu(0)}\varepsilon^{3/2}d \varepsilon=\frac{3N}{5}\mu(0)
$$<h2 id="">
<a class="header-anchor" href="#"></a>

</h2><h3 id="electron-average-energy">
<a class="header-anchor" href="#electron-average-energy"></a>
Electron Average Energy
</h3>$$
u(0)=\frac{U(0)}{N}=\frac{3}{5}\mu(0)
$$<h2 id="">
<a class="header-anchor" href="#"></a>

</h2><h3 id="electron-gas-pressure">
<a class="header-anchor" href="#electron-gas-pressure"></a>
Electron Gas Pressure
</h3>$$
p(0)=\frac{2}{3}\frac{U(0)}{V}=\frac{2}{5}n \mu(0)
$$<h3 id="comparison-of-bose-and-fermi-gases-at">
<a class="header-anchor" href="#comparison-of-bose-and-fermi-gases-at"></a>
Comparison of Bose and Fermi Gases at $T=0$
</h3><p>In stark contrast to an ideal Bose gas, where all particles occupy the zero-energy, zero-momentum ground state with zero pressure at absolute zero temperature, a Fermi gas possesses extremely high average energy and momentum, generating substantial pressure at absolute zero.</p>

        
        <hr><p>Published on 2024-06-06 at <a href='https://www.guzhengsvt.cn/'>孤筝の温暖小家</a>, last modified on 2024-06-06</p><p>All articles on this blog are licensed under the BY-NC-SA license agreement unless otherwise stated. Please indicate the source when reprinting!</p>]]></description><category>Physics</category></item></channel></rss>