JVASP-30314_Si3Sb2O9
JARVIS-ID:JVASP-30314 Functional:optB88-vdW Primitive cell Primitive cell Conventional cell Conventional cell
Chemical formula:Si3Sb2O9 Formation energy/atom (eV):-2.39 a 6.92 Å α:90.0 ° a 6.92 Å α:90.0 °
Space-group :P6_3/m, 176 Relaxed energy/atom (eV):-5.6513 b 6.92 Å β:90.0 ° b 6.92 Å β:90.0 °
Calculation type:Bulk SCF bandgap (eV):4.23 c 9.324 Å γ:120.0 ° c 9.324 Å γ:120.0 °
Crystal system:hexagonal Point group:6/m Density (gcm-3):4.05 Volume (3):386.7 nAtoms_prim:28 nAtoms_conv:28
Download input files

Convergence [Reference]

Calculations are done using VASP software [Source-code]. Convergence on KPOINTS [Source-code] and ENCUT [Source-code] is done with respect to total energy of the system within 0.001 eV tolerance. Please note convergence on KPOINTS and ENCUT is generally done for target properties, but here we assume energy-convergence with 0.001 eV should be sufficient for other properties also. The points on the curves are obtained with single-point calculation (number of ionic steps, NSW=1 ). However, for very accurate calculations, NSW>1 might be needed.


Structural analysis [Reference]

The following shows the X-ray diffraction (XRD)[Source-code] pattern and the Radial distribution function (RDF) plots [Source-code]. XRD peaks should be comparable to experiments for bulk structures. Relative intensities may differ. For mono- and multi-layer structures , we take the z-dimension during DFT calculation for XRD calculations, which may differ from the experimental set-up.


Electronic structure [Reference]

The following shows the electronic density of states and bandstructure [Source-code]. DFT is generally predicted to underestimate bandgap of materials. Accurate band-gaps are obtained with higher level methods (with high computational requirement) such as HSE, GW , which are under progress. If available, MBJ data should be comparable to experiments also. Total DOS, Orbital DOS and Element dos [Source-code] buttons are provided for density of states options. Energy is rescaled to make Fermi-energy zero. In the bandstructure plot [Source-code], spin up is shown with blue lines while spin down are shown with red lines. Non-degenerate spin-up and spin-down states (if applicable) would imply a net orbital magnetic moment in the system. Fermi-occupation tolerance for bandgap calculation is chosen as 0.001.

High-symmetry kpoints based bandgap (eV): 4.23I


Electrostatic potential [Reference]

The following plot shows the plane averaged electrostatic potential (ionic+Hartree) along x, y and z-directions. The red line shows the Fermi-energy while the green line shows the maximum value of the electrostatic potential. For slab structures (with vacuum along z-direction), the difference in these two values can be used to calculate work-function of the material.


Optoelectronic properties Semi-local [Reference]

Incident photon energy dependence of optical is shown below [Source-code]. Only interband optical transitions are taken into account.Please note the underestimatation of band-gap problem with DFT will reflect in the spectra as well. For very accurate optical properties GW/BSE calculation would be needed, which is yet to be done because of their very high computational cost. Optical properties for mono-/multi-layer materials were rescaled with the actual thickness to simulation z-box ratio. Absorption coeffiecient is in cm-1 unit. Also, ionic contributions were neglected.

Dense k-mesh based bandgap is : 4.2296 eV

Static real-parts of dielectric function in x,y,z: 3.62,3.62,3.31


Optoelectronic properties METAGGA-MBJ [Reference]

Single point DFT calculation was carried out with meta-gga MBJ potential [Source-code]. This should give reasonable bandgap, and optical properties assuming the calculation was properly converged. Incident photon energy dependence of optical is shown below. Only interband optical transitions are taken into account. Also, ionic contributions were neglected.

MBJ bandgap is : 5.7226 eV

Static real-parts of dielectric function in x,y,z: 2.79,2.79,2.55


Finite-difference: elastic tensor and derived phonon properties [Reference]

Elastic tensor calculated for the conventional cell of the system with finite-difference method [Source-code]. For bulk structures, elastic constants are given in GPa unit . For layered materials, the elastic constants are rescaled with respect to vacuum padding (see the input files) and the units for elastic coefficients are in N/m . Phonons obtained [Source-code] from this calculation are also shown.

