pith. sign in

arxiv: 1907.00556 · v2 · pith:QY3AK4GEnew · submitted 2019-07-01 · ❄️ cond-mat.mtrl-sci

Highly accurate prediction of material optical properties based on density functional theory

Mitsutoshi Nishiwaki , Hiroyuki Fujiwara This is my paper

Pith reviewed 2026-05-25 12:21 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords density functional theoryoptical absorption spectrasolar cell materialsGGA calculationshybrid functionalssum ruleabsorption coefficient
0
0 comments X

The pith

A uniform energy shift and amplitude adjustment applied to standard DFT spectra reproduces experimental absorption coefficients for solar cell materials.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper presents a method that starts from a generalized gradient approximation calculation performed on a very dense k-point mesh and then applies two corrections to bring the result into close agreement with measured optical absorption spectra. One correction shifts the energy scale using the result of a hybrid functional calculation, while the other adjusts the overall amplitude using the f-sum rule. The same two corrections are used without further tuning for five different solar-cell compounds. A reader would care because this produces spectra that match experiment even on a logarithmic scale, at far lower computational cost than GW or full hybrid-functional approaches.

Core claim

The PHS method obtains highly accurate absorption-coefficient spectra by taking GGA results computed on a high-density k mesh, applying an energy-scale correction derived from a hybrid functional, and applying an amplitude correction derived from the sum rule, yielding spectra that agree closely with experiment for GaAs, InP, CdTe, CuInSe2 and Cu2ZnGeSe4.

What carries the argument

The PHS correction procedure that combines a high-density GGA spectrum with a hybrid-functional energy shift and a sum-rule amplitude adjustment.

If this is right

  • The corrected spectra are more accurate than those from conventional GGA, hybrid functionals, or GW methods for the materials tested.
  • The computational cost remains lower than that of the more advanced methods while still reaching experimental agreement.
  • The same uniform corrections can be applied to other solar-cell or semiconductor materials without additional experimental input.
  • Optical functions beyond absorption coefficient become accessible once the corrected spectra are available.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the corrections remain uniform, the method could be applied to predict optical response in alloys or doped variants of the same base materials.
  • The approach might be tested on optical properties at finite temperature by combining it with existing phonon or molecular-dynamics calculations.
  • Success on absorption spectra suggests the possibility of similar post-processing corrections for other response functions such as dielectric constants or reflectivity.

Load-bearing premise

The same energy shift taken from a hybrid functional and the same amplitude adjustment taken from the sum rule can be applied uniformly to every material without material-specific re-fitting.

What would settle it

Compute the corrected spectrum for an additional compound using only the uniform corrections and compare it directly to a new experimental absorption measurement; large systematic deviation would falsify the claim.

Figures

Figures reproduced from arXiv: 1907.00556 by Hiroyuki Fujiwara, Mitsutoshi Nishiwaki.

