IONIC AND ELECTRONIC-STRUCTURES OF LIQUID ALUMINUM FROM THE QUANTAL HYPERNETTED-CHAIN EQUATIONS COMBINED WITH THE MOLECULAR-DYNAMICS METHOD

被引:16
作者
CHIHARA, J [1 ]
KAMBAYASHI, S [1 ]
机构
[1] JAPAN ATOM ENERGY RES INST,CTR COMP & INFORMAT SYST,TOKAI,IBARAKI 31911,JAPAN
关键词
D O I
10.1088/0953-8984/6/47/005
中图分类号
O469 [凝聚态物理学];
学科分类号
070205 ;
摘要
Liquid aluminium is a simple metal which can be modelled as an electron-ion mixture with pairwise interparticle interactions. On the basis of the density functional theory, the quantal hypernetted-chain (QHNC) formulation for the electron-ion mixture provides exact expressions for the electron-ion and ion-ion radial distribution-functions (RDF), and shows that this electron-ion mixture can be treated as a quasi-one-component liquid interacting only through a pairwise interaction without many-body forces. Using the approximation that the exchange-correlation effect of the conduction electrons in the mixture is represented by the local-field correction (LFC) of the jellium model, the QHNC formulation offers a procedure for performing a first-principles molecular dynamics (MD) simulation where the effective ion-ion interaction is determined self-consistently with its liquid structure. Using the Geldart-Vosko (GV) LFC, the structure factor of liquid aluminium near its melting point is calculated by this procedure (QHNC-MD) and is in fairly good agreement with experiment, in conjunction with the bridge function and the direct correlation function. The calculated ion-ion RDF. exhibits a significant dependence on which kind of LFC is adopted in the QHNC-MD method in contrast with what is found in the case of liquid alkali metals where the LFC of the cv type and of the local-density-approximation type yield almost the same structure factors. The effective ion-ion potential from the QHNC-MD method with the GV LPC becomes a deep negative well where the potential of Dagens, Geldart and Taylor has a positive minimum. The electron-ion RDF is obtained in a consistent way with the density profile rho(r) of a neutral pseudoatom. The Ashcroft model potential with core radius r(c) = 1.12a(B) produces an electron density distribution rho(r) and ion-ion RDF almost identical with the QHNC-MD results.
引用
收藏
页码:10221 / 10236
页数:16
相关论文
共 30 条