WARNING: Please note we provide finite-size cell phonons only. At least 1.2 nm x1.2 nm x1.2 nm size cell or more is generally needed for obtaining reliable phonon spectrum, but we take conventional cell of the structure only. For systems having primitive-cell phonon representation tables, I denotes infrared activity and R denotes Raman active modes (where applicabale). Selection of particular q-point mesh can give rise to unphysical negative modes in phonon density of states and phonon bandstructre. The minimum thermal conductivity was calculated using elastic tensor information following Clarke and Cahill formalism.

Voigt-bulk modulus (KV): 47.63 GPa, Voigt-shear modulus (GV): 41.56 GPa

Reuss-bulk modulus (KR): 45.5 GPa, Reuss-shear modulus (GR): 36.81 GPa

Poisson's ratio: 0.17, Elastic anisotropy parameter: 0.69

Clarke's lower limit of thermal conductivity (W/(m.K)): 0.99

Cahill's lower limit of thermal conductivity (W/(m.K)): 1.08

Elastic tensor
123.5 6.7 25.7 -0.0 0.0 0.0
6.7 123.5 25.7 -0.0 0.0 -0.0
25.7 25.7 65.5 0.0 -0.0 0.0
-0.0 -0.0 0.0 58.4 0.0 0.0
0.0 0.0 0.0 0.0 32.3 -0.0
0.0 -0.0 0.0 0.0 -0.0 32.3

Phonon mode (cm-1)
-0.08
-0.05
-0.01
79.78
79.78
81.98
85.21
85.21
96.98
96.98
118.67
118.67
129.4
129.4
139.12
150.02
158.59
160.56
185.47
185.47
193.49
193.49
196.77
196.77
207.08
207.08
219.74
241.69
247.7
266.23
277.5
277.5
281.66
281.66
301.38
319.38
334.3
334.3
345.49
357.81
357.81
359.11
359.11
363.1
363.1
369.22
414.09
414.09
427.32
427.32
451.7
471.8
494.87
494.87
496.4
496.4
544.41
548.38
548.59
550.2
628.74
636.12
722.81
722.81
732.83
732.83
897.61
897.61
906.8
906.8
910.09
910.09
921.58
921.58
949.79
958.91
960.84
996.89
1008.82
1008.82
1011.96
1011.96
1016.45
1050.42

Point group

point_group_type: 2/m

Visualize Phonons here
Phonon mode (cm-1) Representation
-0.08
-0.0781789372
-0.05
-0.0503220088
-0.01
-0.0125676505
79.78
79.780367609
81.98
81.9767819275
85.21
85.2085187259
96.98
96.9847022395
118.67
118.665636535
129.4
129.403888504
139.12
139.117737546
150.02
150.021074568
158.59
158.594706553
160.56
160.564305095
185.47
185.467484872
193.49
193.490942119
196.77
196.766608689
207.08
207.076464085
219.74
219.735023448
241.69
241.694792493
247.7
247.69993689
266.23
266.232045488
277.5
277.500414663
281.66
281.658227274
301.38
301.379553648
319.38
319.377274668
334.3
334.301337138
345.49
345.486580502
357.81
357.811193851
359.11
359.109004675
363.1
363.10352331
369.22
369.2183787
414.09
414.094081908
427.32
427.320433553
451.7
451.701413146
471.8
471.801397277
494.87
494.868169406
496.4
496.401209959
544.41
544.410339658
548.38
548.377927425
548.59
548.59245956
550.2
550.204914066
628.74
628.737471318
636.12
636.121904725
722.81
722.810395459
732.83
732.829591067
897.61
897.611440671
906.8
906.800185947
910.09
910.085169941
921.58
921.580721611
949.79
949.790002774
958.91
958.914717659
960.84
960.838690537
996.89
996.888760397
1008.82
1008.82219037
1011.96
1011.95954435
1016.45
1016.44811399
1050.42
1050.41915404

See also

Links to other databases or papers are provided below


mp-769028

MP link
mp-769028

NIST Disclaimer