Difference between revisions of "Silicon Slab Structural Relaxation"

From William Parker Wiki
Jump to: navigation, search
(Single Layer)
(Relaxation Quality Parameters)
 
(89 intermediate revisions by 2 users not shown)
Line 2: Line 2:
  
 
===Single Layer===
 
===Single Layer===
 +
====PBE====
  
 +
=====Pair Distribution Functions=====
 +
{| border="1"
 +
|-
 +
! PDF
 +
|-
 +
|[[File:Si.PBE.2 ML.Fd-3m.PDF.png]]
 +
|-
 +
|}
 +
 +
 +
=====Phonon Bands and Density of States=====
 +
{| border ="1"
 +
|-
 +
! Phonons
 +
|-
 +
| [[File:Si.PBE.2_ML.Fd-3m.phonons.png]]
 +
|-
 +
|}
 +
 +
====PBEsol====
 +
=====Structural Visualizations=====
 +
{| border="1"
 +
|-
 +
!
 +
! Unrelaxed Bulk Cutout Slab
 +
! Relaxed Bulk Cutout Slab (0)
 +
! Displacements Built-in Slab (16)
 +
|-
 +
|
 +
| (image)
 +
| (image)
 +
| (image)
 +
|-
 +
| '''Total Energy (eV/atom)'''
 +
|
 +
| -90.40389544
 +
| -90.44063275
 +
|}
 +
=====Relaxation Quality Parameters=====
 
{| border="1"
 
{| border="1"
 
|-
 
|-
Line 11: Line 51:
 
! Total Force (eV/Å/atom)
 
! Total Force (eV/Å/atom)
 
! Cell Pressure (GPa)
 
! Cell Pressure (GPa)
! Resulting Mean &omega;<sub>1</sub> (THz)
+
! Resulting Mean &omega;<sub>1</sub> (cm<sup>-1</sup>)
 +
! Lowest &omega; (cm<sup>-1</sup>)
 
|-
 
|-
 
|0
 
|0
Line 20: Line 61:
 
| 0.010
 
| 0.010
 
| -0.493665
 
| -0.493665
 +
| -63.177
 
|-
 
|-
 
|1
 
|1
| -0.510   0.224   0.000
+
| 9 (-0.526   0.254   0.000)
| 2.0
+
| 1.0
| -90.40383384
+
| -90.33443035
| 1.85e-03
+
| 4.50e-05
| 0.03
+
| 0.000
| -0.603469
+
| -0.718166
 +
| -59.031
 
|-
 
|-
 
|2
 
|2
| -0.510  0.224  0.000
+
| 3 (-0.001 -0.513 0.000)
| 3.0
+
| 1.0
| -90.32054431
+
| -90.40389384
| 1.98e-03
+
| 6.00e-05
| -2.96
+
| -0.000
| -1.18809
+
| -0.597389
 +
| -64.901
 
|-
 
|-
 
|3
 
|3
| 0.241  -0.475   0.000
+
| 9 (-0.510   0.223 0.000)
| 2.0
+
| 1.0
| -90.40388551
+
| -90.39834846
| 4.90e-05
+
| 4.50e-05
| 0.000
+
| -0.010
| -0.721381
+
| 0.839677
 +
| -13.822794
 
|-
 
|-
 
|4
 
|4
| -0.510   0.224   0.000
+
| 9 (-0.477   0.352   0.000)
 
| 1.0
 
| 1.0
| -90.40388552
+
| -90.37901909
| 7.80e-05
+
| 5.20e-05
| -0.01
+
| -0.010
| -0.727365
+
| -0.181455
 +
| -42.396192
 
|-
 
|-
 
|5
 
|5
|
+
| 10 (-0.492  -0.709  0.000)
|
+
| 1.0
|
+
| -90.47275885
|
+
| 4.10e-05
|
+
| -0.050
|
+
| -0.615533
 +
| -91.254540
 
|-
 
|-
 +
|6
 +
| 10 (-0.610  -0.836  0.000)
 +
| 1.0
 +
| -90.47274746
 +
| 4.00e-05
 +
| -0.050
 +
| -0.638118
 +
| -90.184907
 +
|-
 +
|7
 +
| 10 (-0.610  -0.836  0.000)
 +
| 1.0
 +
| -90.40487079
 +
| 4.30e-05
 +
| -0.030
 +
| 0.869615
 +
| -17.865773
 +
|-
 +
|8
 +
| 1 (0.000 0.000 0.000)
 +
| 1.0
 +
| -90.40486887
 +
| 4.80e-05
 +
| -0.020
 +
| 0.901969
 +
| -17.314176
 +
|-
 +
|9
 +
| 4 (0.236  0.0193 0.000)
 +
| 1.0
 +
| -90.41832411
 +
| 4.90e-05
 +
| -0.130
 +
| 0.187767
 +
| -27.658783
 +
|-
 +
|10
 +
| 4 (0.273  0.0425 0.000)
 +
| 1.0
 +
| -90.41833588
 +
| 4.70e-05
 +
| -0.110
 +
| 0.260652
 +
| -26.044145
 +
|-
 +
|11
 +
| 1 (0.000  0.000 0.000)
 +
| 1.0
 +
| -90.41835466
 +
| 4.40e-05
 +
| -0.110
 +
| 0.536324
 +
| -21.519492
 +
|-
 +
|12
 +
| 4 (0.273  0.0427  0.000)
 +
| 1.0
 +
| -90.41833514
 +
| 5.20e-05
 +
| -0.130
 +
| 0.486934
 +
| -20.474298
 +
|-
 +
|13
 +
| 4 (0.273  0.0427  0.000)
 +
| 1.0
 +
| -90.41833714
 +
| 4.10e-05
 +
| -0.150
 +
| 0.296044
 +
| -25.548234
 +
|-
 +
|14
 +
| 4 (0.273 0.0428 0.000)
 +
| 1.0
 +
| -90.41833302
 +
| 5.40e-05
 +
| -0.140
 +
| 0.287938
 +
| -25.461523
 +
|-
 +
|15
 +
| 4 (0.273 0.0428 0.000)
 +
| 1.0
 +
| -90.41833791
 +
| 5.40e-05
 +
| -0.120
 +
| 0.291452
 +
| -23.679543
 +
|-
 +
|16
 +
| 5 (0.353  0.429 0.000)
 +
| 1.0
 +
| -90.44063275
 +
| 4.90e-05
 +
| -0.090
 +
| 1.35723
 +
| -21.043018
 +
|-
 +
|fix_ab (i.)
 +
| -
 +
| -
 +
| -86.42339154
 +
| 6.99e-03
 +
| 0.130
 +
|
 +
|
 +
|-
 +
|fix_ab (f.)
 +
| -
 +
| -
 +
| -90.66327410
 +
| 2.09e-04
 +
| -0.000
 +
| -0.369114
 +
| -74.550002
 
|}
 
|}
  
