No Arabic abstract
We perform a comparative study of the free energies and the density distributions in hard sphere crystals using Monte Carlo simulations and density functional theory (employing Fundamental Measure functionals). Using a recently introduced technique (Schilling and Schmid, J. Chem. Phys 131, 231102 (2009)) we obtain crystal free energies to a high precision. The free energies from Fundamental Measure theory are in good agreement with the simulation results and demonstrate the applicability of these functionals to the treatment of other problems involving crystallization. The agreement between FMT and simulations on the level of the free energies is also reflected in the density distributions around single lattice sites. Overall, the peak widths and anisotropy signs for different lattice directions agree, however, it is found that Fundamental Measure theory gives slightly narrower peaks with more anisotropy than seen in the simulations. Among the three types of Fundamental Measure functionals studied, only the White Bear II functional (Hansen-Goos and Roth, J. Phys.: Condens. Matter 18, 8413 (2006)) exhibits sensible results for the equilibrium vacancy concentration and a physical behavior of the chemical potential in crystals constrained by a fixed vacancy concentration.
In materials science the phase field crystal approach has become popular to model crystallization processes. Phase field crystal models are in essence Landau-Ginzburg-type models, which should be derivable from the underlying microscopic description of the system in question. We present a study on classical density functional theory in three stages of approximation leading to a specific phase field crystal model, and we discuss the limits of applicability of the models that result from these approximations. As a test system we have chosen the three--dimensional suspension of monodisperse hard spheres. The levels of density functional theory that we discuss are fundamental measure theory, a second-order Taylor expansion thereof, and a minimal phase-field crystal model. We have computed coexistence densities, vacancy concentrations in the crystalline phase, interfacial tensions and interfacial order parameter profiles, and we compare these quantities to simulation results. We also suggest a procedure to fit the free parameters of the phase field crystal model.
We report numerical calculations of the concentration of interstitials in hard-sphere crystals. We find that, in a three-dimensional fcc hard-sphere crystal at the melting point, the concentration of interstitials is 2 * 10^-8. This is some three orders of magnitude lower than the concentration of vacancies. A simple, analytical estimate yields a value that is in fair agreement with the numerical results.
The study of zinc oxide, within the homogeneous electron gas approximation, results in overhybridization of zinc $3d$ shell with oxygen $2p$ shell, a problem shown for most transition metal chalcogenides. This problem can be partially overcome by using LDA+$U$ (or, GGA+$U$) methodology. However, in contrast to the zinc $3d$ orbital, Hubbard type correction is typically excluded for the oxygen $2p$ orbital. In this work, we provide results of electronic structure calculations of an oxygen vacancy in ZnO supercell from ab initio perspective, with two Hubbard type corrections, $U_{mathrm{Zn}-3d}$ and $U_{mathrm{O}-2p}$. The results of our numerical simulations clearly reveal that the account of $U_{mathrm{O}-2p}$ has a significant impact on the properties of bulk ZnO, in particular the relaxed lattice constants, effective mass of charge carriers as well as the bandgap. For a set of validated values of $U_{mathrm{Zn}-3d}$ and $U_{mathrm{O}-2p}$ we demonstrate the appearance of a localized state associated with the oxygen vacancy positioned in the bandgap of the ZnO supercell. Our numerical findings suggest that the defect state is characterized by the highest overlap with the conduction band states as obtained in the calculations with no Hubbard-type correction included. We argue that the electronic density of the defect state is primarily determined by Zn atoms closest to the vacancy.
Potassium intercalation in graphite is investigated by first-principles theory. The bonding in the potassium-graphite compound is reasonably well accounted for by traditional semilocal density functional theory (DFT) calculations. However, to investi
The effect of lithium vacancies in the hexagonal structure of $alpha-$Li$_3$N, is studied within the framework of density functional theory. Vacancies ($square$) substituting for lithium in $alpha-$Li$_2$(Li$_{1-x}square_x$)N are treated within the coherent potential approximation as alloy components. According to our results long range N($p$)-ferromagnetism ($sim 1$ $mu_B$) sets in for vacancy substitution within the [Li$_2$N] layers ($x ge 0.7$) with no significant change in unit cell dimensions. By total energies differences we established that in-plane exchange couplings are dominant. Vacancies substituting inter-plane Li, leads to a considerable structural collapse ($c/a approx 0.7$) and no magnetic moment formation.