{"id":"b3c7aec0-cacc-4c53-9076-e348083ef5e6","arxiv_id":"2506.05075","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"An effective-mass model with open-boundary quasi-bound states estimates exciton binding energies of roughly 0.5 eV and room-temperature radiative lifetimes of 0.13 to 0.43 ns for monolayer TMDs, matching published PL decay data within about a factor of two.","lead":"This paper computes exciton binding energies and radiative lifetimes for four monolayer TMD semiconductors (WS2, WSe2, MoS2, MoSe2) using an effective-mass quantum well model with open boundary conditions. The authors compare their results with photoluminescence measurements and argue the model can guide optoelectronic device design.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (8) uses the envelope-function product ψ_e ψ_h instead of the probability densities |ψ_e|^2|ψ_h|^2 in the Coulomb kernel; this changes both the binding energy and the lifetime overlap factor, so the central agreement claim is unverified until the kernel is corrected and recomputed.","rationale":"The reader's weakest_assumption correctly identifies Eq. (8) and, in particular, the use of ψ_e ψ_h instead of probability densities in the Coulomb kernel. I agree that this is the most load-bearing premise: it enters both the binding-energy calculation (through V_C) and the lifetime calculation (through the envelope overlap in the oscillator strength), so an error here invalidates the two central quantitative outputs. The concern is an internal inconsistency, not merely a disagreement with an external consensus: the printed equation does not give the correct expectation value of the Coulomb interaction for the stated wavefunction ansatz. A direct recomputation with the corrected kernel is a concrete and decisive test. The paper also has other issues (garbled Table II, no code, overstated agreement), but those are secondary to the kernel problem. Because the flaw is formal and fixable, and the numerical consequences are not yet quantified, the appropriate disposition remains conditional rather than outright rejection; the current conclusions are unverified until the correction is made and the tables are regenerated. I therefore leave the reader's verdict unchanged.","tokens_in":11270,"tokens_out":9388,"duration_ms":117411,"concrete_test":"Recompute the four monolayer entries of Table II after replacing ψ_eψ_h by |ψ_e|^2|ψ_h|^2 in Eq. (8), keeping every other numerical setting identical (QTBM subband states, finite-difference grid with 300 points, integration limits, all material parameters from Table I). Report the new E_X, E_b, and τ_eff for WS2, WSe2, MoS2, and MoSe2, and compare them with the current Table II and with the cited PL decay data (Ref. 8). A useful internal check is that the corrected kernel satisfies ∫∫ |ψ_e|^2|ψ_h|^2 dx_e dx_h = 1, while the printed kernel does not. If the shifts in E_b or τ_eff exceed ~10–20% for any material, the reported agreement is not robust to this correction; if the shifts are negligible, the printed kernel can be regarded as a typographical error and the numerical results stand.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing problem is Eq. (8). The Coulomb potential entering the Wannier equation is written as V_C(ρ) = −(e^2/ε_b) ∫∫ ψ_e(x_e)ψ_h(x_h) / √((x_e−x_h)^2+ρ^2) dx_e dx_h. For the potential that should appear in the relative-motion equation, the correct matrix element of 1/|r_e−r_h| between the confinement subbands is ∫∫ |ψ_e(x_e)|^2 |ψ_h(x_h)|^2 / √((x_e−x_h)^2+ρ^2) dx_e dx_h (for real envelope functions). The printed form uses amplitudes rather than probability densities; because ψ_e and ψ_h are normalized with ∫|ψ|^2 dx = 1 but ∫ψ dx is not 1, this is not a normalization convention but a different, incorrect kernel. This same kernel sets the in-plane potential in Eqs. (4)–(8) and, through the envelope overlap in Eqs. (10)–(11), the oscillator strength and radiative lifetime. If the code actually used the printed amplitudes, both E_b and τ_eff in Table II are not the solution of the stated Hamiltonian; if the code used densities, then the equation is a misleading typo and the numerical claims cannot be reproduced from the manuscript. Either way, the central claim that the model uses only physical parameters and reproduces PL lifetimes and DFT binding energies cannot be assessed as written.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a self-contained effective-mass model for Wannier-Mott excitons in monolayer TMDs (WS2, WSe2, MoS2, MoSe2). It solves the one-dimensional Schrödinger equation with open boundary conditions (QTBM) for electron and hole quasi-bound states, then solves a radial Wannier equation whose Coulomb potential is computed from the 1D envelope functions, and finally estimates radiative lifetimes from the oscillator strength with a thermalization factor. The central claim is that this model, using only independently sourced physical parameters, gives exciton energies, binding energies, and radiative lifetimes in good agreement with TRPL measurements and DFT estimates. The paper validates the method on a III-V double quantum well and reports binding energies around 0.47-0.54 eV and effective lifetimes of 0.13-0.43 ns at 300 K.","tokens_in":11570,"tokens_out":7159,"duration_ms":78423,"significance":"The proposed approach is attractive as a design-oriented estimator: the input parameters are not fitted to the TMD PL data, and the III-V benchmark is an independent check. If the central claims survive correction, the method would provide a simple route to estimate binding energies and radiative lifetimes for TMD monolayers. However, the current manuscript does not establish those claims: the printed Coulomb kernel in Eq. (8) is not the correct matrix element, the lifetime comparison is only factor-of-two accurate for WS2 and WSe2, and several comparisons (notably the 4 K lifetimes) are not numerically documented. No code is shipped, so reproducibility rests entirely on the equations, which makes the Eq. (8) issue a substantive obstacle.","major_comments":[{"comment":"The Coulomb kernel is written as V_C(rho) = -(e^2/epsilon_b) integral integral psi_e(x_e) psi_h(x_h) / sqrt((x_e-x_h)^2+rho^2) dx_e dx_h. For real envelope functions normalized as integral |psi|^2 dx = 1, the correct matrix element of 1/|r_e-r_h| between the confinement subbands is integral integral |psi_e(x_e)|^2 |psi_h(x_h)|^2 / sqrt((x_e-x_h)^2+rho^2) dx_e dx_h. The printed form is not a normalization convention but a different kernel; it can even change sign if the envelopes have nodes. Since this kernel enters both the in-plane potential in Eqs. (4)-(8) and the oscillator-strength overlap in Eqs. (10)-(11), the reported E_b and tau_eff in Table II are not the solution of the stated Hamiltonian. Please correct the kernel, recompute Table II, and state whether the numerical results change.","section":"II, Eq. (8)"},{"comment":"The abstract and Section IV claim 'good agreement' with TRPL measurements, but the computed tau_eff values are 0.128 ns (WS2), 0.208 ns (WSe2), 0.354 ns (MoS2), and 0.427 ns (MoSe2) against the cited PL averages 0.22, 0.38, 0.42, and 0.36 ns, respectively. The first two disagree by a factor of about 1.7-1.8 and no uncertainties are reported on the computed values. Please provide an uncertainty estimate propagated from epsilon_b, d_cv, effective masses, and band offsets, and either define a quantitative agreement criterion or temper the claim of good agreement.","section":"IV, Table II and Fig. 6"},{"comment":"The text states that the low-temperature (4 K) lifetimes of 2.3-5 ps 'nearly match' the model, but Table II contains only 300 K values and the computed 4 K lifetimes are not given numerically. Please report the 4 K tau_eff values used in Fig. 9 so this comparison can be verified.","section":"IV, Fig. 9 and low-temperature discussion"},{"comment":"The III-V double-quantum-well validation is described as matching 'quantitatively', but the figure is not accompanied by the model parameters, the numerical values being compared, or error bars. Since this is the only independent benchmark in the paper, please include a table or text giving the parameters and the quantitative comparison.","section":"III.A, Fig. 2"},{"comment":"The central results depend sensitively on epsilon_b and d_cv, but no sensitivity analysis is presented. Both the binding energy and the radiative lifetime scale with these inputs, and epsilon_b is taken from zero-strain DFT without considering the dielectric environment (e.g., hBN encapsulation or substrate). Please add a sensitivity study over a physically reasonable range of epsilon_b and d_cv, and discuss how the claimed agreement would be affected.","section":"Tables I-II and Section IV"}],"minor_comments":[{"comment":"Table II is typeset with interleaved citation markers and data (e.g., '17 PL Decay 8' and '710 +/- 1032'), making the column meanings ambiguous; please retypeset it with clearly separated parameter, this-work, and literature columns.","section":"IV, Table II"},{"comment":"The text says 'ground-state wavefunctions shown in Fig. 5 in the device direction (x)', but Fig. 5 displays the in-plane exciton wavefunctions while the device-direction wavefunctions are shown in Fig. 3; please correct the cross-reference.","section":"IV, first paragraph"},{"comment":"The phrase 'As seen in Fig. ref fig4' is an unresolved cross-reference and should read 'Fig. 4'.","section":"IV, first paragraph"},{"comment":"The excitonic wavefunction in Eq. (3) contains an in-plane factor psi(rho)/rho, which is unconventional and potentially singular at rho=0; please define the limiting behavior and normalization used when Eq. (11) evaluates the wavefunction at rho=0.","section":"II, Eq. (3) and Eq. (11)"},{"comment":"There are incomplete sentences in the introduction, such as 'WSe2 for and MoSe2 for and other information technology platforms21' and 'while monolayer MoS2 and MoS2 has been utilized in spintronic 20'; these should be rewritten.","section":"I, Introduction"}],"recommendation":"major_revision","confidential_remarks":"This is a methods-style paper within the scope of cond-mat.mes-hall. The main obstacle is not the physical picture but the verifiability of the numerics: the printed Coulomb kernel in Eq. (8) appears incorrect, the central lifetime comparison is only factor-of-two accurate in two cases, and the 4 K comparison is not documented. I would not recommend rejection if the authors can correct the kernel, rerun the calculations, and add uncertainty and sensitivity analyses. If the corrected results change materially, a further round of review may be needed before acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Honest take: the paper has a useful idea and a clear independent validation, but as written the central Coulomb kernel in Eq. (8) is wrong or misprinted, so the numbers in Table II can't be trusted until that's fixed. The QTBM/open-boundary EMA treatment of monolayer TMDs is new in this combination, and the III-V double quantum well benchmark is a genuine check that the machinery works for confined systems. Credit where due: the model uses only physical inputs, the computed lifetimes land in the same range as PL data, and the authors are candid about the MoS2 discrepancy and non-radiative processes.\n\nWhere it gets soft: Eq. (8) writes the Coulomb matrix element as ∫∫ ψ_e ψ_h / r, with no absolute squares. For a two-band envelope calculation the correct kernel is ∫∫ |ψ_e|² |ψ_h|² / r. Because the envelope functions are normalized but can change sign, the printed form is not a convention issue; it will give a different potential and a different overlap factor in the oscillator strength. If the code actually used densities, the equation is a critical typo and the manuscript can't be reproduced; if it used amplitudes, the stated Hamiltonian isn't the one solved. Either way, the central 'good agreement' claim is unverified. The thermalization formula (Eq. 17) is imported from the literature without derivation, and the d_cv values in Table II appear inconsistent (same value for WS2 and WSe2). Table II itself is garbled: column headers misaligned, experimental values jumbled with superscripts. No code, no uncertainty propagation. These are fixable, but as submitted the paper doesn't support its headline.\n\nThe stress-test note about Eq. (8) holds up on reading. It's the load-bearing issue. The authors should correct the kernel, recompute, and provide a clean table and data/code. If the corrected results remain in the same ballpark, the paper is a reasonable contribution for device screening. For now I'd send it to review rather than desk reject, because the approach is legitimate and the validation is real, but the authors need to address the kernel before publication.","headline":"A useful, parameter-light EMA/QTBM model for TMD exciton lifetimes with a real III-V benchmark, but Eq. (8) has a wrong/typo Coulomb kernel that invalidates the numbers as printed.","tokens_in":12145,"tokens_out":2053,"would_cite":false,"duration_ms":23823,"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":"A Wannier-exciton quantum-well model reproduces monolayer TMD binding energies and lifetimes from physical inputs alone.","keywords":["exciton binding energy","radiative lifetime","transition metal dichalcogenides","monolayer TMDs","Wannier-Mott exciton","effective mass approximation","quantum transmitting boundary method","thermalization"],"falsifier":"Repeat the same calculation for a single monolayer with a screened Coulomb interaction computed from an explicitly nonlocal dielectric function, for example a GW-BSE treatment, instead of the single static $\\epsilon_b$ of Eq. (8); if the binding energies and radiative lifetimes shift by more than the experimental spread relative to the PL data cited in Table II, the static-dielectric kernel is the failing assumption. A complementary experiment is to measure the $1s$ exciton energy and decay lifetime of one monolayer on substrates with different dielectric constants and check whether the model's fixed-$\\epsilon_b$ prediction holds for both.","tokens_in":11017,"feed_emoji":"💡","tokens_out":11251,"duration_ms":125431,"temperature":0.7,"pith_summary":"This paper claims that a Wannier–Mott exciton in a monolayer transition metal dichalcogenide can be described by a one-dimensional quantum-well model whose only inputs are physical parameters: band offsets, effective masses, the bandgap, the dielectric constant, and the transition dipole moment. Solving the open-boundary Schrödinger equation for electron and hole quasi-bound states and then the Wannier equation for the in-plane $1s$ state yields exciton energies, binding energies, and radiative decay widths by Fermi's golden rule. After averaging the decay width over a thermal distribution of in-plane momenta, the model gives room-temperature effective lifetimes of about $0.13$–$0.43$ ns for WS$_2$, WSe$_2$, MoS$_2$, and MoSe$_2$, with binding energies in the range of roughly $0.47$–$0.59$ eV. The authors report that these values match time-resolved photoluminescence measurements and DFT estimates, which matters because exciton binding energy and lifetime are the two parameters that most directly set the usefulness of a TMD monolayer in light-emitting and photodetector devices.","feed_headline":"One model predicts exciton binding and lifetimes in 2D semiconductors","feed_subtitle":"It needs only band offsets, masses, and dielectric constants — and it matches photoluminescence decay and DFT.","key_machinery":"The load-bearing object is the Coulomb kernel $V_C(\\rho)$ in Eq. (8), formed by convolving the $1/r$ electron–hole attraction with the product of the confined electron and hole wavefunctions from Eq. (2). This single kernel does double duty: it sets the in-plane potential in the Wannier equation and therefore the binding energy, and through the value of the in-plane wavefunction at $\\rho=0$ it sets the envelope overlap that enters the oscillator strength and the radiative decay width. The numerical solver is the quantum transmitting boundary method, which turns the open-boundary quantum well into an eigenvalue problem and supplies the complex quasi-bound energies. The final step is the thermalization average of Eq. (17), which weights the decay width over a Maxwell–Boltzmann distribution of exciton momenta.","core_discovery":"The central claim is that a single parameter-light model, using only physical inputs and no fitted exciton parameters, can reproduce the exciton properties of four monolayer TMDs. The electron and hole are treated as quasi-bound states of a finite quantum well with open boundary conditions; their complex eigenenergies give the confinement part of the exciton. The in-plane relative motion is governed by the Wannier equation with a Coulomb kernel built from the envelope-function overlap of Eq. (8), and the radiative decay width follows from Fermi's golden rule for the transition dipole. The final step, Eq. (17), thermally averages the decay width over a Maxwell–Boltzmann distribution of exciton momenta, which is what converts the picosecond intrinsic radiative lifetimes into the hundreds-of-picosecond lifetimes seen at room temperature. On that basis the paper reports exciton energies, binding energies of roughly $0.47$–$0.59$ eV, and effective lifetimes of $0.13$–$0.43$ ns for WS$_2$, WSe$_2$, MoS$_2$, and MoSe$_2$, in the range of measured photoluminescence decay and DFT data.","pith_inferences":["The paper fixes the dielectric constant from zero-strain DFT; an easy extension is to treat $\\epsilon_b$ as a function of the surrounding dielectric environment, which would predict that encapsulated monolayers have smaller binding energies and longer effective lifetimes than suspended ones.","Eq. (8) uses the product $\\psi_e(x_e)\\psi_h(x_h)$, not the product of probability densities; if an independent many-body calculation shows that the relative phase between the envelope functions matters, the oscillator strength and hence the lifetime would need revision.","The thermal average in Eq. (17) assumes equal population of all exciton spin states; valley-selective excitation or a magnetic field should break this and produce a polarization-dependent lifetime that the model does not yet describe."],"forward_implications":["The computed values in Table II give device designers a ready set of exciton binding energies and room-temperature lifetimes for WS$_2$, WSe$_2$, MoS$_2$, and MoSe$_2$ without needing a many-body calculation.","Because the inputs are physical parameters rather than fitted exciton parameters, the same procedure can be rerun for other TMD compositions, layer thicknesses, or dielectric environments by updating the input table.","The thermalization step explains why room-temperature PL lifetimes are orders of magnitude longer than low-temperature radiative lifetimes while remaining radiative in origin.","The successful validation against a III-V double-quantum-well photoluminescence measurement indicates the model transfers beyond TMDs to other heterostructures.","The authors state that the model can be extended to multilayer structures, interlayer excitons, and trions by modifying the oscillator-strength expression."],"supporting_citations":[{"why":"Supplies the dielectric constant $\\epsilon_b$ used in the Coulomb kernel and the radiative linewidth expression.","marker":"40"},{"why":"Supplies band gaps and band offsets that define the quantum-well potential for each monolayer.","marker":"22"},{"why":"Supplies effective masses and the DFT/SVM binding-energy values against which the model is compared.","marker":"62"},{"why":"Provides the room-temperature time-resolved photoluminescence decay times used as the experimental comparison for effective lifetimes.","marker":"8"},{"why":"Provides the h-BN/WS$_2$ coupled-quantum-well photoluminescence data used for the WS$_2$ exciton-energy comparison.","marker":"30"},{"why":"Provides DFT-based exciton radiative lifetimes and low-temperature data that the model is compared with.","marker":"23"},{"why":"Supplies the thermalization procedure that converts intrinsic radiative lifetimes into the room-temperature effective lifetimes.","marker":"44"},{"why":"Introduces the quantum transmitting boundary method used to solve the open-boundary Schrödinger equation for quasi-bound states.","marker":"46"},{"why":"Provides the GaAs/AlGaAs double-quantum-well photoluminescence kinetics used to validate the model before applying it to TMDs.","marker":"48"},{"why":"Provides the Fermi-golden-rule oscillator-strength and lifetime relations used for the decay rate.","marker":"41"}],"fun_headline_variants":["One model spans four monolayer TMD exciton lifetimes","Exciton lifetimes from three physical inputs in 2D TMDs","No fitting: exciton lifetimes match measured decay","Exciton binding and lifetime from a single model"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the electron–hole attraction can be represented by a single static dielectric constant from zero-strain DFT and by the product of the electron and hole envelope functions in Eq. (8); if real screening by the environment or strain, or the correct wavefunction overlap, differs from that, both the predicted binding energy and the predicted lifetime move.","fun_headline_variants_meta":{"raw":{"variants":["One model spans four monolayer TMD exciton lifetimes","Exciton lifetimes from three physical inputs in 2D TMDs","No fitting: exciton lifetimes match measured decay","Exciton binding and lifetime from a single model"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001895,"raw_usage":{"total_tokens":7444,"prompt_tokens":976,"completion_tokens":6468,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":592,"completion_tokens_details":{"reasoning_tokens":6403}},"tokens_in":592,"tokens_out":6468,"duration_ms":53368,"temperature":1.0,"reasoning_tokens":6403,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T10:25:53.475113+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Repeat the same calculation for a single monolayer with a screened Coulomb interaction computed from an explicitly nonlocal dielectric function, for example a GW-BSE treatment, instead of the single static $\\epsilon_b$ of Eq. (8); if the binding energies and radiative lifetimes shift by more than the experimental spread relative to the PL data cited in Table II, the static-dielectric kernel is the failing assumption. A complementary experiment is to measure the $1s$ exciton energy and decay lifetime of one monolayer on substrates with different dielectric constants and check whether the model's fixed-$\\epsilon_b$ prediction holds for both.","supporting_citations":[{"cited_title":"Deng , author L","cited_arxiv_id":null,"evidence_quote":"Supplies band gaps and band offsets that define the quantum-well potential for each monolayer."},{"cited_title":"a np \\\"a \\","cited_arxiv_id":null,"evidence_quote":"Supplies effective masses and the DFT/SVM binding-energy values against which the model is compared."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the room-temperature time-resolved photoluminescence decay times used as the experimental comparison for effective lifetimes."},{"cited_title":"\\ Lee , author D","cited_arxiv_id":null,"evidence_quote":"Provides the h-BN/WS$_2$ coupled-quantum-well photoluminescence data used for the WS$_2$ exciton-energy comparison."},{"cited_title":"Robert , author D","cited_arxiv_id":null,"evidence_quote":"Supplies the thermalization procedure that converts intrinsic radiative lifetimes into the room-temperature effective lifetimes."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the quantum transmitting boundary method used to solve the open-boundary Schrödinger equation for quasi-bound states."},{"cited_title":"Butov , author A","cited_arxiv_id":null,"evidence_quote":"Provides the GaAs/AlGaAs double-quantum-well photoluminescence kinetics used to validate the model before applying it to TMDs."},{"cited_title":"Burstein \\ and\\ author C","cited_arxiv_id":null,"evidence_quote":"Provides the Fermi-golden-rule oscillator-strength and lifetime relations used for the decay rate."}],"review_version":1}