一般域中分子动力学的稳定吸收边界条件
计算物理
2017-10-27 v1
摘要
提出了一种用于分子动力学模拟的新型吸收边界条件。在谐波近似下，晶体固体的精确边界条件表示为动态的 Dirichlet-to-Neumann（DtN）映射。它将边界处原子的位移与这些原子上的牵引力联系起来。该 DtN 映射对具有一般几何的域均有效。为避免计算动态 DtN 映射的时间卷积，我们在 Laplace 域中用有理函数逼近相关核函数。逼近中的参数由插值确定。将给出零阶、一阶与二阶逼近的显式形式。建立了补充这些吸收边界条件的分子动力学模型的稳定性。通过两个数值模拟验证了方法的有效性。
引用
@article{arxiv.1710.09520,
title = {Stable absorbing boundary conditions for molecular dynamics in general domains},
author = {Xiaojie Wu and Xiantao Li},
journal= {arXiv preprint arXiv:1710.09520},
year = {2017}
}
备注
25 pages, 4 figures