====Pair Distribution Functions====
+
=====Pair Distribution Functions=====
 +
''a''<sub>PBEsol</sub> = 5.495 &Aring; ; ''r''<sub>NN, bulk, PBEsol</sub> = 2.379 &Aring;
 
{| border="1"
 
{| border="1"
 
|-
 
|-
! max7
+
! Step
! buildin1
+
! PDF
! buildin2
+
! buildin3
+
! buildin4
+
 
|-
 
|-
! celldm(1) = 10.383624
+
|0
! celldm(1) = 10.383624
+
|[[File:Si.PBEsol.1_ML.Fd-3m.max7.PDF.png]]
! celldm(1) = 10.383624
+
!
+
!
+
 
|-
 
|-
| -1.88312005, 0.3279543, 10.18850117
+
| 1
| -1.74902282  0.34899568  10.26696606
+
| [[File:Si.PBEsol.1_ML.Fd-3m.buildin1.PDF.png]]
| -1.21307207  -1.74352054  11.04130821
+
|
+
|
+
 
|-
 
|-
| 0.21413732  5.63198628  13.25787654
+
| 2
| 0.32032442  5.68006765  13.40047589
+
| [[File:Si.PBEsol.1_ML.Fd-3m.buildin2.PDF.png]]
| 2.10358308  5.56847173  11.49515876
+
|
+
|
+
|
+
 
|-
 
|-
| 5.41505161  1.15781353  13.59935032
+
| 3
| 5.54935756  1.20568723  13.67679658
+
| [[File:Si.PBEsol.1_ML.Fd-3m.buildin3.PDF.png]]
| 5.5521819    3.74684661  14.47954002
+
|
+
|
+
|
+
 
|-
 
|-
| 5.36943116  2.93765287  8.15649656
+
| 4
| 5.47879045  2.96313117  8.29400897
+
| [[File:Si.PBEsol.1 ML.Fd-3m.buildin4.PDF.png]]
| 5.87775732  0.49488975  8.78971383
+
|
+
|
+
|
+
 
|-
 
|-
| 1.96708931  2.27280708  11.02079432
+
| 5
| 2.06776996  2.33190651  11.14231906
+
| [[File:Si.PBEsol.1_ML.Fd-3m.buildin5.PDF.png]]
| 2.19363818  1.12970506  11.28067358
+
|
+
|
+
|
+
|-  
+
| 4.05816739  7.58542112  14.09855031
+
| 4.1872431    7.6226651  14.17736514
+
| 5.11364823  9.94778073  13.8750615
+
|
+
|
+
|
+
|-
+
| 6.82128504  4.59206945  16.12137268
+
| 6.9197841    4.62039702  16.25711261
+
| 2.52782577  6.79555857  15.87718204
+
|
+
|
+
|
+
|
+
|
+
 
|-
 
|-
| 6.77202928  6.36786956  10.68751114
+
| 6
| 6.90963514  6.403177    10.78673497
+
| [[File:Si.PBEsol.1_ML.Fd-3m.buildin6.PDF.png]]
| 6.35621187  4.73753987  10.21102013
+
|-
|  
+
| 7
|  
+
| [[File:Si.PBEsol.1_ML.Fd-3m.buildin7.PDF.png]]
|  
+
|-
|  
+
| 8
 +
| [[File:Si.PBEsol.1_ML.Fd-3m.buildin8.PDF.png]]
 +
|-
 +
| 9
 +
| [[File:Si.PBEsol.1_ML.Fd-3m.buildin9.PDF.png]]
 +
|-
 +
| 10
 +
| [[File:Si.PBEsol.1_ML.Fd-3m.buildin10.PDF.png]]
 +
|-
 +
| 11
 +
| [[File:Si.PBEsol.1_ML.Fd-3m.buildin11.PDF.png]]
 +
|-
 +
| 12
 +
| [[File:Si.PBEsol.1 ML.Fd-3m.max7.buildin12.PDF.png]]
 +
|-
 +
| 13
 +
| [[File:Si.PBEsol.1 ML.Fd-3m.buildin13.PDF.png]]
 +
|-
 +
| 14
 +
| [[File:Si.PBEsol.1 ML.Fd-3m.buildin14.PDF.png]]
 +
|-
 +
| 15
 +
| [[File:Si.PBEsol.1 ML.Fd-3m.buildin15.PDF.png]]
 +
|-
 +
| 16
 +
| [[File:Si.PBEsol.1 ML.Fd-3m.buildin16.PDF.png]]
 +
|-
 +
| fix_ab1
 +
| [[File:Si.PBEsol.1 ML.Fd-3m.fix_ab1.PDF.png]]
 
|}
 
|}
 +
 +
