JVASP-21551_K2BeF4
JARVIS-ID:JVASP-21551 Functional:optB88-vdW Primitive cell Primitive cell Conventional cell Conventional cell
Chemical formula:K2BeF4 Formation energy/atom (eV):-3.287 a 5.656 Å α:90.0 ° a 5.656 Å α:90.0 °
Space-group :Pnma, 62 Relaxed energy/atom (eV):-3.1679 b 7.316 Å β:90.0 ° b 7.316 Å β:90.0 °
Calculation type:Bulk SCF bandgap (eV):7.491 c 9.867 Å γ:90.0 ° c 9.867 Å γ:90.0 °
Crystal system:orthorhombic Point group:mmm Density (gcm-3):2.66 Volume (3):408.29 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): 7.491D


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 : 7.491 eV

Static real-parts of dielectric function in x,y,z: 2.02,2.04,2.02


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): 36.04 GPa, Voigt-shear modulus (GV): 19.25 GPa

Reuss-bulk modulus (KR): 35.34 GPa, Reuss-shear modulus (GR): 18.25 GPa

Poisson's ratio: 0.28, Elastic anisotropy parameter: 0.29

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

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

Elastic tensor
57.5 24.0 18.6 -0.0 -0.0 -0.0
24.0 73.6 23.9 -0.0 0.0 0.0
18.6 23.9 60.3 -0.0 0.0 0.0
-0.0 -0.0 -0.0 20.9 -0.0 -0.0
-0.0 0.0 -0.0 -0.0 21.6 0.0
-0.0 0.0 0.0 -0.0 0.0 12.1

Phonon mode (cm-1)
-0.53
-0.17
-0.06
45.48
47.48
58.2
58.69
63.15
72.91
73.45
75.47
90.68
91.54
94.92
97.84
99.83
102.89
104.1
105.9
105.93
106.0
111.9
117.55
123.32
124.88
126.61
126.76
127.64
129.42
133.45
135.28
139.43
145.18
146.43
147.15
147.59
156.1
157.95
159.73
160.5
162.19
168.6
174.83
182.56
188.54
194.5
195.48
197.06
244.39
247.86
254.8
254.83
255.89
260.96
265.54
271.69
370.36
370.75
371.17
371.71
372.03
372.66
376.97
377.92
382.81
384.88
385.5
389.97
552.87
554.58
555.3
555.53
794.42
799.57
809.97
810.83
811.93
812.88
817.11
825.38
840.49
845.91
861.25
867.67

Point group

point_group_type: mm2

Visualize Phonons here
Phonon mode (cm-1) Representation
-0.53
-0.5267612588
-0.17
-0.1732131426
-0.06
-0.0565014555
45.48
45.4752395168
47.48
47.482642258
58.2
58.2043216301
58.69
58.6853048237
63.15
63.1471928007
72.91
72.9145484274
73.45
73.4452438374
75.47
75.4713994066
90.68
90.6771557902
91.54
91.5434224546
94.92
94.9233826022
97.84
97.8425367543
99.83
99.832698005
102.89
102.894022038
104.1
104.096309091
105.9
105.902283688
105.93
105.925060626
106.0
105.998354401
111.9
111.899977725
117.55
117.546458543
123.32
123.324425001
124.88
124.878273574
126.61
126.614932621
126.76
126.756326074
127.64
127.636036325
129.42
129.417408453
133.45
133.45306285
135.28
135.277503211
139.43
139.430828811
145.18
145.1782008
146.43
146.432909958
147.15
147.151774051
147.59
147.590573568
156.1
156.100585549
157.95
157.954345472
159.73
159.734826036
160.5
160.499086249
162.19
162.190646974
168.6
168.598367802
174.83
174.832016893
182.56
182.563947248
188.54
188.542277384
194.5
194.501862953
195.48
195.482399484
197.06
197.058507061
244.39
244.390445265
247.86
247.856583877
254.8
254.798280428
254.83
254.830707164
255.89
255.887552102
260.96
260.959824453
265.54
265.538164972
271.69
271.694943611
370.36
370.361845884
370.75
370.754618575
371.17
371.166509617
371.71
371.707287493
372.03
372.03263664
372.66
372.65736749
376.97
376.971721017
377.92
377.918311785
382.81
382.813089802
384.88
384.878633715
385.5
385.495510343
389.97
389.965176195
552.87
552.871934853
554.58
554.577537761
555.3
555.296769176
555.53
555.534484644
794.42
794.421436064
799.57
799.572825083
809.97
809.97100065
810.83
810.829652609
811.93
811.925873026
812.88
812.879774637
817.11
817.114052809
825.38
825.377810675
840.49
840.486880808
845.91
845.912151665
861.25
861.251519682
867.67
867.670318667

