JVASP-11775_Na3MoNO3
JARVIS-ID:JVASP-11775 Functional:optB88-vdW Primitive cell Primitive cell Conventional cell Conventional cell
Chemical formula:Na3MoNO3 Formation energy/atom (eV):-1.743 a 5.595 Å α:90.0 ° a 5.595 Å α:90.0 °
Space-group :Pmn2_1, 31 Relaxed energy/atom (eV):-4.4474 b 6.233 Å β:90.0 ° b 6.233 Å β:90.0 °
Calculation type:Bulk SCF bandgap (eV):3.81 c 7.306 Å γ:90.0 ° c 7.306 Å γ:90.0 °
Crystal system:orthorhombic Point group:mm2 Density (gcm-3):2.96 Volume (3):254.76 nAtoms_prim:16 nAtoms_conv:16
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): 3.778I


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

Static real-parts of dielectric function in x,y,z: 4.22,3.41,3.39


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

Static real-parts of dielectric function in x,y,z: 5.1,2.66,2.67


Solar-cell SLME [Reference]

Theoretical solar-cell efficiency (in %) was calculated using spectroscopy limited maximum efficiency (SLME) and TBmBJ for the material with 500 nm thickness and at 300 K. Note that generally there are many factors that contribute towards the efficiency, such as carrier effective mass etc.

SLME is: 0.0


DFPT: IR-intensity, Piezoelecric and Dielectric tensors [Reference]

Calculations are done using density functional perturbation theory (DFPT) method for non-metallic systems for conventional cell and at Gamma-point in phonon BZ.

Static dielecric-tensor

5.73 -0.0 0.0
-0.0 4.86 -0.0
0.0 -0.0 5.02

Piezoelectric-stress-tensor (C/m2)

-0.87 0.15 0.2 0.0 0.0 0.0
0.0 0.0 -0.0 0.21 0.0 0.0
-0.0 -0.0 0.0 0.0 0.0 0.17

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

Reuss-bulk modulus (KR): 43.18 GPa, Reuss-shear modulus (GR): 15.61 GPa

Poisson's ratio: 0.34, Elastic anisotropy parameter: 0.32

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

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

Elastic tensor
48.9 32.6 35.3 0.0 -0.0 -0.0
32.6 73.3 43.2 0.0 0.0 -0.0
35.3 43.2 74.1 0.0 0.0 -0.0
0.0 0.0 0.0 15.8 -0.0 -0.0
-0.0 0.0 -0.0 -0.0 15.7 -0.0
-0.0 -0.0 -0.0 0.0 -0.0 22.0

Phonon mode (cm-1)
-0.05
0.01
0.05
81.37
99.56
104.23
104.36
125.0
140.46
153.51
163.95
177.68
180.67
180.86
192.69
199.94
202.25
204.66
216.57
219.38
227.36
233.15
233.8
251.09
251.19
256.93
267.39
267.53
272.37
279.29
331.55
345.62
347.61
349.19
355.92
366.38
402.93
408.54
422.52
424.47
697.64
706.79
729.9
733.82
745.68
760.86
941.97
950.31

Point group

point_group_type: mm2

Visualize Phonons here
Phonon mode (cm-1) Representation
-0.05
-0.0488009985
0.01
0.0122967493
0.05
0.0486774723
81.37
81.3743804952
99.56
99.5605688222
104.23
104.226777515
104.36
104.364423178
125.0
125.000872661
140.46
140.461358354
153.51
153.507827982
163.95
163.950176051
177.68
177.679709393
180.67
180.671538992
180.86
180.855807152
192.69
192.691259221
199.94
199.94325224
202.25
202.248479401
204.66
204.661979902
216.57
216.574346893
219.38
219.377516877
227.36
227.355257934
233.15
233.152515883
233.8
233.804080337
251.09
251.09063459
251.19
251.194807442
256.93
256.929474147
267.39
267.389709369
267.53
267.529704503
272.37
272.373190126
279.29
279.285451625
331.55
331.546443061
345.62
345.618683933
347.61
347.607864378
349.19
349.191787752
355.92
355.917348639
366.38
366.384960721
402.93
402.929727481
408.54
408.538722077
422.52
422.520661695
424.47
424.472933047
697.64
697.63504688
706.79
706.791455389
729.9
729.902753009
733.82
733.823370627
745.68
745.681375804
760.86
760.855705298
941.97
941.96974417
950.31
950.313174259

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)

4.2 0.0 0.0
0.0 3.51 0.0
0.0 0.0 4.93

Hole mass tensor (me unit)

8.19 -0.0 0.0
-0.0 17.05 0.0
0.0 0.0 14.37

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 -425.57 -408.88 -394.59
n-PowerFactor 480.73 623.25 816.98
n-Conductivity 2875.41 4002.96 4510.93
n-ZT 0.28 0.37 0.48
p-Seebeck 425.28 442.69 517.0
p-PowerFactor 154.67 163.62 924.2
p-Conductivity 834.92 855.19 3457.63
p-ZT 0.09 0.1 0.55

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
Na-0.0-0.00.0-0.0
Na-0.0-0.00.0-0.0
Na-0.0-0.00.0-0.0
Na-0.0-0.00.0-0.0
Na-0.0-0.00.0-0.0
Na-0.0-0.00.0-0.0
Mo0.00.00.00.0
Mo0.00.00.00.0
N-0.0-0.00.0-0.0
N-0.0-0.00.0-0.0
O-0.0-0.00.0-0.0
O-0.0-0.00.0-0.0
O-0.0-0.00.0-0.0
O-0.0-0.00.0-0.0
O-0.0-0.00.0-0.0
O-0.0-0.00.0-0.0

See also

Links to other databases or papers are provided below


mp-19105

ICSD-ID: 98673

AFLOW link

MP link
mp-19105

NIST Disclaimer