=====Phonon Bands and Density of States=====
 +
{| border ="1"
 +
|-
 +
! Step
 +
! Phonons
 +
|-
 +
| 0
 +
| [[File:Si.PBEsol.1_ML.Fd-3m.max7.buildin0.png]]
 +
|-
 +
| 1
 +
| [[File:Si.PBEsol.1_ML.Fd-3m.buildin1.phonons.png]]
 +
|-
 +
| 2
 +
| [[File:Si.PBEsol.1_ML.Fd-3m.buildin2.phonons.png]]
 +
|-
 +
| 3
 +
| [[File:Si.PBEsol.1_ML.Fd-3m.buildin3.phonons.png]]
 +
|-
 +
| 4
 +
| [[File:Si.PBEsol.1_ML.Fd-3m.buildin4.phonons.png]]
 +
|-
 +
| 5
 +
| [[File:Si.PBEsol.1_ML.Fd-3m.buildin5.phonons.png]]
 +
|-
 +
| 6
 +
| [[File:Si.PBEsol.1_ML.Fd-3m.buildin6.phonons.png]]
 +
|-
 +
| 7
 +
| [[File:Si.PBEsol.1_ML.Fd-3m.buildin7.phonons.png]]
 +
|-
 +
| 8
 +
| [[File:Si.PBEsol.1_ML.Fd-3m.buildin8.phonons.png]]
 +
|-
 +
| 9
 +
| [[File:Si.PBEsol.1_ML.Fd-3m.buildin9.phonons.png]]
 +
|-
 +
| 10
 +
| [[File:Si.PBEsol.1_ML.Fd-3m.buildin10.phonons.png]]
 +
|-
 +
| 11
 +
| [[File:Si.PBEsol.1_ML.Fd-3m.buildin11.phonons.png]]
 +
|-
 +
| 12
 +
| [[File:Si.PBEsol.1 ML.Fd-3m.max7.buildin12.phonons.png]]
 +
|-
 +
| 13
 +
| [[File:Si.PBEsol.1 ML.Fd-3m.buildin13.phonons.png]]
 +
|-
 +
| 14
 +
| [[File:Si.PBEsol.1 ML.Fd-3m.buildin14.phonons.png]]
 +
|-
 +
| 15
 +
| [[File:Si.PBEsol.1 ML.Fd-3m.buildin15.phonons.png]]
 +
|-
 +
| 16
 +
| [[File:Si.PBEsol.1 ML.Fd-3m.buildin16.phonons.png]]
 +
|-
 +
| fix_ab1
 +
| [[File:Si.PBEsol.1 ML.Fd-3m.fix_ab1.phonons.png]]
 +
|}
 +
 +