Thermoelectric properties [Reference]

Thermoelectric properties are calculated using BoltzTrap code [Source-code]. Electron and hole mass tensors (useful for semiconductors and insulators mainly)are given at 300 K [Source-code]. Following plots show the Seebeck coefficient and ZT factor (eigenvalues of the tensor shown) at 300 K along three different crystallographic directions. Seebeck coefficient and ZT plots can be compared for three different temperatures available through the buttons given below. Generally very high Kpoints are needed for obtaining thermoelectric properties. We assume the Kpoints obtained from above convergence were sufficient [Source-code].

WARNING: Constant relaxation time approximation (10-14 s) and only electronic contribution to thermal conductivity were utilized for calculating ZT.

Electron mass tensor (me unit)

0.85 0.0 0.0
0.0 0.43 -0.0
0.0 -0.0 0.27

Hole mass tensor (me unit)

13.82 0.0 0.0
0.0 17.4 0.0
0.0 0.0 24.07

n-& p-type Seebeck coeff. (µV/K), power-factor (µW/(mK2)), conductivity (1/(*m)), zT (assuming lattice part of thermal conductivity as 1 W/(mK)) at 600K and 1020 cm-3 doping. For mono/multi-layer materials consider Seebeck-coeff only.)

Property xx yy zz
n-Seebeck -132.83 -107.09 -105.72
n-PowerFactor 607.42 725.16 1030.06
n-Conductivity 41102.1 54348.18 89821.64
n-ZT 0.22 0.28 0.3
p-Seebeck 439.31 466.22 482.1
p-PowerFactor 80.66 182.79 229.08
p-Conductivity 417.93 840.95 985.67
p-ZT 0.05 0.11 0.14

Magnetic moment [Reference]

The orbital magnetic moment was obtained after SCF run. This is not a DFT+U calculation, hence the data could be used to predict zero or non-zero magnetic moment nature of the material only.

Total magnetic moment: 0.0 μB

Magnetic moment per atom: 0.0 μB

Magnetization
Elementsspdtot
K-0.00.0-0.00.0
K-0.0-0.0-0.0-0.0
K-0.0-0.0-0.0-0.0
K-0.00.0-0.00.0
K-0.0-0.00.0-0.0
K-0.0-0.00.0-0.0
K-0.0-0.00.0-0.0
K-0.0-0.00.0-0.0
Be0.0-0.00.00.0
Be0.0-0.00.0-0.0
Be-0.00.00.0-0.0
Be-0.0-0.00.0-0.0
F-0.0-0.00.0-0.0
F-0.0-0.00.0-0.0
F-0.0-0.00.0-0.0
F-0.0-0.00.0-0.0
F-0.0-0.00.0-0.0
F-0.0-0.00.0-0.0
F-0.0-0.00.0-0.0
F-0.0-0.00.0-0.0
F0.00.00.00.0
F0.00.00.00.0
F0.00.00.00.0
F0.00.00.00.0
F0.00.00.00.0
F0.00.00.00.0
F0.00.00.00.0
F0.00.00.00.0

See also

Links to other databases or papers are provided below


mp-5185

ICSD-ID: 153080

AFLOW link

MP link
mp-5185

NIST Disclaimer