对称性可以理解为一种算符，当作用于系统时，使其基本特征保持不变。例如，晶格具有平移对称性；沿晶格矢量移动其位置并不会改变其结构。在量子力学中，这一概念与能级的简并密切相关。简并意味着不同的量子态可以具有相同的能量。虽然对称性有时意味着简并，但能级的简并总是系统中某种潜在对称性的标志。 让我们从旋转对称性出发。设想一个在圆周上运动的粒子，这是一个在旋转下不变的系统。其状态由波函数 \(\psi(\theta)\) 描述。如果我们将系统逆时针旋转一个无穷小角度 \(\epsilon\)，波函数变换为 \(\psi(\theta) \rightarrow \psi(\theta - \epsilon)\)。波函数的变化 \(\Delta\psi\) 可写为： \[ \Delta\psi \approx -\epsilon \frac{\partial\psi}{\partial\theta} = -\frac{i\epsilon}{\hbar} \left(-i\hbar\frac{\partial\psi}{\partial\theta}\right) \] 注意到角动量算符是 \(L_z = -i\hbar\frac{\partial}{\partial\theta}\)，于是有： \[ \Delta\psi = -\frac{i\epsilon}{\hbar}L_z\psi \] 在此语境下，\(L_z\) 是旋转的生成元。将其作用于波函数，刻画一次微小旋转的效果。同样的原理适用于线性动量，其中 \(p_x = -i\hbar\frac{\partial}{\partial x}\) 是沿 x 轴平移的生成元。需注意，为了简化，计算中有时将 \(\hbar\) 取为 1，但这里将其显式保留。 角动量算符的本征值和本征向量是什么？本征方程为 \(L_z|\psi\rangle = m|\psi\rangle\)，其中 \(m\) 是本征值。这给出微分方程： \[ -i\hbar\frac{\partial\psi}{\partial\theta} = m\psi(\theta) \] 其解具有形式 \(\psi(\theta) = C e^{im\theta/\hbar}\)。一个物理约束是波函数必须是单值的，即 \(\psi(\theta) = \psi(\theta + 2\pi)\)。这意味着经过一周旋转后，波函数必须回到原状。这导致条件： \[ e^{im2\pi/\hbar} = 1 \] 当且仅当 \(m/\hbar\) 为整数时成立。传统上，我们将该整数记作 \(m\)，于是角动量被量子化，其允许取值为 \(L_z = m\hbar\)。 ...

在量子力学中，物理系统的状态由一个态矢量描述，使用狄拉克的括号记号（bra-ket）记作 \(|\psi\rangle\)。这些态矢量是称为希尔伯特空间的复向量空间中的元素。 给定两个态矢量 \(|\psi\rangle\) 和 \(|\phi\rangle\)，它们的内积是一个复数，记为 \(\langle\phi|\psi\rangle\)。如果两个不同状态的内积为零，即 \(\langle\phi|\psi\rangle = 0\)，则称这些状态正交，表示完全独立的物理情形。 可以测量的物理量，如位置、动量和能量，称为可观测量。在量子力学中，可观测量由厄米算符表示。算符 \(A\) 若等于其自身的共轭转置 \(A = A^\dagger\)，则为厄米的。厄米算符的一个关键性质是其本征值总是实数，这是必要的，因为它们对应物理测量的结果。 当测量一个可观测量 \(A\) 时，唯一可能的结果是该算符的本征值。本征方程为： $$ A |\alpha\rangle = \alpha |\alpha\rangle $$ 其中，\(\alpha\) 是一个本征值，\(|\alpha\rangle\) 是对应的本征矢。如果我们对 \(A\) 进行测量并得到结果 \(\alpha\)，系统的状态将坍缩到相应的本征态 \(|\alpha\rangle\)。 为了描述空间中的粒子，我们常在位置基底中工作。基底矢量记作 \(|x\rangle\)，表示具有确定位置 \(x\) 的状态。态矢量 \(|\psi\rangle\) 在该基底中的表象是波函数 \(\psi(x)\)，其由内积定义： $$ \psi(x) = \langle x|\psi\rangle $$ 记号 \(\langle x|\psi\rangle\) 可能令人困惑。将其与熟悉的离散向量空间（如三维欧几里得空间）类比会有所帮助。 在离散基底中： 设想一个标准的三维向量 \(\vec{V}\)。我们可以在一个基底中表示它，该基底是一组三个正交的单位矢量 \(B = \{\hat{e}_1, \hat{e}_2, \hat{e}_3\}\)（例如 \(\hat{i}, \hat{j}, \hat{k}\)）。 基底是完整集合 \(B\)。 基底矢量是该集合中的一个成员，例如 \(\hat{e}_1\)。其中的 1 是一个离散指标。 要找到 \(\vec{V}\) 沿 \(\hat{e}_1\) 方向的分量，我们取点积（即内积）：\(V_1 = \hat{e}_1 \cdot \vec{V}\)。 该向量可以通过对其分量求和来重建：\(\vec{V} = V_1 \hat{e}_1 + V_2 \hat{e}_2 + V_3 \hat{e}_3 = \sum_{i=1}^3 (\hat{e}_i \cdot \vec{V}) \hat{e}_i\)。 在连续的位置基底中： ...

