懒猫学术 - GARFIELDTOM
type: channel
link: https://www.gt.ac.cn/author/1/
description: GARFIELDTOM
items: 能源电化学基础：液相传质动力学概述
items: 能源电化学基础：电极的极化过程概述
items: 能源电化学基础：电极/溶液界面结构与性质
items: 加密算法笔记：从椭圆曲线数学到国密SM2算法
items: 加密算法笔记：国密SM4算法与AES算法的对比
items: The Evolution of Materials Science Models From Static to Dynamic
items: The Application of Machine Learning in Materials Science and Engineering
items: 量子力学：近似方法微扰理论
items: 量子力学：矩阵语言本征值与本征函数
items: 量子力学：对易、不确定性与测量

能源电化学基础：液相传质动力学概述
type: item
date: 2025-11-27T00:43:40+08:00
description: 理论基础液相传质的三种不同方式传质机制示意在电解液中，离子的运动并非只有一种方式。我们可以将其想象为三种不同动力驱动的“竞跑者”。要研究纯粹的电化学反应，我们必须学会控制它们。$$ J_{total} = J_{migration} + J_{diffusion} + J_{convection} $$电迁移由电场（电位梯度）驱动。这是电解池中通过电流的必要条件，但在研究电极表面反应动力学时，它是一个干扰项。因为它让计算变得极为复杂（电场与浓度场耦合）。控制方法：加入大量“局外电解质”（支持电解质）。它们承担了绝大部分的离子导电任务，屏蔽了电场对反应离子的作用力，从而“剔除”电迁移项。对流由机械力（搅拌、旋转）或密度差（自然对流）驱动。它将物质从溶液本体“快递”到电极附近，但无法送达“最后一公里”（电极表面极薄的液层）。控制方法：在稳态技术中（如旋转圆盘电极RDE），利用精确控制的对流来维持一个稳定的供给层；在暂态技术中（如静止循环伏安），则需严格避免对流，以免破坏扩散层。扩散由化学位梯度（浓度差）驱动。这是物质抵达电极表面的最后一步，也是必经之路。核心逻辑：电极反应消耗了表面的粒子 $\rightarrow$ 表面浓度降低 $\rightarrow$ 与本体溶液形成浓度差 $\rightarrow$ 驱动本体粒子向表面移动。菲克第一定律与稳态扩散模型当我们通过搅拌或旋转电极建立起一个稳定的流体环境时，浓度分布不再随时间变化，这就是“稳态”。流量与梯度的线性关系菲克第一定律是所有扩散理论的基石。它告诉我们，流得有多快，取决于浓度坡度有多陡。$$ J = -D \frac{dC}{dx} $$能斯特扩散层假说在真实流体中，流速靠近表面逐渐减为0。能斯特（Nernst）提出了一个简化的物理模型：假设在电极表面存在一层静止不动的液层，厚度为 $\delta$。在此层内只发生扩散，层外浓度维持为本体浓度 $C^0$。能斯特扩散层模型极限扩散电流的物理推导电流 $i$ 直接正比于流量 $J$。当电位极负，反应极快，电极表面的反应物来一个“吃”一个，导致表面浓度 $C^s$ 降为 $0$。此时浓度梯度达到最大值，电流也达到极限。$$ \text{Gradient}_{max} = \frac{C^0 - 0}{\delta} $$$$ i_d = nFA J_{max} = n F A D \frac{C^0}{\delta} $$图像与概念的联系这个公式揭示了为什么搅拌越快（对流越强），电流越大。物理图像：搅拌压缩了停滞层厚度 $\delta$。数学关系：$\delta$ 变小 $\rightarrow$ 分母变小 $\rightarrow$ 梯度变陡 $\rightarrow$ $i_d$ 变大。这就是旋转圆盘电极（RDE）的Levich方程的物理本质：$\delta$ 被精确控制为与转速 $\omega$ 相关。浓差极化的热力学解释为什么电流还没到极限扩散电流时，电位就已经偏离了平衡电位？这是传质速度与表面浓度之间的博弈。表面浓度与电位的纽带电极电位并不直接控制电流，它控制的是电极表面的浓度比。只要是可逆反应，表面必须时刻遵守能斯特方程。$$ E = E^{0'} + \frac{RT}{nF} \ln \frac{C_O^s}{C_R^s} $$电流与表面浓度的定量关系对于产物不溶（如金属沉积）的体系，我们假设 $C_R=1$。任何时刻的电流 $i$ 对应着某个表面浓度 $C_O^s$：$$ i = n F A D \frac{C^0 - C_O^s}{\delta} $$极限电流 $i_d$ 对应着表面浓度为0：$$ i_d = n F A D \frac{C^0}{\delta} $$联立上述两式，消去常数项，我们得到了表面浓度与电流比值的无量纲关系：$$ \frac{C_O^s}{C^0} = 1 - \frac{i}{i_d} $$浓差极化公式的诞生将上述表面浓度的表达式代回能斯特方程，就得到了教材中那个经典的极化曲线公式。它描述了$E$与$i$呈现出的S形曲线关系。$$ E = E^{0'} + \frac{RT}{nF} \ln C^0 + \frac{RT}{nF} \ln \left( 1 - \frac{i}{i_d} \right) $$图像理解浓差极化曲线当 $i$ 接近 $i_d$ 时，括号内趋近于0，对数项趋向负无穷，电位 $E$ 急剧负移。这就是极化曲线中出现“平台”的原因——即便电位再负，受扩散限制，电流也上不去了。菲克第二定律与暂态扩散（时间的介入）当我们停止搅拌，并在静止溶液中施加一个电信号时，扩散层不再是一个固定的厚度，它会随着时间向溶液深处“生长”。非稳态的核心浓度不仅随位置 $x$ 变化，还…

