{"id":"fc3a049c-ca8f-465f-8f6d-3cfa31a38d44","arxiv_id":"2607.19003","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"DFT simulations of oxygen-vacancy trap levels under stress show mechanoluminescence arises from piezoelectric fields or from stress-induced atomic rearrangement around defects, with crystal dimensionality deciding which response dominates.","lead":"This paper uses quantum mechanical simulations to show how mechanoluminescent crystals respond differently to squeezing versus shearing, depending on crystal shape and defect type. It offers a design rule for brighter stress-sensing materials by matching host symmetry, dimensionality, and stress-sensitive defects.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Shear protocol and missing polarization check leave the structural-reorganization / dimensionality mechanism under-supported; a stress-controlled shear calculation with Berry-phase polarization would settle it.","rationale":"I read the paper as making a mechanistic, design-oriented claim: stress-induced trap-depth changes in mechanoluminescent crystals come from either a piezoelectric field or defect-site structural reorganization, and dimensionality determines which response dominates. The strongest evidence is the Ba4Si6O16 calculation showing a correlation between Si–Si reconstruction distance and trap depth under stress, supported by the negative-U hole-trap mechanism and by experimental consistency on PersL trap depths. The weakest load-bearing link is the shear simulation. The paper uses a single Δb displacement instead of a full shear stress tensor, borrows the shear modulus from the glass-ceramic, and asserts the absence of a piezoelectric field without presenting computed polarizations. Each of these is a surrogate for a mechanical state that the central comparison must faithfully represent. The reader's weakest assumption identifies exactly this cluster of issues, and I agree with that assessment. I do not think the concern rises to rejection: the qualitative mechanism and the dimensionality argument are plausible, and the hydrostatic results are better grounded. But the conditional status is appropriate because the shear leg of the comparison is not yet directly supported. No code or data is deposited, so independent verification is not currently possible, but the requested check is well defined and would resolve the concern.","tokens_in":20813,"tokens_out":4352,"duration_ms":40109,"concrete_test":"Recompute the shear response of Ba4Si6O16 for all symmetry-distinct shear modes as explicit affine strain states (γ ≈ 0.01, 0.03, 0.06) with full cell-shape relaxation; for each strain, evaluate both (a) the ε(2+/0) trap depths of VO3 and VO4 via Eq. (2) and (b) the Berry-phase polarization difference between the pristine and defect-containing supercells. If any realistic shear strain lowers a VO3/VO4 trap depth by ≳0.1 eV—the scale of the hydrostatic shifts in §3.4—or yields a polarization above ~0.1 µC/cm², the conclusion that shear is ineffective in quasi-1D hosts and the structural-reorganization-only interpretation would need to be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—two contributions, piezoelectric vs. structural reorganization, with mechanoluminescence intensity set by the dominant term—rests on the negative result that shear leaves trap depths in the quasi-1D host Ba4Si6O16 essentially unchanged. Section 3.4 implements shear as a single plane displacement Δb along b (the g23 mode), converts it to stress using the glass-ceramic shear modulus (33 ± 2 GPa) because the crystal modulus is unknown, and states that defect-containing Ba4Si6O16 supercells develop no piezoelectric field without reporting computed polarizations. The trap-depth response in Figure 7(b,c) is therefore a property of one imposed kinematic mode, not of the shear stress tensor delivered in a torsion or compression experiment. If another symmetry-allowed shear component couples to the VO3/VO4 environment, or if a nonzero local polarization develops under the imposed displacement, the conclusion that quasi-1D topology leaves traps shear-insensitive—and the dimensionality design rule built on it—would fail. This is load-bearing because the comparative statements about SrAl2O4 vs. Ba4Si6O16 and the separation of the two mechanisms all depend on this shear response. The manuscript itself flags the missing polarization support by asserting it is unambiguous without showing the computed values.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies mechanoluminescence mechanisms in SrAl2O4:Eu2+,Dy3+ and Ba4Si6O16:Eu2+,Ho3+ using DFT(+U) defect calculations. It computes formation energies and charge-transition levels for oxygen vacancies, identifies bridging-oxygen vacancies in Ba4Si6O16 as negative-U hole traps stabilized by Si–Si bond formation, and evaluates how hydrostatic and shear-type deformations shift the trap depths. The central conclusion is that stress-induced trap-depth changes arise from two distinct contributions: a piezoelectric field in non-centrosymmetric hosts such as SrAl2O4, and local structural reorganization around point defects, which is claimed to be the dominant mechanism in the centrosymmetric quasi-1D host Ba4Si6O16. This leads to a design rule in which crystal dimensionality, symmetry, and stress-sensitive point defects are matched to optimize mechanoluminescence.","tokens_in":21104,"tokens_out":7473,"duration_ms":70645,"significance":"If the central claim is correct, the paper provides a valuable mechanistic framework by separating piezoelectric and structural-reorganization contributions and translating them into a falsifiable dimensionality/symmetry design rule. The computational core is standard and is cross-validated where it matters: absorption energies agree with experiment to within 1%, the computed trap-depth distribution matches TSL data, and the Si–Si rebonding picture is supported by ICOBI analysis. Importantly, no mechanoluminescence intensity was used to fit the computed quantities, so the comparison with prior phenomenological models is not circular. The main risk is the under-supported shear representation and the unquantified polarization assertion, which are addressable with additional calculations rather than invalidating the derivations.","major_comments":[{"comment":"The central negative result that shear leaves trap depths in quasi-1D Ba4Si6O16 essentially unchanged is computed with a single kinematic mode: a displacement Δb of the (a,b) plane along b (the g23 distortion), converted to a stress using the shear modulus of the glass-ceramic (33 ± 2 GPa) rather than a measured or calculated crystal modulus. This is one of the six shear components, and the manuscript does not show that the other symmetry-allowed shear modes, or a general shear stress tensor, couple to the VO3/VO4 environment in the same way. Because the dimensionality design rule and the separation of the two mechanisms depend on this negative result, the authors should either perform stress-controlled shear calculations spanning the independent shear components or explicitly justify why g23 is the only relevant mode under torsion and compression. As written, Figure 7(b,c) demonstrates the response to one imposed displacement, not to shear stress generally.","section":"§3.4, Fig. 7(b,c)"},{"comment":"The manuscript states that defect-containing Ba4Si6O16 supercells develop no piezoelectric field under hydrostatic or shear stress (“our calculations show unambiguously”), but no computed polarization values are reported. This matters because the piezoelectric mechanism is one of the two contributions and because the authors themselves note that local symmetry breaking around defects can produce a local piezoelectric response even in centrosymmetric hosts. A report of Berry-phase (or equivalent) polarization for defect-free and defect-containing supercells as a function of stress, or at minimum a quantitative bound on any induced polarization, is needed to support the claim and to make the proposed mechanism separation falsifiable. Without it, the comparison with SrAl2O4 rests on an unverified premise.","section":"§3.4, piezoelectric contribution"},{"comment":"The cDFT excited-state results show emission-energy deviations up to 27% for Ba4Si6O16 and a strong dependence on the choice of target Kohn-Sham state (LUMO vs LUMO+1), with only averaging bringing agreement to within 7%. The authors report this transparently, but the text should not present the optical-transition calculations as quantitative validation beyond what Table 2 supports. More importantly, the HSE functional, which produces substantially different charge-transition levels and an overestimated gap, was abandoned after convergence difficulties. Since the stress response is computed at the PBE(+U) level, the authors should provide at least one cross-check, for example a single defect under one stress state with PBEsol or HSE, showing that the stress-induced trap-depth shifts, not just the absolute levels, are robust to the functional choice.","section":"§3.3, Table 2"}],"minor_comments":[{"comment":"Several cross-references appear as “Erreur ! Source du renvoi introuvable” (e.g., §2.4 shear schematic and §3.2 structural reorganization), and Figure 1 has duplicate panel labels; these need to be fixed in the final version.","section":"Throughout"},{"comment":"There is a typographical error in “These depths are markedly larder than the ~0.5–0.7 eV energy window”—“larder” should be “larger”—and similar typos (“consumming”, “time-consumming”) appear in the Introduction.","section":"§3.2"},{"comment":"The P/M group labels introduced when the 21 screw axis is broken by the g23 displacement should be defined at first use in the text and in the Figure 7 caption; as written, the reader must infer which sites belong to each group from the symbol colors.","section":"§2.4 / Fig. 7"},{"comment":"The conclusions describe the host as “quasi-1D silicate-chain topology” while earlier text calls it “1D”; one consistent terminology should be used, or the distinction defined, since the dimensionality argument is central to the design rule.","section":"§4 / Conclusions"},{"comment":"The data availability statement says data are available from the corresponding author upon reasonable request; providing a repository with input structures, pseudopotential settings, and defect configurations would materially improve reproducibility.","section":"Data availability"}],"recommendation":"major_revision","confidential_remarks":"The paper is part of a continuing series that leans heavily on the authors' own experimental papers (refs. 8, 28, 29) for the torsion-compression phenomenology. I do not see a circularity problem—trap depths are compared with external TSL data and prior models, and no mechanoluminescence intensity is used to fit computed quantities—but independent replication of the no-torsion-ML result for Ba4Si6O16 would materially strengthen the dimensionality claim. The shear-modulus assumption and the unverified polarization assertion are the main technical risks, and both are addressable within the scope of a major revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThe paper is worth your attention: it gives a concrete mechanism for mechanoluminescence in Ba4Si6O16:Eu,Ho via stress-induced structural reorganization around oxygen-vacancy hole traps, and contrasts it with the piezoelectric mechanism in SrAl2O4:Eu,Dy. The best part is the defect analysis—identifying bridging-oxygen vacancies as negative-U hole traps with Si–Si bond formation, and matching the TSL trap-depth distribution. That is a solid, citable result.\n\nWhat is new is the calculation of stress-dependent charge-transition levels for eight oxygen vacancies under hydrostatic and shear perturbations. Hydrostatic stress systematically shallows traps, consistent with experiment and previous phenomenological models. The correlation between the Si–Si distance and the trap depth is a satisfying mechanistic link.\n\nThe soft spots are in the shear treatment. Shear is modeled as a single plane displacement Δb along the b axis, converted to stress using the glass-ceramic shear modulus because the crystal modulus is unknown. Only one of six possible shear modes is considered. The paper states that defect-containing Ba4Si6O16 supercells develop no piezoelectric field under stress, but no computed polarizations are shown. For a centrosymmetric host that assertion is plausible, but it is load-bearing: the dimensionality design rule—quasi-1D hosts are shear-insensitive—depends on it. If another shear component couples to the vacancy environment, or if a local polarization appears, the conclusion weakens.\n\nAlso, HSE calculations were abandoned after convergence issues, so the PBE results carry that uncertainty, and the emission energies deviate by up to 27% (absorption within 1%, which is good). The 'universal' structural-reorganization mechanism is extrapolated from two systems. And no data or code are deposited, which makes independent verification harder.\n\nAll that said, the central mechanism for Ba4Si6O16 under hydrostatic stress is well-supported. The shear conclusion needs more work. I would send this to peer review, with a request for a proper shear stress tensor calculation, at least the dominant shear mode, and explicit polarization values for the defect supercells. If those hold, this becomes a genuinely useful design rule.\n\nWho is this for? People working on mechanoluminescence, persistent phosphors, and defect engineering in silicates. A serious referee should engage with it.","headline":"A useful defect-chemistry study with a plausible two-mechanism picture, but the shear protocol and unverified no-polarization claim leave the dimensionality design rule under-supported.","tokens_in":21588,"tokens_out":2670,"would_cite":true,"duration_ms":24610,"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":"Mechanoluminescence intensity is set by two independent stress-driven effects—a piezoelectric field in non-centrosymmetric hosts and stress-induced structural rearrangement at the trapping defect—with the host's structural dimensionality…","keywords":["mechanoluminescence","persistent luminescence","trap depth","oxygen vacancy","charge-transition level","density functional theory","structural dimensionality","piezoelectric field"],"falsifier":"Measure the stress-induced electric polarization of defect-containing Ba4Si6O16 supercells under shear; if a net polarization appears, the claim that shear acts only through local reorganization in this material fails. Alternatively, run a torsion test on a centrosymmetric, three-dimensionally connected phosphor with stress-sensitive traps; emission there would contradict the dimensionality rule.","tokens_in":20657,"feed_emoji":"💡","tokens_out":7083,"duration_ms":63968,"temperature":0.7,"pith_summary":"Mechanoluminescent crystals emit light when mechanical stress releases charge carriers trapped at point defects, but it has been unclear why some crystals respond to pressure and others to twisting. This paper uses first-principles calculations on two representative phosphors to argue that two distinct stress-driven effects change the depth of the carrier traps: a piezoelectric field, active in non-centrosymmetric hosts like SrAl2O4, and a structural reorganization around the defect site, active in any host. In the quasi-1D host Ba4Si6O16, hydrostatic compression shallows every oxygen-vacancy hole trap, while shear leaves most trap depths nearly unchanged—explaining why this material glows under isostatic pressure but not under torsion. If correct, the result turns mechanoluminescence into a design problem of matching host dimensionality and symmetry to stress-sensitive point defects.","feed_headline":"Two stress effects control mechanoluminescence brightness","feed_subtitle":"A DFT study isolates piezoelectric fields and defect-site reorganization as the two levers for designing brighter crystals.","key_machinery":"The load-bearing object is the thermodynamic charge-transition level ε(q/q′) of a point defect—the Fermi-level position at which the defect changes charge state, equivalently its trap depth—computed in strained supercells. Around an oxygen vacancy, the key structural coordinate is the Si(1)–Si(2) distance across the vacancy, whose contraction to about 2.4 Å marks Si–Si bond formation and whose expansion to about 4.5 Å breaks that bond; the trap depth varies monotonically with this distance. These quantities connect the applied mechanical load to the emission response, and the paper uses this machinery to separate piezoelectric from structural-reorganization contributions.","core_discovery":"The paper's central claim is that two microscopic contributions govern stress-induced changes in trap depth in mechanoluminescent crystals: generation of a piezoelectric field, as in SrAl2O4:Eu2+,Dy3+, and structural reorganization around the carrier-trapping point defects, as in Ba4Si6O16:Eu2+,Ho3+, with mechanoluminescence intensity scaling with the dominant contribution. In Ba4Si6O16, the active traps are oxygen vacancies at bridging oxygen sites of the silicate chains; hole capture at a neutral vacancy drives formation of a Si–Si bond, and the depth of the ε(2+/0) charge-transition level tracks the Si–Si distance. Under hydrostatic stress every trap depth decreases, releasing holes and producing light; under shear, the quasi-1D chains glide as relatively rigid units, only weakly perturbing the vacancy environment, so the luminescence response is weak or absent. In a 3D non-centrosymmetric host, shear produces a piezoelectric polarization that bends band edges and lowers trap barriers, giving a second route to bright emission.","pith_inferences":["If the structural-reorganization mechanism is universal, then engineering the local stiffness around a trap—through chemical pressure, ligand chemistry, or strain—could tune mechanoluminescence even in non-piezoelectric hosts, a route the paper sketches but does not quantify.","The monotonic trap-depth-versus-Si–Si-distance correlation suggests a testable scaling law: hosts whose vacancy reconstruction has a large configuration-coordinate change should show larger hydrostatic mechanoluminescence responses.","The same negative-U hole-trap picture may transfer to other corner-sharing silicate or phosphate phosphors, and to glass-ceramics, where local rebonding at vacancies could supply mechanoluminescence without long-range order.","A direct calculation of the stress-induced polarization in defect-containing supercells would sharpen the boundary between the piezoelectric and reorganization regimes; the paper asserts the absence of the former in Ba4Si6O16 but does not report computed polarizations."],"forward_implications":["Hydrostatic compression is the dominant mechanical driving force for mechanoluminescence in quasi-1D hosts: it shallows every relevant oxygen-vacancy trap, consistent with observed emission under isostatic pressure.","Shear stress is ineffective in quasi-1D hosts like Ba4Si6O16 because silicate chains glide as rigid units, leaving vacancy environments mostly unchanged; this explains the absence of mechanoluminescence under torsion.","In non-centrosymmetric 3D hosts, shear can generate piezoelectric fields that lower trap barriers, so torsional mechanoluminescence is expected there.","The structural-reorganization contribution is apparently universal and may operate in centrosymmetric crystals and disordered systems, since it does not require bulk piezoelectricity.","Bright mechanoluminescent materials can be designed by combining a host that supports piezoelectric coupling with point defects whose trap environment is sensitive to the applied stress mode."],"supporting_citations":[{"why":"Supplies the first-principles method for defect formation energies and charge-transition levels used throughout the paper.","marker":"20"},{"why":"Provides the prior defect-identification study of SrAl2O4:Eu2+ that the paper builds on for the piezoelectric comparison.","marker":"21"},{"why":"Advances the piezoelectric trap-depth-reduction model that the paper extends and contrasts with structural reorganization.","marker":"23"},{"why":"Reports the experimental hydrostatic mechanoluminescence response of SrAl2O4:Eu2+,Dy3+ that the calculations are designed to explain.","marker":"24"},{"why":"Provides experimental trap-depth distributions and mechanoluminescence behavior of Ba4Si6O16:Eu2+,RE phosphors.","marker":"28"},{"why":"Documents mechanoluminescence of the Ba4Si6O16:Eu2+,Ho3+-containing glass-ceramic under different loading modes, the key experimental counterpart.","marker":"29"},{"why":"Supplies the crystal-orbital bond index used to quantify Si–Si bond formation across the oxygen vacancy.","marker":"39"},{"why":"Documents analogous hole trapping and Si–Si rebonding at oxygen vacancies in amorphous SiO2, supporting the universality of the reorganization mechanism.","marker":"50"}],"fun_headline_variants":["Two mechanisms tune mechanoluminescence","Stress type and structure set emission mechanism","Piezoelectric and defect relaxation drive glow","Hydrostatic vs shear: how traps release light","Dimensionality governs stress-induced luminescence"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole design rule rests on the assumption that a single rigid in-plane displacement, with stress values estimated from the glass-ceramic's shear modulus, faithfully represents how real shear loads change trap depths in the crystal.","fun_headline_variants_meta":{"raw":{"variants":["Two mechanisms tune mechanoluminescence","Stress type and structure set emission mechanism","Piezoelectric and defect relaxation drive glow","Hydrostatic vs shear: how traps release light","Dimensionality governs stress-induced luminescence"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00016,"raw_usage":{"total_tokens":1227,"prompt_tokens":936,"completion_tokens":291,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":552,"completion_tokens_details":{"reasoning_tokens":227}},"tokens_in":552,"tokens_out":291,"duration_ms":4574,"temperature":1.0,"reasoning_tokens":227,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T15:32:57.709979+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the stress-induced electric polarization of defect-containing Ba4Si6O16 supercells under shear; if a net polarization appears, the claim that shear acts only through local reorganization in this material fails. Alternatively, run a torsion test on a centrosymmetric, three-dimensionally connected phosphor with stress-sensitive traps; emission there would contradict the dimensionality rule.","supporting_citations":[],"review_version":1}