遵循 Susskind 的 Statistical Mechanics 课程，第 9 讲。 伊辛模型中单个自旋的能量为： \[E = -J\sigma\]其中： \(E\) 是该自旋的能量 \(J\) 是耦合常数（相互作用强度） \(\sigma\) 是自旋取值，可以是 +1（自旋向上）或 -1（自旋向下） 聚焦于单个自旋。 它的配分函数 \(Z\)，对两个构型求和（\(\sigma = +1, \sigma = -1\)），为： \[ Z = \sum_{\sigma \in \{+1, -1\}} e^{-\beta(-J\sigma)} = \sum_{\sigma \in \{+1, -1\}} e^{\beta J\sigma} = e^{\beta J(+1)} + e^{\beta J(-1)} = e^{\beta J} + e^{-\beta J} = 2 \cosh(\beta J) \]（对于由 \(N\) 个独立、无相互作用的自旋组成的系统，总配分函数的对数是可加的：\(\ln Z_{total} = N \ln Z\)。） 让我们计算该自旋的平均能量 \(\langle E \rangle\)： 平均能量由 \(\langle E \rangle = -\frac{1}{Z} \frac{\partial Z}{\partial \beta}\) 给出。 由于 \(Z = 2 \cosh(\beta J)\)，它对 \(\beta\) 的导数为 \(\frac{\partial Z}{\partial \beta} = 2J \sinh(\beta J)\)。 将这些代入关于 \(\langle E \rangle\) 的表达式中： ...

