{"id":"1f63cd4a-2ece-483d-ad43-f146dd9dca6d","arxiv_id":"2608.09883","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A new unbiased Monte Carlo method shows that a polaron coupled to a double-well lattice potential yields a concave temperature-dependent mobility, qualitatively matching the anomalous transport and Mott-Ioffe-Regel violation in SrTiO3.","lead":"This paper computes how an electron moving through a lattice with a double-well-shaped distortion scatters, and it shows that this simple anharmonicity can produce the puzzling temperature dependence of conductivity seen in doped strontium titanate. The message for a generalist: a nonlinear lattice mode, rather than exotic scattering, may explain decades of anomalous transport in dilute polar metals.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central SrTiO3 claim rests on an untested assertion that the 1D, dispersionless, fixed-Ω double-well result survives in 3D with a dispersive, temperature-dependent soft mode; no calculation or scaling argument supports this transfer.","rationale":"Good-faith reading: the X-propagator extension is a serious methodological contribution, and the internal benchmark against DMC for Holstein and quadratic models lends credibility to the computed Matsubara correlators. I agree with the reader's identification of the weakest assumption: the mapping to SrTiO3. The three-regime concave mobility is the central physical claim; the experimental comparison is its main significance. The authors go further than a caveat by asserting that dimensionality, dispersion, and Ω(T) 'cannot alter' the trend, yet they provide no dimensional analysis, no effective-model reduction, and no higher-dimensional or dispersive-mode calculation. This is precisely the kind of robustness check that would convert a plausible toy-model mechanism into a supported explanation of the SrTiO3 anomaly. The missing error bars on the analytic continuation are a secondary concern; the benchmark of Matsubara correlators and the maturity of the SOM make it less likely that the concave shape is a numerical artifact, although adding continuation uncertainty would strengthen the paper. The concrete test proposed above directly targets the transition from 1D single-site to 3D dispersive reality. If the test is passed, the conditional verdict can be upgraded; if it fails, the experimental claim must be downgraded. Therefore the reader's CONDITIONAL verdict remains appropriate, with the condition being explicit robustness to the neglected physics.","tokens_in":11903,"tokens_out":14421,"duration_ms":137657,"concrete_test":"Perform the same unbiased calculation on a 2D or 3D cubic lattice with a dispersive soft phonon, e.g., ω(q)=Ω(1+ξΣ_i(1−cos q_i a)), and with Ω(T) tuned to mimic the measured Γ15 mode in SrTiO3, for g2=−0.96 and −1.4, and recompute μ(T) and λ_MFP. If the concave three-regime shape and the MIR-onset coincidence with the steep high-T drop do not survive, the proposed mechanism is not robust to the approximations required for SrTiO3.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's claim that its double-well model 'qualitatively reproduces' SrTiO3 transport (Abstract and section 'Anomalous transport in SrTiO3') requires that the 1D, single-site, dispersionless, temperature-independent Ω model (Eqs. 1-2) remains representative of a 3D system with a dispersive Γ15 soft mode whose frequency rises with T. The authors explicitly acknowledge these differences but assert, without calculation, that 'none of the above differences can alter the general trend of concave behavior' because the effect arises from the interplay of phonon energy and barrier height. This is the load-bearing step: the three-regime mobility curve and the MIR-onset coincidence are quantitative features that could easily shift with the phonon density of states and with a temperature-dependent Ω that rescales the regime boundaries. The parameter selection g2∈[-0.96,-1.4] is likewise tied to polaron-mass estimates from the 1D calculation in [27], whereas the experimental m*/m is a thermodynamic mass. No test of robustness to dimensionality, dispersion, or Ω(T) is presented. The theoretical model result may be correct, but the experimental relevance claim is unsupported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":12151,"tokens_out":7378,"duration_ms":69837,"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":[{"comment":"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.","section":"Anomalous transport in SrTiO3"},{"comment":"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.","section":"Method"},{"comment":"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.\"","section":"Model and parameters"},{"comment":"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.","section":"Figure 3"}],"minor_comments":[{"comment":"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.","section":"Model and parameters"},{"comment":"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.","section":"Method"},{"comment":"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.","section":"Eq. (5)"},{"comment":"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.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid methodological contribution, and the concave mobility in the double-well model is a genuine theoretical finding. My main reservation is the SrTiO3 claim: the extrapolation from the 1D dispersionless fixed-Omega model to 3D dispersive and temperature-dependent soft modes is unsupported, and this extrapolation carries the experimental significance. If the authors either provide supporting calculations or explicitly reframe the SrTiO3 comparison as a conjecture, the paper could become acceptable. In its current form, the central claim outstrips the evidence."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The real news here is the method: the X-propagator approach extended to finite-temperature current–current correlations gives unbiased optical conductivity and mobility for arbitrary electron–phonon coupling, benchmarked against Diagrammatic Monte Carlo for Holstein and quadratic models. That is a genuine, useful step. The three-regime concave mobility in the double-well model is also new as a theoretical result—it is absent from the linear-coupling polaron literature and from the purely quadratic models. I buy the physics of the model: low-T confinement in one well, a slow middle regime set by the barrier, and a barrier-insensitive high-T drop. The spectral-weight redistribution at strong coupling is a nice consistency check.\n\nWhere I wince is the SrTiO3 section. The authors correctly list the differences between the model and the real material—1D vs 3D, dispersionless vs dispersive soft mode, temperature-independent vs rising phonon frequency—and then simply assert that none of these can change the concave trend. That is an assertion, not an argument. It could easily fail: the three-regime boundaries depend on the phonon density of states and on how the barrier is renormalized by the temperature-dependent mode. The parameter selection g2 in [-0.96,-1.4] is based on mass enhancement from the same 1D model, so the bracket around the experimental m*/m is partially self-referential. The claim of reproducing the Mott–Ioffe–Regel onset is suggestive but again relies on that transfer. This does not sink the theoretical paper, but the abstract's phrase \"supporting nonlinear coupling ... as a microscopic mechanism\" is stronger than the evidence.\n\nMinor points: the analytic continuation uses the Stochastic Optimization Method, which is reputable, but no error bars are shown on the spectra, and I would have liked to see the raw Matsubara data or at least a statement about the continuation quality. No code or data are provided, which is a shame but common for this community.\n\nBottom line: the numerical method and the model result are solid; the experimental interpretation is plausible but not load-bearing. A referee should ask the authors to either soften the SrTiO3 language or supply a scaling/robustness check (e.g., a dispersive phonon or temperature-dependent Ω) before publication. This deserves full peer review, not desk rejection. I would bring it to a reading group interested in anharmonic polarons, and I would cite the method if I were working in that area.","headline":"Solid unbiased numerics for a nonlinear polaron model, with a plausible but unproven bridge to SrTiO3.","tokens_in":12693,"tokens_out":1548,"would_cite":true,"duration_ms":17576,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["71.38.-k","72.10.Di"],"model":"deepseek-v4-flash","headline":"A nonlinear double-well coupling makes polaron mobility concave in temperature and reproduces SrTiO3's anomalous transport.","keywords":["polaron","nonlinear electron-phonon interaction","double-well potential","anharmonic lattice","mobility","optical conductivity","Mott-Ioffe-Regel limit","SrTiO3"],"falsifier":"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.","tokens_in":11729,"feed_emoji":"⚛️","tokens_out":6588,"duration_ms":59418,"temperature":0.7,"pith_summary":"An unbiased numerical method, the X-propagator approach, is extended to compute the finite-temperature optical conductivity and dc mobility for arbitrary nonlinear electron-phonon interactions. Applied to a polaron in a double-well lattice potential, the calculation shows that at moderate coupling the mobility's temperature dependence is concave, with three regimes: a fast low-temperature drop, a slower intermediate drop when thermal energy is comparable to the barrier, and a fast, barrier-insensitive drop at high temperatures. This concave shape is absent in conventional linear-coupling polaron models and matches the anomalous mobility of doped SrTiO3. The paper argues that nonlinear coupling to a soft anharmonic lattice mode provides a microscopic mechanism for the anomalous transport and the onset of Mott-Ioffe-Regel violation in dilute polar metals.","feed_headline":"Double-well lattice coupling turns polaron mobility concave","feed_subtitle":"Three temperature regimes reproduce SrTiO3's anomalous transport and Mott-Ioffe-Regel violation.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the X-propagator formalism and the double-well model parameters that the present method extends to current-current correlations.","marker":"[27]"},{"why":"Provides the benchmark linear-coupling Holstein polaron mobility and the convex behavior that the double-well result is contrasted with.","marker":"[13]"},{"why":"Provides the purely quadratic electron-phonon model used as an additional benchmark for the correlators.","marker":"[26]"},{"why":"Furnishes the quadratic-phonon scattering mechanism and the $T^{-2}$ to $T^{-3}$ crossover background for SrTiO3.","marker":"[9]"},{"why":"Gives the experimental review of mobility, resistivity, and Mott-Ioffe-Regel violation in doped SrTiO3.","marker":"[1]"},{"why":"Reports the high-temperature metallic transport and the violation of the Mott-Ioffe-Regel limit in SrTiO3 that the model reproduces.","marker":"[4]"}],"fun_headline_variants":["Double-well polaron model yields concave mobility","Three temperature regimes for polaron in double well","Anharmonic double well reproduces SrTiO3 transport anomaly","Polaron mobility in double well: anomalous and concave","Nonlinear coupling to soft mode explains SrTiO3 mobility"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Double-well polaron model yields concave mobility","Three temperature regimes for polaron in double well","Anharmonic double well reproduces SrTiO3 transport anomaly","Polaron mobility in double well: anomalous and concave","Nonlinear coupling to soft mode explains SrTiO3 mobility"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000235,"raw_usage":{"total_tokens":1498,"prompt_tokens":939,"completion_tokens":559,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":555,"completion_tokens_details":{"reasoning_tokens":483}},"tokens_in":555,"tokens_out":559,"duration_ms":5451,"temperature":1.0,"reasoning_tokens":483,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T04:58:15.455065+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Ragni, T","cited_arxiv_id":null,"evidence_quote":"Supplies the X-propagator formalism and the double-well model parameters that the present method extends to current-current correlations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the benchmark linear-coupling Holstein polaron mobility and the convex behavior that the double-well result is contrasted with."},{"cited_title":"Kumar, V","cited_arxiv_id":null,"evidence_quote":"Furnishes the quadratic-phonon scattering mechanism and the $T^{-2}$ to $T^{-3}$ crossover background for SrTiO3."},{"cited_title":"Collignon, X","cited_arxiv_id":null,"evidence_quote":"Gives the experimental review of mobility, resistivity, and Mott-Ioffe-Regel violation in doped SrTiO3."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports the high-temperature metallic transport and the violation of the Mott-Ioffe-Regel limit in SrTiO3 that the model reproduces."}],"review_version":1}