{"id":"2ad06397-6f8c-43bf-8b10-b54cf7442084","arxiv_id":"2509.02682","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"First-principles many-body calculations predict strongly anisotropic exciton diffusion in phosphorene, with long-range exchange interactions boosting transport along the armchair direction at low temperatures.","lead":"This paper calculates how excitons, bound electron-hole pairs, move through phosphorene, a two-dimensional form of black phosphorus. It predicts strongly directional energy flow, with faster motion along the armchair direction, and says the effect is especially strong at low temperatures.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The low-temperature exchange-induced boost in D_AC is likely controlled by the unspecified Lorentzian width in Eq. (4), because the thermal population sits at the Q=0 cusp where phonon dephasing vanishes; the paper reports no sensitivity to this width.","rationale":"The reader identified the long-range exchange model as the weak point of the low-temperature enhancement, and mentioned the unspecified Lorentzian width only in passing. My stress-test focuses on a different but closely related breach: even taking the exchange model at face value, the low-T diffusion constant is not a phonon-limited quantity because the thermal population accumulates at the Q=0 cusp of the linear dispersion, where the phonon dephasing rate vanishes unless a residual broadening or an intervalley channel is invoked. This makes the central low-temperature claim parameter-controlled rather than physically robust. The paper is otherwise internally consistent: the room-temperature anisotropy is supported by the effective-mass hierarchy and a plausible exchange model, and the comparison to experiment, while imperfect, is not load-bearing. I therefore keep the reader's CONDITIONAL verdict rather than moving to reject; the condition should be a concrete sensitivity analysis of the broadening and of the intervalley scattering contribution at low temperature.","tokens_in":13163,"tokens_out":8280,"duration_ms":88324,"concrete_test":"Recompute D_AC at 50 K (and, if feasible, at 10 K) using Eq. (9), varying the Lorentzian full width at half maximum in Eq. (4) over at least 0.1, 1, and 5 meV, and separately with all intervalley ΓV scattering channels disabled. Also report the center-of-mass Q-grid spacing near Q=0. If D_AC changes by more than roughly a factor of two across these choices, or increases without bound as the Q grid is refined, then the low-temperature exchange enhancement is not a well-defined prediction of the model.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The emphasized result, that long-range exchange sharply enhances D_AC at low temperature, depends on evaluating the diffusion constant Eq. (9) with a dephasing rate Γ^η_Q obtained from Eq. (5). At low T, the Boltzmann factor concentrates the diffusion sum at Q≈0. With exchange, Eq. (7) is linear in Q along AC, so the group velocity is finite at Q=0; the numerator in Eq. (9) is maximal exactly where the phonon contribution to Γ is smallest. In the ΓΓ-only transport treatment, the lowest bright 1s state at Q=0 cannot absorb a phonon as T→0, and phonon emission from the band minimum requires an intervalley channel that is neither specified nor tested. Consequently, D_AC(T→0) is set by the residual width of the Lorentzian that replaces the δ-function in Eq. (4), a width the paper never states. If that width is numerical regularization, the low-temperature 'soaring' of D_AC in Fig. 4(a), and the claim that exchange significantly enhances AC diffusion at low T, are not quantitatively established. The 300 K values (D_AC=7.40, D_ZZ=0.48 cm2/s) are less affected because phonon occupancy dominates the dephasing, so the room-temperature anisotropy may survive; but the paper's emphasized low-temperature regime is unconstrained without a specified Γ(Q≈0) or a sensitivity analysis.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript develops a first-principles-informed many-body description of exciton optics, relaxation, and transport in monolayer phosphorene. The electronic band structure and interband dipoles come from HSE06 DFT, phonons and electron-phonon matrix elements from PBE DFT, and exciton states from the Wannier equation with a Keldysh-screened Coulomb potential. Exciton-phonon scattering rates are then used to compute optical absorption, relaxation cascades, and the anisotropic diffusion constants D_AC and D_ZZ. The central claim is that the long-range exchange interaction strongly enhances the already large anisotropy of exciton diffusion, particularly at low temperatures, where D_AC is reported to 'soar'; at 300 K the reported values are D_AC = 7.40 cm2/s and D_ZZ = 0.48 cm2/s.","tokens_in":13403,"tokens_out":9695,"duration_ms":91315,"significance":"If correct, the paper provides a concrete, largely parameter-free microscopic prediction for directional exciton energy transport in an anisotropic two-dimensional semiconductor, with practical relevance for exciton-based circuitry and light harvesting. The strengths are the ab initio parametrization with no transport parameters fitted to the target diffusion data, the good match of the 1s exciton binding energy to the experimental value (0.29 eV versus about 0.3 eV), and the explicit treatment of exciton-phonon scattering beyond effective-mass arguments. The low-temperature anisotropy enhancement is a falsifiable prediction that would make the paper influential, which is why the numerical robustness of that regime needs to be demonstrated.","major_comments":[{"comment":"The Lorentzian width used to regularize the energy-conserving delta function in Eq. (4) is never stated, and no sensitivity analysis is given. This is load-bearing for the low-temperature claim: in Eq. (9) the diffusion sum is weighted by 1/Gamma_Q and by the Boltzmann factor, so as T goes to zero the population is concentrated at Q approximately equal to 0. With exchange, the AC dispersion is linear in Q at small Q (the pink dashed curve in Fig. 3(a)), so the group velocity is finite at Q equals 0, whereas Gamma_Q at Q equals 0 is smallest; if no phonon-emission channel exists from the lowest bright Gamma-Gamma 1s state, then Gamma(0) is controlled by the Lorentzian tail and D_AC(T -> 0) diverges as the width goes to zero. The paper says only that a Lorentzian is used and justifies this by higher-order phonon processes, but gives no value and no convergence test. Please provide the width and show that D_AC in Fig. 4(a) is insensitive to it, or otherwise establish that the low-temperature 'soaring' is physical and not regularization-controlled.","section":"Sec. III.C.3, Eqs. (4), (5), (9)"},{"comment":"The low-temperature enhancement of D_AC is generated entirely by the linear small-Q term in the long-range exchange K_Q, yet the implementation of Eq. (8) rests on three assumptions that are stated but not tested: long-range exchange is included only for Gamma-Gamma excitons, short-range exchange is neglected, and the dipole factor is taken as Q dot d with d along AC. Agreement of the 1s binding energy with experiment validates the Wannier potential, not these exchange-model choices. Because the headline result is the exchange-induced low-temperature anisotropy, I ask for a sensitivity test: for example, include a screened short-range exchange term or a finite-Q correction to the dipole matrix element and recompute D_AC(T). If the result is robust, this should be stated; if not, the central claim needs to be qualified.","section":"Sec. III.C.1, Eq. (8)"},{"comment":"The exciton-phonon scattering rates are computed using only a selected subset of phonon modes, with a constant deformation potential for optical modes, g proportional to sqrt(q) for acoustic modes, and linear interpolation between high-symmetry directions. These choices affect the absolute values of Gamma_Q and hence both the linewidths and the diffusion constants. The good agreement of the 300 K average with experiment is encouraging, but the arithmetic mean of the reported D_AC = 7.40 and D_ZZ = 0.48 cm2/s is 3.94 cm2/s, and the experimental value of 5.0 cm2/s is for a 2 nm amorphous sample rather than a crystalline monolayer, so the comparison is not a sharp validation. Please provide an uncertainty estimate from the omitted modes and the interpolation scheme, or justify that the selected modes dominate the transport-relevant scattering.","section":"Sec. II.A, Tables I-II, Sec. III.C.3"}],"minor_comments":[{"comment":"The text refers to the 'Elliot formula'; this should be the Elliott formula, and 'interband diople' should read 'interband dipole'.","section":"Sec. III.B"},{"comment":"The notation for direction-resolved coupling coefficients, e.g. [g_x/sqrt(q), g_xy/sqrt(q), g_y/sqrt(q)], and the pairing of a coupling coefficient with a phonon energy in a single cell should be defined explicitly in the caption; as printed, the two numbers per cell with mixed units are easy to misread.","section":"Table I"},{"comment":"The comparison with experiment should be stated more carefully: the average of the reported 300 K diffusion constants is 3.94 cm2/s, not 5.0 cm2/s, and the experimental system is amorphous few-layer rather than crystalline monolayer phosphorene.","section":"Sec. III.C.3"},{"comment":"The caption uses 'Diff. Constant' where 'Diffusion constant' would be clearer, and the order in which panels (a)-(f) are described in the text should match the layout in the figure.","section":"Fig. 4"}],"recommendation":"major_revision","confidential_remarks":"I do not see a circularity problem: essentially all material parameters come from first-principles calculations. The main risk is the unstated Lorentzian width and the sensitivity of the low-temperature D_AC result to the exchange model. Both are fixable in revision, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The thing to know: this is a competent, mostly first-principles calculation of exciton diffusion in phosphorene, and the central qualitative claim—that the armchair direction carries excitons far faster than zigzag and that long-range exchange sharpens that anisotropy at low temperature—is plausible. The paper earns credit for parametrizing everything from DFT and the Wannier equation, with no transport data fitted. The 1s binding energy matching experiment is a good sanity check, and the 300 K values (D_AC = 7.40, D_ZZ = 0.48 cm2/s) are concrete, falsifiable predictions.\n\nThe main soft spot is the low-temperature regime, and I think the stress-test note lands here. The diffusion constant in Eq. (9) divides by the dephasing rate Gamma(Q). At T→0, the population sits at Q≈0, where phonon absorption is frozen out and emission is forbidden if the exchange-induced linear dispersion makes Q=0 the minimum. That leaves the Lorentzian width in Eq. (4) as the effective regulator, and the paper never states that width or tests sensitivity to it. The “soaring” D_AC in Fig. 4(a) is therefore not quantitatively established. The room-temperature result is less exposed, since phonon occupancy dominates the dephasing there, so the anisotropy at 300 K may survive even if the low-temperature boost is fragile.\n\nOther soft spots are proportionate and mostly quantitative: the electron-phonon coupling uses a selected subset of phonon modes, a constant deformation-potential approximation for optical modes, and linear interpolation in q-space; the dielectric constant is averaged over in-plane directions. These approximations are reasonable for a first pass but mean the absolute diffusion constants should be taken with some slack. The comparison to the amorphous-sample experiment is acknowledged as approximate, so I don't treat it as a claim of quantitative agreement.\n\nWho should read this: people working on exciton transport in anisotropic 2D materials, especially phosphorene and black phosphorus, and theorists doing many-body parametrization of TMD-like systems. It advances the subfield without reorganizing it.\n\nMy recommendation: this deserves a serious referee rather than a desk reject. But the referees should be told to focus on the low-T sensitivity and to require the authors to report the Lorentzian width, justify it, or provide a sensitivity scan. With that addressed, I'd be comfortable with it in a good specialty journal.","headline":"Solid first-principles study of anisotropic exciton transport in phosphorene, but the headline low-temperature exchange enhancement is quantitatively underdetermined by an unspecified Lorentzian dephasing width.","tokens_in":13982,"tokens_out":1846,"would_cite":true,"duration_ms":21394,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["71.35.-y","72.20.-i"],"model":"deepseek-v4-flash","headline":"Phosphorene channels exciton energy along one crystal direction, with long-range exchange interactions sharpening the directional contrast at low temperatures.","keywords":["exciton transport","phosphorene","anisotropic diffusion","long-range exchange interaction","exciton-phonon coupling","two-dimensional semiconductors","first-principles many-body theory","armchair direction"],"falsifier":"Measure the two in-plane diffusion constants of an exciton population in a single phosphorene monolayer (on SiO₂ or suspended) from 10 K to 300 K using spatially resolved photoluminescence or transient absorption microscopy after a focused excitation spot, and compare the ratio $D_{\\mathrm{AC}}/D_{\\mathrm{ZZ}}$ with the predicted strong increase at low temperature; a plateau or decline of the ratio at low temperature would indicate that the exchange-induced linear dispersion is not the relevant transport channel.","tokens_in":12908,"feed_emoji":"⚡","tokens_out":6768,"duration_ms":56733,"temperature":0.7,"pith_summary":"Phosphorene is an atomically thin semiconductor whose crystal structure is strongly directional, and this paper argues that the directionality carries over to the motion of excitons—bound electron-hole pairs—making energy transport highly anisotropic. Combining a microscopic many-body description of excitons with first-principles electronic structure and phonon data, the authors compute optical absorption, phonon-mediated relaxation, and diffusion constants along the two in-plane axes. Their central finding is that both the anisotropic effective masses and, crucially, the long-range exchange interaction make the armchair direction far more conductive than the zigzag direction, with computed diffusion constants of $D_{\\mathrm{AC}}=7.40$ cm$^2$/s and $D_{\\mathrm{ZZ}}=0.48$ cm$^2$/s at 300 K. They also find that the exchange interaction's linear-in-momentum term strongly enhances the anisotropy at low temperature, because exciton population concentrates at small center-of-mass momentum where the exchange-induced dispersion is steep. If correct, phosphorene could serve as a directional channel for exciton energy in ultrathin optoelectronic devices.","feed_headline":"Excitons race 15x faster along one phosphorene axis","feed_subtitle":"First-principles calculation predicts anisotropy that grows when the material cools.","key_machinery":"The load-bearing objects are the Wannier equation for the exciton wavefunction, the Keldysh potential for the screened Coulomb interaction, and a microscopic exciton-phonon scattering rate built from first-principles deformation potentials. The directional contrast is carried by the long-range exchange interaction $K_Q = \\frac{V_Q}{A}\\left(\\sum_k \\phi(k)\\,\\mathbf Q\\cdot\\mathbf d_{k+Q}\\right)\\left(\\sum_{k'}\\phi(k')\\,\\mathbf Q\\cdot\\mathbf d_{k'-Q}\\right)$, in which the interband dipole $\\mathbf d$ is oriented along the armchair direction, so the scalar product $\\mathbf Q\\cdot\\mathbf d$ vanishes for zigzag momentum and the small-$Q$ dispersion is linear only along armchair. Because the diffusion constant $D_i = \\sum_{Q\\eta} \\frac{(v^{\\eta,i}_Q)^2}{2\\Gamma^\\eta_Q}\\, e^{-E_{Q,\\eta}/k_B T}/Z$ weights states by their group velocity and inverse dephasing, the exchange-distorted dispersion directly translates into enhanced AC mobility that becomes most pronounced at low temperature, where population piles up at small $Q$.","core_discovery":"The paper's central claim is that the intrinsic in-plane anisotropy of phosphorene suffices to produce ultrafast, highly directional exciton transport, and that the long-range exchange interaction is a primary driver of the directional contrast. In the authors' calculation, the bright 1s exciton has an elliptical momentum-space wavefunction tied to the anisotropic band masses, the interband dipole points along the armchair direction, and the exciton energy dispersion acquires a linear-in-momentum contribution from exchange only along that dipole direction. Since the diffusion constant is a thermal average of group-velocity squared divided by the phonon-induced dephasing rate, the combination of a light effective mass and a steep exchange dispersion along armchair yields $D_{\\mathrm{AC}}=7.40$ cm$^2$/s versus $D_{\\mathrm{ZZ}}=0.48$ cm$^2$/s at 300 K, with the AC value growing sharply as temperature drops. The exchange term also suppresses the zero-momentum linewidth by reducing the low-energy density of states, which further helps mobility at low temperature.","pith_inferences":["The same exchange-enhancement mechanism should operate in other anisotropic two-dimensional semiconductors with an in-plane dipole, such as black arsenic or ReS₂, where the magnitude of the effect would depend on dipole strength and dielectric screening.","Strain engineering that rotates or renormalizes the interband dipole could switch the fast-transport axis, a control knob the present calculation does not explore but directly follows from Eq. (8).","The calculations neglect short-range exchange and intervalley exchange for the dark ΓV excitons; including those terms could reduce the low-temperature boost, so the sharp rise of $D_{\\mathrm{AC}}$ below roughly 100 K is the most fragile quantitative prediction.","Time-resolved spatial imaging of the exciton cloud after a focused pump pulse should reveal a characteristic elliptical profile whose aspect ratio varies with temperature, providing a direct test of the exchange-enhanced anisotropy."],"forward_implications":["Phosphorene monolayers on SiO₂ should show exciton diffusion lengths several times larger along the armchair axis than along the zigzag axis at room temperature, with the contrast growing at cryogenic temperatures.","The computed 300 K average diffusion constant, about 4 cm²/s, is consistent with the experimental value reported for amorphous phosphorene, so the claimed anisotropy should be observable in a crystalline monolayer.","The reduction of the zero-momentum linewidth by exchange means that the bright exciton emission line should narrow more than a mass-only model predicts, and the narrowing should be anisotropic.","Exciton energy can be channeled along the armchair direction over ultrafast timescales, supporting designs such as directional exciton highways and linearly polarized light emission."],"supporting_citations":[{"why":"Supplies the microscopic exciton-phonon scattering and relaxation-cascade formalism used for dephasing and population dynamics.","marker":"[39]"},{"why":"Provides the long-range (dipole-dipole) exchange interaction form used in Eq. (8) to modify the exciton dispersion.","marker":"[62]"},{"why":"Establishes the anisotropic Fickian diffusion framework and the diffusion-constant expression used for the transport calculation.","marker":"[34]"},{"why":"Provides the experimental amorphous-phosphorene diffusion value used to benchmark the computed room-temperature average.","marker":"[65]"},{"why":"Underlies the Wannier equation, Elliot formula, and semiconductor Bloch equations from which the optical and transport quantities are derived.","marker":"[50]"},{"why":"Justifies neglecting the short-range exchange interaction in two-dimensional materials, a premise of the exchange model.","marker":"[63]"},{"why":"Provides the measured 1s exciton binding energy in phosphorene on SiO₂ used to validate the Wannier-equation calculation.","marker":"[59]"}],"fun_headline_variants":["Phosphorene excitons travel 15x faster along one axis","Cooling sharpens phosphorene's directional exciton transport","Exchange interactions steer anisotropic exciton flow in phosphorene","Phosphorene excitons pick a fast axis at low temperature"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire low-temperature enhancement of the armchair diffusion constant rests on the long-range exchange term of Eq. (8), evaluated with the interband dipole fixed along armchair, applied only to the bright ΓΓ excitons, and with short-range exchange neglected; if that model is quantitatively incomplete, the predicted sharp rise in $D_{\\mathrm{AC}}$ at low temperature is not guaranteed.","fun_headline_variants_meta":{"raw":{"variants":["Phosphorene excitons travel 15x faster along one axis","Cooling sharpens phosphorene's directional exciton transport","Exchange interactions steer anisotropic exciton flow in phosphorene","Phosphorene excitons pick a fast axis at low temperature"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001043,"raw_usage":{"total_tokens":4361,"prompt_tokens":899,"completion_tokens":3462,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":515,"completion_tokens_details":{"reasoning_tokens":3389}},"tokens_in":515,"tokens_out":3462,"duration_ms":24025,"temperature":1.0,"reasoning_tokens":3389,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T16:35:59.312457+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the two in-plane diffusion constants of an exciton population in a single phosphorene monolayer (on SiO₂ or suspended) from 10 K to 300 K using spatially resolved photoluminescence or transient absorption microscopy after a focused excitation spot, and compare the ratio $D_{\\mathrm{AC}}/D_{\\mathrm{ZZ}}$ with the predicted strong increase at low temperature; a plateau or decline of the ratio at low temperature would indicate that the exchange-induced linear dispersion is not the relevant transport channel.","supporting_citations":[{"cited_title":"Anisotropic plasmons, excitons, and electron energy loss spectroscopy of phosphorene,","cited_arxiv_id":null,"evidence_quote":"Supplies the microscopic exciton-phonon scattering and relaxation-cascade formalism used for dephasing and population dynamics."},{"cited_title":"Op- 9 tical tuning of exciton and trion emissions in monolayer phos- phorene,","cited_arxiv_id":null,"evidence_quote":"Provides the long-range (dipole-dipole) exchange interaction form used in Eq. (8) to modify the exciton dispersion."},{"cited_title":"2d materials for optical modulation: challenges and op- portunities,","cited_arxiv_id":null,"evidence_quote":"Establishes the anisotropic Fickian diffusion framework and the diffusion-constant expression used for the transport calculation."},{"cited_title":"Valley- exchange coupling probed by angle-resolved photolumines- cence,","cited_arxiv_id":null,"evidence_quote":"Provides the experimental amorphous-phosphorene diffusion value used to benchmark the computed room-temperature average."},{"cited_title":"Dielectric properties of hexagonal boron nitride and transition metal dichalcogenides: from monolayer to bulk,","cited_arxiv_id":null,"evidence_quote":"Underlies the Wannier equation, Elliot formula, and semiconductor Bloch equations from which the optical and transport quantities are derived."},{"cited_title":"Semiconduc- tor excitons in new light,","cited_arxiv_id":null,"evidence_quote":"Justifies neglecting the short-range exchange interaction in two-dimensional materials, a premise of the exchange model."},{"cited_title":"Ex- citonic thermalization bottleneck in twisted tmd heterostruc- tures,","cited_arxiv_id":null,"evidence_quote":"Provides the measured 1s exciton binding energy in phosphorene on SiO₂ used to validate the Wannier-equation calculation."}],"review_version":2}