Temperature- and charge carrier density-dependent electronic response in methylammonium lead iodide

Hai I. Wang; Heejae Kim; Jaco J. Geuchies; Jiacheng Wang Jungmin Park; Lei Gao; Li-Lin Wu; Lucia Di Virgilio; Mischa Bonn; Sheng Qu; Wen Zeng

REVIEW 3 major objections 6 minor 6 references

Temperature- and charge carrier density-dependent electronic response in methylammonium lead iodide

T0 review · 3 major / 6 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read This paper maps the temperature-dependent electronic response of methylammonium lead iodide across carrier densities from $10^{14}$ to $10^{19}\,\mathrm{cm}^{-3}$, identifying three distinct regimes.

desk verdict A well-executed but calibration-dependent extension of the authors' earlier Mott-polaron work; the temperature-dependent Mott-density map is useful, but the density axis needs an independent high-density check before the phase diagram is trusted quantitatively. read the letter →

arxiv 2505.18887 v1 pith:U7JQ4MBQ submitted 2025-05-24 cond-mat.mtrl-sci physics.chem-ph

Jiacheng Wang Jungmin Park , Lei Gao , Lucia Di Virgilio , Sheng Qu , Heejae Kim , Hai I. Wang , Li-Lin Wu , Wen Zeng

show 3 more authors

Mischa Bonn Zefeng Ren Jaco J. Geuchies

This is my paper · ORCID

classification cond-mat.mtrl-sciphysics.chem-ph

keywords metal-halideperovskitemethylammoniumleadiodideMottdensitycarriertrappingpolaronannihilationoptical-pump/THzprobespectroscopytransientabsorptionelectronicphasediagram

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper aims to establish that in methylammonium lead iodide (MAPI), the fate of photoexcited carriers is set by where the excitation density falls on a temperature-dependent map, not by a single intrinsic recombination time. By combining optical-pump/THz-probe spectroscopy with highly sensitive transient absorption, the authors cover five orders of magnitude in carrier density and temperatures from 78 K to 315 K. They identify a low-density regime (below roughly $10^{15}\,\mathrm{cm}^{-3}$) where shallow traps remove carriers within picoseconds, an intermediate regime (up to about $10^{18}\,\mathrm{cm}^{-3}$) where direct bimolecular recombination dominates and signals persist for hundreds of picoseconds, and a high-density regime above the Mott density where overlapping polaron wavefunctions annihilate via an Auger-type process within tens of picoseconds. The result is a quantified electronic phase diagram that places the Mott density between $1.34\times10^{18}\,\mathrm{cm}^{-3}$ at 78 K and $4.5\times10^{18}\,\mathrm{cm}^{-3}$ at 315 K, and it shows that these Mott densities sit at or above the densities needed for population inversion. If correct, this gives device designers a direct way to anticipate whether a given operating condition will be trap-limited, recombination-limited, or annihilation-limited.

What carries the argument

The load-bearing object is the electronic phase diagram, assembled from two complementary spectroscopies. Optical-pump/THz-probe (OPTP) spectroscopy supplies the high-density branch, in which the THz photoconductivity (proportional to carrier density times mobility) is calibrated against fluence; the critical Mott density is obtained by extrapolating the late-time plateau of the transient back to zero delay and reading where it crosses the linear low-fluence calibration. Transient absorption spectroscopy (TAS) with roughly $10^{-7}$ sensitivity supplies the low-density branch, in which the ratio of the signal at 10 ps to the instantaneous signal tracks the fraction of carriers lost to shallow traps. The two datasets share a carrier-density axis computed from fluence, absorption, reflection, and a photon-to-carrier quantum yield determined independently from Drude-Smith fits to THz conductivity spectra. The Mott-density values also feed a geometric estimate of polaron radii through the filling condition $N_{\mathrm{Mott}}\,\frac{4}{3}\pi r^3 = 0.74$.

What would settle it

