Abstract

We report the study of the mechanical structure and optical properties of lead halide perovskite CH3NH3PbI3 using ab initio methods. The ground state energy calculations were performed within density functional theory and generalized gradient approximation using the pseudopotential method with plane-wave basis sets. The norm conserving pseudopotential was used. The ground state properties of the electronic structure of the perovskite were used and elastic parameters such as bulk modulus B, Young’s modulus E, shear modulus G, and Poisson’s ratio were determined and found to be in good agreement with experimental values. The ratio obtained was found to be greater than 1.75. Poisson’s ratio () was obtained as 0.25 implying that CH3NH3PbI3 is a ductile material. The absorption coefficient within the energy range of 0 to 6 eV was found to be 5.76 × 105 cm−1 indicating maximum absorption. The absorption coefficient compares well with the available experimental and computed values.

1. Introduction

The desire to achieve high efficiency solar cells has led to research in perovskites that are lead-based, in particular lead halide perovskites. Since lead is not environmentally friendly, there has been a need to have organic lead halide perovskites with less lead, but organic lead-free perovskites are less stable [1]. Saliba et al. [2] reported that lead halide perovskites have become popular as photovoltaic materials due to their high power conversion efficiency (PCE) which stands at over 22%. The optical properties of these new materials are important not only to device design but also because of the insight they provide into less accessible properties such as energy band structures and binding energies, among other properties.

A study by Feng (2014), investigated the mechanical properties of hybrid organic-inorganic CH3NH3BX3 (X = Br, I; B = Sn, Pb) perovskites used in photovoltaic cell absorbers. The structure consists of a network of corner-sharing BX6 octahedra, where the B atom is a metal cation (Sn2+ or Pb2+) and X is the hybrid anion (Br-, Cl or I). The dependent octahedra give room for a broad readjustment of the B-X-B bond angle, for this case, Pb-I-Pb. The various sets of joint rotations are called tilt transitions [3]. This develops symmetry, which shows different structures at varied temperatures. For instance, the three structural phases in which the compound CH3NH3PbI3 exists include pseudocubic at increasing temperatures (above 327 K), tetragonal at room temperature, and orthorhombic at reducing temperatures (below 162.2 K) [3].

According to Zhu et al. 2021, the tetragonal phase changes to a cubic phase at 0.3 GPa. In addition to the 2H-phase, Young’s modulus (E) reduces as pressure increases. This implies that the stiffness reduces with increasing pressure. Nonetheless, in orthorhombic, tetragonal, and cubic phases, as the pressure increases, E first increases and then reduces. In particular, the cubic phase has only three independent elastic constants: C11, C12, and C44 which were determined in this work [3, 4].

The crystal structures, elastic, and anisotropy properties of the CH3NH3PbI3 were studied using the ab initio calculations.

The rationale of employing first-principles (ab initio) calculations was premised on the assertion that the mechanical properties of the compounds are relatively difficult to measure through an experimental approach. It was revealed that the absorption performance of perovskite solar cells greatly depends on the crystalline and stress state of the perovskite layer. The study by Bretschneider et al. [5] concurred with previously conducted work, and it indicated that the mechanical properties of perovskite are crucial for practical applications.

The study conducted by Sun et al. [6] analyzed factors that influence the mechanical properties of formamidinium lead halides and associated perovskites. The study formulated a systematic way of probing the mechanical properties of hybrid perovskite single crystals under nanoindentation. It was revealed that the shape, size, and hydrogen bonding resulting from the organic cations have a significant influence on their mechanical properties. Similarly, it was found that bonding in the inorganic framework and hydrogen bonding play a vital role in determining elastic stiffness.

The perovskite structure with the formula ABX3 is used in many oxide compounds. In this study, A represents an organic positively charged ion (CH3NH3+), X represents n halide specifically Iodide (I), and B is a divalent metal ion, in this case Pb2+. This structure is a cubic unit cell which contains an A atom (CH3NH3+) in the center of the cube, B atoms Pb2+ at the corners, and X atoms (I) at the center of the cell edges. The modelled pseudocubic structure is shown in Figure 1.