能源电化学基础：电极的极化过程概述
type: item
date: 2025-11-24T10:23:00+08:00
description: 核心概念归纳核心在于研究电极过程的速率及其影响因素。这一理论体系的逻辑起点是电极的平衡状态，通过引入电流，观察其偏离平衡的现象（极化），然后探究现象背后的原因（基本历程和RDS），最终形成一套研究方法（极化曲线）。电极的极化现象1. 平衡态 (Equilibrium State)定义： 处于热力学平衡状态的电极体系（可逆电极）。特征：氧化反应和还原反应速度相等（动态平衡）。电荷交换和物质交换处于动态平衡。净反应速度为零，外电路电流 $I=0$。电极电位为平衡电位 ($\phi_{eq}$ 或 $\phi_{平}$)。2. 极化现象的引入 (Introduction to Polarization)电极极化曲线示意图极化定义： 当电极上有电流通过时，电极失去原有平衡状态，其电极电位偏离平衡电位的现象。极化规律：阴极极化 (Cathodic Polarization)： 发生还原反应，电极电位 $\phi_c$ 偏离平衡电位 $\phi_{eq,c}$ 向更负方向移动 ($\phi_c < \phi_{eq,c}$)。阳极极化 (Anodic Polarization)： 发生氧化反应，电极电位 $\phi_a$ 偏离平衡电位 $\phi_{eq,a}$ 向更正方向移动 ($\phi_a > \phi_{eq,a}$)。3. 极化程度的定量描述过电位 (Overpotential, $\eta$) 是衡量极化程度的直接参数，是电极实际工作电位与平衡电位的差值。$$\eta_c = \phi_{eq,c} - \phi_c$$$$\eta_a = \phi_a - \phi_{eq,a}$$注意： 教材中规定 $\eta$ 习惯取正值，故阳极极化时 $\eta_a = \phi_a - \phi_{eq,a}$，阴极极化时取 $\eta_c = \phi_{eq,c} - \phi_c$。极化的实质与作用机制1. 极化的实质 (Essence of Polarization)电极极化的根本原因是：电子运动速度 ($v_{电子}$) 与电极反应速度 ($v_{反}$) 之间的矛盾。逻辑递进： 电子在外电路中的传输速度极快，远大于电化学反应的化学变化速度（$v_{电子} \gg v_{反}$）。当外加电流时，电极反应速度跟不上电子的供给或消耗速度，导致电荷在电极/溶液界面积累。极化类型电子与反应速度关系净电荷积累极化方向阴极极化电子流入电极速度 > 还原反应消耗电子速度负电荷积累$\phi$ 向负移阳极极化氧化反应产生电子速度 > 电子流出电极速度正电荷积累$\phi$ 向正移2. 作用机制 (Polarization vs. Depolarization)电极性质的变化取决于两种相互矛盾作用的对立统一：极化作用 (Polarization Effect)： 电子流动导致电极表面电荷积累，使电极电位偏离平衡状态。去极化作用 (Depolarization Effect)： 缓慢的电极反应吸收或释放电荷，使电极电位恢复平衡状态。极化曲线极化曲线是电化学动力学研究的基础工具。1. 电极反应速度与电流密度的关系在稳态条件下，电极反应的净速度 ($v$) 可以用外电路测得的净电流密度 ($j$) 来表示。对于电极反应 $O + ne \rightleftharpoons R$：非均相化学反应速度：$$v = - \frac{1}{S} \frac{dc_O}{dt} = \frac{1}{S} \frac{dc_R}{dt}$$法拉第定律关联：$$j = nFv$$结论： 当电极反应达到稳定状态时，实验测得的外电流密度 ($j$) 值就代表了电极反应的净速度 ($v$)。因此，极化曲线 ($\phi$ vs $j$ 或 $\eta$ vs $j$) 实际上反映了电极反应速度与电极电位之间的特征关系。2. 极化曲线图像与概念关系图2_锌在0.38mol/L $ ZnCl_2 $ 溶液（1）和氯化镀锌溶液（2）中的阴极极化曲线图像关系： 极化曲线的形状（如斜率）直接反映了电极过程的动力学特性。极化度 (Polarization Rate, $P$):$$P = \frac{d\phi}{dj} \quad \text{或} \quad P = \frac{d\eta}{dj}$$极化度在数值上具有电阻的量纲（反应电阻），斜率越大（即极化度越大），表明该电极过程进行越困难，受到的阻力越大。电极过程的基本历程与速度控制步骤电极过程不是一个简单反应，而是由一系列性质不同的单元步骤串联（有时并联）组成的复杂过程。1. 电极过程的基本历程 (Basic Steps)电极过程的基本历程示意图序号步骤名称核心行为动力学特点极化类型必需步骤(1)液相传质反应粒子从…

