Calculations are done using VASP software. Convergence on KPOINTS and ENCUT 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 (nuber of ionic steps,NSW=1). However, for very accurate calculations, NSW>1 might be needed.
The following shows the X-ray diffraction (XRD) pattern and the Radial distribution function (RDF) plots. XRD peaks should be comparable to experiments for bulk structures. Relative intensities may differ.
Incident photon energy dependence of optical is shown below. 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 layered materials needs to be rescaled with the actual thickness to simulation z-box ratio. Absorption coeffiecient is in cm-1 unit.
Dense k-mesh based bandgap is : 0.9809 eV
Elastic tensor calculated for the conventional cell of the system with finite-difference method. 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 from this calcuation are also shown.
WARNING: Please note this may not be the exact phonon modes of the system as we did not test the cell-size dependence of phonons yet. At least 1.2 nm x1.2 nm x1.2 nm or more is needed for obtaining reliable phonon spectrum. For systems having primitive-cell phonon representation tables, I denotes infrared activity and R denotes Raman active modes (where applicabale). The minimum thermal conductivity was calculated using elastic tensor information following Clarke and Cahill formalism.
Bulk Modulus BV 10.356 GPa
Shear Modulus GV 6.093 GPa
|
|
|||||
|---|---|---|---|---|---|
| 41.3 | 8.6 | -2.3 | -0.0 | -0.3 | -0.0 |
| 8.6 | 41.3 | -2.3 | -0.0 | 0.3 | -0.0 |
| -0.8 | -0.8 | -0.4 | -0.0 | -0.0 | -0.0 |
| -0.0 | -0.0 | -0.0 | 5.5 | 0.0 | -0.1 |
| -0.1 | 0.1 | -0.0 | 0.0 | -0.3 | -0.0 |
| -0.0 | -0.0 | -0.0 | -0.1 | -0.0 | -0.3 |
| Phonon mode (cm-1) |
|---|
| -17.6032476313 |
| -0.3041137814 |
| -0.3029098758 |
| 33.9856389929 |
| 33.9856452459 |
| 86.3302590565 |
| 86.3302617484 |
| 93.9403733061 |
| 104.236085802 |
| 119.503604525 |
| 134.939045348 |
| 185.044017829 |
| 185.044022669 |
| 186.017109891 |
| 186.017113846 |
| 208.990021537 |
| 208.990025326 |
| 216.092907748 |
| 216.092913629 |
| 286.783354815 |
| 340.212618143 |
point_group_type: 3m
| Phonon mode (cm-1) | Visualize Phonons hereRepresentation |
|---|---|
| -17.6032476313 | A1 I+R |
| -0.3041137814 | None |
| -0.3029098759 | None |
| 33.9856389929 | E I+R |
| 86.3302590565 | E I+R |
| 93.9403733061 | A1 I+R |
| 104.236085802 | A1 I+R |
| 119.503604525 | A1 I+R |
| 134.939045348 | A1 I+R |
| 185.044017829 | E I+R |
| 186.017109891 | E I+R |
| 208.990021537 | E I+R |
| 216.092907748 | E I+R |
| 286.783354815 | A1 I+R |
| 340.212618143 | A1 I+R |
The orbital magnetic moment was obtained after SCF run. Please note no DFT+U parameters were taken into account.
-0.001 μB
ICSD-ID: None