Kress J D, Cohen James S, Horner D A, Lambert F, Collins L A
Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA.
Phys Rev E Stat Nonlin Soft Matter Phys. 2010 Sep;82(3 Pt 2):036404. doi: 10.1103/PhysRevE.82.036404. Epub 2010 Sep 14.
We have calculated viscosity and mutual diffusion of deuterium-tritium (DT) in the warm, dense matter regime for densities from 5 to 20 g/cm{3} and temperatures from 2 to 10 eV, using both finite-temperature Kohn-Sham density-functional theory molecular dynamics (QMD) and orbital-free molecular dynamics (OFMD). The OFMD simulations are in generally good agreement with the benchmark QMD results, and we conclude that the simpler OFMD method can be used with confidence in this regime. For low temperatures (3 eV and below), one-component plasma (OCP) model simulations for diffusion agree with the QMD and OFMD calculations, but deviate by 30% at 10 eV. In comparison with the QMD and OFMD results, the OCP viscosities are not as good as for diffusion, especially for 5 g/cm{3} where the temperature dependence is significantly different. The QMD and OFMD reduced diffusion and viscosity coefficients are found to depend largely, though not completely, only on the Coulomb coupling parameter Γ , with a minimum in the reduced viscosity at Γ≈25 , approximately the same position found in the OCP simulations. The QMD and OFMD equations of state (pressure) are also compared with the hydrogen two-component plasma model.
我们利用有限温度的科恩-沈密度泛函理论分子动力学(QMD)和无轨道分子动力学(OFMD),计算了氘-氚(DT)在密度为5至20 g/cm³、温度为2至10 eV的温稠密物质状态下的粘度和互扩散系数。OFMD模拟结果与基准QMD结果总体上吻合良好,我们得出结论,在该状态下可以放心使用更简单的OFMD方法。对于低温(3 eV及以下),单组分等离子体(OCP)模型的扩散模拟结果与QMD和OFMD计算结果相符,但在10 eV时偏差达30%。与QMD和OFMD结果相比,OCP粘度的模拟效果不如扩散,特别是在5 g/cm³时,其温度依赖性明显不同。发现QMD和OFMD的约化扩散系数和粘度系数在很大程度上(尽管不是完全)仅取决于库仑耦合参数Γ,在Γ≈25时约化粘度出现最小值,这与OCP模拟中的位置大致相同。还将QMD和OFMD的状态方程(压力)与氢双组分等离子体模型进行了比较。