能源电化学基础：电极/溶液界面结构与性质
type: item
date: 2025-11-22T23:11:00+08:00
description: 理想极化电极与界面参数理想极化电极 (IPE) 是研究电极/溶液界面结构和性质的理论前提。它被定义为在研究电位范围内不发生任何电极反应的电化学体系。$R_f \to \infty$此时，通过电极的所有电流都用于改变界面电荷和电位，其等效电路是一个纯电容 $C$。在 IPE 体系中，界面的状态可以由三个核心参数随电极电位 $\phi$ 变化的函数关系来描述：界面张力 $\sigma$表面剩余电荷密度 $q$微分电容 $C_d$Lippmann 公式与零电荷电位研究界面性质首先需要建立热力学关系，即 Lippmann 公式。它将宏观可测量的界面张力 $\sigma$ 与界面的电荷状态 $q$ 联系起来。根据 Gibbs 等温吸附方程（忽略化学反应和组分变化）：$$\text{d}\sigma = - \sum \Gamma_i \text{d}\mu_i$$对于 IPE，溶液相组分 $\mu_i$ 不变 ($\text{d}\mu_i=0$)，但考虑到电极上的电子吸附（即 $q$ 的存在）：$$\text{d}\sigma = - \Gamma_e \text{d}\mu_e$$其中，电子的表面吸附量 $\Gamma_e$ 与表面电荷密度 $q$ 的关系为 $\Gamma_e = -q/F$，电子化学位变化 $\text{d}\mu_e$ 与电极电位 $\phi$ 的关系为 $\text{d}\mu_e = -F\text{d}\phi$。代入并简化可得著名的 Lippmann 公式：$$\left(\frac{\partial \sigma}{\partial \phi}\right)_{\text{T}, \text{p}, \mu_i} = -q$$该公式是连接宏观测量（电毛细曲线的斜率）和界面微观状态（表面电荷密度）的桥梁。汞电极的电毛细曲线与表面剩余电荷密度-电位曲线零电荷电位 ($\phi_0$) 的概念正是由 Lippmann 公式导出：当界面张力 $\sigma$ 达到最大值时，其微分值 $\left(\frac{\partial \sigma}{\partial \phi}\right)$ 为零。$$\left(\frac{\partial \sigma}{\partial \phi}\right)_{\phi=\phi_0} = 0$$此时，电极表面剩余电荷密度 $q$ 为零，该电位即为 $\phi_0$。微分电容与积分电容微分电容 $C_d$ 定义为表面电荷密度 $q$ 对电极电位 $\phi$ 的一阶导数，表征了界面储存电荷的能力：$$C_d = \left(\frac{\partial q}{\partial \phi}\right)_{\text{T}, \text{p}, \mu_i}$$微分电容 $C_d$ 与界面张力 $\sigma$ 的关系，即 $C_d$ 是 $\sigma$ 对 $\phi$ 的二阶微分：$$C_d = -\left(\frac{\partial^2 \sigma}{\partial \phi^2}\right)_{\text{T}, \text{p}, \mu_i}$$积分电容 $C_i$ 定义为从零电荷电位 $\phi_0$ 到某一电位 $\phi$ 之间的平均电容：$$C_i = \frac{q}{\phi - \phi_0} = \frac{1}{\phi - \phi_0} \int_{\phi_0}^{\phi} C_d \text{d}\phi$$微分电容法与电毛细曲线法的对比根据上述关系，微分电容 $C_d$ 是 $q$ 的微分函数，而 $q$ 是 $\sigma$ 的微分函数，$\sigma$ 是 $q$ 的积分函数。$C_d$ 是 $\sigma$ 的二阶微分。微分函数比积分函数更灵敏地反映原函数的细微变化。因此，通过 微分电容法 测量 $C_d$ 进而计算 $q$ ($q = \int C_d \text{d}\phi$) 比通过 电毛细曲线法 测量 $\sigma$ 进而计算 $q$ ($q = -\text{d}\sigma/\text{d}\phi$) 更精确、更灵敏。双电层模型定量推导与对比双电层模型从平板电容器模型（Helmholtz）发展到考虑热运动影响的分散层模型（Gouy-Chapman），最终到 Stern 模型，形成了紧密层和分散层串联的基本结构。1. Helmholtz (H) 模型（紧密层模型）紧密层模型示意图将双电层简化为平行板电容器，忽略热运动和电位在紧密层内的非均匀性。$$C_{\text{H}} = \frac{\epsilon_0 \epsilon_r}{l}$$其中，$\epsilon_0$ 是真空介电常数，$\eps…