The calculated elastic parameters in this work include C11, C12, C44, bulk modulus B, Young’s modulus E, shear modulus G, Poisson’s ratio , and anisotropy which describe the mechanical stability and ductility of CH3NH3PbI3 [7, 8]. The study also noted that the measured Young’s moduli (9.7–12.3 GPa) and hardness (0.36–0.45 GPa) reflected good mechanical flexibility and ductility. The study indicated that the mechanical properties of lead halide and related perovskites are important in device fabrication and performance.

The objective of this work is to establish the ideal mechanical and optical properties for photovoltaic applications from the first principle method. The rest of this paper is organized such that Section 2 is computational details, Section 3 is the results and discussions, and Section 4 deals with conclusions.

2. Computational Details

In this work, scalar relativistic electronic-structure calculations on this material were carried out based on density-functional theory (DFT) [9], plane waves, and the pseudopotential approach as implemented in the Quantum Espresso computer code [10]. In the calculation of total energy, the exchange-correlation potential is treated with the generalized gradient approximation of Perdew–Burke–Ernzerhof (GGA). The lattice constant of CH3NH3PbI3 was well optimized [11, 12]. The total energy convergence in the iterative solution of the Kohn–Sham equations was fixed at 2 × 10−8 Rydberg (Ry), and self-consistency was achieved [13, 14]. All calculations were carried out under ground state conditions. The energy cut off was obtained as 30 Ry and 4 × 4 × 4 k-point grid or more was found to be sufficient for this material.

Core electrons were replaced by ab initio norm conserving pseudopotentials, generated using the Troullier–Martins scheme [15]. In the Kleinman–Bylander fully nonlocal separable representation [16]. According to this arrangement, s1, 2s22p, 2s22p3, 5s25p5, and 5d106s26p2 were used as valence electrons for H, C, N, I, and Pb, respectively. The large overlap between the semicore and valence states makes it possible for the semicore 5d electrons of Pb to be treated as valence electrons and distinctly included in the simulations, as explained in [17]. In this study, the electronic density, Hartree, and exchange correlation potentials were calculated in a uniform real-space grid with an equivalent plane-wave cutoff of 55 Ry in the representation of charge density. To integrate the Brillouin zone, we employed a Monkhorst–Pack sampling [18] equivalent to 13 × 13 × 13 in a twelve-atom CH3NH3PbI3 unit cell [19]. The ground state energies were determined by running the files in Appendices A and B. The determination of e-cut energy, alat, k-points, and the pseudocubic structure are the preliminary calculations determined for the optical and mechanical properties of CH3NH3PbI3.

The approximations used in the calculation include the following: the exchange correlation GGA and the Voigt–Reuss–Hill approximation in calculating bulk, shear, and Young’s modulus. Other types of approximations are beyond the objective of this study.

3. Results and Discussions

3.1. Mechanical Properties of CH3NH3PbI3

Mechanical properties of solids include elasticity, strength, abrasion, hardness, ductility, brittleness, malleability, and toughness. The elastic constants give the relation between dynamic properties and mechanical properties with respect to forces existing in materials.

The elastic constants of CH3NH3PbI3 determined in this work are given in Table 1.

Table 1 shows the mechanical properties calculation of CH3NH3PbI3 using the Quantum Espresso code, and results were compared with other available data. The average of the Voigt–Reuss–Hill approximations was used in calculating bulk (B), shear (G), and Young’s modulus (E).

The results show that our calculated elastic constants for CH3NH3PbI3 compare well with other work. These results confirm that pseudocubic CH3NH3PbI3 is stable mechanically because it fulfills the Born stability criteria of cubic crystals whereby C44 > 0, C11 > C12, and C11 + 2C 12 > 0 [20]. The ratio obtained was found to be greater than 1.75 implying that CH3NH3PbI3 is a ductile material. These results confirm the ductility property of the perovskite [14, 2123].

The computed anisotropic ratio which was obtained as 0.6 indicated that cubic CH3NH3PbI3 was elastically anisotropic because of substituting it in the equation . The value of A was obtained to be smaller than one (A <1).