Measure the THz conductivity spectrum (not just the peak transient) as a function of pump fluence above $10^{18}\,\mathrm{cm}^{-3}$ at a fixed temperature: if the extracted mobility decreases with density or the late-time photoconductivity keeps growing with fluence rather than saturating, the assignment of the plateau to a density-independent Mott value would be contradicted.

Watch

Extended reading notes

Core claim

The central claim is that the photoexcited-carrier response of MAPI is organized by three density regimes whose boundaries move with temperature, and that the high-density cutoff is a genuine material property: the Mott density. Below roughly $10^{15}\,\mathrm{cm}^{-3}$, a fast picosecond decay of the band-edge bleach reflects shallow-trap capture, and the fraction of trapped carriers grows as temperature drops. From about $10^{15}$ to $10^{18}\,\mathrm{cm}^{-3}$ the shallow traps saturate and the carrier population decays slowly through bimolecular recombination, surviving for hundreds of picoseconds. Above roughly $10^{18}\,\mathrm{cm}^{-3}$, the photoconductivity peak rises sublinearly with fluence and relaxes within tens of picoseconds to a plateau; the authors interpret this as polaron-polaron annihilation until the density settles at the Mott density. They extract Mott densities from the crossing of the low-fluence linear photoconductivity calibration with the extrapolated late-time plateau, obtaining values from $(1.34\pm0.06)\times10^{18}\,\mathrm{cm}^{-3}$ at 78 K to $(4.5\pm0.1)\times10^{18}\,\mathrm{cm}^{-3}$ at 315 K. These numbers imply polaron radii of about 6.5 nm at 78 K shrinking to about 4 nm at 315 K, consistent with large polarons spanning multiple lattice constants.

Load-bearing premise

Every reported carrier density, including all Mott densities, assumes that the THz photoconductivity signal is directly proportional to carrier density with a mobility that does not change with density, so that a low-fluence linear calibration can be extrapolated to the high-density plateau; if mobility falls at high density, all density boundaries shift.

Editorial extensions

If this is right

Below about $10^{15}\,\mathrm{cm}^{-3}$, as in solar illumination, shallow-trap capture dominates, so the fast decay reports sample quality rather than intrinsic recombination.
Between $10^{15}$ and roughly $10^{18}\,\mathrm{cm}^{-3}$, carrier populations remain stable for hundreds of picoseconds, meaning bimolecular recombination governs the response.
Above the Mott density (from $1.3\times10^{18}$ to $4.5\times10^{18}\,\mathrm{cm}^{-3}$ depending on temperature), excess carriers annihilate within tens of picoseconds, capping the achievable free-carrier density.
The measured Mott densities lie at or above the calculated population-inversion thresholds, so MAPI can in principle support optical gain before annihilation truncates the carrier density.
The phase diagram implies that transient signals can change sign or timescale with fluence, so ultrafast studies must report carrier density before assigning a decay to a specific physical effect.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

A testable extension is that single-crystal MAPI, with far lower trap densities, should show the trap-dominated regime pushed to lower fluences, sharpening the low-density boundary of the phase diagram.
The same combined OPTP/TAS protocol applied to other halide perovskites would reveal whether the temperature dependence of the Mott density tracks lattice stiffness or the structural phase, a comparison the present data cannot settle.
Because Auger-type annihilation is expected to scale with the cube of carrier density, fitting the fast-decay amplitude versus initial density would directly test the annihilation mechanism assumed above the Mott density.
The proximity of the Mott density to the gain threshold suggests a practical design rule: perovskite lasers should operate just below $N_{\mathrm{Mott}}$, where bimolecular recombination still preserves a long-lived carrier population.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