加密算法笔记：从椭圆曲线数学到国密SM2算法
type: item
date: 2025-09-30T14:21:00+08:00
description: 1. 椭圆曲线上的 Abel 群结构设 $\mathbb{F}_q$ 为有限域， $q = p$ 为奇素数（密码学中最常见情形）。一条椭圆曲线 $E$ 的短 Weierstrass 形式为$$ E : \; y^2 = x^3 + A x + B, \quad A, B \in \mathbb{F}_p, $$其中判别式$$ \Delta = -16(4A^3 + 27B^2) \not\equiv 0 \pmod{p}. $$曲线上的仿射点集合连同无穷远点 $\mathcal{O}$ 构成 $E(\mathbb{F}_p)$。由 Hasse 定理，$$ |E(\mathbb{F}_p)| = p + 1 - t, \quad |t| \le 2\sqrt{p}. $$1.1 弦切群律令 $P = (x_1, y_1), Q = (x_2, y_2) \in E(\mathbb{F}_p) \setminus \{\mathcal{O}\}$ ，且 $P \neq -Q$。a) 点加法 ($P \neq Q$)直线 $L$ 通过 $P, Q$ 的斜率$$ \lambda = \frac{y_2 - y_1}{x_2 - x_1} \pmod{p} $$将直线方程 $y = \lambda(x - x_1) + y_1$ 代入曲线方程，解得第三个交点 $R = (x_3, y_3)$ 的坐标：$$ x_3 = \lambda^2 - x_1 - x_2, $$$$ y_3 = \lambda(x_1 - x_3) - y_1. $$群加法定义为 $R$ 关于 $x$ 轴的反射点：$$ P + Q = (x_3, -y_3). $$图1 点加法：弦切群律 $P+Q=R_{reflection}$椭圆曲线上的点加法 $P+Q$ 的几何过程。作一条直线 $L$ 穿过 $P$ 和 $Q$，它与曲线交于第三点 $R$。 $P+Q$ 定义为 $R$ 关于 $x$ 轴的对称点（反射）。这个几何操作对应于群律的代数公式 $P+Q=(x_3, -y_3)$。b) 点倍乘 ($P = Q$)在 $P$ 处作曲线的切线 $L$。其斜率：$$ \lambda = \frac{3x_1^2 + A}{2y_1} \pmod{p} \quad \text{if } y_1 \neq 0. $$$2P$ 的计算公式与点加法相同，即：$$ x_3 = \lambda^2 - 2x_1, $$$$ y_3 = \lambda(x_1 - x_3) - y_1. $$群加法定义为 $R=(x_3, y_3)$ 关于 $x$ 轴的反射点：$$ 2P = (x_3, -y_3). $$图2 点倍乘：弦切群律 $2P=R_{reflection}$椭圆曲线上的点倍乘 $2P$ 的几何过程。在点 $P$ 处作曲线的切线 $L$，它与曲线交于第二点 $R$。 $2P$ 定义为 $R$ 关于 $x$ 轴的对称点（反射）。这个几何操作对应于代数公式中的 $P=Q$ 特殊情形，即 $2P=(x_3, -y_3)$。特殊情形：逆元： 任意点 $P=(x_1, y_1)$ 的逆元为 $-P = (x_1, -y_1)$。若 $P = -Q$，则 $P + Q = \mathcal{O}$。单位元： $\mathcal{O}$ 为单位元， $P + \mathcal{O} = P$。结合律的严格代数证明见 Silverman [1, Theorem III.2.3]，共需验证九类情形，此处从略。最终得到 $(E(\mathbb{F}_p), +)$ 为有限交换群。1.2 点乘与阶对 $k \in \mathbb{Z}$，记 $[k]P = P + \cdots + P$（ $k$ 次）。最优算法为二进制 double-and-add，复杂度 $O(\log k)$ 次群运算。选取基点 $G \in E(\mathbb{F}_p)$ 使得其阶$$ n = ord(G) = \min\{ m > 0 : [m]G = \mathcal{O} \} $$为大素数，子群 $\langle G \rangle$ 循环且有 $n$ 个元素。1.3 椭圆曲线离散对数问题（ECDLP）给定 $Q = [k]G$，求 $k \in \mathbb{Z}_n$。当前最快通用算法为 Pollard-ρ，期望复杂度 $\sqrt{\pi n}/2$ 次群运算。对于 256-bit 阶，约 $2^{128}$ 次运算，构成现代 ECC 全部安全性的基础。2. 标准椭圆曲线密码算法图3 ECDH密钥交换流程：基于椭圆曲线离散对数问题 (ECDLP)椭圆曲线 Diffie-Hellman (ECDH) 密钥交换流程…

