钠修饰多环芳烃上D2/H2吸附选择性（I）：精确建模方法与低压选择性

《Physical Chemistry Chemical Physics》：On the D2/H2 adsorption selectivity of sodium decorated polycyclic aromatic hydrocarbons. I. Accurate modeling approaches and low pressure selectivity

　　

在寻找经济可行的氘(D)与普通氢(H)同位素分离方法的进程中，研究人员通过理论计算研究了H2与D2在钠(Na)修饰多环芳烃——质子化石墨烯(coronene, [H–Cor–Na]+)及作为石墨烯原型的circumcircumcoronene([CCCor–Na])上的平衡吸附行为。研究人员构建了H2/D2与上述基底相互作用的精确势能面（采用类原子近似PsAt与刚体近似RB），并求解吸附态H2和D2的薛定谔方程以获得相应的配分函数，进而估算吸附等温线与D2/H2吸附选择性，完整表征了低压条件下的行为。结果表明，类原子模型(PsAt)不足以正确预测上述性质，且忽略热激发态贡献会得到精确选择性的下界。定量方面，两种材料的氘选择性均远高于工业应用所需阈值6，同时H2自旋异构体（正氢ortho-H2与仲氢para-H2）的分离亦显示有效性。此类基底有望作为色谱固定相，值得进一步研究其在更高压力下的分离性能，相关理论框架与结果详见系列第二部分。

论文解读：钠修饰多环芳烃上D₂/H₂吸附选择性（I）——精确建模与低压行为

研究背景与意义

重氢同位素（氘D、氚T）在科技领域及聚变堆燃料中具有重要价值，但其天然丰度极低，需从微量H/D混合气中分离。由于同位素理化性质极为相近，分离须依赖二者量子力学运动性质的微小差异（零点能不同），通常需在低温下操作以放大此效应。化学亲和力量子筛分（chemical affinity quantum sieving）利用物理吸附（physisorption）或化学吸附（chemisorption）的强度差异进行分离，其中物理吸附可避免表面催化同位素扰拌（scrambling）。碳质材料经碱金属（如Na、K）修饰可增强与H₂的相互作用（physisorption），增大同位素吸附能差，是潜在的经济型分离介质。本文（发表于Physical Chemistry Chemical Physics）针对Na修饰多环芳烃（Polycyclic Aromatic Hydrocarbons, PAHs）在低温低压下分离D₂/H₂的能力开展理论研究，旨在验证其工业应用潜力并阐明精确建模的必要性。

主要技术方法

研究人员选取Na修饰质子化晕苯（[H–Cor–Na]⁺）与Na修饰circumcircumcoronene（[CCCor–Na]，模拟中性石墨烯片层）为模型体系。首先基于高精度电子结构计算校准X₂（X=H,D）与基底的相互作用势：X₂分别处理为球对称点电荷伪原子（pseudo-atom, PsAt）及考虑取向与转动的刚性转子（rigid-body, RB），后者含库仑与改进Lennard-Jones（Improved Lennard-Jones, ILJ）项。随后在独立吸附位点假设下，采用离散变量表象（Discrete Variable Representation, DVR）基组与球谐函数Y_JM(θ,φ)展开吸附分子哈密顿量，数值对角化求解束缚态能谱E_l。由吸附配分函数q_ads(T,X₂)=Σ_lexp[-βE_l(X₂)]及气相化学势推导Langmuir型低压吸附等温线θ_X2=p_X2χ_r(T,X₂)q_ads(T,X₂)/(1+p_X2χ_r(T,X₂)q_ads(T,X₂))，并定义D₂/H₂选择性S_D2/H2=(θ_D2/θ_H2)/(x_D2/x_H2)，扩展至多组分（考虑ortho/para自旋异构体统计权重）情形。

研究结果

3.1 The case of [H–Cor–Na]⁺

通过求解RB与PsAt模型下吸附H₂/D₂的束缚态，得到基态能量E₀（相对无转动J=0无穷远分离态）：RB-even-J H₂为-1116.3 K，D₂为-1250.7 K（ΔE=134.4 K）；PsAt模型下H₂为-871.1 K，D₂为-940.3 K（ΔE=69.2 K）。RB模型因考虑分子取向各向异性（H₂四极矩与Na正电荷静电作用），基态结合更深。激发态谱显示低能级对应平面内绕Na原子的角向激发（类粒子在环模型），高能级含z或ρ方向激发。D₂的求和over-states大于H₂，故忽略热激发会低估选择性。吸附等温线表明RB模型使半覆盖压力显著左移，且D₂较H₂更易吸附（20 K时间距约300倍）。选择性方面：PsAt模型仅20 K时略超6；RB模型考虑正常混合物（n-H₂: 75% ortho-J=1, 25% para-J=0；n-D₂: 67% para-J=0, 33% ortho-J=1）下20–40 K选择性均高于工业标准（20 K达~15）；若只考虑even-J物种（假设自旋转化平衡），选择性进一步提升（20 K约~100）。此外，[H–Cor–Na]⁺对H₂的ortho/para异构体具分离能力（S_ortho/para(H₂) 20 K约30）。

3.2 The case of [CCCor–Na]

中性更大PAH [CCCor–Na]结合能略弱于[H–Cor–Na]⁺（RB-even-J H₂E₀=-1097.9 K，D₂-1239.2 K，ΔE=141.3 K），但D₂/H₂基态能差更大，且激发态能隙更窄致热激发贡献更显著。波函数分析显示H₂质心在Na上方呈C₆对称环形离域，J=odd块M_J分量分裂较大。吸附等温线与[H–Cor–Na]⁺相似，但RB-normal混合物的D₂/H₂选择性略高于前者（20 K约18，仍远超6）。Ortho/para-H₂分离选择性亦更高（20 K约59）。表明中性延伸石墨烯模型同样具备优良同位素筛分性能，且未必需基底带净电荷。

讨论与结论翻译

研究表明，钠修饰PAHs可用于氢同位素（D₂/H₂及ortho/para-H₂/D₂）的吸附分离。通过电子结构校准的RB与PsAt势能面及求解吸附分子薛定谔方程获得束缚态，发现PsAt模型因忽略各向异性不足以准确预测等温线与选择性；忽略热激发态得选择性下界。在单分子/金属位点假设下，两种Na修饰PAHs对正常混合D₂/H₂的低压（20–40 K）选择性均高于工业要求值6（circumcircumcoronene-Na于20 K可达约48）。更大PAH显示更高选择性，说明有效色谱固定相不必带净电荷。若自旋异构体能达平衡（odd-J→even-J），选择性可进一步提升。此类材料是优于未修饰碳材料的氘氢分离候选介质，后续将研究多分子共吸附（operando条件）对选择性的影响（系列第二部分）。