在这篇文章中，我想介绍现代宇宙学的基本假设，即FLRW模型，同时也作为爱因斯坦场方程的一个应用。 思路很简单：给定一个FLRW度规，代入场方程，然后我们就能推导出尺度因子\(a\)（它是度规的一个参数）如何随\(t\)坐标变化。这里我用一个Python程序来表达这些计算。 import sympy as sp from sympy.diffgeom import Manifold, Patch, CoordSystem, TensorProduct from sympy.diffgeom.diffgeom import metric_to_Ricci_components, metric_to_Christoffel_2nd from sympy import Symbol, Function, sin, simplify, Matrix # Define spacetime dimensions dim = 4 # Define manifold and coordinate system (using spherical coordinates) M = Manifold('M', dim) patch = Patch('P', M) coords = CoordSystem('coords', patch, [ Symbol('t', real=True), Symbol('r', real=True, positive=True), Symbol('theta', real=True, positive=True), Symbol('phi', real=True) ]) # Get coordinate functions and basic one-forms t, r, theta, phi = coords.coord_functions() dt, dr, dtheta, dphi = coords.base_oneforms() # Define scale factor and curvature parameter (using natural units c=1) a = Function('a')(t) # Scale factor as a function of time k = Symbol('k') # Spatial curvature parameter # Construct FLRW metric # ds² = -dt² + a(t)²[dr²/(1-kr²) + r²(dθ² + sin²θ dφ²)] g = -TensorProduct(dt, dt) + \ a**2 * (TensorProduct(dr, dr)/(1-k*r**2) + \ r**2 * TensorProduct(dtheta, dtheta) + \ r**2 * sin(theta)**2 * TensorProduct(dphi, dphi)) # Matrix form of the metric (for Einstein tensor calculation) g_matrix = Matrix([ [-1, 0, 0, 0], [0, a**2/(1-k*r**2), 0, 0], [0, 0, a**2*r**2, 0], [0, 0, 0, a**2*r**2*sin(theta)**2] ]) # Calculate Christoffel symbols Christoffel = metric_to_Christoffel_2nd(g) # Calculate Ricci tensor Ricci_tensor = metric_to_Ricci_components(g) # Define coordinate index names coord_names = ['t', 'r', 'θ', 'φ'] # Calculate inverse metric tensor g_inverse = Matrix([ [-1, 0, 0, 0], [0, (1-k*r**2)/a**2, 0, 0], [0, 0, 1/(a**2*r**2), 0], [0, 0, 0, 1/(a**2*r**2*sin(theta)**2)] ]) # Calculate Ricci scalar R = g^{μν} * R_{μν} Ricci_scalar = 0 for i in range(dim): for j in range(dim): Ricci_scalar += g_inverse[i, j] * Ricci_tensor[i, j] Ricci_scalar = simplify(Ricci_scalar) # Calculate Einstein tensor G_μν = R_μν - (1/2) * g_μν * R Einstein_tensor = Matrix([[0 for _ in range(dim)] for _ in range(dim)]) for i in range(dim): for j in range(dim): Einstein_tensor[i, j] = Ricci_tensor[i, j] - (1/2) * g_matrix[i, j] * Ricci_scalar Einstein_tensor[i, j] = simplify(Einstein_tensor[i, j]) # Print Christoffel symbols print("Christoffel symbols of FLRW metric (Γ^μ_νρ):") for i in range(dim): for j in range(dim): for k in range(dim): if Christoffel[i, j, k] != 0: print(f"Γ^{coord_names[i]}_{coord_names[j]}{coord_names[k]} = {simplify(Christoffel[i, j, k])}") # Print Ricci tensor print("\nRicci tensor of FLRW metric (R_μν):") for i in range(dim): for j in range(dim): if Ricci_tensor[i, j] != 0: print(f"R_{coord_names[i]}{coord_names[j]} = {simplify(Ricci_tensor[i, j])}") # Print Ricci scalar print("\nRicci scalar of FLRW metric (R):") print(f"R = {Ricci_scalar}") # Print Einstein tensor print("\nEinstein tensor of FLRW metric (G_μν):") for i in range(dim): for j in range(dim): if Einstein_tensor[i, j] != 0: print(f"G_{coord_names[i]}{coord_names[j]} = {Einstein_tensor[i, j]}") Christoffel symbols of FLRW metric (Γ^μ_νρ): Γ^t_rr = -a(t)*Subs(Derivative(a(_xi), _xi), _xi, t)/(k*r**2 - 1) Γ^t_θθ = a(t)*r**2*Subs(Derivative(a(_xi), _xi), _xi, t) Γ^t_φφ = a(t)*sin(theta)**2*r**2*Subs(Derivative(a(_xi), _xi), _xi, t) Γ^r_tr = Subs(Derivative(a(_xi), _xi), _xi, t)/a(t) Γ^r_rt = Subs(Derivative(a(_xi), _xi), _xi, t)/a(t) Γ^r_rr = -k*r/(k*r**2 - 1) Γ^r_θθ = k*r**3 - r Γ^r_φφ = (k*r**2 - 1)*sin(theta)**2*r Γ^θ_tθ = Subs(Derivative(a(_xi), _xi), _xi, t)/a(t) Γ^θ_rθ = 1/r Γ^θ_θt = Subs(Derivative(a(_xi), _xi), _xi, t)/a(t) Γ^θ_θr = 1/r Γ^θ_φφ = -sin(2*theta)/2 Γ^φ_tφ = Subs(Derivative(a(_xi), _xi), _xi, t)/a(t) Γ^φ_rφ = 1/r Γ^φ_θφ = 1/tan(theta) Γ^φ_φt = Subs(Derivative(a(_xi), _xi), _xi, t)/a(t) Γ^φ_φr = 1/r Γ^φ_φθ = 1/tan(theta) Ricci tensor of FLRW metric (R_μν): R_tt = -3*Subs(Derivative(a(_xi), (_xi, 2)), _xi, t)/a(t) R_rr = (-2*k - a(t)*Subs(Derivative(a(_xi), (_xi, 2)), _xi, t) - 2*Subs(Derivative(a(_xi), _xi), _xi, t)**2)/(k*r**2 - 1) R_θθ = (2*k + a(t)*Subs(Derivative(a(_xi), (_xi, 2)), _xi, t) + 2*Subs(Derivative(a(_xi), _xi), _xi, t)**2)*r**2 R_φφ = (2*k + a(t)*Subs(Derivative(a(_xi), (_xi, 2)), _xi, t) + 2*Subs(Derivative(a(_xi), _xi), _xi, t)**2)*sin(theta)**2*r**2 Ricci scalar of FLRW metric (R): R = 6*(k + a(t)*Subs(Derivative(a(_xi), (_xi, 2)), _xi, t) + Subs(Derivative(a(_xi), _xi), _xi, t)**2)/a(t)**2 Einstein tensor of FLRW metric (G_μν): G_tt = 3.0*(k + Subs(Derivative(a(_xi), _xi), _xi, t)**2)/a(t)**2 G_rr = (1.0*k + 2.0*a(t)*Subs(Derivative(a(_xi), (_xi, 2)), _xi, t) + 1.0*Subs(Derivative(a(_xi), _xi), _xi, t)**2)/(k*r**2 - 1) G_θθ = (-1.0*k - 2.0*a(t)*Subs(Derivative(a(_xi), (_xi, 2)), _xi, t) - 1.0*Subs(Derivative(a(_xi), _xi), _xi, t)**2)*r**2 G_φφ = (-1.0*k - 2.0*a(t)*Subs(Derivative(a(_xi), (_xi, 2)), _xi, t) - 1.0*Subs(Derivative(a(_xi), _xi), _xi, t)**2)*sin(theta)**2*r**2 将最终的爱因斯坦张量代入爱因斯坦方程。以\( G_{tt} \)为例： ...

