3D plot of the charge density of diamond using ELK, Octave an VESTA
In this project we use ELK to simulate the spatial distribution of electronic charge density in diamond. 3D distribution can then be plotted using Matlab/Octave.
Semiconductor material often described by diamond crystal lattice (or the related zincblende lattice). Diamond lattice can be visualized as a FCC cubic lattice with 4 additional atoms quarter way along the 4 body diagonals. These 4 additional atoms can actually be seen as part of another FCC cubic lattice. Therefore, both diamond lattice as well as zincblende lattice can also be visualized as two interpenetrating FCC lattices.
VESTA is a 3D visualization program for structural models, volumetric data such as electron/nuclear densities, and crystal morphologies. VESTA is a successor to two 3D visualization programs, VICS and VEND, in the VENUS (Visualization of Electron/NUclear and Structures) software package.
For more details see VESTA .
GNU Octave is software featuring a high-level programming language, primarily intended for numerical computations. Octave helps in solving linear and nonlinear problems numerically, and for performing other numerical experiments using a language that is mostly compatible with Matlab. It may also be used as a batch-oriented language. Since it is part of the GNU Project, it is free software under the terms of the GNU General Public License.
For more details see Octave .
Elk is an all-electron full-potential linearized augmented-plane-wave (FP-LAPW) code for determining the properties of crystalline solids. It was developed originally at the Karl-Franzens-Universit¨at Graz as part of the EXCITING EU Research and Training Network project. The guiding philosophy during the implementation of the code was to keep it as simple as possible for both users and developers without compromising on its capabilities.
Using Elk on Kogence
For more details see Elk FP-LAPW .
Model File Tree
task 33 is to calculate the charge density. Charge density is output in the file
RHO3D.OUT. This file is a simple 4 column file. First three columns are coordinates while the fourth column is the charge density. We use Matlab to extract an isosurface and then plot this data in 3D. Script
ChargeDensityPlot.m is used to massage the data. We have to first interpolate the data on a regular grid.
Following are some related models available for cloning/copying by anyone:
- Fermi surface plot of aluminium (Al) - Copy
- Electronic band structure of AlAs using ELK test
- 3D plot of the charge density of diamond using ELK, Octave an VESTA
- Octave Charting Capability Demonstration (2)
- 3D plot of the charge density of diamond using ELK and OpenDX sdClone
- A 3D plot of the charge density of diamond using ELK and OpenDX sdClone
- Fermi surface plot of aluminium (Al)
- 3D plot of the charge density of diamond using ELK and OpenDX
- Octave Charting Capability Demonstration
- Electronic band structure of AlAs using ELK
- Calculation of ground state of Aluminum using ELK
- Fermi surface of MgB2
- Fast converging ELK simulation of electronic states in the ferromagnetic bcc iron.
- ELK simulation of electron-phonon coupling parameters and superconducting transition temperature T c in Niobium.
Click on the category links at the bottom of this page to navigate to a full list of simulation models in similar subject area or similar computational methodology.