700eV cutoff energy and 7×7×7 kpoint set is used and the optimized value is a=5.468 Å. The calculation starts from the experimentally reported result with a=5.381 A. The ‘on the fly’ generated ultrasoft pseudopotential for Si has a core radius of 1.8 Bohr(0.95 Angstroms) and was generated with 4 electrons in the valence panel with (3s2 3p2) as the electronic configuration.įirst, the optimization of bulk material is done to make sure there is no artificial stress in the model. The exchange and correlation function were calculated using the Perdew, Burke and Ernzerhof(PBE) functional described within the generalized gradient approximation(GGA). The CASTEP package is used to carry out the DFT calculation. To explain and predict such phenomena, we used plane-wave basis sets with ultrasoft pseudopotentials to perform DFT methods to calculate the energy of the surface with and without reconstruction. Such phenomenon have been broadly investigated and confirmed through Scanning Tunneling Microscope. Due to the recombination of dangling bonds, usually there will be reconstruction patterns formed at the surface. However, the real surface of Si doesn’t resemble the plane directly cut from bulk material. Silicon(Si) is a widely used substrate for molecule beam epitaxy growth due to its high quality and low cost. Reconstruction pattern of Si(001) surface This entry was posted in 2nd Post 2020 on Apby njk22. Proceedings of the National Academy of Sciences of the United States of America, 114(44), E9188–E9196. Properties of real metallic surfaces: Effects of density functional semilocality and van der Waals nonlocality. Density Functional Theory: A Practical Introduction.“First principles methods using CASTEP” Zeitschrift fuer Kristallographie 220(5-6) pp.With better computational ability, these results should be checked against varying layers for the surface calculations as well as a higher ENCUT for both the bulk calculation and the surface calculations. This result shows that the (111) surface is the preferred face for crystal growth. In previous comparisons against multiple functionals, the percentage difference was lower. The energy of the (111) surface is found to be 15% lower than the (100) surface. The results are in agreement with the known result that says that the (111) surface is energetically favorable to the (100) surface of Ag. This gave the following results for the surface energy densities: The change in energy from the previous k-point energy calculation. For all the rest of the options, the defaults were used. The convergence tolerances were chosen somewhat arbitrarily to be small: energy at 2.0e-5 eV/atom, force at 0.05 eV/Å, stress at 0.1 GPa, and displacement at 0.002 Å. Ag has the electron configuration of 1s2 2s2 2p6 3s2 3p6 3d10 4s2 4p6 4d10 5s1, and the pseudopotential treats 4s2 4p6 4d10 5s1 as the valence electrons. Using CASTEP with the GGA PBE functional and OTFG ultrasoft pseudopotentials, we run through different energy cutoffs to find an appropriate energy cutoff for energy convergence. So for both 100 and 111 surfaces, we need to calculate the energy of bulk Ag. Surface Energy Calculationįirst, in order to calculate the surface energy for a particular surface, we need to refer to the following equation found in “Density Functional Theory : A Practical Introduction”: Note that silver is an FCC crystal with a lattice constant of 4.09 Å and that this experimental result will be used in the calculations. These properties can reveal the preferred behavior in surface Ag formation.
#Materials studio dispersion correction software
In this post, we will use DFT calculations using CASTEP software to calculate the energy of bulk silver as well as relaxing a cleaved surface. We will investigate some properties of crystalline silver: the relative surface energies for the 100 and 111 cleaved faces.