Poisson’s ratio obtained was 0.25 which when rounded off to one decimal place is closely equal to 0.3 [24] which confirms that CH3NH3PbI3 is a ductile material.

3.2. Optical Properties of CH3NH3PbI3

Optical properties of a material include reflectance, conductivity, absorption, refraction, transmittance, dispersion, diffraction, and so on. Solar panels, for example, need to absorb sufficient light energy. Therefore, a lower solar reflectance index rating is required of the material for high absorption. The solar transmittance of a surface is the fraction of the sun’s radiation that is transmitted through the surface of the solar material. What is utilized by solar materials to produce electricity is the amount of solar radiation absorbed. This therefore prompted the study of absorption other than reflectance and transmittance as the main optical property in this work with respect to its suitability as a photovoltaic material. The absorption coefficient of CH3NH3PbI3 has therefore been computed in this work.

The absorption coefficient shows how deep light rays will penetrate into a layer before it is absorbed. Light is rarely absorbed in a material with a low absorption coefficient. For a direct gap semiconductor like CH3NH3PbI3, the absorption coefficient can reach elevated values.

A material that has a high optical absorption is suitable to be used in photovoltaic applications [25].

The high photovoltaic performance of CH3NH3PbI3 is associated with optically high absorption characteristics. It is clear from Figure 2 that the material has a high absorption ability. Its absorption coefficient is about 1.4 × 106 cm−1 for energy between 0 and 40 eV which contributes to the effective utilization of solar radiation.

This means that, with a high absorption coefficient, CH3NH3PbI3 readily absorbs photons, which excite electrons from the valence band into the conduction band. In the solar spectrum, visible light is found within the range of wavelength λ ≈ 380–780 nm which translates to energy of approximately 1.6–3.5 eV. Figure 3 indicates that at about 3.5 eV there is a peak in the absorption of the visible spectrum by CH3NH3PbI3. The peak in Figure 3 gives the absorption coefficient within the given range of energy and it was found to be 5.76 × 105 cm−1 indicating maximum absorption. This value can be compared with the derived absorption coefficient of 1.0978 × 105 cm−1 [26] and 8.1448 × 104 cm−1 [27] which greatly depends on the range at which the energy or wavelength was measured from. The data plotted in Figure 2 are in Appendix C.

The absorption coefficient determines how far into a material light of a particular wavelength or energy can pass through a material before it is absorbed. Absorption of solar energy will occur within the near UV region, visible region, and near infrared region between energy 0.886 eV and 3.98 eV. This is equivalent to a wavelength of 1400 nm–380 nm. This means that the material has a relatively wide absorption range, as shown by the wavelength. The high absorption coefficient of CH3NH3PbI3 indicates that the material readily absorbs photons with energy equivalent to the band gap energy which excites electrons from the valence band into the conduction band.

The imaginary part of Figure 4 indicates a maximum absorption of photon energy between 0 and 5.0 eV, which is within the region where solar energy is utilized in a solar cell. CH3NH3PbI3 has a band gap of 1.58 eV [28]. This means that only photons of energy equal to 1.58 eV will be absorbed by the material. Photons of higher or lower energy value than the band gap are wasted.

Figure 4 shows that the appearance of a sharp peak of the imaginary part of the dielectric constant of the CH3NH3PbI3 implies the occurrence of strong absorption in this spectral region. The real part is related to the refractive index, while the imaginary part gives the absorption coefficient. The data plotted in the figure are in Appendix D.

Generally, a material with a high dielectric constant has a relatively less charge carrier recombination rate. As a result, the overall performance of optoelectronic devices is enhanced.

4. Conclusion

The result from the elastic constants done compares well with previously carried out work. The elastic parameters determined include bulk modulus B, Young’s modulus E, shear modulus G, and Poisson’s ratio . The ratio and Poisson’s ratio obtained implied that CH3NH3PbI3 is a ductile material. The evaluated elastic parameters enable us to conclude that the investigated compound is mechanically stable. This means that it can be molded into different shapes and sizes; therefore, suitable for the fabrication and modelling of solar cells.