加密算法笔记：国密SM4算法与AES算法的对比
type: item
date: 2025-09-25T15:55:00+08:00
description: SM4 分组密码算法推导SM4是一种分组密码算法，其分组长度和密钥长度均为128位。它采用32轮迭代的非平衡Feistel结构。1. 初始设定明文分组长度为128位，密钥长度为128位。首先，将128位的输入明文 $P$ 分为4个32位的字 $(X_0, X_1, X_2, X_3)$。$$ P = (X_0, X_1, X_2, X_3) \quad \text{其中 } X_i \in \{0,1\}^{32} $$加密过程共进行 $i=0, 1, \dots, 31$ 共32轮迭代。2. 轮函数 FSM4的精髓在于其轮函数 $F$。每一轮的输入为4个32位的字 $(X_i, X_{i+1}, X_{i+2}, X_{i+3})$ 和一个32位的轮密钥 $rk_i$。$$ X_{i+4} = F(X_i, X_{i+1}, X_{i+2}, X_{i+3}, rk_i) $$这个函数具体定义为：$$ X_{i+4} = X_i \oplus T(X_{i+1} \oplus X_{i+2} \oplus X_{i+3} \oplus rk_i) $$3. T 变换轮函数中的核心是合成变换 $T$，$T$ 是一个可逆的32位输入到32位输出的变换，由一个非线性变换 $\tau$ 和一个线性变换 $L$ 复合而成。$$ T(A) = L(\tau(A)) $$非线性变换 $\tau$:这个变换将一个32位的字 $A$ 分为4个8位的字节 $(a_0, a_1, a_2, a_3)$，然后分别对每个字节进行S盒（S-box）代换，最后将结果合并成一个32位的字 $B$。$$ A = (a_0, a_1, a_2, a_3) \quad \text{其中 } a_j \in \{0,1\}^8 $$$$ B = \tau(A) = (Sbox(a_0), Sbox(a_1), Sbox(a_2), Sbox(a_3)) $$S盒是一个预定义的8位输入到8位输出的非线性置换表，提供了抗击差分分析和线性分析的能力。线性变换 $L$:这个变换对32位的输入 $B$ 进行线性混合，以增强扩散性。它通过异或和循环移位操作实现。$$ L(B) = B \oplus (B \lll 2) \oplus (B \lll 10) \oplus (B \lll 18) \oplus (B \lll 24) $$其中 $\lll k$ 表示循环左移 $k$ 位。4. 输出经过32轮迭代后，得到 $(X_{32}, X_{33}, X_{34}, X_{35})$。最后将这4个字进行一次反序变换 $R$，得到密文输出 $C$。$$ C = R(X_{32}, X_{33}, X_{34}, X_{35}) = (X_{35}, X_{34}, X_{33}, X_{32}) $$5. 密钥扩展SM4的128位主密钥 $MK$ 也被分为4个32位的字:$$ MK_{0}, MK_{1}, MK_{2}, MK_{3} $$轮密钥 $rk_i$ 由主密钥通过一个类似的迭代过程生成，同样也使用了S盒。与AES的不同之处SM4和AES（高级加密标准）虽然都是128位分组的对称密码算法，但它们的内部结构和设计哲学有显著不同。特性SM4 (GM/T 0002-2012)AES (Rijndael)核心结构Feistel 网络SPN (Substitution-Permutation Network)分组长度128 位128 位密钥长度128 位128, 192, 或 256 位迭代轮数32 轮10, 12, 或 14 轮 (取决于密钥长度)加密/解密结构相同，仅轮密钥顺序相反解密是加密的逆过程，需要逆向操作（如InvMixColumns）1. 结构差异 (Feistel vs. SPN)SM4 (Feistel网络):Feistel结构将数据块分成两半（在SM4中是分成4个字，但本质上是非平衡Feistel），在每一轮中，一块数据通过一个轮函数 $F$ 进行处理，其结果与另一块数据进行异或。这种结构的一个巨大优势是加解密过程非常相似，解密时只需按相反的顺序使用轮密钥即可，这在硬件实现上可以节省资源。AES (SPN网络):SPN结构在每一轮中都对整个数据块进行操作。AES的一轮由四个阶段组成：SubBytes: 类似于SM4的S盒变换，进行非线性替换。$$ S'_{i,j} = \text{S-box}(S_{i,j}) $$ShiftRows: 对状态矩阵的行进行循环移位，实现置换。$$ S'_{i,j} = S_{i, (j + \text{shift}(i, N_b)) \pmod{N_b}} $$MixColumns: 对状态矩阵的列进行线性混合，提供扩散。…

The Evolution of Materials Science Models From Static to Dynamic
type: item
date: 2025-09-21T14:26:00+08:00
description: The Foundation - Constructing the High-Dimensional Element Space1. Introduction: Translating Chemistry into GeometryThe primary challenge in applying machine learning to materials science is one of translation. A model does not understand the abstract concept of an element like "Iron" (Fe) with its associated chemical intuition. Models understand mathematics: vectors, matrices, and the geometric relationships between them.The foundational step, therefore, is to create a robust, numerical representation for every element. We achieve this by constructing a high-dimensional "Materials Property Space". In this space, every element is defined not by its name, but by a unique set of coordinates—a vector—where each coordinate corresponds to a fundamental physical or chemical property.This act of "featurization" transforms the periodic table from a chart into a geometric map. The model can then learn by measuring distances, identifying clusters, and finding directional trends within this space, effectively discovering the rules of chemistry and physics through the language of geometry.2. Mathematical Formulation: The Element and Compound VectorsThe Element VectorWe define each element as a vector $\mathbf{v}$ in an $n$-dimensional property space. Each component of the vector corresponds to a quantifiable property.$$ \mathbf{v}_{\text{element}} = [ p_1, p_2, p_3, \dots, p_n ] $$Where:$\mathbf{v}_{\text{element}}$ is the feature vector for a single element.$p_k$ is the value of the $k$-th property. These properties can include:Positional Properties: Group, Period on the periodic table.Electronic Properties: Pauling Electronegativity, First Ionization Energy, Electron Affinity.Physical Properties: Atomic Radius, Covalent Radius, Atomic Mass, Melting Point.Orbital Properties: Number of valence s, p, d, and f electrons.The Compound VectorTo represent a compound with a specific stoichiometry (e.g., $A_x B_y$), we cannot simply concatenate the vectors of its constituent elements. …