The paper is a solid, systematic extension of the Bonn group's earlier Mott-polaron work. What is actually new is the temperature dependence of the Mott density (1.34e18 at 78 K to 4.5e18 at 315 K) and the three-regime T–N “phase diagram” built from combined TAS and OPTP. That is genuinely useful for the halide perovskite subfield. The measurements are extensive: densities from 1e14 to 1e19 cm–3, nine to ten temperatures, both probes, and a gain-threshold model benchmarked against Suárez et al. Credit where it is due: the low-density trap saturation data are clean, the separation of trap-dominated, direct-recombination, and Mott-dominated regimes is clearly presented, and the polaron radii derived from N_Mott are consistent with large polaron formation.

The soft spot is the density calibration. Every N_Mott value is read off by taking the late-time photoconductivity plateau and inserting it into a linear low-fluence calibration. That assumes THz photoconductivity is proportional to carrier density up to ~1e19 with a density-independent mobility and a fluence-independent photon-to-carrier quantum yield. The quantum yields (0.55 at 78 K, 0.30 at 292 K) come from low-fluence Drude-Smith fits. If mobility decreases at high densities, the sublinear peak and the plateau could be largely a mobility artifact, and the reported Mott densities would be shifted. The paper does not show high-density conductivity spectra, only low-fluence ones, and the TAS and OPTP density ranges barely overlap, so there is no independent cross-check of the absolute density scale. This is not an internal contradiction—the authors are transparent about the procedure—but it means the central phase diagram is interpretation-heavy rather than directly demonstrated. The “direct recombination-dominated” label is also inferred from the absence of trapping and Mott effects, not from a measured bimolecular coefficient, which is a minor point.

My take: this deserves a serious referee. The dataset is worth publishing, and the phase diagram is a reasonable synthesis, but the Mott-density values should be treated as provisional until the linearity assumption is checked—for example with high-density THz conductivity spectra or an overlapping TAS/OPTP density point. With that check, this would be a useful reference; without it, it is a well-documented map with a calibrated-but-unvalidated density axis. For a journal, I would send it to review with a request for those controls.

Referee Report

3 major / 6 minor

Summary. The paper combines optical-pump/THz probe (OPTP) spectroscopy and highly sensitive transient absorption spectroscopy (TAS) on polycrystalline methylammonium lead iodide (MAPbI3) films to map carrier dynamics over a nominal density range of 10^14–10^19 cm^-3 and temperatures from 78 to 315 K. The authors identify three density regimes—trap-dominated dynamics below ~10^15 cm^-3, direct/bimolecular recombination from ~10^15 to ~10^18 cm^-3, and a Mott-dominated regime above ~10^18 cm^-3 with fast polaron–polaron annihilation—and combine these with literature deep-trap densities and a gain-threshold model into a temperature-dependent electronic 'phase diagram'. The central quantitative output is the temperature-dependent Mott density extracted from OPTP transients, which increases from 1.34±0.06×10^18 cm^-3 at 78 K to 4.5±0.1×10^18 cm^-3 at 315 K.

Significance. If the Mott-density interpretation is correct, the paper provides a useful density–temperature map of carrier dynamics in a prototypical perovskite, with direct relevance to photovoltaic (low-density) and lasing (high-density) operating conditions. The experimental dataset is extensive, spanning five orders of magnitude in density and bridging TAS and OPTP, and the low-density TAS results on shallow-trap saturation are presented with clear density-dependent trends. The gain-threshold model is benchmarked against the prior work of Suárez et al., which strengthens the optical-gain comparison. However, the central quantitative claim—the absolute Mott densities and the phase boundaries in Fig. 4—rests on an assumption of density-independent THz photoconductivity that is not demonstrated in the high-density regime, and the intermediate 'direct recombination-dominated' regime is labeled without a quantitative recombination analysis. These gaps are significant but appear addressable with additional measurements and analysis.

major comments (3)

