Jr. Chelikowsky et al., Ab initio simulations of liquid semiconductors using the pseudopotential-density functional method, J PHYS-COND, 13(41), 2001, pp. R817-R854
One of the most difficult problems in condensed matter physics is describin
g the microscopic nature of the liquid state. Owing to the dynamical nature
of the liquid state, it is not possible to discuss a particular microscopi
c structure; only ensemble averages can be specified. Such averages can be
performed via well crafted molecular dynamics simulations: the length of th
e simulation, the size of the ensemble and the nature of the interatomic fo
rces must all be carefully analysed. Historically, a problematic issue in d
oing such simulations is that of how to describe the interatomic forces in
the liquid state. This matter is especially challenging for the melt of sem
iconductors, such as silicon or gallium arsenide, where the chemical bond c
ontains a strong covalent component. It is difficult to use pairwise intera
tomic potentials in such cases. Although many-body potentials can be utiliz
ed for simulations of these materials, one must map quantum phenomena such
as hybridization onto classical interatomic potentials. This mapping is com
plex and difficult. In this review, we illustrate how one can avoid this pr
oblem by utilizing quantum forces to simulate liquids. Our focus is on the
pseudopotential-density functional method. Within the pseudopotential metho
d, only the valence electrons are explicitly treated and within the density
functional theory, exchange and correlation terms are mapped onto an effec
tive one-electron potential. These two approximations allow one to extract
quantum forces at every time step of the simulation. The pseudopotential-de
nsity functional method is highly accurate and well tested for semiconducto
rs in the solid state, but has only recently been applied to liquids. In th
is review, we illustrate this approach for a number of semiconducting liqui
ds such as liquid Si, Ge, GaAs, CdTe and GeTe. For these liquids, we will p
resent results for the microstructure, the dynamical properties such as the
diffusion constants and the electronic properties such as the conductivity
.