==Relaxation Procedure==
 +
 +
# Create slab using copies of bulk structure coordinates
 +
# Relax slab using loose total force convergence (~10<sup>-3</sup> Ry bohr)
 +
# Vary celldm(1) by fixed perentage around loosely relaxed slab and relax each new structure to tighter total force convergence (~10<sup>-4</sup> Ry bohr)
 +
# Relax the minimum-energy structure from that step to the tightest force convergence (~10<sup>-5</sup> Ry bohr)
 +
# Vary celldm(3) around the relaxed minimum energy structure and relax the in-plane (x & y) coordinates only
 +
# Fit the equation of state to these energies and volumes

Latest revision as of 20:27, 2 May 2022

Diamond

Single Layer

PBE

Pair Distribution Functions
PDF
Si.PBE.2 ML.Fd-3m.PDF.png


Phonon Bands and Density of States
Phonons
Si.PBE.2 ML.Fd-3m.phonons.png

PBEsol

Structural Visualizations
Unrelaxed Bulk Cutout Slab Relaxed Bulk Cutout Slab (0) Displacements Built-in Slab (16)
(image) (image) (image)
Total Energy (eV/atom) -90.40389544 -90.44063275
Relaxation Quality Parameters
Step Displacement q (b1, b2, b3) Build-in Factor Total Energy (eV/atom) Total Force (eV/Å/atom) Cell Pressure (GPa) Resulting Mean ω1 (cm-1) Lowest ω (cm-1)
0 - - -90.40389544 6.90e-05 0.010 -0.493665 -63.177
1 9 (-0.526 0.254 0.000) 1.0 -90.33443035 4.50e-05 0.000 -0.718166 -59.031
2 3 (-0.001 -0.513 0.000) 1.0 -90.40389384 6.00e-05 -0.000 -0.597389 -64.901
3 9 (-0.510 0.223 0.000) 1.0 -90.39834846 4.50e-05 -0.010 0.839677 -13.822794
4 9 (-0.477 0.352 0.000) 1.0 -90.37901909 5.20e-05 -0.010 -0.181455 -42.396192
5 10 (-0.492 -0.709 0.000) 1.0 -90.47275885 4.10e-05 -0.050 -0.615533 -91.254540
6 10 (-0.610 -0.836 0.000) 1.0 -90.47274746 4.00e-05 -0.050 -0.638118 -90.184907
7 10 (-0.610 -0.836 0.000) 1.0 -90.40487079 4.30e-05 -0.030 0.869615 -17.865773
8 1 (0.000 0.000 0.000) 1.0 -90.40486887 4.80e-05 -0.020 0.901969 -17.314176
9 4 (0.236 0.0193 0.000) 1.0 -90.41832411 4.90e-05 -0.130 0.187767 -27.658783
10 4 (0.273 0.0425 0.000) 1.0 -90.41833588 4.70e-05 -0.110 0.260652 -26.044145
11 1 (0.000 0.000 0.000) 1.0 -90.41835466 4.40e-05 -0.110 0.536324 -21.519492
12 4 (0.273 0.0427 0.000) 1.0 -90.41833514 5.20e-05 -0.130 0.486934 -20.474298
13 4 (0.273 0.0427 0.000) 1.0 -90.41833714 4.10e-05 -0.150 0.296044 -25.548234
14 4 (0.273 0.0428 0.000) 1.0 -90.41833302 5.40e-05 -0.140 0.287938 -25.461523
15 4 (0.273 0.0428 0.000) 1.0 -90.41833791 5.40e-05 -0.120 0.291452 -23.679543
16 5 (0.353 0.429 0.000) 1.0 -90.44063275 4.90e-05 -0.090 1.35723 -21.043018
fix_ab (i.) - - -86.42339154 6.99e-03 0.130
fix_ab (f.) - - -90.66327410 2.09e-04 -0.000 -0.369114 -74.550002
Pair Distribution Functions

