晶格系统多尺度计算的一种新方法

【摘要】：本文提出了一种基于PML算法的多尺度模拟方法。在本文的计算工作中,将被研究的区域分成三个部分:MD(Molecular Dynamics)区域、MC(MacrosCopic)区域、以及PML区域。粗网格方程用于整个区域中,而MD计算仅限于MD+PML区域。PML区域用于吸收MD区域与MC区域之间数值边界上的人为反射。对于粗网格计算,我们使用一种由匹配微分算子方法得到的方程作为粗网格的控制方程,理论上说,该粗网格方程可以达到直至谱精度任意精度。本文使用该算法模拟几个一维晶格系统,取得了比较满意的结果。