Here, the optical constants which include absorption coefficients and real and imaginary components of the dielectric constant are reported. The value of the absorption coefficient within the spectral range of 0–6 eV within which the visible spectrum exists is 5.76 × 105 cm−1. There was, however, inadequate literature from experimental and computational for the absorption coefficient to compare these results with, but from the results obtained, it can be concluded that CH3NH3PbI3 has a high absorption coefficient of solar radiation and is highly suitable to be used in the fabrication and modelling of solar cells.

Appendix

A. Input file for pwscf code: Ecut.in, k-points.in and alat.in

&CONTROL   title = “cubic_CH3NH3PbI3”,  calculation = “scf”,restart_mode = “from_scratch”,  pseudo_dir = “. /”,  outdir = “./tmp”,   prefix = “pseudo”,   tstress = .true.,  tprnfor = .true.,/&SYSTEM   ibrav = 0,   celldm(1) = 12.0772387957,   nat = 12,   ntyp = 5,   nbnd = 50,  ecutwfc = 40.0,  ecutrho = 500.0,occupations = “smearing”,  smearing = “marzari-vanderbilt”,  degauss = 0.0075,  !London_s6 = 0.750,! London_rcut = 50.0,! noncolin = .true.,! lda_plus_u = .false.,! lspinorb = .true.,! nosym = .true.,/&ELECTRONS   conv_thr = 1.0D-7,/ATOMIC_SPECIESPb 207.200Pb.pbe-mt_fhi.UPFI 126.900 I.pbe-mt_fhi.UPFH 1.008 H.pbe-mt_fhi.UPFC 12.011C.pbe-mt_fhi.UPFN 14.007N.pbe-mt_fhi.UPFATOMIC_POSITIONS {crystal}C 0.3939882648252961 0.4998165502286085 0.4578307086810938N 0.6090744378334350 0.5000072812331666 0.5477685929486071H 0.4101503344735065 0.4996421926480608 0.2857772518934283H 0.3086756133568471 0.6447043733807547 0.5110306384926346H 0.3088856672784672 0.3549521469503247 0.5113014134206253H 0.6962691547534092 0.6357485658144455 0.5015129819664779H 0.6964106005070789 0.3642097749426654 0.5017817996551273H 0.6063401459460280 0.5001870434389417 0.7112184893396645Pb 0.9448835519577230 0.9999795897298682 0.9776228940296079I 0.9094873486598019 0.9999820990246278 0.4751454839169895I 0.8992708791277906 0.4999580210910821 0.0236644354684188I 0.4323428768805755 0.9999948727175436 0.9224794036872126K_POINTS {automatic}13 13 13 0 0 0

B. Input file for pwscf code: Pseudo-Cubic structure

&CONTROLtitle = “cubic_CH3NH3PbI3”,calculation = “scf”,restart_mode = “from_scratch”,  pseudo_dir = “./”,   outdir = “./tmp”,   prefix = “pseudo”,   tstress = .true.,   tprnfor = .true.,/&SYSTEM   ibrav = 0,celldm(1) = 12.0772387957,   nat = 12,   ntyp = 5,   nbnd = 50,  ecutwfc = 50.0,  ecutrho = 500.0,occupations = “smearing”,  smearing = “marzari-vanderbilt”,   degauss = 0.0075,  !London_s6 = 0.750,”! London_rcut = 50.0,! noncolin = .true.,! lda_plus_u = .false.,! lspinorb = .true.,! nosym = .true.,/&ELECTRONS   conv_thr = 1.0D-7,/ATOMIC_SPECIESPb 207.200Pb.pbe-mt_fhi.UPFI 126.900 I.pbe-mt_fhi.UPFH 1.008 H.pbe-mt_fhi.UPFC 12.011C.pbe-mt_fhi.UPFN 14.007N.pbe-mt_fhi.UPFATOMIC_POSITIONS {crystal}C 0.3939882648252961 0.4998165502286085 0.4578307086810938N 0.6090744378334350 0.5000072812331666 0.5477685929486071H 0.4101503344735065 0.4996421926480608 0.2857772518934283H 0.3086756133568471 0.6447043733807547 0.5110306384926346H 0.3088856672784672 0.3549521469503247 0.5113014134206253H 0.6962691547534092 0.6357485658144455 0.5015129819664779H 0.6964106005070789 0.3642097749426654 0.5017817996551273H 0.6063401459460280 0.5001870434389417 0.7112184893396645Pb 0.9448835519577230 0.9999795897298682 0.9776228940296079I 0.9094873486598019 0.9999820990246278 0.4751454839169895I 0.8992708791277906 0.4999580210910821 0.0236644354684188I 0.4323428768805755 0.9999948727175436 0.9224794036872126K_POINTS {automatic}13 13 13 0 0 0CELL_PARAMETERS {alat}1.0 0.0 0.00.0 1.0 0.0

