{"id":"540a01e9-79fd-4ebf-a591-a6b3dec882f3","arxiv_id":"2501.01410","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"Biaxial strain activates symmetry-forbidden A' phonon-assisted recombination in monolayer MoS2, predicting a peak in indirect photoluminescence at +1% strain.","lead":"A computational study shows that stretching a single layer of MoS2 can switch on a normally forbidden light-emission pathway, with the effect strongest at 1% strain. The finding suggests strain engineering could make molybdenum-based 2D semiconductors brighter emitters and explains strain-dependent photoluminescence seen in WSe2.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Phase/gauge mismatch admitted by the authors between DFPT g-elements and BSE A-coefficients can reverse computed intensity trends, so the +1% strain PL peak is not quantitatively secure.","rationale":"The reader's weakest_assumption names the same issue: the gauge/phase mismatch between DFPT g-elements and BSE A-coefficients. I agree this is the most load-bearing concern because the paper's headline quantitative output — a sharp PL intensity maximum at +1% strain and assignment of the lower-energy replica to mode 8 — is exactly what Eq. (2) produces after coherent summation, and the authors explicitly concede that phase mismatch can change intensity levels and even reverse trends. The comparison against Chan et al. supports the methodology for unstrained linewidths, and the direct-PL inset agrees with experiment, so I would not recommend rejection; these anchors give the qualitative activation mechanism some credibility. But the strained PL enhancement is not a symmetry-forced result: symmetry only decides which modes are allowed, not whether the allowed mode-8 coupling dominates and peaks at +1%. Without a gauge-consistent recalculation or strained-MoS2 PL data, the central quantitative claim cannot be distinguished from a phase artifact. The abstract's 'no phonon-assisted PL' phrasing is also too strong relative to the body's 'faint indirect emission' under unstrained conditions, but that is a presentation issue; the phase issue is the technical load-bearing one. Therefore the verdict should remain CONDITIONAL, pending the concrete test.","tokens_in":15820,"tokens_out":4211,"duration_ms":41858,"concrete_test":"Recompute the lowest-exciton |Gβλ,ν(Q,q)|^2 and PL spectra at -0.5%, 0%, +0.75%, +1%, +2%, +3% using a gauge-consistent workflow: propagate the exact QE Kohn-Sham phases into the Yambo A-coefficients (or use the Yambo-PERTURBO interface used successfully for hBN), then compare mode-resolved couplings and Fig. 5. If the +1% intensity maximum and the mode-8 hole-phonon replica survive unchanged, the phase-mismatch concern is retired; if the peak shifts in strain, reverses, or mode 8 drops to noise level, the central claim must be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central strain-activation prediction rests on combining electron-phonon matrix elements from Quantum Espresso with BSE exciton coefficients from Yambo through Eq. (2). The authors themselves state in the Limitations paragraph that this combination is only valid if the same Kohn-Sham wavefunction phases enter both steps, that current QE/Yambo interfaces cannot enforce this, that \"there could be a phase-mismatch which can show an opposite trend in the intensity levels,\" and that their own Fig. 2 contains an \"anomalous emergence of high coupling strength near the Γ point\" attributed to phase mismatch. Since the +1% peak intensity and the identification of mode 8 as the dominant hole-phonon scattering channel (Fig. 2c and Fig. 6d) are precisely the quantities most sensitive to such coherent-summation/cancellation errors, the central claim that biaxial strain activates a maximum phonon-assisted PL at +1% is not currently established quantitatively. The qualitative symmetry selection rule (A' allowed, A'' forbidden) is plausible and supported by group theory, but symmetry alone does not fix the magnitude or strain dependence of the allowed couplings. The unstrained case is also described inconsistently: the abstract says no phonon-assisted PL, while the body reports faint indirect emission near 1.73 eV under unstrained conditions, so the \"forbidden\" baseline is overstated.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports first-principles calculations of exciton-phonon coupling and phonon-assisted photoluminescence in monolayer MoS2 under biaxial strain. Using finite-momentum BSE exciton coefficients from Yambo and DFPT electron-phonon matrix elements from Quantum Espresso, the authors compute the exciton-phonon coupling strengths, scattering rates, linewidths, and PL spectra through the van Roosbroeck-Shockley relation. Group-theoretic analysis in D3h and its subgroups identifies A' phonon modes as symmetry-allowed and A'' modes as forbidden for intervalley exciton scattering. The central claim is that unstrained MoS2 has no phonon-assisted PL at cryogenic temperature, while tensile biaxial strain activates A'-mediated channels, mainly hole-phonon scattering by an out-of-plane optical mode (mode 8, about 42 meV), producing a maximum indirect PL near +1% strain. The unstrained linewidth and exciton dispersion are benchmarked against Chan et al. and Wu et al., and the computed direct PL is compared with experimental data.","tokens_in":16162,"tokens_out":4435,"duration_ms":44815,"significance":"If the central result holds, the work would establish a symmetry-based strain-engineering route to brighten dark intervalley excitons in Mo-based TMD monolayers and would offer a plausible explanation for the non-monotonic strain dependence of PL in WSe2. The paper's strengths are its explicit group-theoretic selection argument, its benchmarking against external unstrained calculations (Chan et al. for linewidths, Wu et al. for the exciton dispersion) and experiment, and its falsifiable prediction of a maximum phonon-assisted PL near +1% biaxial strain. However, the admitted wavefunction-phase mismatch between the DFPT and BSE codes is not a peripheral numerical detail: it enters the central coupling matrix element of Eq. (2) and, by the authors' own statement, can reverse intensity trends. The quantitative strain dependence and the assignment of mode 8 as the dominant hole-phonon channel are therefore not yet secure, although the qualitative symmetry selection argument is plausible.","major_comments":[{"comment":"The authors explicitly state that Eq. (2) requires the same Kohn-Sham wavefunction phases in the DFPT and BSE steps, that current Quantum Espresso/Yambo interfaces cannot enforce this, and that 'there could be a phase-mismatch which can show an opposite trend in the intensity levels.' Because the central prediction of a +1% strain maximum and the assignment of mode 8 as the dominant hole-phonon channel (Fig. 2(c) and Fig. 6(d)) are based on the magnitudes of the gauge-sensitive complex sums in Eq. (2), the admitted phase mismatch directly undermines the quantitative strain-dependence claim. The anomalous strong coupling near the Gamma point in Fig. 2, attributed by the authors to this mismatch, confirms that the contamination is present in the plotted quantity. The authors should either repeat the calculation with a phase-consistent interface (for example, the YAMBO-PERTURBO route they cite) or provide an explicit gauge-invariance check, such as randomizing wavefunction phases and showing that the couplings and PL spectra are stable.","section":"Limitations paragraph (before Summary)"},{"comment":"The abstract claims that unstrained MoS2 exhibits no phonon-assisted PL at cryogenic temperatures, but the body reports that 'faint indirect emission is visible near 1.73 eV' under unstrained conditions in the discussion of Fig. 7, and Fig. 7(a) shows a weak indirect feature. This internal inconsistency makes the baseline 'forbidden' claim stronger than the computed result supports. Please state quantitatively whether the unstrained emission is zero, negligible, or merely weak, and adjust the abstract and summary wording accordingly.","section":"Abstract and Fig. 7 discussion"},{"comment":"The linewidth curves in Fig. 4 are obtained by fitting Toyozawa and linear-plus-Bose models with fitted parameters (S_A, S_O, E_A, E_O, a, b, gamma0), and the PL calculation uses a hand-set 2 meV damping factor. These auxiliary models are not derived from the ab initio calculation and introduce several free parameters. The authors should test and report the sensitivity of the predicted +1% strain maximum to reasonable variations of these parameters, so that the reader can see that the central strain trend is not an artifact of the fitting choices.","section":"Eqs. (4)-(5) and damping in Eq. (6)"}],"minor_comments":[{"comment":"The symbols D^±_{beta lambda, q nu}, R_lambda, and the sign conventions in omega ∓ 2 omega_{q mu} are not fully defined in the text; please define all symbols or refer the reader to the relevant equation in the Supplemental Information.","section":"Eq. (6)"},{"comment":"The legend labels the curves as 'Exciton 1 & 2', 'Exciton 3 & 4', and 'Exciton 5', while the text says the lowest five excitons are considered; please clarify whether excitons 2 and 4 are exactly degenerate with excitons 1 and 3, respectively, and whether they are omitted from the plot or simply degenerate.","section":"Fig. 3(a)"},{"comment":"The infinitesimal eta appearing in the self-energy and PL expressions is never given a numerical value; please specify the broadening parameter used and state whether the results are converged with respect to it.","section":"Eq. (3) and Eq. (6)"},{"comment":"The text refers to 'Lechifflart et al.' in the limitations paragraph; the corresponding reference is listed as P. Lechifflart et al., and the spelling should be made consistent throughout.","section":"General presentation"}],"recommendation":"major_revision","confidential_remarks":"The admitted phase/gauge mismatch between the DFPT and BSE steps is the central obstacle to publication. The authors themselves note that it can reverse intensity trends, and the anomalous coupling near Gamma in Fig. 2 indicates the issue is present in the central quantity. If the authors can demonstrate gauge invariance or rerun with a consistent interface, the paper could be suitable. The abstract overstates the unstrained case relative to the body of the manuscript, and the fitted linewidth parameters add further uncertainty that should be addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThe genuinely new thing here is the symmetry-resolved finite-momentum BSE analysis showing that A′ phonon modes activate phonon-assisted PL in strained MoS2, with a predicted peak near +1% strain. The symmetry argument is coherent, and the unstrained benchmarks against Chan et al. and Wu et al. are solid. But the quantitative prediction rests on a phase-consistency approximation the authors themselves say cannot be enforced with current QE/Yambo interfaces, and which they admit can reverse intensity trends. So treat the +1% peak as plausible, not established.\n\nWhat the paper does well: it fills a real gap—phonon-assisted PL in Mo-based TMDs is much less studied than in W-based ones. The group-theoretic reduction to Cs(m) and the A′/A″ distinction is a clean, testable idea that does not depend on the phase issue. The unstrained linewidths match Chan et al., and the direct PL matches experiment, which shows the machinery is mostly under control. The authors also list their fitting parameters and write an unusually candid limitations paragraph acknowledging the phase mismatch and even attributing an anomalous coupling feature near Γ to it. That is honest scholarship.\n\nThe soft spots are real but proportioned. The load-bearing one is the phase mismatch. The coupling in Eq. (2) is a coherent sum of g-elements from QE and A-coefficients from Yambo; if the wavefunction phases are inconsistent, the sum can go wrong, and the authors themselves flag that this can show \"an opposite trend in the intensity levels.\" Since the strain-dependent intensity and the identification of mode 8 as the dominant channel are precisely what depend on coherent sums, the central prediction is fragile. The abstract overstates the unstrained case, saying there is no phonon-assisted PL, while the body and Fig. 7 show faint indirect emission near 1.73 eV; that inconsistency weakens the \"forbidden to allowed\" narrative. Finally, no code or data is released, and the SI link is a placeholder, which makes independent verification harder.\n\nWho this is for: anyone working on exciton-phonon coupling in TMDs, especially strain engineering of dark excitons. The symmetry selection rule is a useful contribution regardless of the numerical peak. It deserves a serious referee—the fix is a gauge-consistent calculation (e.g., via the PERTURBO-Yambo interface they cite) and softened claims about the unstrained case. My recommendation: send it to review; the referee should ask for that recalculation before the +1% peak is taken as quantitative.\n\nBest","headline":"A clean symmetry argument for strain-activated phonon-assisted PL in MoS2, but the quantitative +1% peak rests on a known phase-mismatch that could reverse the predicted intensity trends.","tokens_in":16721,"tokens_out":3156,"would_cite":true,"duration_ms":29524,"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":"Strain switches on phonon-assisted light emission in MoS2.","keywords":["exciton-phonon coupling","monolayer MoS2","photoluminescence","biaxial strain","dark excitons","symmetry selection rules","Bethe-Salpeter equation","first-principles calculation"],"falsifier":"Measure the cryogenic photoluminescence of a MoS2 monolayer under controlled biaxial tensile strain from 0 to +3% with a piezo stage. The central claim is falsified if no doublet of phonon sidebands appears near 50 and 70 meV below the bright exciton as strain approaches +1%, or if the indirect emission does not peak around +1% but instead rises monotonically or stays absent across the whole window.","tokens_in":15590,"feed_emoji":"💡","tokens_out":7431,"duration_ms":64512,"temperature":0.7,"pith_summary":"Monolayer MoS2 is a bright emitter at room temperature, but its lowest-energy intervalley excitons are momentum-dark and, this paper argues, cannot recombine with the help of phonons at cryogenic temperatures because the only phonons that could supply the momentum belong to a symmetry-forbidden representation (A''). The paper uses finite-momentum Bethe-Salpeter and density-functional perturbation calculations, backed by group-theoretic selection rules, to show that biaxial tensile strain activates the even A' phonon channels. Near +1% strain the phonon-assisted photoluminescence should peak sharply, dominated by a hole-phonon process that scatters the hole from K to K' and recombines there. If correct, this gives a concrete route to making dark excitons radiative in Mo-based dichalcogenides through strain engineering.","feed_headline":"Strain turns on dark-exciton luminescence in MoS2","feed_subtitle":"First-principles calculations show a symmetry-forbidden channel opens under tension, peaking near +1% strain.","key_machinery":"The central object is the exciton-phonon coupling matrix element $G^{\\beta\\lambda,\\nu}(\\mathbf{Q},\\mathbf{q})$ of Eq. (2), which the paper constructs as a quantum superposition of electron and hole scattering events, each weighted by the exciton wavefunction coefficients from the finite-momentum Bethe-Salpeter equation. Symmetry analysis is the engine: along the M-K line the point group reduces to $C_s(m)$, where the $E'$ exciton of $C_{3h}$ at K becomes $A'$; only the phonon modes that are even under the mirror plane ($A'$) can couple, while the $A''$ modes are forbidden, explaining the absence of phonon-assisted emission in the unstrained crystal. Biaxial strain does not change the symmetry but alters the coupling strengths, and the van Roosbroeck-Shockley relation with phonon sidebands converts those couplings into a predicted PL spectrum. The strain window from -0.5% to +3% is chosen because it is experimentally reachable.","core_discovery":"On the paper's own terms, the central discovery is a symmetry switch: in unstrained 2H-MoS2 the lowest intervalley excitons at K and M are dark, and the A'' phonon modes along the M-K route are forbidden to couple to them, so no phonon-assisted photoluminescence appears at cryogenic temperature. Biaxial tensile strain, while leaving the D3h point group intact, redistributes the exciton-phonon coupling strength so that A' modes become active; in particular an out-of-plane optical mode (mode 8, about 42 meV) and an in-plane acoustic shear mode (mode 3) dominate. The calculation further shows that the corresponding emission is carried mainly by hole-phonon scattering from K to K', producing two phonon replicas about 50 and 70 meV below the bright exciton. The indirect emission is non-monotonic in strain and maximized near +1%.","pith_inferences":["If the strain-activated hole-phonon channel is as strong as predicted, it would compete with spin and valley relaxation channels, so strained MoS2 could also show strain-tunable valley depolarization times; the paper does not discuss this consequence.","The phase-mismatch caveat means the quantitative ordering of modes 3 and 8 and the exact optimal strain are the least secure numbers; a consistent gauge between the phonon and exciton calculations might shift the peak by a few tenths of a percent.","A direct experimental check could be done on a piezo-strained MoS2 monolayer: monitor the PL doublet separation under strain; if the two replicas move with the phonon energies rather than with the exciton, the assignment to modes 3 and 8 is confirmed.","The symmetry argument is generic: any TMD monolayer with a $C_{3h}$ valley and $C_s(m)$ M-K line should show the same A'/A'' dichotomy, so the strain-activation mechanism likely transfers to MoSe2, WS2, and WSe2."],"forward_implications":["Cryogenic PL of unstrained MoS2 should show only the bright intravalley line; under roughly +1% biaxial tension, two additional replicas should appear about 50 and 70 meV below it.","The lower-energy replica is predicted to be dominated by hole-phonon scattering from K to K', which can be tested by comparing spectra under different excitation polarizations or by time-resolved measurements.","The indirect PL intensity is predicted to be non-monotonic in strain, peaking near +1%; stretching beyond that should reduce the phonon sidebands even though the direct line redshifts.","The same symmetry-based mechanism offers a possible explanation for the non-monotonic strain dependence of PL observed in WSe2, where an exciton-phonon channel would complement the defect-resonance picture.","For other Mo-based TMDs, the selection rule says no A''-assisted emission should appear at low temperature; only the even A' modes can be switched on by strain, which narrows the search for brightening dark excitons."],"supporting_citations":[{"why":"Provides the concurrent exciton linewidth benchmark for monolayer MoS2 that the present scattering rates are compared against.","marker":"[60]"},{"why":"Supplies the phonon-assisted PL formalism and the warning that phase mismatch between codes can reverse intensity trends, which the paper's limitation discussion builds on.","marker":"[56]"},{"why":"Gives the exciton-phonon coupling matrix-element formulation (quantum superposition of electron and hole scattering) used in Eq. (2).","marker":"[59]"},{"why":"Sets up the exciton-phonon coupling treatment of indirect luminescence in h-BN that the present method adapts to MoS2.","marker":"[50]"},{"why":"Establishes the comparison baseline: phonon-assisted PL in Mo-based TMDs is absent in the unstrained case, motivating the strain activation study.","marker":"[31]"},{"why":"Provides the experimental cryogenic PL and dark-bright splittings in MoS2 that the unstrained direct PL calculation is validated against.","marker":"[27]"},{"why":"Supplies the finite-momentum exciton band structure of monolayer MoS2 used to validate the present exciton dispersion.","marker":"[64]"},{"why":"Is the WSe2 strain experiment whose non-monotonic PL variation the paper offers an alternative explanation for.","marker":"[41]"}],"fun_headline_variants":["Strain flips symmetry switch, lighting up dark MoS2 excitons","Tension activates forbidden phonon route for MoS2 dark excitons","Symmetry-forbidden emission turned on by strain in MoS2","Strain unblocks dark-exciton luminescence in MoS2"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the numerical phase mismatch between the phonon and exciton wavefunctions, which the authors explicitly say exists with current codes, does not change the strain-dependent ordering and sign of the A' coupling strengths that produces the +1% peak.","fun_headline_variants_meta":{"raw":{"variants":["Strain flips symmetry switch, lighting up dark MoS2 excitons","Tension activates forbidden phonon route for MoS2 dark excitons","Symmetry-forbidden emission turned on by strain in MoS2","Strain unblocks dark-exciton luminescence in MoS2"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000469,"raw_usage":{"total_tokens":2349,"prompt_tokens":974,"completion_tokens":1375,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":590,"completion_tokens_details":{"reasoning_tokens":1296}},"tokens_in":590,"tokens_out":1375,"duration_ms":10144,"temperature":1.0,"reasoning_tokens":1296,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T22:28:25.757610+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the cryogenic photoluminescence of a MoS2 monolayer under controlled biaxial tensile strain from 0 to +3% with a piezo stage. The central claim is falsified if no doublet of phonon sidebands appears near 50 and 70 meV below the bright exciton as strain approaches +1%, or if the indirect emission does not peak around +1% but instead rises monotonically or stays absent across the whole window.","supporting_citations":[{"cited_title":"Paleari, First-principles approaches to the description of indirect absorption and luminescence spectroscopy: exciton- 10 phonon coupling in hexagonal boron nitride, Ph.D","cited_arxiv_id":null,"evidence_quote":"Provides the concurrent exciton linewidth benchmark for monolayer MoS2 that the present scattering rates are compared against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the phonon-assisted PL formalism and the warning that phase mismatch between codes can reverse intensity trends, which the paper's limitation discussion builds on."},{"cited_title":"Marini, M","cited_arxiv_id":null,"evidence_quote":"Gives the exciton-phonon coupling matrix-element formulation (quantum superposition of electron and hole scattering) used in Eq. (2)."},{"cited_title":"Cannuccia and A","cited_arxiv_id":null,"evidence_quote":"Sets up the exciton-phonon coupling treatment of indirect luminescence in h-BN that the present method adapts to MoS2."},{"cited_title":"Yu and M","cited_arxiv_id":null,"evidence_quote":"Establishes the comparison baseline: phonon-assisted PL in Mo-based TMDs is absent in the unstrained case, motivating the strain activation study."},{"cited_title":"Zhang, Y","cited_arxiv_id":null,"evidence_quote":"Provides the experimental cryogenic PL and dark-bright splittings in MoS2 that the unstrained direct PL calculation is validated against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the finite-momentum exciton band structure of monolayer MoS2 used to validate the present exciton dispersion."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Is the WSe2 strain experiment whose non-monotonic PL variation the paper offers an alternative explanation for."}],"review_version":1}