为了推导爱因斯坦场方程，我们从涉及能量与动量的物理考量出发，并结合来自时空曲率的几何考量。 爱因斯坦场方程依赖于三个基本假设： 能量守恒 动量守恒 牛顿引力方程 对于假设1和2，能量-动量张量的守恒律成立： \[ \nabla_\mu T^{\mu \nu} = 0 \]能量-动量张量 \(T^{\mu \nu}\) 概括了关键的物理量，其表示为： \[ T^{\mu \nu} = \begin{pmatrix} T^{00} & T^{01} & T^{02} & T^{03} \\ T^{10} & T^{11} & T^{12} & T^{13} \\ T^{20} & T^{21} & T^{22} & T^{23} \\ T^{30} & T^{31} & T^{32} & T^{33} \end{pmatrix} \]其中： \(T^{00} = \rho\)：能量密度 \(T^{i0} = P_i\)：在第\(i\)个空间方向上的动量密度 \(T^{0i} = S_i\)：在第\(i\)个空间方向的能流 \(\sigma_{ij}\)：空间应力，与压力和应力相关 一个关键特征是该张量的对称性 \(T^{\mu \nu} = T^{\nu \mu}\)，这表明能流 \(S_i\) 等于动量密度 \(P_i\)。尽管这看起来似乎违反直觉，但它自然源于爱因斯坦的质能等价： ...