[SI 'Determination of the photoexcitation density'; Fig. 3 insets; Fig. S11] The extraction of N_Mott from the OPTP plateau assumes that the measured -ΔE/E is proportional to carrier density with a density-independent mobility up to ~10^19 cm^-3. The peak photoconductivity at low fluence is linearly fit to N, and the late-time plateau value is inserted into this same linear relation to read off N_Mott. At such high densities, carrier–carrier scattering, hot-carrier effects, and phase-space filling can all reduce the THz mobility; if the mobility decreases with density, the sublinear peak and the plateau could be explained without invoking polaron–polaron annihilation, and the reported Mott densities would be systematically affected. The paper provides no direct evidence for mobility constancy in this regime: the TAS and OPTP density ranges barely overlap, and the Drude-Smith parameters in Table S1 are obtained only at low fluence. This affects every blue datapoint in Fig. 4 and therefore the central phase diagram. Please provide either high-density THz conductivity spectra at several fluences (to extract the density dependence of the mobility) or an explicit cross-check of the absolute density scale in the overlap region.
[Fig. 1(a), Fig. 4, and 'Results and Discussion'] The intermediate regime is labeled 'direct recombination-dominated' even though no bimolecular recombination coefficient is extracted and no quantitative rate-equation fit is presented. The TAS traces in Fig. 2 show near-constant band-edge bleach over hundreds of picoseconds, but this is also consistent with a slow effective decay that is not uniquely assigned to direct recombination. The boundary between regime 2 and regime 3 therefore depends on the same OPTP plateau that is the subject of the previous comment. Fitting the density- and temperature-dependent TAS lifetimes (e.g., from Table S2) to a trap-assisted/direct/Auger model would substantiate the regime assignment and make the phase diagram more robust.
[SI 'Photoconductivity spectra and photon-to-carrier quantum yield'] The photon-to-carrier quantum yield Φ (0.55 at 78 K, 0.30 at 292 K) is determined from low-fluence Drude-Smith fits and then assumed to be fluence-independent in the conversion from absorbed photon density to carrier density for all measurements, including those above 10^18 cm^-3. Any fluence dependence of Φ would rescale every reported density, including the Mott densities on the phase diagram. At a minimum, the manuscript should state this assumption explicitly and discuss its likely magnitude (e.g., whether high-fluence Auger or hot-carrier processes could change the initial quantum yield).

minor comments (6)

[Fig. 1(b) caption] The density ranges for TAS and OPTP are given in the text as 3×10^14–4×10^16 cm^-3 and 5×10^16–10^19 cm^-3, leaving no overlap region between the two techniques; this lack of overlap is important for assessing the absolute density calibration and should be stated clearly in the caption or main text.
[Main text, first use of N_Mott (Fig. 3 inset)] The symbol N_Mott is used in the inset of Fig. 3 without an explicit definition; define it at first appearance and state that it is extracted from the late-time plateau of the OPTP signal.
[SI 'Determination of the photoexcitation density'] The equation for the absorbed carrier density is difficult to parse because of the inline formatting and undefined symbols (r, A, L, and the integral over depth). Please rewrite it with all variables defined and with balanced parentheses in the reflectance factor.
[SI 'Estimation of the temperature-dependent polaron size'] The polaron radius estimate assumes a fixed packing fraction of 0.74 and spherical wavefunctions; the resulting radii of ~4–6.5 nm therefore depend directly on this geometric assumption. This should be stated as a model-dependent estimate rather than a direct measurement.
[Conclusion and Fig. 1(a)] The text claims that the Mott density is 'an intrinsic material property' while trap densities are sample-dependent. This claim is not proven by the data; the extraction depends on the sample thickness, absorption coefficient, and the assumed linear photoconductivity–density relation. Please temper the claim or provide supporting evidence.
[Throughout] Minor typographical issues include 'Feymann's polaron theory' in the conclusion (should be 'Feynman's') and the inconsistent spacing in 'one-to-several ps' (should be 'one to several ps').

Circularity Check

0 steps flagged · score 2.0 of 10

No significant circularity; the central Mott-density values are a calibration-based readout, and the polaron interpretation relies on the authors' prior work only as interpretive support, not as a substitute for the new data.

full rationale

