[论文解读] Van der Waals interactions of parallel and concentric nanotubes

Elsebeth Schröder, Per Hyldgaard|arXiv (Cornell University)|Jul 1, 2003

Graphene research and applications参考文献 11被引用 14

一句话总结

本文提出一种基于第一性原理密度泛函理论（DFT）的方法，通过引入非局域电子响应和屏蔽效应，精确描述平行排列和同轴碳纳米管之间的范德华（vdW）相互作用。该方法采用频率积分、与几何相关的相互作用模型，将原本复杂的6+1维积分简化为可处理的二维或一维数值积分，所得结合能定量结果与渐近极限下的宏观d⁻⁵标度一致。

ABSTRACT

For sparse materials like graphitic systems and carbon nanotubes the standard density functional theory (DFT) faces significant problems because it cannot accurately describe the van der Waals interactions that are essential to the carbon-nanostructure materials behavior. While standard implementations of DFT can describe the strong chemical binding within an isolated, single-walled carbon nanotube, a new and extended DFT implementation is needed to describe the binding between nanotubes. We here provide the first steps to such an extension for parallel and concentric nanotubes through an electron-density based description of the materials coupling to the electrodynamical field. We thus find a consistent description of the (fully screened) van der Waals interactions that bind the nanotubes across the low-electron-density voids between the nanotubes, in bundles and as multiwalled tubes.

研究动机与目标

• 解决标准DFT在描述碳纳米管等低电子密度稀疏材料中长程范德华（vdW）相互作用时的失效问题。
• 开发一种一致的第一性原理方法，通过将电子密度响应与电磁场耦合，计算管间vdW结合能。
• 考虑纳米管中的局域场屏蔽和各向异性介电响应，这些因素对准确预测结合能至关重要。
• 提供一种适用于平行纳米管束和同轴多壁纳米管的框架，实现与实验体系的定量比较。

提出的方法

• 使用第一性原理DFT计算单壁纳米管的径向电子密度和静态极化率。
• 利用具有频率依赖性的极化率χ₀(n(r), u, u₀)的等离子体极点近似，建模局域动态电子响应。
• 通过关系式χ_eff(u)E_applied = −χ₀(u)∇φ(s, u)推导包含局域场屏蔽效应的有效极化率张量χ_eff。
• 应用薄壁近似，将纳米管表示为径向狄拉克δ函数，将问题简化为一维径向响应函数。
• 将单位长度的vdW相互作用能表示为频率积分J_tot(R₁, R₂, u₀⁽¹⁾, u₀⁽²⁾)与角坐标空间积分的乘积。
• 利用圆柱对称性，将二维角积分简化为一个超几何函数积分，实现高效数值计算。

实验结果

研究问题

• RQ1 如何将标准DFT扩展以精确描述稀疏碳纳米结构之间非局域范德华相互作用？
• RQ2 局域场屏蔽和各向异性介电响应在决定管间结合能的强度与距离依赖性方面起什么作用？
• RQ3 几何因素（半径、间距、同轴性）如何影响纳米管体系中的vdW相互作用能？
• RQ4 相互作用能能否以一种形式表达，既能恢复已知的宏观极限（如d⁻⁵标度），又保持完全第一性原理特性？
• RQ5 对称性与解析简化在原本高维的管间vdW能量积分中能发挥多大作用？

主要发现

• 该方法仅基于第一性原理DFT输入和一致的电磁模型，成功计算了平行和同轴纳米管的单位长度vdW相互作用能。
• 对于相同平行纳米管，相互作用能可表示为包含超几何函数的单一积分，显著提升计算效率。
• 相互作用能的渐近极限遵循d⁻⁵依赖关系，与宏观伦敦理论一致，但因中空圆柱几何结构引入了高阶修正。
• 对于同轴纳米管，空间积分可解析表示为完全椭圆积分K(k)和E(k)的形式，从而获得几何因子的闭式表达。
• 计算的单位长度结合能表明，在相同中心间距下，半径更大的纳米管（曲率更低）因接触面积更大而表现出更强的吸引力。
• 该模型通过引入屏蔽和各向异性响应，优于传统Hamaker估算，减少了体相介电模型近似带来的误差。