The Application of Machine Learning in Materials Science and Engineering
type: item
date: 2025-09-17T15:45:00+08:00
description: Advances in Machine Learning and Their Application to Materials ScienceMachine learning is rapidly transforming materials science, shifting the field from a traditional trial-and-error approach to a more predictive and accelerated paradigm of materials discovery and design. This evolution is driven by the growing availability of large materials databases, increased computational power, and the development of sophisticated machine learning algorithms capable of uncovering complex relationships between a material's structure and its properties. Current developments focus on using machine learning for rapid screening of candidate materials, predicting their performance, and even generating novel materials with desired characteristics. For instance, machine learning models are being employed to discover new materials for solar panels by predicting how structural changes will affect their efficiency in converting sunlight to electricity.A significant area of progress lies in the application of advanced AI techniques to solve fundamental challenges in materials science. Researchers at Caltech, for example, have recently developed an AI-based method to dramatically speed up complex calculations of quantum interactions, such as those among atomic vibrations (phonons), which are critical for determining a material's thermal properties. This approach uses machine learning to compress the enormous tensors that describe these interactions, reducing computational complexity by orders of magnitude without sacrificing accuracy. Furthermore, the integration of machine learning is poised to enhance every stage of the scientific process, from generating hypotheses and planning experiments to analyzing results, fostering a human-machine partnership that promises to accelerate research at an unprecedented rate. This data-driven approach is significantly reducing the cost, risk, and time associated with materials research and development.An Introduction to the Most Advanced Models in Ma…