The paper's central quantitative outputs are the temperature-dependent Mott densities. These are obtained by fitting a line to the low-fluence peak photoconductivity versus nominal carrier density and then using that line to convert the late-time plateau photoconductivity into a density (main text: 'We utilize the linear relationship established for the photoconductivity at low pump fluences to determine NMott'). This is a calibration, not a reduction-by-construction: the low-fluence slope is fitted to data below roughly 10^18 cm^-3, and the Mott density is read off from a separate high-fluence observable. The result is only as secure as the assumption of a density-independent THz mobility and a fluence-independent quantum yield, but those are stated assumptions (SI: 'the photon-to-carrier quantum yield Φ ... determined independently' from Drude-Smith fits; 'we have used the room temperature refractive index' at the pump wavelength), not circular definitions. A density-dependent mobility would indeed invalidate the conversion, but that is a falsifiable experimental assumption, not a logical equivalence. The gain-threshold model is benchmarked against Suárez et al. and uses literature effective masses, so the claim that gain thresholds align with Mott densities is an independent comparison rather than a fitted coincidence. The polaron radius is a geometric conversion of the measured N_Mott (SI: 'N_Mott * 4/3 pi r_polaron^3 = 0.74'), so it is not an independent confirmation, but the paper presents it as an estimate rather than as a prediction derived from polaron theory. The high-density interpretation as polaron-polaron annihilation does cite the authors' own earlier work (refs. 32 and 35), but the present OPTP transients provide new experimental evidence—sublinear peak scaling, a late-time plateau, and fast decay—so the self-citation is interpretive support, not a logical substitution for the data. Overall, no equation in the derivation chain is equivalent to its own input by construction, and the central phase-diagram claim has independent experimental content.

Assumptions & free parameters 3 free parameters · 6 assumptions · 0 invented entities

The central numbers depend on the density calibration (Phi from Drude-Smith fits), the assumed proportionality between THz photoconductivity and carrier density, literature effective masses for the gain model, and the packing model for polaron radii. No new entities are introduced. One free parameter set (Phi) scales every reported density, and the polaron packing fraction is a hand-chosen geometric constant.

free parameters (3)

Photon-to-carrier quantum yield Phi = 0.30 (tetragonal), 0.55 (orthorhombic)
Fitted from Drude-Smith plasma frequency using an assumed effective mass; directly multiplies all reported carrier densities, including every N_Mott value.
Drude-Smith parameters (tau0, C, omega_p) = 78K: tau0=53+-10 fs, C=-0.55; 292K: tau0=72+-4 fs, C=-0.67
Fitted to low-fluence THz conductivity spectra to obtain the plasma frequency and hence the quantum yield; anchors the density calibration.
Polaron packing fraction = 0.74
Chosen by hand from the densest packing of spheres; used to translate N_Mott into polaron radii.

assumptions (6)

domain assumption Carrier density is calibrated using photon-to-carrier quantum yields Phi=0.30 (tetragonal) and Phi=0.55 (orthorhombic), obtained from Drude-Smith fits to low-fluence THz spectra.
Systematic error in Phi scales all densities; extracting Phi requires an assumed effective mass and the Drude-Smith model.
domain assumption The measured THz photoconductivity (-DE/E) remains proportional to carrier density with approximately constant mobility up to 1e19 cm-3.
Used to convert the late-time OPTP plateau into N_Mott; no high-density Drude-Smith calibration is provided.
domain assumption The fast few-picosecond TAS bleach decay is attributed to shallow-trap localization rather than to hot-carrier cooling or exciton formation.
Interpretation follows refs 3 and 29; no independent trap spectroscopy is performed on this film.
domain assumption Effective masses for electrons and holes are taken from cited literature (0.11/0.13, 0.15/0.18, 0.23/0.29 m0) for the gain-threshold model.
The calculated gain thresholds scale with these inputs; the paper benchmarks the code but does not measure masses in this sample.
domain assumption Polaron radius is derived from N_Mott via spherical packing with 0.74 volume fraction.
Geometric conversion; if polarons are not hard spheres or packing differs, the extracted radii change.
domain assumption The refractive index at 515 nm is assumed nearly temperature-independent for density calculations.
Stated with citation; small errors here shift the whole density scale.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Temperature- and charge carrier density-dependent electronic response in methylammonium lead iodide." pith.science (2026). https://pith.science/paper/U7JQ4MBQ