在这篇文章中，我希望学习关于\( \nabla \)符号的一些基础知识，以及一个与黎曼曲率张量相关的重要公式。 I. 联络 定义：光滑流形\( (M, \mathcal{O}) \)上的联络\( \nabla \)是一个映射，它把由向量场\( X \)和\( (p, q) \)-张量场\( T \)组成的有序对映射为\( (p, q) \)-张量场\( \nabla_X T \)，满足： \( \nabla_X f = X f \), \(\forall f \in C^\infty(M) \) \( \nabla_X (T + S) = \nabla_X T + \nabla_X S \) \( \nabla_X (T(\omega, Y)) = (\nabla_X T)(\omega, Y) + T(\nabla_X \omega, Y) + T(\omega, \nabla_X Y) \) （对于\( (1, 1) \)-张量场\( T \)，但对任意\( (p, q) \)-张量场\( T \)同理）。 （“莱布尼茨”） ...

在本文中，我们将更深入探讨微分几何的概念，主要聚焦于黎曼张量、里奇张量以及标量曲率的概念。 一、黎曼张量 A. 定义 我们从一个\(d\)维流形\(M\)出发，其上赋予了（伪）黎曼度规\(g_{\mu\nu}\)。与该度规相容的联络为Levi-Civita联络\(\nabla\)。黎曼曲率张量（常简称为黎曼张量）\(R^\rho_{\ \sigma\mu\nu}\)定义为作用在向量场\(V^\rho\)上的协变导数的对易子： \[ \bigl[\nabla_\mu, \nabla_\nu\bigr] V^\rho \;=\; R^\rho_{\ \sigma\mu\nu} \, V^\sigma. \] 用指标记号，另一种常见的写法是： \[ R^\rho_{\ \sigma\mu\nu} \;=\; \partial_\mu \Gamma^\rho_{\nu\sigma} \;-\; \partial_\nu \Gamma^\rho_{\mu\sigma} \;+\; \Gamma^\rho_{\mu\lambda}\,\Gamma^\lambda_{\nu\sigma} \;-\; \Gamma^\rho_{\nu\lambda}\,\Gamma^\lambda_{\mu\sigma}, \] 其中\(\Gamma^\rho_{\mu\nu}\)是与Levi-Civita联络相关的Christoffel符号。 记号约定 希腊指标\(\mu, \nu, \rho, \sigma, \ldots\)从0到\(d-1\)（若流形是\(d\)维）。 通过度规\(g_{\mu\nu}\)升降指标，例如\(V_\mu = g_{\mu\nu}V^\nu\)。 \(\Gamma^\rho_{\mu\nu}\)为Christoffel符号，定义为 \[ \Gamma^\rho_{\mu\nu} \;=\; \frac{1}{2}\,g^{\rho\lambda} \Bigl( \partial_\mu g_{\lambda\nu} +\partial_\nu g_{\lambda\mu} -\partial_\lambda g_{\mu\nu} \Bigr). \] B. 通过沿闭合回路的平行移动推导 1. 平行移动 考虑流形上的一个向量场\(V^\rho\)。沿曲线\(\gamma(\tau)\)对\(V^\rho\)进行的“平行移动”由条件给出：\(V^\rho\)沿该曲线的切向方向的协变导数为零。数学上表示为： \[ \frac{dV^\rho}{d\tau} = \frac{d x^\mu}{d\tau} \nabla_\mu V^\rho = 0. \] 展开协变导数，得到： ...

在这篇文章中，我希望建立一些微分几何的基础，以便理解爱因斯坦场方程。 1. 黎曼曲率张量及其性质 黎曼曲率张量的定义： 黎曼曲率张量 \( R^i_{jkl} \) 是微分几何中的一个基本对象，表示流形的内禀曲率。其定义为： \[ R^i_{jkl} = \partial_k \Gamma^i_{jl} - \partial_l \Gamma^i_{jk} + \Gamma^i_{km} \Gamma^m_{jl} - \Gamma^i_{lm} \Gamma^m_{jk} \] 其中 \( \Gamma^i_{jk} \) 是克里斯托费尔符号。 黎曼曲率张量的性质： 比安基恒等式： \[ \nabla_m R^i_{jkl} + \nabla_k R^i_{jlm} + \nabla_l R^i_{jmk} = 0 \] 对称性： 在最后两个指标上反对称： \[ R^i_{jkl} = -R^i_{jlk} \] 在交换第一对与第二对指标时反对称： \[ R^i_{jkl} = -R^j_{ikl} \] 应用： 黎曼张量刻画当沿闭合回路进行平行输运时，向量被旋转或改变的程度。 2. 测地线及其推导原理 测地线方程： 测地线是在弯曲空间中连接两点的最短路径，其方程可通过极小化作用量得到： \[ \int ds = \int \sqrt{g_{\mu \nu} \frac{dx^\mu}{d\tau} \frac{dx^\nu}{d\tau}} \, d\tau \] 其中 \( g_{\mu \nu} \) 是度量张量，\( \tau \) 是仿射参数。 ...

