REVIEW 4 major objections 4 minor 39 references
Anomalous temperature dependence of polaron mobility in a nonlinear double-well potential: unbiased X-propagator approach
T0 review · 4 major / 4 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read A nonlinear double-well coupling makes polaron mobility concave in temperature and reproduces SrTiO3's anomalous transport.
desk verdict Solid unbiased numerics for a nonlinear polaron model, with a plausible but unproven bridge to SrTiO3. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The method is an extension of the numerical X-propagator technique, which samples partition-function diagrams in the electron-site and oscillator-displacement representation. Free-oscillator propagators $U$ and electron-occupied propagators $\tilde U$ carry the weight of the diagrams, and the current-current correlation function $C_{JJ}(\tau)$ is measured directly by inserting current operators at pairs of hopping vertices, giving the Matsubara correlator from the Monte Carlo average over kink-antikink configurations. The mobility spectrum $\mu(\omega)$ is then obtained by inverting the relation $$C_{JJ}(i\omega_n) = \frac{2}{\pi}\int_0^\infty d\omega\, \frac{\$omega^{2}$}{\$omega^{2}$+\$omega_n^{2}$}\,\mu(\omega)$$ with the Stochastic Optimization Method. The double-well potential defines two energy scales—the barrier height $W$ and the intra-well frequency $\tilde\omega$—which set the temperature regimes that produce the concave mobility.
What would settle it
Measure the optical conductivity spectrum of a dilute polar metal such as doped SrTiO3 across the temperature range from roughly 0.1 to 2 times the soft-mode energy. The model predicts that in the intermediate-temperature regime the spectrum has a single peak at nonzero frequency and no zero-frequency Drude peak; if the measured spectrum retains a Drude peak in that temperature range, the three-regime mechanism is falsified. Equally decisive would be a mobility measurement showing convex (saturating) temperature dependence in a clean sample where the model predicts concavity.
Extended reading notes
Core claim
The central discovery is that a polaron coupled to a double-well potential—a minimal model of a strongly anharmonic soft mode—predicts a mobility that decreases with temperature in a concave fashion at moderate coupling. In the model Hamiltonian, the electron-phonon coupling contains a negative quadratic term $g_2(2\Omega)x^2$ and a positive quartic term $g_4(2\Omega)^2 x^4$, which for $g_2<0$ produce a double-well potential with barrier height $W$. The mobility divides into three temperature regimes governed by the relation between thermal energy and $W$, and the optical conductivity $\mu(\omega)$ shows distinct spectral shapes in each: a double peak at low $T$, a single finite-frequency peak at intermediate $T$, and a broad zero-frequency peak at high $T$. For the experimentally relevant range $g_2$ between $-0.96$ and $-1.4$, the model reproduces both the concave mobility and the mean free path dropping below the lattice spacing, i.e., the Mott-Ioffe-Regel violation, in SrTiO3.
Load-bearing premise
The claim relies on the assumption that the one-dimensional, single-site, dispersionless double-well model, with parameters chosen so that its effective mass brackets the experimental $m^*/m \approx 3$ of SrTiO3, preserves the qualitative temperature dependence of the mobility despite neglecting dimensionality, phonon dispersion, and the temperature dependence of the soft mode.
Editorial extensions
If this is right
- For a range of couplings, the model predicts that the mean free path drops below one lattice spacing exactly in the temperature regime where the mobility falls fastest, explaining the simultaneous onset of Mott-Ioffe-Regel violation in SrTiO3.
- The three regimes have distinct signatures in the optical conductivity, so frequency-resolved measurements can be used to locate the crossover temperatures and infer the barrier height.
- The method is applicable to any polynomial electron-phonon interaction, opening the way to calculations for other anharmonic materials beyond the double-well example.
- The concave mobility is a unique fingerprint of nonlinear double-well coupling; seeing it in a material would indicate that a soft anharmonic mode dominates polaron scattering.
Reading between the lines
- Extending beyond the paper, the same mechanism should apply to other quantum paraelectrics such as KTaO3, where a similar concave mobility is observed, suggesting a general explanation for anomalous transport in dilute polar metals.
- If the model is correct, the temperature at which the mobility switches from slow to fast drop provides a direct estimate of the double-well barrier height $W$, turning transport data into a probe of lattice anharmonicity.
- A testable extension would be to compute the mobility for a two- or three-dimensional version of the model with a dispersive soft mode; if the concave profile survives, the mapping to real materials is more secure than the paper's single-mode assumption.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper extends the numerical X-propagator Monte Carlo method of Ref. [27] to the finite-temperature current-current correlation function for polarons with arbitrary nonlinear electron-phonon coupling. The method is applied to a one-dimensional double-well potential model (Eqs. 1-2) with quartic coupling g4>0 and negative quadratic coupling g2. The authors report a concave temperature dependence of the dc mobility at moderate coupling, with three regimes associated with the double-well barrier, in contrast to the convex dependence of linear-coupling polaron models. They also report a strong-coupling regime with nonmonotonic mobility and a correspondence between the temperature domains and the shape of the optical conductivity. Finally, they argue that for g2 between -0.96 and -1.4, selected by bracketing the experimental effective mass of SrTiO3, the model qualitatively reproduces the anomalous concave mobility and the onset of Mott-Ioffe-Regel violation in SrTiO3.
Significance. If the results are correct, this is a valuable methodological and conceptual contribution. The unbiased Monte Carlo sampling of partition-function diagrams with direct measurement of Matsubara correlators is a principled extension of the X-propagator approach and is benchmarked against independent Diagrammatic Monte Carlo for Holstein and purely quadratic models. The predicted three-regime mobility structure is a falsifiable qualitative signature of double-well electron-phonon coupling that is absent in linear-coupling models, and the connection to the SrTiO3 data is suggestive. The main weaknesses are the unsupported extrapolation from the 1D dispersionless fixed-frequency model to 3D SrTiO3 and the reliance on stochastic analytic continuation for the central mobility values; these points need to be addressed before the experimental claim can be accepted.
major comments (4)
- [Anomalous transport in SrTiO3] The claim that "none of the above differences can alter the general trend of concave behavior" is an unsupported assertion. The paper correctly acknowledges that the model is 1D, dispersionless, and has a fixed phonon frequency, whereas SrTiO3 has a dispersive soft mode whose frequency increases with temperature. The three-regime mobility structure and the coincidence between the onset of accelerated mobility suppression and MIR-limit violation are quantitative features that could shift with the phonon density of states and with a temperature-dependent Omega that rescales the phonon and barrier energies. Since the experimental relevance of the model is the central motivation of the paper, this transfer step needs support: either a calculation for a model with dispersive or temperature-dependent phonons, a scaling argument, or an explicit downgrading of the SrTiO3 comparison to a conjecture.
- [Method] The paper repeatedly describes the calculation as "unbiased" and "approximation-free" (Abstract, Section "Method", Conclusions), but the mobility spectra and the dc mobility are reconstructed from Matsubara data via the Stochastic Optimization Method, whose systematic errors are not quantified. The benchmark in Section "Method" only verifies that the Matsubara correlators agree with Diagrammatic Monte Carlo; it does not verify the continued spectra or the dc mobility. Given that the central claim of a concave temperature dependence rests on the values of mu(omega=0) obtained by this continuation, the paper should state this limitation explicitly, provide reproducibility tests across independent SOM runs, and ideally compare the continued dc mobility with an independent calculation or sum-rule constraint.
- [Model and parameters] The "experimentally relevant" range g2 in [-0.96, -1.4] is selected by matching the polaron mass enhancement from the 1D calculation of Ref. [27] (m*/m about 1.6 and 5.8) to the experimental m*/m less than or about 3 in SrTiO3. However, the experimental mass is a 3D thermodynamic mass that depends on carrier density, while the theoretical mass is a 1D quasiparticle mass; no argument establishes that these masses are comparable. The selection of g2 therefore rests on an unverified identity, and the bracket is sensitive to this assumption. This weakens the claim that the model parameters are "close to their experimental values."
- [Figure 3] The mean free path estimate lambda_MFP approximately equal to mu sqrt(C_JJ(0)/<K_xx>) is introduced without derivation and without specifying whether the polaron renormalization of the mass and velocity is included. Since the MIR-violation claim depends on lambda_MFP dropping below the lattice spacing, an order-of-magnitude error in this estimate could shift the onset relative to experiment. Please clarify the derivation and the renormalization convention, or cite a source that establishes this formula for the polaron case.
minor comments (4)
- [Model and parameters] The condition for a double well, g2 < -Omega/4, is stated without derivation; a one-line derivation from the curvature of V(x) at x=0 would improve readability.
- [Method] The term "X-propagator" is not defined at first use in the Abstract; define it or explicitly refer to Ref. [27] at the first occurrence.
- [Eq. (5)] For Eq. (5), it would be helpful to state that the Monte Carlo average is over partition-function diagrams normalized by the partition function, since the factor 1/beta alone does not make this explicit.
- [General] The paper would benefit from a brief discussion of the energy scale in the comparison with SrTiO3: the model units set t=1 and Omega=0.25, but the experimental temperatures (about 50 K) are not mapped to these units, so the qualitative nature of the comparison should be reinforced.
Circularity Check
No significant circularity: the mobility is a genuine unbiased computation; the only self-referential element is a benign calibration of the g2 range to prior mass estimates, not a fitted prediction.
full rationale
The central claim—concave temperature dependence of polaron mobility in the double-well model—is the output of an unbiased Monte Carlo calculation of the current–current correlator (Eqs. 3–6), benchmarked against independent DMC results for linear Holstein and quadratic models. No fitting to SrTiO3 mobility data is involved. The only self-referential element is the selection of the 'experimentally relevant' g2∈[−0.96,−1.4] range using mass-enhancement estimates m*/m≈1.6–5.8 from the same authors' prior work [27]; this is a calibration to a separate observable (effective mass), not an input that forces the concave mobility shape. The paper's own admission that dimensionality, dispersion, and Ω(T) differ in SrTiO3, followed by the assertion that 'none of the above differences can alter the general trend,' is a robustness assumption, not a circular reduction. Accordingly, no load-bearing argument reduces to its own input by construction.
Assumptions & free parameters
free parameters (3)
- g2 (quadratic electron-phonon coupling) =
-0.2, -0.45, -0.8, -0.96, -1.4, -1.8; experimentally relevant range -0.96 to -1.4
- g4 (quartic electron-phonon coupling) =
0.1
- Omega (bare phonon frequency) =
0.25 in units of t = 1
assumptions (4)
- domain assumption The single-site dispersionless Einstein oscillator with Hamiltonian (1)-(2) is an adequate minimal model for the soft anharmonic mode in SrTiO3 and KTaO3.
- domain assumption The 1D model's temperature dependence of mobility is representative of 3D systems; dimensionality only changes density of states and does not alter the concave trend.
- domain assumption The X-propagator Monte Carlo estimates of Matsubara correlators are unbiased and converged, and the Stochastic Optimization Method analytic continuation provides reliable mu(omega) and mu(omega->0).
- domain assumption Mott-Ioffe-Regel mean free path estimate lambda_MFP approximately mu sqrt(CJJ(0)/<Kxx>) is a valid estimator of the carrier mean free path in this model.
Cite this review
Pith. "Pith review of Anomalous temperature dependence of polaron mobility in a nonlinear double-well potential: unbiased X-propagator approach." pith.science (2026). https://pith.science/paper/K6WLHTIR
@misc{pith2026260809883,
author = {Pith},
title = {Pith review of: Anomalous temperature dependence of polaron mobility in a nonlinear double-well potential: unbiased X-propagator approach},
year = {2026},
howpublished = {\url{https://pith.science/paper/K6WLHTIR}},
note = {Machine review of arXiv:2608.09883}
}
abstract
We develop an unbiased X-propagator method for calculating finite-temperature optical conductivity $\sigma(\omega)$ and dc mobility $\mu$ for arbitrary nonlinear electron-phonon interaction. We apply it to a polaron coupled to a double-well lattice potential, a minimal model for strongly anharmonic polar materials. At moderate coupling, the mobility exhibits three temperature regimes associated with confinement within one well, thermal competition with the barrier, and barrier-insensitive high-temperature dynamics. This sequence produces a concave temperature dependence of the mobility that is absent in conventional linear-coupling polaron models. At strong coupling, the mobility becomes nonmonotonic, and its temperature evolution is reflected in a characteristic redistribution of optical spectral weight. For parameters relevant to SrTiO$_3$, our results qualitatively reproduce both the anomalous concave mobility and the onset of violation of the Mott-Ioffe-Regel limit, thereby supporting nonlinear coupling to a soft anharmonic lattice mode as a microscopic mechanism for anomalous transport in dilute polar metals.
Figures
Reference graph
Works this paper leans on
- [27]
-
[1]
C. Collignon, X. Lin, C. W. Rischau, B. Fauqu´ e, and K. Behnia, Metallicity and superconductivity in doped strontium titanate, Annual Review of Condensed Matter Physics10, 25 (2019)
work page 2019
-
[2]
M. J. Schilcher, P. J. Robinson, D. J. Abramovitch, L. Z. Tan, A. M. Rappe, D. R. Reichman, and D. A. Eg- ger, The significance of polarons and dynamic disorder in halide perovskites, ACS Energy Letters6, 2162 (2021)
work page 2021
-
[3]
O. N. Tufte and P. W. Chapman, Electron mobility in semiconducting strontium titanate, Phys. Rev.155, 796 (1967)
work page 1967
-
[4]
X. Lin, C. W. Rischau, L. Buchauer, A. Jaoui, B. Fauqu´ e, and K. Behnia, Metallicity without quasi-particles in 6 room-temperature strontium titanate, npj Quantum Ma- terials2, 41 (2017)
work page 2017
-
[5]
X. Lin, B. Fauqu´ e, and K. Behnia, ScalableT2 resistivity in a small single-component Fermi surface, Science349, 945 (2015)
work page 2015
-
[6]
C. Collignon, P. Bourges, B. Fauqu´ e, and K. Behnia, Heavy nondegenerate electrons in doped strontium ti- tanate, Phys. Rev. X10, 031025 (2020)
work page 2020
-
[7]
S. H. Wemple, Some transport properties of oxygen- deficient single-crystal potassium tantalate (KTaO 3), Phys. Rev.137, A1575 (1965)
work page 1965
Show all 39 references
-
[8]
Sakai, T
A. Sakai, T. Kanno, S. Yotsuhashi, H. Adachi, and Y. Tokura, Thermoelectric properties of electron-doped KTaO3, Japanese Journal of Applied Physics48, 097002 (2009)
2009
-
[9]
Kumar, V
A. Kumar, V. I. Yudson, and D. L. Maslov, Quasipar- ticle and nonquasiparticle transport in doped quantum paraelectrics, Phys. Rev. Lett.126, 076601 (2021)
2021
-
[10]
J.-J. Zhou, O. Hellman, and M. Bernardi, Electron- phonon scattering in the presence of soft modes and elec- tron mobility in SrTiO 3 perovskite from first principles, Phys. Rev. Lett.121, 226603 (2018)
2018
-
[11]
Y. Luo, J. Park, and M. Bernardi, First-principles dia- grammatic Monte Carlo for electron–phonon interactions and polaron, Nature Physics21, 1275 (2025)
2025
-
[12]
Gunnarsson, M
O. Gunnarsson, M. Calandra, and J. E. Han, Colloquium: Saturation of electrical resistivity, Rev. Mod. Phys.75, 1085 (2003)
2003
-
[13]
A. S. Mishchenko, N. Nagaosa, G. De Filippis, A. de Can- dia, and V. Cataudella, Mobility of Holstein polaron at finite temperature: An unbiased approach, Phys. Rev. Lett.114, 146401 (2015)
2015
-
[14]
beyond quasiparticles
A. S. Mishchenko, L. Pollet, N. V. Prokof’ev, A. Kumar, D. L. Maslov, and N. Nagaosa, Polaron mobility in the “beyond quasiparticles” regime, Phys. Rev. Lett.123, 076601 (2019)
2019
-
[15]
Esposito, M
V. Esposito, M. Fechner, R. Mankowsky, H. Lemke, M. Chollet, J. M. Glownia, M. Nakamura, M. Kawasaki, Y. Tokura, U. Staub, P. Beaud, and M. F¨ orst, Nonlin- ear electron-phonon coupling in doped manganites, Phys. Rev. Lett.118, 247601 (2017)
2017
-
[16]
Hoesch, P
M. Hoesch, P. Piekarz, A. Bosak, M. Le Tacon, M. Krisch, A. Koz lowski, A. M. Ole´ s, and K. Parlin- ski, Anharmonicity due to electron-phonon coupling in magnetite, Phys. Rev. Lett.110, 207204 (2013)
2013
-
[17]
W. A. Saidi, S. Ponc´ e, and B. Monserrat, Temperature dependence of the energy levels of methylammonium lead iodide perovskite from first-principles, The Journal of Physical Chemistry Letters7, 5247 (2016)
2016
-
[18]
Zacharias, G
M. Zacharias, G. Volonakis, F. Giustino, and J. Even, Anharmonic electron-phonon coupling in ultrasoft and locally disordered perovskites, npj Computational Mate- rials9, 153 (2023)
2023
-
[19]
Migoni, H
R. Migoni, H. Bilz, and D. B¨ auerle, Origin of Raman scattering and ferroelectricity in oxidic perovskites, Phys. Rev. Lett.37, 1155 (1976)
1976
-
[20]
K. G. Nazaryan and M. V. Feigel’man, Conductivity and thermoelectric coefficients of doped SrTiO 3 at high tem- peratures, Phys. Rev. B104, 115201 (2021)
2021
-
[21]
Ranalli, C
L. Ranalli, C. Verdi, M. Zacharias, J. Even, F. Giustino, and C. Franchini, Electron mobilities in SrTiO 3 and KTaO3: Role of phonon anharmonicity, mass renor- malization, and disorder, Phys. Rev. Mater.8, 104603 (2024)
2024
-
[22]
Esswein and N
T. Esswein and N. A. Spaldin, Ferroelectric, quantum paraelectric, or paraelectric? calculating the evolution from BaTiO3 to SrTiO3 to KTaO3 using a single-particle quantum mechanical description of the ions, Phys. Rev. Res.4, 033020 (2022)
2022
-
[23]
C. P. J. Adolphs and M. Berciu, Single-polaron properties for double-well electron-phonon coupling, Phys. Rev. B 89, 035122 (2014)
2014
-
[24]
C. P. J. Adolphs and M. Berciu, Going beyond the linear approximation in describing electron-phonon coupling: Relevance for the Holstein model, Europhysics Letters 102, 47003 (2013)
2013
-
[25]
C. P. J. Adolphs and M. Berciu, Strongly bound yet light bipolarons for double-well electron-phonon cou- pling, Phys. Rev. B90, 085149 (2014)
2014
-
[26]
Ragni, T
S. Ragni, T. Hahn, Z. Zhang, N. Prokof’ev, A. Kuklov, S. Klimin, M. Houtput, B. Svistunov, J. Tempere, N. Na- gaosa, C. Franchini, and A. S. Mishchenko, Polaron with quadratic electron-phonon interaction, Phys. Rev. B107, L121109 (2023)
2023
-
[28]
Supplemental material (2026)
2026
-
[29]
B. A. A. Martin and J. M. Frost, Multiple phonon modes in Feynman path-integral variational polaron mobility, Phys. Rev. B107, 115203 (2023)
2023
-
[30]
R. A. Cowley, Temperature dependence of a transverse optic mode in strontium titanate, Phys. Rev. Lett.9, 159 (1962)
1962
-
[31]
W. G. Stirling, Neutron inelastic scattering study of the lattice dynamics of strontium titanate: harmonic models, Journal of Physics C: Solid State Physics5, 2711 (1972)
1972
-
[32]
J. L. M. van Mechelen, D. van der Marel, C. Grimaldi, A. B. Kuzmenko, N. P. Armitage, N. Reyren, H. Hage- mann, and I. I. Mazin, Electron-phonon interaction and charge carrier mass enhancement in SrTiO 3, Phys. Rev. Lett.100, 226403 (2008)
2008
-
[33]
McCalla, M
E. McCalla, M. N. Gastiasoro, G. Cassuto, R. M. Fer- nandes, and C. Leighton, Low-temperature specific heat of doped SrTiO 3: Doping dependence of the effective mass and Kadowaki-Woods scaling violation, Phys. Rev. Mater.3, 022001(R) (2019)
2019
-
[34]
Mori, Transport, collective motion, and brownian mo- tion, Progress of Theoretical Physics33, 423 (1965)
H. Mori, Transport, collective motion, and brownian mo- tion, Progress of Theoretical Physics33, 423 (1965)
1965
-
[35]
Mori, A continued-fraction representation of the time- correlation functions, Progress of Theoretical Physics34, 399 (1965)
H. Mori, A continued-fraction representation of the time- correlation functions, Progress of Theoretical Physics34, 399 (1965)
1965
-
[36]
Goleˇ z, J
D. Goleˇ z, J. Bonˇ ca, L. Vidmar, and S. A. Trugman, Re- laxation dynamics of the Holstein polaron, Phys. Rev. Lett.109, 236402 (2012)
2012
-
[37]
A. S. Mishchenko, N. V. Prokof’ev, A. Sakamoto, and B. V. Svistunov, Diagrammatic quantum Monte Carlo study of the Fr¨ ohlich polaron, Phys. Rev. B62, 6317 (2000)
2000
-
[38]
Goulko, A
O. Goulko, A. S. Mishchenko, L. Pollet, N. Prokof’ev, and B. Svistunov, Numerical analytic continuation: An- swers to well-posed questions, Phys. Rev. B95, 014102 (2017)
2017
-
[39]
Anomalous temperature dependence of polaron mobility in a nonlinear double-well potential: unbiased X-propagator approach
I. Krivenko and A. S. Mishchenko, TRIQS/SOM 2.0: Implementation of the stochastic optimization with con- sistent constraints for analytic continuation, Computer Physics Communications280, 108491 (2022). Supplemental material for “Anomalous temperature dependence of polaron mob...
2022
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.