{"id":"2b281ec0-bcca-4151-ae37-97252fbebd0b","arxiv_id":"2606.13833","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Machine-learned nanosecond simulations show small polarons on reduced rutile TiO2(110) hop almost exclusively within subsurface rows, suppressing surface mobility by orders of magnitude versus bulk.","lead":"Using a machine-learned potential trained on quantum mechanical simulations, the authors followed individual electron polarons on an oxygen-deficient titanium dioxide surface for nanoseconds. They find that surface polaron mobility is suppressed by orders of magnitude relative to the bulk crystal, matching measurements on porous samples.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Room-temperature mobility and 'orders of magnitude' suppression depend on a 400–700 K Arrhenius extrapolation from trajectories with zero hops at 300 K; the fit itself may be biased by a dominant subsurface channel and unverified Arrhenius curvature.","rationale":"The reader's weakest assumption is exactly the Arrhenius extrapolation from 400–700 K to 300 K, and my analysis of the manuscript confirms that this is the most load-bearing point. The strongest claim — quantitative room-temperature mobility and orders-of-magnitude suppression — depends on that extrapolation. The concern is not manufactured: the paper explicitly states 'no hopping event was observed over a total simulation time of ~40 ns' at 300 K (Sec. III.A) and the Arrhenius fit is the only bridge to room temperature. The paper also contains independent evidence that the low-T transport regime may differ from the high-T diffusive regime: Sec. III.B shows polaron–vacancy trapping is strongly temperature-dependent, with 400 K already showing pronounced localization near the vacancy, which suggests the 400 K mobility may mix drift/trapping with diffusion. The activation barrier (386 meV) is compared to infrared spectroscopy estimates, but the mobility prefactor and its temperature dependence are not independently verified. The absence of error bars on the mobility and the small number of fit points further weaken the extrapolation, but these are secondary to the regime-change risk. My recommendation is CONDITIONAL, matching the reader's verdict, because the concern is addressable (run 350 K simulations or analyze 300 K trajectories) and does not invalidate the qualitative story of suppressed surface mobility at high T or the topological explanation. I agree with the reader's identification of the weakest assumption, and no additional load-bearing concern was identified that would change the verdict.","tokens_in":15126,"tokens_out":1763,"duration_ms":16813,"concrete_test":"Compute the polaron mobility at 350 K (and, if feasible, 320 K) with LEOPOLD using long trajectories; if the 350 K mobility falls on the same Arrhenius line as 400–700 K within statistical uncertainty, the extrapolation is supported. Alternatively, re-analyze the existing 300 K trajectories for any localization-length/trapping signature (e.g., distribution of polaron–vacancy distances, waiting-time distribution, or MSD subdiffusive exponent) and compare with the 400 K behavior; if the 300 K trajectory shows increasing localization or subdiffusion, report the 300 K value as an upper bound rather than an Arrhenius-based number.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The headline quantitative result — a room-temperature surface mobility of 1.9×10^-5 cm2/Vs and a 'several orders of magnitude' suppression relative to bulk — is obtained by extrapolating an Arrhenius fit to LEOPOLD mobilities computed at 400–700 K (Fig. 4), because the authors state in Sec. III.A that 'no hopping event was observed over a total simulation time of ~40 ns' at 300 K. This is the load-bearing step: if the Arrhenius form is not valid down to 300 K (e.g., due to polaron–vacancy binding that becomes increasingly dominant at low T, sub-Arrhenius quantum tunneling, or a crossover to trap-limited transport), then the reported room-temperature value and the inferred suppression relative to bulk are both unsupported. The paper itself provides evidence that the trapping regime is temperature-dependent: Sec. III.B shows the polaron–vacancy distance distribution at 400 K is strongly peaked near the vacancy when initialized nearby, implying the diffusive regime sampled by the mobility fit at 400 K may already be contaminated by drift/trapping. The Arrhenius fit is also performed on only four temperatures (400–700 K, with only 2–3 degrees of freedom), with no uncertainty estimates and no stated check of fit residuals or curvature. Because the central comparison to experiment is quantitative ('very good agreement' with porous-TiO2 range), the argument's weight rests on an extrapolation whose validity is neither demonstrated nor stress-tested.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript extends the LEOPOLD machine-learned polaron-dynamics framework from bulk rutile TiO2 to the reduced rutile TiO2(110) surface. The authors construct a polaron-aware NequIP-style potential trained on ~100 ps of DFT+U FPMD with active learning, then perform multiple nanosecond MLMD trajectories at 300–700 K. From mean-square-displacement fits they extract surface polaron mobilities and report a suppression of several orders of magnitude relative to bulk, a dominant subsurface [001] hopping channel, rare interlayer hops, and a temperature-dependent polaron–vacancy attraction. An Arrhenius extrapolation gives a room-temperature mobility of 1.9×10^-5 cm2 V^-1 s^-1, which the authors state is in very good agreement with porous-TiO2 experiments.","tokens_in":15503,"tokens_out":6422,"duration_ms":76861,"significance":"If the quantitative claims hold, this work is significant: it demonstrates a transferable machine-learning strategy for simulating polaron dynamics at oxide surfaces on nanosecond time scales and provides a microscopic rationalization for the much lower electron mobilities observed in porous versus single-crystal rutile TiO2. The qualitative pathway analysis — preferred subsurface [001] hopping and rare interlayer events — is consistent with earlier short FPMD studies and is a useful advance in its own right. The model is not circular in the objectionable sense: it is trained on DFT data, and the porous-TiO2 comparison is an external experimental benchmark. However, the headline room-temperature number and the 'orders of magnitude' claim at 300 K rest on an Arrhenius extrapolation from temperatures at which hopping is actually observed, with no uncertainty quantification. This limits the current support for the central quantitative conclusion and needs to be addressed before the paper can be accepted.","major_comments":[{"comment":"The room-temperature value μ = 1.9×10^-5 cm2 V^-1 s^-1 and the claimed 'very good agreement' with porous-TiO2 experiments rest entirely on extrapolating the Arrhenius fit of mobilities computed at 400, 500, 600 and 700 K down to 300 K, because the authors state that 'no hopping event was observed over a total simulation time of ~40 ns' at 300 K. This extrapolation is load-bearing. The paper gives no uncertainty for E_a (386 meV) or μ, no fit residuals, and no check of Arrhenius curvature. I request confidence intervals and a discussion of whether sub-Arrhenius behavior, quantum tunneling, or trap-limited transport could set in below 400 K. If longer or accelerated sampling at 300 K is infeasible, the conclusions should distinguish the directly simulated 400–700 K suppression from the extrapolated room-temperature value, and the 'very good agreement' wording should be softened accordingly","section":"Sec. III.A and Fig. 4"},{"comment":"Figure 5(b) shows that at 400 K the polaron–vacancy distance distribution initialized near V_O is strongly peaked at short separations, indicating trap-dominated local exploration rather than homogeneous diffusion. If the 400 K MSD includes this trapping/drift component, the diffusion coefficient entering the Arrhenius fit is not the free-hopping mobility, and the extracted 386 meV barrier mixes trapping and hopping contributions. The claim that the suppression arises from the 'loss of favourable hopping pathways' rather than from the vacancy itself requires a control calculation (e.g., a surface without the vacancy, or trajectories initialized far from V_O with demonstrated linear MSD over the fitted window) or a quantitative estimate of the trapping contribution. As written, the microscopic interpretation is not fully disentangled from the defect potential.","section":"Sec. III.B"},{"comment":"The quantitative mobility estimates lack statistical uncertainties. Each temperature is the average of four independent runs, but Fig. 4 shows no error bars, and the text does not report the number of hopping events per run, the duration of the fitted diffusive regime, or the variance of D across runs. Given the rare-event character of the dynamics (interlayer hops are ~1% of events, and only six such hops occur in a 3 ns trajectory at 700 K), Poisson counting errors are non-negligible. Reporting event counts and error bars is necessary to support the stated orders-of-magnitude comparison and the fitted activation energy.","section":"Sec. III.A, Fig. 4 and Methods"},{"comment":"The slab model keeps the two deepest layers fixed, while small-polaron formation requires local lattice distortion. The statement that 'transport to deeper layers was never observed' is therefore partly a consequence of the model constraint: a polaron cannot stabilize in a fixed layer because the surrounding atoms cannot relax. This could bias the surface/bulk comparison and the conclusion that surface confinement controls transport. I ask the authors to test with a thicker relaxed slab, or at least to discuss explicitly how the fixed bottom layers affect the confinement claim and the inferred mobility suppression.","section":"Sec. II.A and Sec. III.A"}],"minor_comments":[{"comment":"Typo: 'hooping trajectories' should be 'hopping trajectories'.","section":"Sec. III.A"},{"comment":"The activation barrier of 386 meV is described as within 'approximately 15%' of infrared estimates of 300–330 meV; the actual discrepancy is 17–29%. Please recalculate or reword.","section":"Conclusions / Sec. IV"},{"comment":"Grammar: 'straightforward its use' should be 'streamline its use'.","section":"Supplemental Material, Sec. I"},{"comment":"The statement that data and code 'will be made available upon publication' does not allow reviewers to reproduce the results. Consider depositing the LEOPOLD version, configuration files, and representative trajectories in a public repository at the revision stage.","section":"Data Availability"}],"recommendation":"major_revision","confidential_remarks":"The qualitative finding — that surface polaron mobility is strongly suppressed and dominated by subsurface [001] hops — appears likely to be correct and is a useful extension of the LEOPOLD methodology. The obstacle to acceptance is the headline quantitative claim: the room-temperature mobility and the 'several orders of magnitude' suppression at 300 K rest on a four-point Arrhenius extrapolation from a temperature regime where the polaron may already be trap-influenced, with no uncertainty estimates. The fixed-layer issue compounds the concern. I would be willing to accept a revision that either provides error bars and a low-temperature validation (or explicit trap correction) or rephrases the central claim to focus on the directly simulated 400–700 K regime. The self-reliance on Ref. [49] is acceptable given that this is a direct extension, but a brief discussion of shared systematic errors in the bulk/surface comparison would strengthen the paper."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper extends the LEOPOLD machine-learning scheme to polaron dynamics at a reduced rutile TiO2(110) surface and is the first to get nanosecond-time-scale surface polaron trajectories. The main qualitative result—that surface mobility is suppressed by several orders of magnitude relative to bulk at all simulated temperatures—looks solid. It is directly supported by multiple independent runs at 400–700 K, matches the known preference for [001] subsurface inter-row hopping, and is consistent with earlier picosecond FPMD. The active-learning training strategy and the direct prediction of on-site magnetizations are sensible technical improvements, and the validation errors are reasonable for a surface with a polaron. The activation barrier from the Arrhenius fit (386 meV) also lands close to the infrared-spectroscopy estimate. The authors deserve credit for making this work and for a clearly written Letter.\n\nThe soft spot is the quantitative room-temperature claim. No hop is seen in ~40 ns at 300 K, so the paper extrapolates an Arrhenius fit over 400–700 K to obtain µ = 1.9×10^-5 cm2/Vs and calls it very good agreement with porous-TiO2 experiments. That extrapolation is load-bearing for the abstract's central comparison, and there are three reasons to be cautious. First, four temperatures with no uncertainty estimates is a thin fit, and the paper does not report residuals or curvature checks. Second, Sec. III.B shows that at 400 K the polaron–vacancy distance distribution is strongly peaked near the vacancy when initialized close to it, which suggests the lowest fitted temperature may already mix trapping/drift with true diffusion. Third, the bulk comparison is pristine bulk versus a surface that contains an oxygen vacancy, so part of the suppression could be vacancy trapping rather than surface topology. The qualitative suppression likely survives these concerns, but the specific room-temperature value and the claimed 'orders of magnitude' relative to bulk at 300 K should be read with a grain of salt.\n\nThis is a paper worth engaging with, and a serious referee should see it. The authors should be asked to provide error bars, a bulk-with-vacancy control, and either longer 300 K runs or a less assumption-dependent way to estimate the low-temperature mobility. If those are addressed, the quantitative claim would become much stronger.\n\nFor a reading group, it is a good example of extending ML potentials to charge dynamics, and I would probably cite it if I worked on polarons in oxides.","headline":"Strong qualitative result on surface polaron suppression, but the room-temperature mobility is an extrapolation from zero observed hops and should be treated as provisional.","tokens_in":15946,"tokens_out":1670,"would_cite":true,"duration_ms":20840,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper argues that small polarons—electrons trapped by lattice distortion—hop several orders of magnitude more slowly on the reduced rutile TiO2(110) surface than in bulk rutile, and that this surface bottleneck, not trapping defects al","keywords":["small polaron","TiO2(110)","surface","mobility","machine learning interatomic potential","oxygen vacancy","hopping transport","porous TiO2"],"falsifier":"Run the machine-learned polaron dynamics at 300 K for tens of nanoseconds and count hops: if the observed hopping rate is statistically incompatible with the Arrhenius-extrapolated rate (equivalently, mobility at 300 K differs from about 2×10^-5 cm²/Vs by more than an order of magnitude), the central claim fails. Alternatively, a terahertz or time-resolved conductivity measurement on porous rutile films with controlled particle size that yields a mobility much higher than 10^-4 cm²/Vs and independent of surface fraction would contradict the surface-suppression picture.","tokens_in":15026,"feed_emoji":"⚡","tokens_out":4181,"duration_ms":46042,"temperature":0.7,"pith_summary":"The paper tries to establish that small polarons—electrons localized at titanium sites and coupled to lattice distortion—move far more slowly on the oxygen-deficient rutile TiO2(110) surface than in bulk rutile, by roughly four orders of magnitude, and that this explains why porous TiO2 conducts electrons much more poorly than single crystals. It reaches this conclusion by using a machine-learned interatomic potential, trained on first-principles molecular dynamics, to simulate several nanoseconds of polaron hopping at 400–700 K, then extrapolating to room temperature. The underlying cause is geometric: surface hopping is confined mostly to flat rows in the subsurface layer, while jumps between layers are rare. A reader should care because the result gives a microscopic, transferable picture of how surface structure, not just defects, controls charge transport in redox-active oxides used in catalysis, photovoltaics, and energy storage.","feed_headline":"Polaron hops on TiO2 surface run four orders slower than bulk","feed_subtitle":"Nanosecond machine-learned simulations trace the slowdown to missing interlayer pathways and match porous-electrode data.","key_machinery":"The carrying mechanism is a polaron-aware machine-learned interatomic potential: a graph neural network that augments each atom's chemical identity with a one-hot charge-state flag identifying which Ti site currently hosts the polaron, and that predicts both energies/forces and site-resolved spin magnetizations. At each molecular-dynamics step the polaron is reassigned to the Ti atom with the largest predicted magnetization, so lattice motion and charge position evolve together. This enables nanosecond-scale trajectories with explicit hopping events, from which diffusion coefficients and mobilities are obtained via mean-square displacement fits and the Einstein relation.","core_discovery":"The central claim is that small-polaron mobility on the reduced rutile TiO2(110) surface is suppressed by several orders of magnitude relative to bulk rutile across the investigated temperature range, with an Arrhenius extrapolation giving approximately 1.9×10^-5 cm² V^-1 s^-1 at 300 K, matching the experimental range for porous TiO2 (7×10^-6 to 5×10^-4). The suppression is attributed to loss of favorable hopping pathways: about 87% of hops occur along [001] Ti rows in the subsurface S-1 layer, about 12% along surface rows, and only about 1% are interlayer hops; inter-row and deeper-layer hops were not observed. The oxygen vacancy acts as an attractive center that biases polaron positions to","pith_inferences":["If the Arrhenius extrapolation holds, then at device operating temperatures surface polaron transport must rely on a sparse network of subsurface-row hops, so experimental probes that distinguish subsurface from surface charge carriers could directly test this picture.","The results imply a design rule: nanostructuring that increases exposure of (110)-like facets will suppress carrier mobility, so porous electrodes might be optimized by selecting facets that preserve bulk-like hopping rows or by introducing dopants that deepen favorable subsurface pathways.","Because the model simulates a single polaron with one vacancy, the strong polaron–vacancy attraction suggests that at realistic carrier densities multi-polaron interactions and vacancy clustering will modify mobilities; extending the approach to two or more excess charges is a natural next test.","The paper's claim predicts that measured electron mobility in porous TiO2 should depend systematically on particle size and surface-to-volume ratio; a controlled experimental series varying those parameters could quantitatively validate the surface-confined transport scenario."],"forward_implications":["Porous TiO2's low electron mobility is explained by surface-restricted hopping topology, not solely by trapping defects, so bulk single-crystal mobility values should not be used directly in models of nanostructured or porous electrodes.","The machine-learning strategy transfers to other reducible oxide surfaces and defect/adsorbate environments, opening a route to predictive simulations of excess-charge dynamics in catalytic and energy-conversion materials.","Because the oxygen vacancy reshapes the polaron free-energy landscape, vacancy concentration and spatial distribution directly control surface transport; engineering vacancy placement could tune conductivity.","The computed activation barrier (~386 meV) provides a parameter for device-level transport models of TiO2-based photoelectrodes and sensors.","The strong anisotropy of hopping—dominant subsurface row hops, rare interlayer hops—implies that surface transport is effectively quasi-one-dimensional along [001] rows in the second layer."],"fun_headline_variants":["Surface polarons on TiO2 hop 4 orders slower than bulk","Rare interlayer hops stall surface polarons on TiO2","87% of polaron hops stay in subsurface layer: mobility drops","Oxygen vacancy traps polarons in planar paths on TiO2","ML reveals polarons avoid interlayer hops on rutile surface"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The room-temperature mobility and the 'orders of magnitude slower' claim rest on the assumption that polaron hopping follows the same Arrhenius law measured at 400–700 K all the way down to 300 K, where no hop was observed in about 40 ns of simulation; a regime change or increased trapping below 400 K would invalidate the extrapolation.","fun_headline_variants_meta":{"raw":{"variants":["Surface polarons on TiO2 hop 4 orders slower than bulk","Rare interlayer hops stall surface polarons on TiO2","87% of polaron hops stay in subsurface layer: mobility drops","Oxygen vacancy traps polarons in planar paths on TiO2","ML reveals polarons avoid interlayer hops on rutile surface"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00013,"raw_usage":{"total_tokens":989,"prompt_tokens":798,"completion_tokens":191,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":542,"completion_tokens_details":{"reasoning_tokens":103}},"tokens_in":542,"tokens_out":191,"duration_ms":2731,"temperature":1.0,"reasoning_tokens":103,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T11:34:09.108346+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the machine-learned polaron dynamics at 300 K for tens of nanoseconds and count hops: if the observed hopping rate is statistically incompatible with the Arrhenius-extrapolated rate (equivalently, mobility at 300 K differs from about 2×10^-5 cm²/Vs by more than an order of magnitude), the central claim fails. Alternatively, a terahertz or time-resolved conductivity measurement on porous rutile films with controlled particle size that yields a mobility much higher than 10^-4 cm²/Vs and independent of surface fraction would contradict the surface-suppression picture.","supporting_citations":[],"review_version":1}