Figure 2
Figure 2. Figure 2: (a) ε2 spectra and (b) ε1 spectra of GaAs obtained from experiment (open circles) and theoretical DFT calculations (solid lines). For the calculations, the results determined by PBE, HSE06, and the developed PHS method are shown. The experimental data were taken from Ref. [17]. and the optical transition energies observed in the experimental spectrum are reproduced well, as reported previously [9]. When th… view at source ↗
Figure 3
Figure 3. Figure 3: (a) GaAs α spectra calculated by the PHS method using different k-mesh densities and (b) comparison of GaAs α spectra calculated by the PHS, HSE06 and GW methods. The experimental result (open circles) is taken from Ref. [17] and the GW spectrum in (b) is adopted from Ref. [8]. The ∆Eg represents the energy-shift value of the PBE spectra (∆Eg = Eg,HSE – Eg,PBE) [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: α spectra of (a) CdTe, (b) InP, (c) CISe and (d) CZGSe, obtained from experiment (open circles) and the calculations using the PHS method (solid lines). The ∆Eg values of each spectrum are also indicated. The experimental data were adopted from Refs. [14], [17] and [29]. universality of our PHS approach. To justify our method further, we have calculated the band structures of the solar cell materials [PIT… view at source ↗
Figure 5
Figure 5. Figure 5: Band structures of (a) CdTe, (b) InP, (c) C [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
read the original abstract

Theoretical material investigation based on density functional theory (DFT) has been a breakthrough in the last century. Nevertheless, the optical properties calculated by DFT generally show poor agreement with experimental results particularly when the absorption-coefficient ({\alpha}) spectra in logarithmic scale are compared. In this study, we have established an alternative DFT approach (PHS method) that calculates highly accurate {\alpha} spectra, which show remarkable agreement with experimental spectra even in logarithmic scale. In the developed method, the optical function estimated from generalized gradient approximation (GGA) using very high-density k mesh is blue-shifted by incorporating the energy-scale correction by a hybrid functional and the amplitude correction by sum rule. Our simple approach enables high-precision prediction of the experimental {\alpha} spectra of all solar-cell materials (GaAs, InP, CdTe, CuInSe2 and Cu2ZnGeSe4) investigated here. The developed method is superior to conventional GGA, hybrid functional and GW methods and has clear advantages in accuracy and computational cost.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 1 minor

Summary. The manuscript introduces the PHS method, which computes optical absorption coefficient (α) spectra via GGA on a high-density k-mesh, followed by a uniform blue-shift derived from the hybrid-functional bandgap difference and an amplitude rescaling via the sum rule. The central claim is that this yields highly accurate agreement with experimental log-scale α spectra for GaAs, InP, CdTe, CuInSe₂ and Cu₂ZnGeSe₄, outperforming standard GGA, hybrid functionals and GW calculations while remaining computationally cheaper.

Significance. If the rigid-shift construction proves robust and parameter-free, the approach would offer a practical, low-cost route to accurate optical spectra for photovoltaic materials, with clear computational advantages over GW. The emphasis on log-scale fidelity for solar-cell compounds is a useful target.

major comments (3)
  1. [Abstract] Abstract and PHS-method description: the claim of 'remarkable agreement' and 'high-precision prediction' of experimental α spectra is unsupported by any quantitative error metric (e.g., RMS deviation, MAE or R² on log α). Without such numbers the superiority statements cannot be evaluated.
  2. [PHS-method description] PHS-method description (energy-scale correction paragraph): the hybrid bandgap difference is applied as a single rigid shift to the entire GGA spectrum. This scissor-operator assumption is load-bearing for the log-scale agreement claim, yet the manuscript supplies no test that the GGA–hybrid discrepancy is constant across all transitions rather than k- or band-character dependent.
  3. [Methods] Methods/Results sections: the procedure for extracting and applying the hybrid energy correction is not described in sufficient detail to rule out post-hoc adjustment or material-specific choices, raising the possibility that the final spectra are not fully independent predictions.
minor comments (1)
  1. [Abstract] The abstract contains raw LaTeX fragments (e.g., {α}); these should be rendered in the published version.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive comments on our manuscript. We address each major point below and have revised the manuscript to incorporate quantitative metrics, expanded methodological details, and additional discussion of the underlying approximations.

read point-by-point responses

  1. Referee: [Abstract] Abstract and PHS-method description: the claim of 'remarkable agreement' and 'high-precision prediction' of experimental α spectra is unsupported by any quantitative error metric (e.g., RMS deviation, MAE or R² on log α). Without such numbers the superiority statements cannot be evaluated.

    Authors: We agree that quantitative metrics are necessary to rigorously support the claims. In the revised manuscript we have added Table II reporting MAE, RMSD and R² values computed on log₁₀(α) for the PHS spectra versus experiment (and versus GGA, hybrid and GW results) over the 0–3 eV range for all five materials. The tabulated RMSD values are typically <0.25, confirming the visual fidelity shown in the figures. revision: yes

  2. Referee: [PHS-method description] PHS-method description (energy-scale correction paragraph): the hybrid bandgap difference is applied as a single rigid shift to the entire GGA spectrum. This scissor-operator assumption is load-bearing for the log-scale agreement claim, yet the manuscript supplies no test that the GGA–hybrid discrepancy is constant across all transitions rather than k- or band-character dependent.

    Authors: The PHS method adopts the conventional scissor-operator approximation based on the difference in fundamental band gaps. While a k-resolved decomposition of the hybrid correction would be desirable, performing such calculations for all materials lies outside the present scope. We have added a paragraph in the revised Discussion section that explicitly states this assumption, notes its common use in the literature, and acknowledges that its validity is ultimately supported by the empirical agreement across the studied compounds rather than by a direct k-dependent validation. revision: partial

  3. Referee: [Methods] Methods/Results sections: the procedure for extracting and applying the hybrid energy correction is not described in sufficient detail to rule out post-hoc adjustment or material-specific choices, raising the possibility that the final spectra are not fully independent predictions.

    Authors: The procedure is already specified in the Methods section: the HSE06 band gap is computed once per material on a Γ-centered 8×8×8 mesh, the difference from the GGA gap is taken as a single scalar, and this scalar is applied uniformly to the dense-k GGA spectrum before sum-rule rescaling. No additional fitting parameters are introduced. To eliminate any ambiguity we have expanded the Methods section with a numbered step-by-step protocol and a short pseudocode block in the revised manuscript. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected in derivation chain

full rationale

The PHS method computes GGA optical spectra on a dense k-mesh, then applies an energy-scale blue-shift taken from the hybrid-functional bandgap and an amplitude rescaling from the f-sum rule. These steps are independent DFT calculations plus a physical constraint; the resulting spectra are validated against external experimental data for multiple materials rather than being forced to match by construction or by self-referential fitting. No quoted equations reduce the output to the input, no self-citation chain bears the central claim, and no ansatz is smuggled via prior work. The approach therefore remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Based solely on the abstract, the central claim rests on the transferability of a hybrid-functional energy shift and a sum-rule amplitude correction to GGA results; no explicit free parameters, new axioms, or invented entities are named, but the corrections themselves may embed fitted choices whose independence cannot be verified from the abstract.

axioms (1)
  • domain assumption Standard GGA and hybrid-functional DFT approximations remain valid starting points for the optical-function calculation.
    The method begins from GGA and hybrid results without questioning their foundational validity.

pith-pipeline@v0.9.0 · 5710 in / 1396 out tokens · 54939 ms · 2026-05-25T12:21:57.988889+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

38 extracted references · 38 canonical work pages

  1. [1]

    The DFT methods for optical-function calculation have already been established [5-9] and a vast amount of optical spectra deduced from DFT calculations have been reported [1 -15]

    Introduction Prediction of material optical properties based o n density functional theory (DFT) has been a revolutionary technique that allows quite ef fective optical-material searches even without forming materials experimentally [1-4]. The DFT methods for optical-function calculation have already been established [5-9] and a vast amount of optical spe...

    work page

  2. [2]

    with HSE06 and sum rule and, from this PHS app roach (PBE+HSE06+Sum rule), 3 the dielectric function ( ε = ε1 – i ε2) and optical constants (refractive index n, extinction coefficient k and α) are readily obtained in a consistent manner. As a result, we find that the α spectra of representative solar cells materials [G aAs, CdTe, InP, CuInSe 2 (CISe), and...

    work page

  3. [3]

    Here, as an example, the calculation of a Ga As dielectric function is shown

    PHS method Figure 1 explains the calculation procedure of th e PHS method developed in this study. Here, as an example, the calculation of a Ga As dielectric function is shown. In our approach, the ε2 spectrum of the dielectric function is calculated first using very high k-mesh density within GGA-PBE. As known well [23], Eg is seriously underestimated wh...

    work page

  4. [4]

    For the calculatio ns of GGA within PBE, the Advance/PHASE software was employed, while the V ASP software was applied for HSE06 calculations

    DFT calculation The DFT calculations were performed using Advance /PHASE and the Vienna Ab initio Simulation Package (V ASP) [24]. For the calculatio ns of GGA within PBE, the Advance/PHASE software was employed, while the V ASP software was applied for HSE06 calculations. For the DFT calculations of zin cblende crystals (GaAs, CdTe, InP), two-atom primit...

    work page

  5. [5]

    For the DFT results, those obtained applying PBE, HSE06 and the PHS meth od are shown

    Results and Discussion Figure 2 shows (a) the ε2 spectra and (b) the ε1 spectra of GaAs obtained from experiment (open circles) and the DFT calculations (solid lines). For the DFT results, those obtained applying PBE, HSE06 and the PHS meth od are shown. The experimental spectrum was taken from Ref. [17]. Since Eg is seriously underestimated in PBE, the w...

    work page

  6. [6]

    The α spectrum estimated from GW is adopted from Ref

    and the PHS method and (b) the GaAs α spectra calculated by the PHS, HSE06 and GW methods. The α spectrum estimated from GW is adopted from Ref. [8 ]. In Fig. 3(a), 7 the α spectra calculated by the PHS method using differe nt k-mesh densities are shown. It can be seen that the DFT calculation using a ver y high mesh density is vital for the accurate calc...

    work page

  7. [7]

    Conclusion We have developed a new DFT approach that can acc urately predict material α spectra in logarithmic scale. In this method, the ε2 spectrum calculated from the PBE functional using very high- k mesh density is blue-shifted and the underestimate d Eg 11 contribution in PBE is corrected using the energy s cale determined by HSE06 calculation, whil...

    work page

  8. [8]

    Hinuma, T

    Y . Hinuma, T. Hatakeyama, Y . Kumagai, L.A. Burton, H. Sato, Y . Muraba, S. Iimura, H. Hiramatsu, I. Tanaka, H. Hosono, F. Oba, Nat. Commun. 7 (2016) 11962

    work page 2016

  9. [9]

    Shaposhnikov, A.V

    V .L. Shaposhnikov, A.V . Krivosheeva, V .E. Boris enko, J.-L. Lazzari, F. Arnaud d’Avitaya, Phys. Rev. B 85 (2012) 205201

    work page 2012

  10. [10]

    W. Meng, B. Saparov, F. Hong, J. Wang, D.B. Mitzi, Y . Yan, Chem. Mater. 28 (2016) 821-829

    work page 2016

  11. [11]

    M. Kato, T. Fujiseki, T. Miyadera, T. Sugita, S . Fujimoto, M. Tamakoshi, M. Chikamatsu, H. Fujiwara, J. Appl. Phys. 121 (2017) 115501

    work page 2017

  12. [12]

    Gajdoš, K

    M. Gajdoš, K. Hummer, G . Kresse, J. Furthmüller , F. Bechstedt, Phys. Rev. B 73 (2006) 045112

    work page 2006

  13. [13]

    Ambrosch-Draxl, J.O

    C. Ambrosch-Draxl, J.O. Sofo, Comput. Phys. 175 (2006) 1-14

    work page 2006

  14. [14]

    Kootstra, P.L

    F. Kootstra, P.L. de Boeij, J.G . Snijders, Phys. Rev. B 62 (2000) 7071-7083

    work page 2000

  15. [15]

    Botti, F

    S. Botti, F. Sottile, N. Vast, V . Olevano, L. R eining, H.-C. Weissker, A. Rubio, G . Onida, R.D. Sole, R.W. Godby, Phys. Rev. B 69 (2004) 155112

    work page 2004

  16. [16]

    Paier, M

    J. Paier, M. Marsman, G . Kresse, Phys. Rev. B 78 (2008) 121201

    work page 2008

  17. [17]

    Fujiwara, M

    H. Fujiwara, M. Kato, M. Tamakoshi, T. Miyader a, M. Chikamatsu, Phys. Status Solidi A 215 (2018) 1700730

    work page 2018

  18. [18]

    Paier, R

    J. Paier, R. Asahi, A. Nagoya G . Kresse, Phys. Rev. B 79 (2009) 115126

    work page 2009

  19. [19]

    Persson, J

    C. Persson, J. Appl. Phys. 107 (2010) 053710

    work page 2010

  20. [20]

    Nishiwaki, K

    M. Nishiwaki, K. Nagaya, M. Kato, S. Fujimoto, H. Tampo, T. Miyadera, M. Chikamatsu, H. Shibata, H. Fujiwara, Phys. Rev. Materials 2 (2018) 085404

    work page 2018

  21. [21]

    Nagaya, S

    K. Nagaya, S. Fujimoto, H. Tampo, S. Kim, M. N ishiwaki, Y . Nishigaki, M. Kato, H. Shibata, H. Fujiwara, Appl. Phys. Lett. 113 (2018) 093901

    work page 2018

  22. [22]

    E.M. Chen, L. Williams, A. Olvera, C. Zhang, M . Zhang, G . Shi, J.T. Heron, L. Qi, L.J. Guo, E. Kioupakis, P.F.P. Poudeu, Chem. Sci. 9 (2018) 5405-5414

    work page 2018

  23. [23]

    Nakane, H

    A. Nakane, H. Tampo, M. Tamakoshi, S. Fujimoto , K.M. Kim, S. Kim, H. Shibata, S. Niki, H. Fujiwara, J. Appl. Phys. 120 (2016) 064505

    work page 2016

  24. [24]

    Fujiwara, R.W

    H. Fujiwara, R.W. Collins, editors, Spectroscopic Ellipsometry for Photovoltaics: Volume 2: Applications and Optical Data of Solar Ce ll Materials , Springer, Cham, 2018

    work page 2018

  25. [25]

    Perdew, K

    J.P. Perdew, K. Burke, M. Ernzerhof, Phys. Rev. Lett. 77 (1996) 3865-3868

    work page 1996

  26. [26]

    J. Wang, H. Chen, S.-H. Wei, W.-J. Yin, Adv. Mater. 31 (2019) 1806593

    work page 2019

  27. [27]

    Heyd, G .E

    J. Heyd, G .E. Scuseria, M. Ernzerhof, J. Chem. Phys. 118 (2003) 8207-8215

    work page 2003

  28. [28]

    Heyd G .E

    J. Heyd G .E. Scuseria, M. Ernzerhof, J. Chem. Phys. 124 (2006) 219906. 13

    work page 2006

  29. [29]

    Marder, Condensed Matter Physics, Wiley, Hoboken, 2010

    M.P. Marder, Condensed Matter Physics, Wiley, Hoboken, 2010

    work page 2010

  30. [30]

    Marsman, J

    M. Marsman, J. Paier, A. Stroppa, G . Kresse, J. Phys.: Condens. Matter 20 (2008), 064201

    work page 2008

  31. [31]

    Kresse, J

    G . Kresse, J. Furthmüller, Phys. Rev. B 54 (1996) 11169-11186

    work page 1996

  32. [32]

    Kageshima, K

    H. Kageshima, K. Shiraishi, Phys. Rev. B 56 (1997) 14985-14992

    work page 1997

  33. [33]

    Zhang, X

    Y . Zhang, X. Yuan, X. Sun, B.-C. Shih, P. Zhan g, W. Zhang, Phys. Rev. B 84 (2011) 075127

    work page 2011

  34. [34]

    M. Kato, M. Nishiwaki, H. Fujiwara, 2019, https://arxiv.org/abs/1906.03005

    work page arXiv 2019

  35. [35]

    Shirayama, H

    M. Shirayama, H. Kadowaki, T. Miyadera, T. Sug ita, M. Tamakoshi, M. Kato, T. Fujiseki, D. Murata, S. Hara, T.N. Murakami, S. Fujimoto, M. Chikamatsu, H. Fujiwara, Phys. Rev. Applied 5 (2016) 014012

    work page 2016

  36. [36]

    Minoura, K

    S. Minoura, K. Kodera, T. Maekawa, K. Miyazaki , S. Niki, H. Fujiwara, J. Appl. Phys. 113 (2013) 063505

    work page 2013

  37. [37]

    F. Gygi, A. Baldereschi, Phys. Rev. Lett. 62 (1989) 2160-2163

    work page 1989

  38. [38]

    Bechstedt, Many-Body Approach to Electronic Excitations: Conce pts and Applications , Springer, Heidelberg, 2015

    F. Bechstedt, Many-Body Approach to Electronic Excitations: Conce pts and Applications , Springer, Heidelberg, 2015. 14 Supplementary Information Highly accurate prediction of material optical prop erties based on density functional theory Mitsutoshi Nishiwaki and Hiroyuki Fujiwara Department of Electrical, Electronic and Computer E ngineering, Gifu Unive...

    work page 2015