我想了解一些关于国债的基础信息。我感兴趣的问题包括： 票面利率如何确定？ 收益率如何确定？ 收益率与价格的关系是什么？ 收益率与利率的关系是什么？ 从宏观经济角度，美国、中国和日本的国债市场处于什么状态？ 我们参考这些讲座。 国债是政府为筹资而发行的债务证券。它们以稳定性著称，是金融市场中的关键工具。以下是与国债相关关键概念的结构化说明。 1. 国债的内在属性 国债具有两个内在属性： 面值（本金）： 到期时支付给债券持有人的金额。通常面值为$1,000，但可能有所不同。 票面利率： 发行时设定的固定年度利率，以面值的百分比表示。 示例： 面值为$1,000、票面利率为5%的债券每年支付$50利息。 2. 国债如何运作 持有国债时，你将收到： 定期票息支付： 这些支付基于债券的票面利率。通常为半年付。 到期偿还面值： 在债券期限结束时，将向持有人返还面值。 示意图：债券支付时间线 |--------------------|--------------------|--------------------| Maturity $25 $25 $25 $1,000 (Coupon Payment) (Coupon Payment) (Coupon Payment) (Face Value + Coupon Payment) 3. 市场交易：价格与收益率 国债在二级市场交易，其价格和收益率会随市场状况变化。 价格： 债券价格是投资者愿意为其支付的金额。由于利率和需求变化，价格可能与面值不同。 如果市场利率低于债券的票面利率，债券价格会上涨（溢价）。 如果市场利率高于债券的票面利率，债券价格会下跌（折价）。 收益率： 投资者以当前价格买入并持有至到期所获得的回报。 反向关系： 价格与收益率反向变动。债券价格上升时，收益率下降，反之亦然。 4. 债券价格与到期收益率的计算 债券价格等于所有未来现金流（票息支付和面值）按到期收益率（YTM）贴现后的现值之和。 公式： \[ P = \sum \frac{C}{(1 + YTM)^t} + \frac{F}{(1 + YTM)^T} \] 其中： ...

在这篇文章中，我想构建一些关于爱因斯坦场方程的基础知识。首先，我将基于上一篇关于宇宙学的文章提出几个问题： 什么是爱因斯坦张量和能量-动量张量？ 爱因斯坦场方程是如何推导出来的？ 在宇宙学背景下，什么是完美流体近似？ 我们如何从爱因斯坦场方程推导出弗里德曼方程和加速度方程？ 我们参考这门课程。 1. 四速度与固有时 四速度 \( u^\mu \) 描述粒子在时空中的运动，定义为： \[ u^\mu = \frac{dx^\mu}{d\tau} \] 其中： \( x^\mu = (ct, x, y, z) \)：时空坐标。 \( \tau \)：固有时，即在粒子静止系中测得的时间。 在广义相对论中，采用 \( (-, +, +, +) \) 度规符号约定时，四速度的归一化条件为： \[ u^\mu u_\mu = -c^2 \]如果我们使用自然单位制（\( c = 1 \)），则可简化为： \[ u^\mu u_\mu = -1 \] 2. 能量与动量 四动量 \( P^\mu \) 将能量 \( E \) 与空间动量 \( \vec{p} \) 组合在一起： \[ P^\mu = (E/c, \vec{p}) \]关键组成： ...