Temperature dependence of electronic conductivity from ab initio thermal simulation
Pith reviewed 2026-05-16 07:41 UTC · model grok-4.3
The pith
Thermally averaging the square of the density of states near the Fermi level estimates how electronic conductivity changes with temperature.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
By thermally averaging the square of the density of states near the Fermi level from finite-temperature ab initio molecular dynamics, the temperature dependence of electronic conductivity is estimated through a direct extension of the Hindley-Mott formula. When applied to crystalline aluminum and grain-boundary aluminum the method reproduces the Bloch-Gruneisen decrease; for semiconducting amorphous silicon and GST it captures the thermally activated increase; and for the aluminum-graphene composite it registers the microstructure-induced conduction channel.
What carries the argument
the thermally-averaged Hindley-Mott (TAHM) method, which replaces the static density of states in the original Hindley-Mott formula with a time average of its square taken near the Fermi level from ab initio MD trajectories
If this is right
- Metallic systems exhibit the Bloch-Gruneisen conductivity decrease with rising temperature.
- Semiconducting amorphous materials show conductivity increase from thermal activation of carriers.
- Microstructure effects, such as grain boundaries or layered interfaces, produce distinct activation signatures in the conductivity trend.
- The method supplies a computationally light route to conductivity trends in disordered or composite materials.
- Temperature-dependent electronic structure fluctuations are directly mapped onto transport without separate scattering-rate calculations.
Where Pith is reading between the lines
- The same averaging procedure could be applied to other response functions that depend on the density of states, such as the Seebeck coefficient.
- In highly disordered systems the method may offer a way to separate intrinsic thermal activation from defect-induced hopping without additional models.
- Because the input is only the time series of the density of states, the approach could be inserted into existing MD workflows for rapid screening of high-temperature conductors.
- Extension to magnetic or spin-polarized cases would require only replacing the total density of states with the appropriate spin channel.
Load-bearing premise
The original Hindley-Mott approximate formula remains sufficiently accurate when its static density-of-states input is replaced by a time average taken from finite-temperature molecular dynamics.
What would settle it
Direct comparison of TAHM-predicted conductivity versus temperature curves against measured data for any one of the five systems (for example, the magnitude and sign of the slope for a-Si) would confirm or refute the estimate.
Figures
read the original abstract
We present a temperature-dependent extension of the approximate electronic conductivity formula of Hindley and Mott that leverages time-averaged fluctuations of the electronic density of states obtained from ab initio molecular dynamics. By thermally averaging the square of the density of states near the Fermi level, we obtain an estimate of the temperature dependence of the conductivity. This approach termed the thermally-averaged Hindley-Mott (TAHM) method was applied to five representative systems: crystalline aluminum (c-Al), aluminum with a grain boundary (AlGB), a four-layer graphene-aluminum composite (Al-Gr), amorphous silicon (a-Si) and amorphous germanium-antimony-telluride (a-GST). The method reproduces the expected Bloch-Gruneisen decrease in conductivity for c-Al and AlGB. Generally, the reduction (increase) in conductivity for metallic (semiconducting) materials are reproduced. It captures microstructure-induced, thermally activated conduction in multilayer Al-Gr, a-Si and a-GST. Overall, the approach provides a computationally efficient link between time-dependent electronic structure and temperature-dependent transport, offering a simple and approximate tool for exploring electronic conductivity trends in complex and disordered materials.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes the thermally-averaged Hindley-Mott (TAHM) method, which extends the approximate Hindley-Mott conductivity formula by replacing its static density-of-states input with a time average of [N(E_F, t)]^2 extracted from ab initio molecular dynamics trajectories. The approach is applied to five systems (crystalline Al, Al grain boundary, Al-graphene multilayer, amorphous Si, and amorphous GST) and is reported to recover the expected Bloch-Grüneisen decrease in conductivity for the metallic cases and the increase for the semiconducting cases, while also capturing microstructure-induced conduction trends.
Significance. If the central approximation holds, the method supplies a parameter-free, computationally lightweight route to temperature-dependent conductivity trends directly from standard AIMD runs, without explicit scattering-rate calculations. This could be useful for screening complex or disordered materials where full transport theory remains expensive.
major comments (2)
- [methods and c-Al results] The application to crystalline Al (reported to recover the Bloch-Grüneisen T-dependence) rests on the assumption that thermal DOS fluctuations alone reproduce the leading electron-phonon scattering effect. The original Hindley-Mott relation was derived for static structural disorder; the manuscript does not demonstrate that the time-averaged N(E_F)^2 is mathematically equivalent to the proper thermal average over the electron-phonon scattering operator (see the derivation in the methods section and the c-Al results paragraph).
- [results and discussion] No quantitative comparison to independent transport calculations (e.g., Boltzmann transport or explicit electron-phonon matrix elements) or to experimental conductivity values with error bars is provided for any system. Without such benchmarks, it is unclear whether the observed trends are physically grounded or arise from incidental broadening of the DOS by MD-induced lattice distortions.
minor comments (2)
- [abstract and computational details] The abstract and introduction should explicitly state the MD timestep, total simulation length, and k-point sampling used to compute the time-averaged DOS for each material.
- [methods] Notation for the thermally averaged quantity (e.g., whether it is <N(E_F)^2> or <N(E_F)>^2) should be defined once in the methods section and used consistently in all figures and equations.
Simulated Author's Rebuttal
We thank the referee for their detailed and constructive review. We address each major comment below, clarifying the scope of the TAHM approximation and adding discussion where needed. The revisions focus on emphasizing the method's approximate nature for trend capture rather than claiming equivalence to full transport theory.
read point-by-point responses
-
Referee: [methods and c-Al results] The application to crystalline Al (reported to recover the Bloch-Grüneisen T-dependence) rests on the assumption that thermal DOS fluctuations alone reproduce the leading electron-phonon scattering effect. The original Hindley-Mott relation was derived for static structural disorder; the manuscript does not demonstrate that the time-averaged N(E_F)^2 is mathematically equivalent to the proper thermal average over the electron-phonon scattering operator (see the derivation in the methods section and the c-Al results paragraph).
Authors: We agree that TAHM is not derived as a direct mathematical equivalent to the thermal average of the electron-phonon scattering operator. The original Hindley-Mott formula applies to static disorder, and our extension uses time-averaged DOS fluctuations from AIMD as a practical proxy to capture temperature-induced changes in conductivity trends. This is presented as an approximation suitable for complex systems where full e-ph calculations are expensive. We have revised the methods section to explicitly state the approximate character of the approach, note that it is not intended to replace scattering-rate methods, and added a brief discussion of its relation to Bloch-Grüneisen behavior as an empirical observation rather than a rigorous derivation. revision: partial
-
Referee: [results and discussion] No quantitative comparison to independent transport calculations (e.g., Boltzmann transport or explicit electron-phonon matrix elements) or to experimental conductivity values with error bars is provided for any system. Without such benchmarks, it is unclear whether the observed trends are physically grounded or arise from incidental broadening of the DOS by MD-induced lattice distortions.
Authors: We acknowledge the value of quantitative benchmarks. The manuscript's focus is on demonstrating that TAHM recovers expected qualitative trends (Bloch-Grüneisen decrease for metals, increase for semiconductors, microstructure effects) across five diverse systems using only standard AIMD data. Full Boltzmann transport or explicit e-ph matrix element calculations for the disordered cases (a-Si, a-GST, AlGB) remain computationally demanding and were outside the scope of this initial method paper. We have added a new paragraph in the discussion comparing the c-Al trend to literature experimental conductivity temperature dependence and note the absence of absolute-value benchmarks as a limitation of the current study. revision: partial
Circularity Check
No significant circularity; TAHM estimate uses independent MD input
full rationale
The derivation replaces the static DOS input of the external Hindley-Mott formula with a time average of N(E_F,t)^2 computed from separate ab initio MD trajectories. This step does not reduce the conductivity estimate to a fitted parameter, a self-definition, or a self-citation chain; the MD data constitute independent dynamical input. No load-bearing uniqueness theorem or ansatz is imported from the authors' prior work, and the reported trends for c-Al, a-Si, etc., are presented as numerical outcomes rather than algebraic identities. The procedure therefore remains self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Ab initio molecular dynamics trajectories provide statistically representative samples of electronic density-of-states fluctuations at finite temperature.
Reference graph
Works this paper leans on
-
[1]
J. M. Ziman, Formal transport theory, in: Electrons and Phonons: The Theory of Transport Phenomena in Solids, Oxford University Press, 2001, Ch. 7.doi:10.1093/acprof:oso/9780198507796. 003.0007
-
[2]
Giustino, Electron-phonon interactions from first principles, Rev
F. Giustino, Electron-phonon interactions from first principles, Rev. Mod. Phys. 89 (2017) 015003.doi:10.1103/RevModPhys.89. 015003
-
[3]
G. S. Rocha, M. N. Ferreira, G. S. Denicol, J. Noronha, Transport co- efficients of quasiparticle models within a new relaxation time approx- imation of the boltzmann equation, Phys. Rev. D 106 (2022) 036022. doi:10.1103/PhysRevD.106.036022
-
[4]
R. Claes, S. Poncé, G.-M. Rignanese, G. Hautier, Phonon-limited elec- tronic transport through first principles, Nature Reviews Physics 7 (2) (2025) 73–90.doi:10.1038/s42254-024-00795-0
-
[5]
G. K. Madsen, D. J. Singh, Boltztrap. a code for calculating band- structure dependent quantities, Computer Physics Communications 175 (1) (2006) 67–71.doi:10.1016/j.cpc.2006.03.007
-
[6]
G. K. Madsen, J. Carrete, M. J. Verstraete, Boltztrap2, a program for in- terpolating band structures and calculating semi-classical transport coef- ficients, Computer Physics Communications 231 (2018) 140–145.doi: 10.1016/j.cpc.2018.05.010
-
[7]
S. Poncé, E. Margine, C. Verdi, F. Giustino, Epw: Electron–phonon cou- pling, transport and superconducting properties using maximally local- ized wannier functions, Computer Physics Communications 209 (2016) 116–133.doi:10.1016/j.cpc.2016.07.028
-
[8]
C. Ringhofer, Numerical methods for the semiconductor boltzmann equa- tion based on spherical harmonics expansions and entropy discretiza- tions, Transport Theory and Statistical Physics 31 (4-6) (2002) 431–452. doi:10.1081/TT-120015508
-
[9]
R. Kubo, M. Yokota, S. Nakajima, Statistical-mechanical theory of ir- reversible processes. ii. response to thermal disturbance, Journal of the Physical Society of Japan 12 (11) (1957) 1203–1211
work page 1957
-
[10]
D. A. Greenwood, The Boltzmann equation in the theory of electrical con- duction in metals, Poceedings of the Physical Society of London 71 (460) (1958) 585–596.doi:10.1088/0370-1328/71/4/306
-
[11]
S. K. Bose, O. Jepsen, O. K. Andersen, Real-space calculation of the electrical resistivity of liquid 3d transition metals using tight-binding linear muffin-tin orbitals, Phys. Rev. B 48 (1993) 4265–4275.doi: 10.1103/PhysRevB.48.4265. URLhttps://link.aps.org/doi/10.1103/PhysRevB.48. 4265
-
[12]
L. Calderín, V . Karasiev, S. Trickey, Kubo–Greenwood electrical con- ductivity formulation and implementation for projector augmented wave datasets, Computer Physics Communications 221 (2017) 118–142.doi: 10.1016/j.cpc.2017.08.008
-
[13]
J. Dufty, J. Wrighton, K. Luo, S. Trickey, On the kubo-greenwood model for electron conductivity, Contributions to Plasma Physics 58 (2-3) (2018) 150–154.doi:10.1002/ctpp.201700102
-
[14]
K. Prasai, K. N. Subedi, K. Ferris, P. Biswas, D. A. Drabold, Spatial projection of electronic conductivity: The example of conducting bridge memory materials, physica status solidi (RRL) – Rapid Research Letters 12 (9) (2018) 1800238.doi:10.1002/pssr.201800238
-
[15]
K. N. Subedi, K. Prasai, D. A. Drabold, Space-projected conductivity and spectral properties of the conduction matrix, physica status solidi (b) 258 (2021) 2000438.doi:10.1002/pssb.202000438
-
[16]
N. Mott, Conduction in glasses containing transition metal ions, Journal of Non-Crystalline Solids 1 (1) (1968) 1–17. doi:https://doi.org/10.1016/0022-3093(68)90002-1. URLhttps://www.sciencedirect.com/science/ article/pii/0022309368900021
-
[17]
N. F. Mott, E. A. Davis, Theory of electrons in a non-crystalline medium, in: Electronic Processes in Non-crystalline Materials, 2nd Edi- tion, Clarendon/Oxford University Press, Oxford, New York, 1979, Ch. 2, pp. 11–15
work page 1979
-
[18]
N. K. Hindley, Random phase model of amorphous semiconduc- tors i. transport and optical properties, Journal of Non-Crystalline Solids 5 (1) (1970) 17–30.doi:https://doi.org/10.1016/ 0022-3093(70)90193-6
work page 1970
-
[19]
K. Nepal, C. Ugwumadu, F. Kraft, Y . Al-Majali, D. Drabold, The ef- fects of crystal orientation and common coal impurities on electronic conductivity in copper–carbon composites, Carbon 231 (2025) 119711. doi:10.1016/j.carbon.2024.119711
-
[20]
C. Ugwumadu, D. A. Drabold, R. M. Tutchton, Effects of galactic irradi- ation on thermal and electronic transport in tungsten, physica status solidi (b) 262 (11) (2025) 2500109.doi:10.1002/pssb.202500109
-
[21]
T. A. Abtew, M. Zhang, D. A. Drabold, Ab initio estimate of tempera- ture dependence of electrical conductivity in a model amorphous mate- rial: Hydrogenated amorphous silicon, Phys. Rev. B 76 (2007) 045212. doi:10.1103/PhysRevB.76.045212
-
[22]
K. N. Subedi, K. Kappagantula, F. Kraft, A. Nittala, D. A. Drabold, Electrical conduction processes in aluminum: Defects and phonons, Phys. Rev. B 105 (2022) 104114.doi:10.1103/PhysRevB.105. 104114
-
[23]
A. Chowdhury, S. Bhattacharjee, Experimental investigation of change in sheet resistance and debye temperatures in metallic thin films due to low-energy ion beam irradiation, Journal of Physics D: Applied Physics 46 (43) (2013) 435304
work page 2013
-
[24]
R. Thapa, C. Ugwumadu, K. Nepal, J. Trembly, D. A. Drabold, Ab initio simulation of amorphous graphite, Phys. Rev. Lett. 128 (2022) 236402. doi:10.1103/PhysRevLett.128.236402
-
[25]
C. Ugwumadu, R. Thapa, K. Nepal, D. A. Drabold, Atomistic nature of amorphous graphite, European Journal of Glass Science and Technology Part B 64 (2023) 16–22.doi:10.13036/17533562.64.1.18
-
[26]
K. Nepal, C. Ugwumadu, A. Gautam, K. Kappagantula, D. A. Drabold, Electronic conductivity in metal-graphene composites: the role of disor- dered carbon structures, defects, and impurities, Journal of Physics: Ma- terials 7 (2) (2024) 025003.doi:10.1088/2515-7639/ad261a
-
[27]
A. Lewis, Evidence for the mott model of hopping conduction in the an- neal stable state of amorphous silicon, Phys. Rev. Lett. 29 (1972) 1555– 1558.doi:10.1103/PhysRevLett.29.1555
-
[28]
S. S. Ashwin, U. V . Waghmare, S. Sastry, Metal-to-semimetal transition in supercooled liquid silicon, Phys. Rev. Lett. 92 (2004) 175701.doi: 10.1103/PhysRevLett.92.175701
- [29]
-
[30]
K. Nepal, A. Gautam, R. Hussein, K. Konstantinou, S. Elliott, C. Ugwu- madu, D. Drabold, Electronic and thermal properties of the phase-change memory material, ge2sb2te5, and results from spatially resolved transport calculations, Solid State Sciences 173 (2026) 108182.doi:https:// doi.org/10.1016/j.solidstatesciences.2025.108182
-
[31]
G. Kresse, J. Furthmüller, Efficient iterative schemes for ab initio total- energy calculations using a plane-wave basis set, Phys. Rev. B 54 (1996) 11169–11186.doi:10.1103/PhysRevB.54.11169
-
[32]
A. Nittala, J. Smith, B. Gwalani, J. Silverstein, F. F. Kraft, K. Kappagan- 7 tula, Simultaneously improved electrical and mechanical performance of hot-extruded bulk scale aluminum-graphene wires, Materials Science and Engineering: B 293 (2023) 116452.doi:10.1016/j.mseb.2023. 116452
- [33]
-
[34]
B. R. Djordjevi ´c, M. F. Thorpe, F. Wooten, Computer model of tetra- hedral amorphous diamond, Phys. Rev. B 52 (1995) 5685–5689.doi: 10.1103/PhysRevB.52.5685
-
[35]
F. Wooten, K. Winer, D. Weaire, Computer generation of structural mod- els of amorphous si and ge, Phys. Rev. Lett. 54 (1985) 1392–1395. doi:10.1103/PhysRevLett.54.1392
-
[36]
D. Igram, B. Bhattarai, P. Biswas, D. Drabold, Large and realistic models of amorphous silicon, Journal of Non-Crystalline Solids 492 (2018) 27– 32.doi:10.1016/j.jnoncrysol.2018.04.011
-
[37]
C. Ugwumadu, A. Gautam, Y . G. Lee, D. A. Drabold, Mapping ther- mal conductivity at the atomic scale: A step toward the thermal design of materials, physica status solidi (b) 2500316doi:10.1002/pssb. 202500316
-
[38]
P. E. Blöchl, Projector augmented-wave method, Phys. Rev. B 50 (1994) 17953–17979.doi:10.1103/PhysRevB.50.17953
-
[39]
J. P. Perdew, K. Burke, M. Ernzerhof, Generalized gradient approxi- mation made simple, Phys. Rev. Lett. 77 (1996) 3865–3868.doi: 10.1103/PhysRevLett.77.3865
-
[40]
S. Nosé, A unified formulation of the constant temperature molecular dy- namics methods, The Journal of Chemical Physics 81 (1) (1984) 511–519. doi:10.1063/1.447334
-
[41]
B. Smyrak, M. Gniełczyk, Study of the effect of graphene content on the electrical and mechanical properties of aluminium–graphene composites, Materials 18 (2025) 590.doi:10.3390/ma18030590
-
[42]
T. Kato, K. Tanaka, Electronic properties of amorphous and crystalline ge2sb2te5 films, Japanese Journal of Applied Physics 44 (10R) (2005) 7340.doi:10.1143/JJAP.44.7340
-
[43]
N. Ashcroft, N. Mermin, Solid State Physics, Holt, Rinehart and Winston, 1976
work page 1976
-
[44]
K. S. Kappagantula, Reza-E-Rabby, A. M. Nittala, H. Das, J. Christu- dasjustus, T. Poudel, S. Wolff, Bulk-scale ultra-conductors via low-c manufacturing pathways, Tech. rep., DOE Vehicles Technology Office Annual Merit Review (2025). URLhttps://www1.eere.energy.gov/ vehiclesandfuels/downloads/2025_AMR/MAT294_ Kappagantula_2025_o.pdf
work page 2025
-
[45]
E. P. Donovan, F. Spaepen, D. Turnbull, J. M. Poate, D. C. Jacobson, Heat of crystallization and melting point of amorphous silicon, Applied Physics Letters 42 (8) (1983) 698–700.doi:10.1063/1.94077
-
[46]
D. A. Drabold, P. A. Fedders, S. Klemm, O. F. Sankey, Finite-temperature properties of amorphous silicon, Phys. Rev. Lett. 67 (1991) 2179–2182. doi:10.1103/PhysRevLett.67.2179
-
[47]
J. Li, D. A. Drabold, Electron hopping between localized states: A simu- lation of the finite-temperature anderson problem using density functional methods, Phys. Rev. B 68 (2003) 033103.doi:10.1103/PhysRevB. 68.033103
-
[48]
H. Sasaki, A. Ikari, K. Terashima, S. K. Shigeyuki Kimura, Temperature dependence of the electrical resistivity of molten silicon, Japanese Journal of Applied Physics 34 (7R) (1995) 3426.doi:10.1143/JJAP.34. 3426
-
[49]
J. Li, D. Drabold, Atomistic simulation of the finite-temperature an- derson localization problem, physica status solidi (b) 233 (1) (2002) 10–17.doi:10.1002/1521-3951(200209)233:1<10:: AID-PSSB10>3.0.CO;2-V
-
[50]
B.-S. Lee, J. R. Abelson, S. G. Bishop, D.-H. Kang, B.-k. Cheong, K.-B. Kim, Investigation of the optical and electronic properties of ge2sb2te5 phase change material in its amorphous, cubic, and hexagonal phases, Journal of Applied Physics 97 (9) (2005) 093509.doi:10.1063/1. 1884248
work page doi:10.1063/1 2005
-
[51]
S. Muneer, J. Scoggin, F. Dirisaglik, L. Adnane, A. Cywar, G. Bakan, K. Cil, C. Lam, H. Silva, A. Gokirmak, Activation energy of metastable amorphous ge2sb2te5 from room temperature to melt, AIP Advances 8 (6) (2018) 065314.doi:10.1063/1.5035085. 8 SUPPLEMENTARY MATERIAL Temperature dependence of electronic conductivity fromab initiothermal simulation R. ...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.