aPBEsol = 5.495 Å ; rNN, bulk, PBEsol = 2.379 Å

Step PDF
0 Si.PBEsol.1 ML.Fd-3m.max7.PDF.png
1 Si.PBEsol.1 ML.Fd-3m.buildin1.PDF.png
2 Si.PBEsol.1 ML.Fd-3m.buildin2.PDF.png
3 Si.PBEsol.1 ML.Fd-3m.buildin3.PDF.png
4 Si.PBEsol.1 ML.Fd-3m.buildin4.PDF.png
5 Si.PBEsol.1 ML.Fd-3m.buildin5.PDF.png
6 Si.PBEsol.1 ML.Fd-3m.buildin6.PDF.png
7 Si.PBEsol.1 ML.Fd-3m.buildin7.PDF.png
8 Si.PBEsol.1 ML.Fd-3m.buildin8.PDF.png
9 Si.PBEsol.1 ML.Fd-3m.buildin9.PDF.png
10 Si.PBEsol.1 ML.Fd-3m.buildin10.PDF.png
11 Si.PBEsol.1 ML.Fd-3m.buildin11.PDF.png
12 Si.PBEsol.1 ML.Fd-3m.max7.buildin12.PDF.png
13 Si.PBEsol.1 ML.Fd-3m.buildin13.PDF.png
14 Si.PBEsol.1 ML.Fd-3m.buildin14.PDF.png
15 Si.PBEsol.1 ML.Fd-3m.buildin15.PDF.png
16 Si.PBEsol.1 ML.Fd-3m.buildin16.PDF.png
fix_ab1 Si.PBEsol.1 ML.Fd-3m.fix ab1.PDF.png
Phonon Bands and Density of States
Step Phonons
0 Si.PBEsol.1 ML.Fd-3m.max7.buildin0.png
1 Si.PBEsol.1 ML.Fd-3m.buildin1.phonons.png
2 Si.PBEsol.1 ML.Fd-3m.buildin2.phonons.png
3 Si.PBEsol.1 ML.Fd-3m.buildin3.phonons.png
4 Si.PBEsol.1 ML.Fd-3m.buildin4.phonons.png
5 Si.PBEsol.1 ML.Fd-3m.buildin5.phonons.png
6 Si.PBEsol.1 ML.Fd-3m.buildin6.phonons.png
7 Si.PBEsol.1 ML.Fd-3m.buildin7.phonons.png
8 Si.PBEsol.1 ML.Fd-3m.buildin8.phonons.png
9 Si.PBEsol.1 ML.Fd-3m.buildin9.phonons.png
10 Si.PBEsol.1 ML.Fd-3m.buildin10.phonons.png
11 Si.PBEsol.1 ML.Fd-3m.buildin11.phonons.png
12 Si.PBEsol.1 ML.Fd-3m.max7.buildin12.phonons.png
13 Si.PBEsol.1 ML.Fd-3m.buildin13.phonons.png
14 Si.PBEsol.1 ML.Fd-3m.buildin14.phonons.png
15 Si.PBEsol.1 ML.Fd-3m.buildin15.phonons.png
16 Si.PBEsol.1 ML.Fd-3m.buildin16.phonons.png
fix_ab1 Si.PBEsol.1 ML.Fd-3m.fix ab1.phonons.png

Relaxation Procedure

  1. Create slab using copies of bulk structure coordinates
  2. Relax slab using loose total force convergence (~10-3 Ry bohr)
  3. Vary celldm(1) by fixed perentage around loosely relaxed slab and relax each new structure to tighter total force convergence (~10-4 Ry bohr)
  4. Relax the minimum-energy structure from that step to the tightest force convergence (~10-5 Ry bohr)
  5. Vary celldm(3) around the relaxed minimum energy structure and relax the in-plane (x & y) coordinates only
  6. Fit the equation of state to these energies and volumes