C. Absorption coefficient. Out

0.0000000000000000 0.00000000000000001.9999978624805612E-002 0.249889288909469853.9999957249611223E-002 0.564337626545359355.9999935874416835E-002 0.938804090359954387.9999914499222446E-002 1.36373633479785489.9999893124028058E-002 1.82492224111180710.11999987174883367 2.30430545492078930.13999985037363930 2.78142195995021790.15999982899844489 3.23505685598979520.17999980762325052 3.64519557310504670.19999978624805612 3.99475838529902740.21999976487286174 4.27114639154359650.23999974349766734 4.46720206740008230.25999972212247296 4.58162543722798610.27999970074727859 4.61873717761052040.29999967937208416 4.58767035858269430.31999965799688979 4.50118512869115680.33999963662169541 4.37420030315634280.35999961524650104 4.22240780187399170.37999959387130661 4.06102220725911070.39999957249611223 3.90404521363536270.41999955112091786 3.76406325673717920.43999952974572348 3.65294555195807020.45999950837052905 3.58360061046790350.47999948699533468 3.57338977994208080.49999946562014030 3.64992739578270340.51999944424494593 3.86066123489111620.53999942286975156 4.28825869751666210.55999940149455718 5.07474714441122380.57999938011936270 6.45852563704788100.59999935874416832 8.82910544744260580.61999933736897395 12.8059962075689260.63999931599377957 19.3466079050963810.65999929461858520 29.8909705410863750.67999927324339082 46.5389395883708940.69999925186819645 72.2722559632764360.71999923049300207 111.173991614910700.73999920911780759 168.706578060926150.75999918774261321 251.827325741331440.77999916636741884 369.297210311687080.79999914499222446 531.283292141323050.81999912361703009 749.891842902998750.83999910224183572 1037.29558294846700.85999908086664134 1407.54411400569850.87999905949144697 1871.1235629043078

D. Epsilon_imaginary data