@misc{pith2026250518887,
  author       = {Pith},
  title        = {Pith review of: Temperature- and charge carrier density-dependent electronic response in methylammonium lead iodide},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/U7JQ4MBQ}},
  note         = {Machine review of arXiv:2505.18887}
}

read the original abstract

Understanding carrier dynamics in photoexcited metal-halide perovskites is key for optoelectronic devices such as solar cells (low carrier densities) and lasers (high carrier densities). Trapping processes at low carrier densities and many-body recombination at high densities can significantly alter the dynamics of photoexcited carriers. Combining optical-pump/THz probe and transient absorption spectroscopy we examine carrier responses over a wide density range (10^14-10^19 cm-3) and temperatures (78-315K) in the prototypical methylammonium lead iodide perovskite. At densities below ~10^15 cm-3 (room temperature, sunlight conditions), fast carrier trapping at shallow trap states occurs within a few picoseconds. As excited carrier densities increase, trapping saturates, and the carrier response stabilizes, lasting up to hundreds of picoseconds at densities around ~10^17 cm-3. Above 10^18 cm-3 a Mott transition sets in: overlapping polaron wavefunctions lead to ultrafast annihilation through an Auger recombination process occurring over a few picoseconds. We map out trap-dominated, direct recombination-dominated, and Mott-dominated density regimes from 78-315 K, ultimately enabling the construction of an electronic phase diagram. These findings clarify carrier behavior across operational conditions, aiding material optimization for optoelectronics operating in the low (e.g. photovoltaics) and high (e.g. laser) carrier density regimes.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

6 extracted references · 5 canonical work pages

[1]

State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, 457 Zhongshan Road, Dalian 116023, P. R. China

work page
[2]

University of Chinese Academy of Sciences, 19A Yuquan Road, Beijing 100049, P.R. China

work page
[3]

Max Planck Institute for Polymer Research, 55128 Mainz, Germany

work page
[4]

Department of Physics, Pohang University of Science and Technology, 37673, Pohang, Korea

work page
[5]

Leiden Institute of Chemistry, Leiden University, Einsteinweg 55, 2333CC, Leiden, the Netherlands

work page
[6]

𝜔#$𝜏"(1−𝑖𝜔𝜏

School of Physics, Xidian University, Xi’an 710071, P. R. China 2 Table of Contents SAMPLE PREPARATION .............................................................................................................. 3 ULTRAFAST SPECTROSCOPIC METHODS ......................................................................................... 4 Optical-pump/THz ...

arXiv

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.

[1] [1]

State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, 457 Zhongshan Road, Dalian 116023, P. R. China

work page

[2] [2]

University of Chinese Academy of Sciences, 19A Yuquan Road, Beijing 100049, P.R. China

work page

[3] [3]

Max Planck Institute for Polymer Research, 55128 Mainz, Germany

work page

[4] [4]

Department of Physics, Pohang University of Science and Technology, 37673, Pohang, Korea

work page

[5] [5]

Leiden Institute of Chemistry, Leiden University, Einsteinweg 55, 2333CC, Leiden, the Netherlands

work page

[6] [6]

𝜔#$𝜏"(1−𝑖𝜔𝜏

School of Physics, Xidian University, Xi’an 710071, P. R. China 2 Table of Contents SAMPLE PREPARATION .............................................................................................................. 3 ULTRAFAST SPECTROSCOPIC METHODS ......................................................................................... 4 Optical-pump/THz ...

arXiv