{"id":"81cac792-ae49-41f5-aa8e-3d3343112f54","arxiv_id":"2505.19999","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Using a dynamic disorder tight-binding model, the authors trace carrier movement through orbital channels and identify the ppπ channel as a transport bottleneck in MAPbBr3 and MAPbI3.","lead":"This paper simulates how electrons and holes move inside two halide perovskite crystals, MAPbBr3 and MAPbI3, by tracking which atomic orbitals they temporarily occupy. The authors identify the ppπ bonding channel as a bottleneck that slows charge transport, and they show that energy gaps between on-site levels and the strength of spin-orbit coupling control how fast charges spread.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The ppπ bottleneck claim is inferred from a small mean hopping and slow spreading of X-pπ-initialized wavepackets, but the full hole wavefunction is never tested against perturbation of t_ppπ itself; a bottleneck in a multi-channel network must be demonstrated by perturbing the channel, not by its…","rationale":"The reader's weakest_assumption focuses on the fidelity of the MD trajectories and the TB mapping rule. That is a legitimate input-validity concern, but the single most load-bearing issue for the paper's headline claim is internal: the paper identifies a causal bottleneck without performing the control experiment that would establish causation. The bottleneck claim is central to the abstract and conclusions, and it is supported only by (a) a parameter comparison and (b) single-orbital-initiation MSDs. Both are necessary but not sufficient. A bottleneck in a network of parallel transport channels must be tested by removing or rescaling that channel and observing the effect on the full wavefunction dynamics. The paper's own methodology already includes exactly this kind of parameter perturbation for on-site energies and halide SOC, so the missing test is natural and inexpensive. If t_ppπ perturbation leaves the full hole MSD essentially unchanged, the central claim would need to be substantially weakened; if it changes the MSD strongly, the claim would be robustly supported. This does not change the overall CONDITIONAL verdict, but it sharpens the condition: the authors should either provide this test or temper the bottleneck language.","tokens_in":13085,"tokens_out":8161,"duration_ms":101179,"concrete_test":"Perform the full h+ wavefunction MSD calculation of Fig. 5 (right) with t_ppπ set to 0, ×0.5, and ×2 in the TB Hamiltonian, leaving all other parameters fixed (same protocol as the ε_Pb^p and γ_X perturbations). Compute the diffusive MSD slope (e.g., 50–100 fs) for each case. A monotonic, sizable (>~10%) change would confirm the bottleneck; negligible change would falsify the 'critical bottleneck' claim. Running the same test for the e− wavefunction (where X-pπ occupancy is small) would further verify that the effect is channel-specific.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central causal claim — 'the ppπ channel is a critical bottleneck for charge transport' — is never directly tested. The evidence offered is (1) the averaged TB parameter \\bar{t}_{ppπ} is smaller than \\bar{t}_{spσ} and \\bar{t}_{ppσ} (Table I), and (2) wave packets initialized on X-pπ orbitals show suppressed MSD (Fig. 2). Neither step establishes that the ppπ channel limits transport of a realistic full hole wave packet. The full h+ state contains a mixture of all three channels (Fig. 7, lower panels), and a slow parallel pathway need not throttle a network with two faster pathways: carriers can reroute through spσ/ppσ. The paper therefore infers the bottleneck rather than demonstrating it. The authors already use parameter-perturbation tests for on-site energies and γ_X in Fig. 5; the decisive missing control is the analogous perturbation of t_ppπ itself. Without it, the attribution of the MAPbI3-vs-MAPbBr3 hole mobility difference specifically to halide SOC, and the general 'critical bottleneck' conclusion, remain correlative rather than causal.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript applies the authors' dynamic disorder (DD) model to MAPbBr3 and MAPbI3, using tight-binding Hamiltonians parameterized from hybrid-DFT (HSE+SOC) calculations and Wannier90, with nuclear trajectories taken from force-field molecular dynamics. The authors compute mean-squared displacements (MSD) of carrier wave packets and time-resolved orbital occupations for both orbital-localized and full conduction/valence-band initializations. Their central claims are that (i) energy gaps between on-site levels and the halide spin-orbit coupling strength govern which real-space transport channels (spσ, ppσ, ppπ) become populated over time, and (ii) the ppπ channel is a critical bottleneck for charge transport, particularly for holes in MAPbI3 where halide SOC is stronger. The paper further proposes that the difference in hole mobility between MAPbBr3 and MAPbI3 is partly attributable to the larger halide SOC in the iodide.","tokens_in":13269,"tokens_out":4393,"duration_ms":49445,"significance":"If the bottleneck and channel-occupation claims are correct, the work offers a useful real-space decomposition of carrier transport in halide perovskites and could inform composition-based design rules (e.g., why MASnI3 is more conductive). The first-principles parameterization of the TB Hamiltonian (HSE, not PBE) is a step forward relative to earlier DD-model studies, and the orbital-occupation analysis is a genuinely informative diagnostic. However, the paper's central causal assertions about the ppπ bottleneck and about SOC controlling hole transport are not directly tested: the evidence is correlative (small mean hopping, suppressed MSD for a ppπ-localized initial state, and occupation changes under arbitrary parameter shifts). The missing control is a perturbation of t_ppπ itself in the full hole dynamics. Given that the authors already perform parameter-perturbation scans for on-site energies and γ_X, adding such a test is feasible and would substantially strengthen the paper. The manuscript is therefore of interest but requires additional work before the bottleneck conclusion is established.","major_comments":[{"comment":"The MSD curves in Fig. 2 (and Fig. 5) are presented without error bars or any statistical measure. The methods section states that mobilities are averaged over different starting positions only later in Sec. II.C, but it is not clear whether the MSD curves in Fig. 2 are averaged over starting positions, over multiple trajectories, or over the 100 snapshots. Given that the curves for MAPbBr3 and MAPbI3 are close (e.g., Pb-s vs Pb-p in Fig. 2), the reader cannot assess whether the claimed differences are statistically meaningful. The authors should report standard errors or confidence intervals, and should state explicitly how many initial positions and trajectory realizations were used for each curve.","section":"§III.A and Fig. 2"}],"minor_comments":[{"comment":"The phrase 'proceeding the initialization' should likely be 'preceding the initialization' or 'following the initialization'; the intended temporal order is unclear.","section":"Fig. 3 caption"},{"comment":"Reference [33] contains a typo: 'Metal-Insulator Tranisitions' should be 'Metal-Insulator Transitions'.","section":"References"},{"comment":"The text says 'when we apply shifts that align the Pb-p on-site levels for the two HaPs more (cf. Tab. I)'. Since the shifts (−1 eV for Br, +1 eV for I) move the Pb-p levels in opposite directions and make the absolute on-site energies less similar (3.23→2.23 eV, 2.58→3.58 eV), 'align' should be clarified: the intended effect is to make the Pb-p to X-p gaps more similar (5.29→4.29 eV for Br, 4.15→5.15 eV for I), not the Pb-p levels themselves.","section":"Sec. III.B"},{"comment":"The relation between on-site energies and PDOS is stated without derivation or reference; a short justification or citation would help readers who are not familiar with the Wannier90 projection scheme.","section":"Eq. (1)"},{"comment":"The sentence 'the extracted hopping elements, on-site energies, and SOC parameters are combined into a first-principles database' could be clearer about whether the hopping mapping is a single exponential fit per bond type or includes a dependence on local environment beyond the Pb-X distance.","section":"Sec. II.A"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of the journal and reports an interesting extension of the dynamic-disorder model with a real-space channel-resolved analysis. The central bottleneck claim is plausible but under-tested; adding a direct perturbation of t_ppπ in the full hole dynamics would move the paper from correlative to causal. The parameter scans are also somewhat arbitrary in magnitude, but that is fixable with a sensitivity analysis. I do not see grounds for rejection, as the methodological framework and the qualitative observations (small ppπ hopping, suppressed ppπ-initialized MSD, occupation shifts under on-site/SOC changes) are internally consistent and the main missing element is an additional control calculation. I would encourage the authors to also report error bars and to clarify the averaging procedure."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper's real contribution is a real-space channel decomposition: separating carrier dynamics into spσ, ppσ, ppπ pathways and tracking orbital occupations over time, with the TB parameters now from HSE+SOC rather than PBE. That is a genuine step beyond the earlier dynamic-disorder papers and yields a concrete, testable picture of why holes lag electrons and why MAPbI3 and MAPbBr3 differ. The mapping from DFT to TB is not fitted to mobilities, so the circularity concern is low — credit where due.\n\nWhere I push back is on the bottleneck claim. The paper identifies ppπ because its average hopping is smallest and because wavepackets initialized on X-pπ spread slowly. That is suggestive but not demonstrative that ppπ throttles a realistic full hole state, which contains a mixture of all three channels. The authors already run perturbation scans for on-site levels and halide SOC; the missing control is scaling t_ppπ itself. Without that, the 'critical bottleneck' language overstates what is actually shown. The difference in hole mobility between the two materials is attributed to halide SOC, but that attribution is correlative given the parameter set; the same scan logic would need to isolate γX.\n\nOther soft spots are minor but real: 100 fs MSD curves with no error bars or extracted diffusion coefficients, arbitrary ±1 eV and 2× shifts, and a placeholder for the supplemental data. None of these sinks the qualitative mechanism, but they cap how much quantitative weight the paper can carry right now. The reliance on ref. [24] trajectories is a reasonable inheritance but worth stating as a limitation.\n\nIf I were the editor, I'd send this to a referee. The channel decomposition and the candidate bottleneck are interesting enough to warrant scrutiny, and the missing perturbation control is exactly the kind of thing a good referee would catch. The paper is honest about its own scaffolding, and the authors have a track record of shipping reproducible parameterizations. A revision that adds the t_ppπ scan, error bars on MSD, and actual diffusion coefficients would make the central claim much stronger.","headline":"Real-space channel decomposition is a genuine step forward, but the 'ppπ bottleneck' is a plausible inference, not a tested causal claim.","tokens_in":13856,"tokens_out":1864,"would_cite":true,"duration_ms":21312,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"In halide perovskites, the smallest hopping channel—the ppπ bond—acts as the bottleneck that suppresses hole transport, and its filling is set by on-site energy gaps and halide spin-orbit coupling.","keywords":["dynamic disorder model","halide perovskites","charge transport","spin-orbit coupling","tight-binding","orbital occupations","ppπ channel","MAPbBr3"],"falsifier":"A direct test would be to recompute the carrier dynamics using molecular dynamics trajectories generated from first principles (for example, HSE-quality ab initio MD) rather than a force field, and check whether the ppπ orbital occupation and the associated hole-mobility suppression are reproduced. Alternatively, an experiment that systematically varies halide spin-orbit coupling—such as a pressurized or alloyed MAPb(Br,I)3 series—while measuring hole mobility could falsify the claim if larger halide SOC does not produce the predicted drop in hole transport.","tokens_in":12838,"feed_emoji":"⚛️","tokens_out":3936,"duration_ms":44726,"temperature":0.7,"pith_summary":"This paper aims to establish which microscopic features of halide perovskites control how fast electrons and holes move, using a dynamic disorder model that couples molecular-dynamics geometries to a time-dependent tight-binding Hamiltonian. The authors find that charge carriers travel through three real-space bonding channels—spσ, ppσ, and ppπ—and that the ppπ channel, which has a markedly weaker hopping parameter, becomes a critical bottleneck. They further show that the energy gaps between on-site orbital levels and the magnitude of halide spin-orbit coupling determine how much each channel is populated over time. If this picture is correct, it explains why MAPbBr3 and MAPbI3 differ in electron versus hole transport, and it offers concrete levers—level alignment and halide composition—for engineering carrier mobility.","feed_headline":"One weak orbital channel throttles hole transport in halide perovskites","feed_subtitle":"On-site energy gaps and halide spin-orbit coupling decide which transport channel fills—and the weakest one traps holes.","key_machinery":"The central object is the dynamic disorder (DD) model: a time-dependent tight-binding Hamiltonian whose hopping elements are mapped from first-principles DFT (HSE functional including SOC) data via an exponential fit to Pb–X distances, with on-site energies from Ewald summations and SOC parameters fixed to their mean values. The transport channels are defined by the three hopping types—tspσ, tppσ, and tppπ—and the analysis tool is the time evolution of orbital occupations, which reveals which channel carries the current and which one traps carriers.","core_discovery":"The central claim is that charge carrier transport in MAPbBr3 and MAPbI3 is governed by the time-dependent filling of three real-space transport channels—spσ, ppσ, and ppπ—where the ppπ channel acts as a bottleneck. The paper demonstrates that the relative occupations of these channels are set by two material-specific parameters: the energy gaps between on-site levels (particularly Pb-p and X-p) and the halide spin-orbit coupling strength γX. A smaller on-site gap enhances transport by allowing more balanced occupation and steadier flow along the ppσ channel, whereas a larger halide SOC increases X-pπ occupation, drawing carriers into the weak ppπ channel and suppressing hole mobility. This mechanism rationalizes why electron transport is stronger in MAPbI3 than in MAPbBr3 while hole transport shows the opposite trend, and it attributes the lower hole mobility of the iodide in part to its roughly twofold larger halide SOC.","pith_inferences":["If the mechanism generalizes, any dynamically disordered semiconductor with a weak π-type hopping channel could show the same bottleneck, meaning the design rule 'minimize the energy gap between channel orbitals' may extend beyond halide perovskites to other heavy-halide or heavy-metal systems.","The Rabi-like proportionality between transition probability and on-site level difference implies that tuning on-site gaps by even a few tenths of an electronvolt could measurably change mobility—this is a testable prediction for alloy series such as MAPb(BrxI1−x)3.","Experimental probes that vary the halide SOC independently of lattice structure, such as high-pressure studies or epitaxial strain, could isolate the SOC contribution and directly test whether the ppπ occupation follows the predicted trend."],"forward_implications":["The ppπ bottleneck implies that hole transport in halide perovskites is intrinsically limited by the occupation of X-pπ orbitals, so reducing halide spin-orbit coupling—for example, by substituting a lighter halide—should raise hole mobility.","Because MAPbI3 exhibits stronger electron transport but weaker hole transport than MAPbBr3, the results provide a microscopic rationale for composition-dependent mobility trends that experiments already observe.","Aligning the Pb-p and X-p on-site energy levels, whether by chemical substitution or strain, should enhance carrier transport by making the ppσ channel more efficient.","The finding that halide SOC acts by repopulating the bottleneck orbital suggests that spin-orbit effects can be understood as a channel-filling mechanism rather than merely a band-structure correction."],"supporting_citations":[{"why":"Supplies the force-field molecular dynamics trajectories for the large supercells and the 300 K reference data used to construct the time-dependent Hamiltonians.","marker":"[24]"},{"why":"Introduces the dynamic disorder model that this paper extends and provides the original framework for computing mobilities from wave-packet evolution.","marker":"[22]"},{"why":"Establishes the transient-localization picture of carrier transport in halide perovskites, which motivates the use of a time-dependent tight-binding description.","marker":"[14]"},{"why":"Provides the experimental carrier-mobility trends for MAPbBr3 and MAPbI3 that the paper compares its computed differences against.","marker":"[60]"},{"why":"Reports experimental evidence supporting dynamic disorder as the dominant scattering mechanism in these materials, underpinning the model's relevance.","marker":"[46]"},{"why":"The Wannier90 code is used to extract the tight-binding parameters (on-site energies, hoppings, SOC) that enter the mapping rule.","marker":"[53]"},{"why":"Provides the connection between crystal orbital bond indices and hopping elements, used to interpret the hopping-parameter variations.","marker":"[54]"}],"fun_headline_variants":["ppπ channel is the bottleneck for halide perovskite transport","Spin-orbit coupling steers carrier flow in halide perovskites","Weak ppπ channel limits hole mobility in perovskites","On-site gaps and SOC dictate transport channels in perovskites"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire analysis rests on the assumption that the force-field molecular dynamics trajectories taken from the earlier work faithfully reproduce the relevant anharmonic nuclear motion of MAPbBr3 and MAPbI3, and that the tight-binding mapping rule—exponentially fitted hoppings, Ewald-averaged on-site energies, and mean SOC parameters—yields a trustworthy time-dependent Hamiltonian.","fun_headline_variants_meta":{"raw":{"variants":["ppπ channel is the bottleneck for halide perovskite transport","Spin-orbit coupling steers carrier flow in halide perovskites","Weak ppπ channel limits hole mobility in perovskites","On-site gaps and SOC dictate transport channels in perovskites"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000604,"raw_usage":{"total_tokens":2795,"prompt_tokens":897,"completion_tokens":1898,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":513,"completion_tokens_details":{"reasoning_tokens":1830}},"tokens_in":513,"tokens_out":1898,"duration_ms":12409,"temperature":1.0,"reasoning_tokens":1830,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T14:02:00.073541+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct test would be to recompute the carrier dynamics using molecular dynamics trajectories generated from first principles (for example, HSE-quality ab initio MD) rather than a force field, and check whether the ppπ orbital occupation and the associated hole-mobility suppression are reproduced. Alternatively, an experiment that systematically varies halide spin-orbit coupling—such as a pressurized or alloyed MAPb(Br,I)3 series—while measuring hole mobility could falsify the claim if larger halide SOC does not produce the predicted drop in hole transport.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the force-field molecular dynamics trajectories for the large supercells and the 300 K reference data used to construct the time-dependent Hamiltonians."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the dynamic disorder model that this paper extends and provides the original framework for computing mobilities from wave-packet evolution."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the transient-localization picture of carrier transport in halide perovskites, which motivates the use of a time-dependent tight-binding description."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the experimental carrier-mobility trends for MAPbBr3 and MAPbI3 that the paper compares its computed differences against."},{"cited_title":"D¨ orflinger, P","cited_arxiv_id":null,"evidence_quote":"Reports experimental evidence supporting dynamic disorder as the dominant scattering mechanism in these materials, underpinning the model's relevance."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the connection between crystal orbital bond indices and hopping elements, used to interpret the hopping-parameter variations."}],"review_version":1}