0.0000000000000000 5.3344163605233434E-0041.9999978624805612E-002 6.1341647492641534E-0043.9999957249611223E-002 6.9270033738413005E-0045.9999935874416835E-002 7.6831501735604530E-0047.9999914499222446E-002 8.3720472518939446E-0049.9999893124028058E-002 8.9646173515729776E-0040.11999987174883367 9.4355558235892508E-0040.13999985037363930 9.7654034997798542E-0040.15999982899844489 9.9421597189845528E-0040.17999980762325052 9.9622441246949539E-0040.19999978624805612 9.8306940473013942E-0040.21999976487286174 9.5605795049326185E-0040.23999974349766734 9.1717153941446345E-0040.25999972212247296 8.6888354097770174E-0040.27999970074727859 8.1394522609164774E-0040.29999967937208416 7.5516561144114915E-0040.31999965799688979 6.9520960676005248E-0040.33999963662169541 6.3643521379392102E-0040.35999961524650104 5.8078475824888121E-0040.37999959387130661 5.2973874225868385E-0040.39999957249611223 4.8433556982319375E-0040.41999955112091786 4.4525793070104391E-0040.43999952974572348 4.1298883291785365E-0040.45999950837052905 3.8804883906635231E-0040.47999948699533468 3.7134244973854934E-0040.49999946562014030 3.6466686038921385E-0040.51999944424494593 3.7147087673180520E-0040.53999942286975156 3.9799477086350474E-0040.55999940149455718 4.5497040612939641E-0040.57999938011936270 5.6010896872353295E-0040.59999935874416832 7.4164052921932842E-0040.61999933736897395 1.0431792142133784E-0030.63999931599377957 1.5301479011784957E-0030.65999929461858520 2.2978804873596687E-0030.67999927324339082 3.4813009216057358E-0030.69999925186819645 5.2657362272757581E-0030.71999923049300207 7.8979530811855176E-0030.73999920911780759 1.1695941614998524E-0020.75999918774261321 1.7055375221479332E-0020.77999916636741884 2.4450168760568370E-0020.79999914499222446 3.4424301272587242E-0020.81999912361703009 4.7572249380705797E-0020.83999910224183572 6.4506167837405079E-0020.85999908086664134 8.5809463239980802E-0020.87999905949144697 0.111978621251519390.89999903811625248 0.143357877767025300.91999901674105811 0.180074183183597330.93999899536586373 0.221982331913983030.95999897399066936 0.268631451897507980.97999895261547498 0.319263636779867300.99999893124028061 0.372852917690719541.0199989098650861 0.428187932758826221.0399988884898919 0.483994958983819561.0599988671146974 0.539090327117887851.0799988457395031 0.592543965683534621.0999988243643086 0.643830439536894671.1199988029891144 0.692941788748749121.1399987816139199 0.740438732483997391.1599987602387254 0.787423660094039861.1799987388635311 0.835429682613328951.1999987174883366 0.886233377059471921.2199986961131424 0.941612553637258041.2399986747379479 1.00308191845986721.2599986533627536 1.07164653648656131.2799986319875591 1.14761378018719841.2999986106123647 1.23049824696742991.3199985892371704 1.31904144914879271.3399985678619759 1.41135067839770861.3599985464867816 1.50514209377986121.3799985251115872 1.59805512413298741.3999985037363929 1.68799203592198601.4199984823611984 1.77343065948103921.4399984609860041 1.85366122089472411.4599984396108097 1.92890982256062671.4799984182356152 2.00032950415458231.4999983968604209 2.06986174193191141.5199983754852264 2.13999263910008251.5399983541100322 2.21344485483218101.5599983327348377 2.29285531280718141.5799983113596434 2.38048830605522581.5999982899844489 2.47802412958529181.6199982686092544 2.58644709278146491.6399982472340602 2.70603735028948971.6599982258588657 2.83645263541913951.6799982044836714 2.97687241285786051.6999981831084769 3.12617057990789251.7199981617332827 3.28308417332088441.7399981403580882 3.44635323496611031.7599981189828939 3.61481829432630471.7799980976076994 3.78747349741217711.7999980762325050 3.96348227676384651.8199980548573107 4.14216680714916931.8399980334821162 4.32298209834173091.8599980121069219 4.50548171719883591.8799979907317275 4.68927706086616601.8999979693565332 4.87398824872023001.9199979479813387 5.05918383857951601.9399979266061442 5.24430924283457231.9599979052309500 5.42860906182762901.9799978838557555 5.61105460373952881.9999978624805612 5.79029222309985952.0199978411053667 5.96462875629020942.0399978197301722 6.13206630366185752.0599977983549782 6.29039037872223532.0799977769797837 6.43730483928595732.0999977556045892 6.57059670314235332.1199977342293947 6.68830667929047442.1399977128542003 6.7888790441831928Epsilon_real data0.0000000000000000 6.19314307148631741.9999978624805612E-002 6.19304070536089933.9999957249611223E-002 6.19385631654395665.9999935874416835E-002 6.19526026507185367.9999914499222446E-002 6.19741887831759369.9999893124028058E-002 6.20016351240672490.11999987174883367 6.20362769784909140.13999985037363930 6.20770479554940740.15999982899844489 6.21249367007555890.17999980762325052 6.21793329140753230.19999978624805612 6.22408731355195140.21999976487286174 6.23093527817928820.23999974349766734 6.23850865827811150.25999972212247296 6.24682057933217210.27999970074727859 6.25587805041387310.29999967937208416 6.26571828315889960.31999965799688979 6.27633471276581380.33999963662169541 6.28777879774458800.35999961524650104 6.30004183372950700.37999959387130661 6.31318173972132610.39999957249611223 6.32719742033901870.41999955112091786 6.34214889422096380.43999952974572348 6.35805001578009770.45999950837052905 6.37496182133869030.47999948699533468 6.39291942788583970.49999946562014030 6.41198513632271890.51999944424494593 6.43222467998139000.53999942286975156 6.45369817575557820.55999940149455718 6.47652288600928470.57999938011936270 6.50073741611663620.59999935874416832 6.52656049180325940.61999933736897395 6.55394587184157280.63999931599377957 6.58332751462393610.65999929461858520 6.61440918025057380.67999927324339082 6.64808123609759160.69999925186819645 6.68342916806222130.71999923049300207 6.72226923596733260.73999920911780759 6.76235278168916490.75999918774261321 6.80725016833468730.77999916636741884 6.85216559272653390.79999914499222446 6.90372659990873140.81999912361703009 6.95282174715095640.83999910224183572 7.01092251127993920.85999908086664134 7.06245306335163470.87999905949144697 7.12577638323430660.89999903811625248 7.17684586190970730.91999901674105811 7.24268680333884340.93999899536586373 7.28974166233883650.95999897399066936 7.35441378984596380.97999895261547498 7.39437237620118370.99999893124028061 7.45435594730161631.0199989098650861 7.48604256188613261.0399988884898919 7.53957050272943091.0599988671146974 7.56470618490611901.0799988457395031 7.61299670792553501.0999988243643086 7.63591976263466511.1199988029891144 7.68312485190378691.1399987816139199 7.70892897190866841.1599987602387254 7.76036873759983961.1799987388635311 7.79211647008061981.1999987174883366 7.85139039944265621.2199986961131424 7.88790426587086291.2399986747379479 7.95441214718522721.2599986533627536 7.99017754679228001.2799986319875591 8.05870874603499131.2999986106123647 8.08638560408285481.3199985892371704 8.14943525802319791.3399985678619759 8.16386487412300891.3599985464867816 8.21567748218483371.3799985251115872 8.21711957484640371.3999985037363929 8.25719614192675341.4199984823611984 8.25172765544369471.4399984609860041 8.28561088466764421.4599984396108097 8.28221473262303581.4799984182356152 8.31890452500757501.4999983968604209 8.32481739935015241.5199983754852264 8.37223818001637101.5399983541100322 8.38922900174486231.5599983327348377 8.45036578691415491.5799983113596434 8.47397472860590281.5999982899844489 8.54576164815933661.6199982686092544 8.56758910347434011.6399982472340602 8.64277441654414781.6599982258588657 8.65408385043583911.6799982044836714 8.72449796805247501.6999981831084769 8.71895808219584011.7199981617332827 8.77813496524354871.7399981403580882 8.75258242124549211.7599981189828939 8.79664229723685761.7799980976076994 8.75023841112260391.7999980762325050 8.77731218720476751.8199980548573107 8.71020009269462391.8399980334821162 8.71933660488774541.8599980121069219 8.63160058201298241.8799979907317275 8.62182152884850341.8999979693565332 8.51290524072811121.9199979479813387 8.48254245182987181.9399979266061442 8.35106767408848331.9599979052309500 8.29742103692090541.9799978838557555 8.14157363686575411.9999978624805612 8.06119250326999292.0199978411053667 7.87989232338982592.0399978197301722 7.76977905211694252.0599977983549782 7.56434876300656712.0799977769797837 7.42379368481350182.0999977556045892 7.19919577590463082.1199977342293947 7.03125406311260102.1399977128542003 6.7959420381638100

Data Availability

The density functional theory (DFT) codes were used in this work to run the input files to obtain the output data used in this work. The code used is available at https://www.quantum-espresso.org/. Pseudopotentials and other resources on the study material are available on the website. The other data supporting this research work are from previously reported studies and datasets, which have been cited. Some of the input files used to process data have been given in the appendices of the document.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors gratefully acknowledge the computer resources, technical expertise, and assistance provided by the Centre for High-Performance Computing (CHPC), Cape Town, South Africa.