量子力学：近似方法微扰理论
type: item
date: 2025-06-29T15:48:00+08:00
description: 在现实世界中，能够被精确求解的量子系统（如一维无限深势阱、氢原子）少之又少。大多数系统都因为存在各种复杂的相互作用而难以求解。微扰理论 (Perturbation Theory) 正是为此而生的，它是一种强大而优雅的近似方法。核心思想：将一个复杂系统的哈密顿量 $\hat{H}$ 分解为一个我们可以精确求解的简单部分 $\hat{H}_0$ 和一个很小的修正部分 $\hat{H}'$（称为微扰）。$$ \hat{H} = \hat{H}_0 + \hat{H}' $$然后，我们从 $\hat{H}_0$ 的已知解出发，系统地计算微扰 $\hat{H}'$ 对能量和波函数造成的各级修正。一、 定态非简并微扰理论当我们要修正的能级 $E_n^{(0)}$ 是非简并的（即只有一个量子态对应这个能量），我们使用非简并微扰理论。1. 能量的修正公式 (必记！)一级能量修正 $E_n^{(1)}$:$$ E_n^{(1)} = \langle\psi_n^{(0)}| \hat{H}' |\psi_n^{(0)}\rangle = H'_{nn} $$物理意义: 一级修正是微扰哈密顿量在未受扰动状态下的平均值。它只考虑了微扰对该能级自身的直接影响。二级能量修正 $E_n^{(2)}$:$$ E_n^{(2)} = \sum_{k \neq n} \frac{|\langle\psi_k^{(0)}| \hat{H}' |\psi_n^{(0)}\rangle|^2}{E_n^{(0)} - E_k^{(0)}} = \sum_{k \neq n} \frac{|H'_{kn}|^2}{E_n^{(0)} - E_k^{(0)}} $$物理意义: 二级修正描述了微扰如何通过将态 $n$ 与所有其他态 $k$ “混合”起来，从而间接地影响能量。这个修正的大小与能级间的耦合强度 $|H'_{kn}|^2$ 成正比，与能级间距 $E_n^{(0)} - E_k^{(0)}$ 成反比。能量的二级近似值:$$ E_n \approx E_n^{(0)} + E_n^{(1)} + E_n^{(2)} $$2. 波函数的一级修正 $ |\psi_n^{(1)}\rangle $$$ |\psi_n^{(1)}\rangle = \sum_{k \neq n} \frac{\langle\psi_k^{(0)}| \hat{H}' |\psi_n^{(0)}\rangle}{E_n^{(0)} - E_k^{(0)}} |\psi_k^{(0)}\rangle = \sum_{k \neq n} \frac{H'_{kn}}{E_n^{(0)} - E_k^{(0)}} |\psi_k^{(0)}\rangle $$物理意义: 修正后的波函数 $|\psi_n\rangle \approx |\psi_n^{(0)}\rangle + |\psi_n^{(1)}\rangle$ 不再是纯粹的 $|\psi_n^{(0)}\rangle$，而是“混入”了其他所有本征态的成分。二、 定态简并微扰理论当我们要修正的能级 $E_n^{(0)}$ 是简并的（即有多个量子态对应同一个能量），非简并微扰理论的公式会因为分母出现零而失效。此时必须使用简并微扰理论。核心思路: 微扰会“打破”或“部分打破”原来的简并，将一个简并的能级分裂成几个不同的能级。我们的任务就是找到这些分裂后的能级。求解步骤:确定简并子空间: 找出所有对应于简并能级 $E_n^{(0)}$ 的本征态 $\{|\psi_{n1}^{(0)}\rangle, |\psi_{n2}^{(0)}\rangle, \dots, |\psi_{ng}^{(0)}\rangle\}$。构建微扰矩阵 $W$: 在这个简并子空间中，计算微扰哈密顿量 $\hat{H}'$ 的矩阵。这是一个 g x g 的矩阵，其矩阵元为：$$ W_{ij} = \langle\psi_{ni}^{(0)}| \hat{H}' |\psi_{nj}^{(0)}\rangle $$对角化矩阵 $W$: 求解矩阵 $W$ 的本征值 $\{E_1^{(1)}, E_2^{(1)}, \dots, E_g^{(1)}\}$。得到一级修正: 这些本征值就是原来简并能级的一级能量修正。分裂后的新能级为：$$ E_{ni} = E_n^{(0)} + E_i^{(1)} $$三、 综合应用例题问题:考虑一个二维的量子系统，其未微扰的哈密顿量 $\hat{H}_0$ 在某组基底下 $\{|1\rangle, |2\rangle, |3\rangle\}$ 的表示是对角的，但能级 $E_1^{(0)}$ 和 $E_2^{(0)}$ 简并。现在给系统施加一…

量子力学：矩阵语言本征值与本征函数
type: item
date: 2025-06-20T15:44:00+08:00
description: 在量子力学中，当系统只有有限个离散状态时（例如电子自旋的“上”和“下”），使用矩阵力学来描述会比使用波函数更简洁、更直观。在这种描述下，物理量算符是矩阵，量子态是列向量。而求解系统的物理性质，就转化为了一个核心的数学问题——求解矩阵的本征值和本征向量。一、 核心概念与求解步骤1. 物理意义本征值 (Eigenvalue): 代表了对物理量进行测量时，可能得到的确定数值。对于哈密顿量矩阵，其本征值就是系统的能量能级。本征函数/本征向量 (Eigenvector): 代表了系统处于该物理量的本征态。当系统处于一个本征态时，测量对应的物理量，结果是确定无疑的。2. 核心方程 (本征值方程)对于一个算符（矩阵）$\hat{A}$，其本征值 $λ$ 和本征向量 $v$ 满足以下关系：$$ \hat{A}v = \lambda v $$这个方程的物理含义是：当算符 $\hat{A}$ 作用于它的本征态 $v$ 时，效果仅仅是给这个态乘以一个常数 $λ$，而没有改变这个态的“方向”。3. 求解流程 (四步法)步骤操作公式1. 构造特征方程目标是求解本征值 $λ$$\det(\hat{A} - \lambda I) = 0$2. 解代数方程解上一步得到的关于 $λ$ 的方程(无通用公式)3. 求解本征向量将每个求出的 $λ$ 代回，解线性方程组$(\hat{A} - \lambda I)v = 0$4. 归一化使本征向量的模长为1，满足概率守恒$v_{\text{norm}} = \frac{v}{\sqrt{v^\dagger v}}$二、 综合应用例题问题:考虑一个与电子自旋相关的哈密顿量，它描述了电子在一个指向 $(1, 0, 1)$ 方向的磁场中的能量。该哈密顿量可以表示为：$$ \hat{H} = \omega_0 (\hat{S}_x + \hat{S}_z) $$其中 $\omega_0$ 是一个具有能量量纲的常数，$\hat{S}_x$ 和 $\hat{S}_z$ 是自旋算符。求:写出该哈密顿量 $\hat{H}$ 的 2x2 矩阵形式。求解该系统的能量本征值。求解与每个能量本征值对应的、归一化的能量本征态。解答:1. 写出 $\hat{H}$ 的矩阵形式我们知道自旋算符的矩阵形式为 $\hat{S}_i = (\hbar/2)\sigma_i$。$$ \hat{S}_x = \frac{\hbar}{2}\begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix}, \quad \hat{S}_z = \frac{\hbar}{2}\begin{pmatrix} 1 & 0 \\ 0 & -1 \end{pmatrix} $$将它们代入哈密顿量的表达式：$$ \hat{H} = \omega_0 \left( \frac{\hbar}{2}\begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix} + \frac{\hbar}{2}\begin{pmatrix} 1 & 0 \\ 0 & -1 \end{pmatrix} \right) $$$$ \boxed{\hat{H} = \frac{\hbar\omega_0}{2} \begin{pmatrix} 1 & 1 \\ 1 & -1 \end{pmatrix}} $$2. 求解能量本征值我们求解特征方程 $\det(\hat{H} - E \cdot I) = 0$。为了方便，我们先令 $\alpha = \frac{\hbar\omega_0}{2}$。$$ \det\left( \alpha \begin{pmatrix} 1 & 1 \\ 1 & -1 \end{pmatrix} - \begin{pmatrix} E & 0 \\ 0 & E \end{pmatrix} \right) = \det\begin{pmatrix} \alpha-E & \alpha \\ \alpha & -\alpha-E \end{pmatrix} = 0 $$计算行列式：$$ (\alpha-E)(-\alpha-E) - \alpha^2 = 0 $$$$ -(\alpha-E)(\alpha+E) - \alpha^2 = 0 $$$$ -( \alpha^2 - E^2) - \alpha^2 = 0 $$$$ E^2 - 2\alpha^2 = 0 $$$$ E^2 = 2\alpha^2 \quad \implies \quad E = \pm \sqrt{2}\alpha $$将 $\alpha = \frac{\hbar\omega_0}{2}$ 代回，得到两个能量本征值：$$ \boxed{E_1 …

量子力学：对易、不确定性与测量
type: item
date: 2025-06-12T15:40:00+08:00
description: 经典世界里，我们可以随心所欲地测量一个物体的任何属性。但在量子世界，测量行为本身会深刻地影响系统，测量的顺序变得至关重要。对易关系正是描述这种顺序效应的数学语言，它直接导向了著名的不确定性原理。一、 对易子：算符的“交换律”1. 定义两个算符 $\hat{A}$ 和 $\hat{B}$ 的对易子 (Commutator) 定义为：$$ [\hat{A}, \hat{B}] = \hat{A}\hat{B} - \hat{B}\hat{A} $$它衡量的是两种操作顺序（先A后B vs 先B后A）之间的差异。2. 核心物理意义如果 $[\hat{A}, \hat{B}] = 0$:称 $\hat{A}$ 与 $\hat{B}$ 对易 (commute)。物理含义: 对应的物理量 $A$ 和 $B$ 是相容的 (compatible)。这意味着：它们可以同时被精确测量。它们拥有共同的本征函数系。测量 $A$ 不会影响 $B$ 的状态，反之亦然。如果 $[\hat{A}, \hat{B}] \neq 0$:称 $\hat{A}$ 与 $\hat{B}$ 不对易 (do not commute)。物理含义: 对应的物理量 $A$ 和 $B$ 是不相容的 (incompatible)。这意味着：它们不能同时被精确测量。测量其中一个，必然会干扰另一个的状态，使其变得不确定。二、 不确定性原理：不对易的定量描述不对易的关系可以用一个精确的数学不等式来量化，这就是海森堡不确定性原理。1. 通用形式对于任意两个物理量 $A$ 和 $B$，其测量结果的标准差（不确定度）$\Delta A$ 和 $\Delta B$ 必须满足：$$ (\Delta A)^2 (\Delta B)^2 \ge \left( \frac{1}{2i} \langle[\hat{A}, \hat{B}]\rangle \right)^2 $$核心思想: 两个算符的对易子越大（离零越远），我们就越不可能同时把它们测准。不确定性是量子世界内禀的属性，而非测量技术不精导致。2. 两个最重要的特例物理量对易关系不确定性关系位置-动量$[\hat{x}, \hat{p}_x] = i\hbar$$\Delta x \cdot \Delta p_x \ge \frac{\hbar}{2}$角动量分量$[\hat{L}_x, \hat{L}_y] = i\hbar \hat{L}_z$$(\Delta L_x)^2 (\Delta L_y)^2 \ge \frac{\hbar^2}{4}\langle\hat{L}_z\rangle^2$三、 综合应用例题问题:一个粒子被限制在一维无限深势阱中（$0 \le x \le a$），并处于基态 $\psi_1(x) = \sqrt{\frac{2}{a}}\sin\left(\frac{\pi x}{a}\right)$。求:计算该状态下粒子位置的平均值 $\langle x \rangle$ 和动量的平均值 $\langle p_x \rangle$。计算该状态下粒子位置的平方的平均值 $\langle x^2 \rangle$ 和动量的平方的平均值 $\langle p_x^2 \rangle$。利用以上结果，计算位置不确定度 $\Delta x$ 和动量不确定度 $\Delta p_x$。验证你的计算结果是否满足海森堡不确定性关系 $\Delta x \cdot \Delta p_x \ge \hbar/2$。解答:1. 计算平均值 $\langle x \rangle$ 和 $\langle p_x \rangle$$\langle x \rangle$:$$ \langle x \rangle = \int_0^a \psi_1^*(x) \hat{x} \psi_1(x) dx = \frac{2}{a} \int_0^a x \sin^2\left(\frac{\pi x}{a}\right) dx $$由于被积函数 $x \sin^2(\dots)$ 是关于中心点 $x=a/2$ 对称的，所以积分结果就是中心点的值。$$ \langle x \rangle = \frac{a}{2} $$$\langle p_x \rangle$:$$ \langle p_x \rangle = \int_0^a \psi_1^*(x) (-i\hbar\frac{d}{dx}) \psi_1(x) dx = \frac{2}{a} \int_0^a \sin\left(\frac{\pi x}{a}\right) (-i\hbar) \left(\frac{\pi}{a}\right)\cos\left(\frac{\pi x}{a}\r…