{"id":"2711f6e9-6839-43e6-9574-f780cef7a863","arxiv_id":"2506.05809","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Angle-resolved Raman spectroscopy and density functional perturbation theory trace unusual four-fold Ag-mode polarization patterns in Ta2Ni3Te5 to anisotropic electron-phonon interactions.","lead":"This paper reports that Raman light-scattering patterns from the layered crystal Ta2Ni3Te5 show unexpected four-lobe shapes in vibration modes that normally have two lobes, and attributes the effect to direction-dependent electron-phonon interactions. If correct, it adds a new material platform for studying anisotropic quantum behavior in low-dimensional crystals.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central claim rests on monolayer A1 DFPT tensors substituting for measured few-layer bulk-like Ag modes; identical tensor form fixes allowed components, not the complex phases that generate the four-fold pattern.","rationale":"The reader's weakest-assumption analysis identifies exactly the load-bearing point: the measured few-layer Ag modes are bulk-like (9-12 layers, D2h), while the DFPT Raman tensors are computed for a monolayer (C2v). I agree that the shared Raman tensor form is necessary but not sufficient. Equation 5 shows that the four-fold angular pattern is produced by the phase δ between the complex b and c tensor elements; group theory only tells you which elements are nonzero, not their relative phase. In Equation 13 those phases are determined by sums over intermediate states with energy denominators; monolayer and bulk have different symmetries and band structures (the monolayer Γ-point gap is ~50 meV), so the computed complex tensors need not transfer. The paper's own admission that bulk calculations exceed current resources (Section 2.3) makes this a real gap rather than a stylistic choice. The structural characterization, group-theory analysis, and use of an open-source code (QERaman) are independent strengths, and the temperature-dependent Raman analysis is internally consistent; nothing in my reading questions those parts. The concern is falsifiable: a bulk DFPT Raman calculation would settle whether the monolayer result is representative. Because the reader already returned CONDITIONAL on essentially this basis, my read does not move the verdict.","tokens_in":34778,"tokens_out":9681,"duration_ms":103309,"concrete_test":"Run the same QERaman DFPT Raman-tensor calculation for bulk Ta2Ni3Te5 (D2h, 40-atom conventional cell) at the same 532 nm excitation, with a k-mesh comparable to the monolayer calculation (e.g., 1×8×4 or 1×12×6), and compute the Ag-mode polar plots for the modes corresponding to the measured peaks. Overlay calculated and measured 9L/12L polar plots and compare the fitted b/c amplitude ratio and phase δ from Equation 5. If bulk DFPT reproduces the measured four-fold patterns, the monolayer-to-bulk transferability assumption is validated; if the bulk polar plots revert to two-fold or give substantially different phases, the central claim is unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 2.3 states that bulk DFPT 'exceeds our current resources' and therefore all Raman intensity calculations are for a monolayer (C2v), justified by the statement that 'the monolayer A1 mode retains the same Raman tensor form as the bulk Ag mode.' That group-theoretic statement fixes the allowed nonzero tensor components, but it does not fix the magnitudes and phases of those components. The four-fold Ag response in Equation 5 is controlled by the phase difference δ between the b and c elements; in the quantum formula (Equation 13) these phases arise from sums over intermediate states with complex resonant denominators. Monolayer Ta2Ni3Te5 has C2v symmetry and a Γ-point gap of only ~50 meV (Figure S16a), while the measured 9-12L flakes have D2h symmetry and a different band structure, so the B2u and B1u intermediate-state spectra can differ substantially. The paper provides no bulk DFPT Raman tensor calculation and no layer-dependent validation (e.g., 1L/3L vs 9L/bulk). Figure 4d-f is a qualitative shape comparison, not a quantitative fit to the experimental polar plots. The central attribution to anisotropic electron-phonon interactions is therefore plausible but not quantitatively established; the four-fold pattern could be a monolayer-specific or electron-photon effect whose separation from EPI anisotropy is asserted rather than demonstrated.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports angle-resolved polarized Raman spectroscopy, temperature-dependent Raman measurements, and DFPT-based Raman tensor calculations for few-layer Ta2Ni3Te5, aiming to explain an unusual four-fold angular response of Ag phonon modes. The authors attribute this response to anisotropic electron-photon and electron-phonon interactions, supported by QERaman calculations of complex Raman tensors for a monolayer model. They also analyze temperature-dependent peak shifts with three- and four-phonon decay models and report dominant four-phonon contributions for some modes.","tokens_in":35055,"tokens_out":5875,"duration_ms":62947,"significance":"If the central attribution is quantitatively established, the paper would offer a useful demonstration of combining angle-resolved polarized Raman spectroscopy with first-principles quantum Raman calculations to probe anisotropic electron-phonon interactions in low-symmetry 2D materials. The experimental data set is extensive, including structural imaging, polar plots for many modes, and temperature-dependent measurements. A notable strength is that the QERaman tensors are computed ab initio rather than fitted to the polar data, and the derivation showing the identical angular forms of absorption and birefringence in the parallel configuration is a useful contribution. However, the load-bearing connection between the monolayer DFPT result and the measured few-layer flakes is asserted rather than quantitatively demonstrated, and the comparison between calculated and measured polar plots is qualitative.","major_comments":[{"comment":"The step from monolayer DFPT to few-layer experiment is not quantitatively justified. The manuscript states in Section 2.3 that bulk DFPT \"exceeds our current resources\" and therefore uses a monolayer C2v calculation, justified by the statement that \"the monolayer A1 mode retains the same Raman tensor form as the bulk Ag mode.\" That group-theoretic statement fixes the allowed nonzero components but not their magnitudes and phases. Equation 5 shows that the four-fold contribution in the Ag response is controlled by the phase difference between the b and c tensor elements, while Equation 13 shows that these phases arise from sums over intermediate states with complex resonant denominators. The monolayer has C2v symmetry and a ~50 meV gap at the Gamma point (Figure S16a), whereas the measured 9L/12L flakes have D2h symmetry and a different band structure, so the intermediate-state spectra can differ substantially. The paper provides no bulk DFPT Raman tensor and no layer-dependent validation such as comparing 1L, 3L, and 9L polar plots. I recommend either providing bulk or few-layer QERaman tensors or quantitatively testing the monolayer tensors against measured few-layer polar data.","section":"Section 2.3 and Supporting Section S5"},{"comment":"The comparison between calculated and measured polar plots is qualitative only. The text says the DFPT intensities \"successfully reproduce\" the complex angular dependence, but no numerical metric is reported, such as fitted phase differences, amplitude ratios, or a chi-squared/overlap measure between the calculated and experimental polar plots. Since the experimental fits in Figure S9 already assign per-mode phase differences and amplitude ratios, these fitted parameters should be tabulated and compared with the corresponding values extracted from the calculated complex tensors. Without such a quantitative comparison, the central claim that the four-fold Ag response originates from the computed anisotropic electron-phonon interactions remains plausible but not established.","section":"Figure 4d-f and Figure S9"},{"comment":"The inference that a dominant four-phonon process in several modes is \"a possible manifestation of strong anisotropic electron-phonon interactions\" is not supported by the presented analysis. The fitted A and B coefficients in Equation 16 are empirical anharmonic parameters; no calculation or model connects their relative magnitude to electron-phonon matrix elements or to their anisotropy. This statement appears in the abstract and conclusion as a substantive finding, but the paper offers no independent evidence for that link. I suggest either adding a microscopic calculation of anharmonic phonon decay weighted by electron-phonon coupling or explicitly labeling this as a speculative remark.","section":"Section 2.4 and Table S13"}],"minor_comments":[{"comment":"Several equations in Supporting Section S3 contain corrupted or unreadable mathematical symbols, particularly Equations S2, S3, S8, and S9. These should be regenerated so that the derivations are fully legible.","section":"Supporting Information S3"},{"comment":"In the crystal growth description, the text lists \"selenium (Alfa Aesar, 99.99% purity)\" as a starting material, but Ta2Ni3Te5 contains tellurium. This is presumably a typo and should be corrected.","section":"Section 4, Experimental Section"},{"comment":"The caption of Figure 3 states the flake is 12L, while the text in Section 2.2 describes initial measurements on a 9L flake and complementary measurements on a 12L flake. The manuscript should clarify which flake is used for the polar plots in Figures 3 and S9-S10.","section":"Figure 3 caption and Section 2.2"},{"comment":"The statement that \"the four-phonon process is not only significant but even dominant\" for several modes is somewhat stronger than the table indicates: only one listed mode has B/A clearly greater than unity, and several modes have B/A below 0.1. The wording should be aligned with the fitted values.","section":"Table S13"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The core experiment is solid and the DFPT work is real: the paper reports a careful growth and characterization of Ta2Ni3Te5, a full angle-resolved polarized Raman dataset, and a genuine ab initio calculation of complex Raman tensors that is not fitted to the polar data. That is more than many Raman papers do. The generalization of Kranert's equivalence between absorption and birefringence in the parallel configuration is a useful, clean observation and appears to be new. The SI tables of computed complex tensors for all A1 modes are a good faith effort at reproducibility.\n\nThe soft spot is exactly where the stress-test note lands. The central claim—that the four-fold Ag patterns come from anisotropic electron-phonon interactions—rests on substituting monolayer C2v DFPT tensors for the measured few-layer D2h flakes. The paper's justification, that the monolayer A1 mode retains the same Raman tensor form as the bulk Ag mode, fixes the allowed nonzero components but not their magnitudes or phases. Those phases are what produce the four-fold pattern in Equation 5, and they depend on intermediate-state spectra that a 50 meV-gap monolayer and a 9-12 layer flake are not guaranteed to share. The comparison in Figure 4d-f is qualitative; there is no numerical measure of how well the computed and measured polar plots agree. The authors say bulk DFPT exceeds their resources, which is honest, but that does not make the substitution valid. This is a load-bearing gap, though not a fatal one: the measurements stand, and the attribution is plausible.\n\nMinor points: the temperature-dependent four-phonon analysis fits the data, but the link between four-phonon dominance and EPI anisotropy is speculative and could be flagged as such. No data or code are provided, and the experimental polar plots have no error bars. The absence of error bars matters more for the fitted phase differences than for the raw intensities.\n\nOverall, this is a worthwhile paper for the subfield: a new anisotropic 2D material with clear quasi-1D structure, a thorough Raman investigation, and an independent DFPT prediction that is at least presented in full. It deserves peer review, but the referees should push for either bulk or few-layer DFPT validation, or a quantitative comparison that does not rely on the asserted tensor-form equivalence.","headline":"A careful Raman study with a real DFPT calculation, but the monolayer-for-few-layer substitution leaves the central EPI-anisotropy claim plausible rather than proven.","tokens_in":35618,"tokens_out":1114,"would_cite":false,"duration_ms":14839,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["78.30.-j","71.38.-k"],"model":"deepseek-v4-flash","headline":"Anisotropic electron-phonon coupling explains the four-fold Raman patterns measured in Ta2Ni3Te5 flakes.","keywords":["anisotropic electron-phonon interaction","Ta2Ni3Te5","quasi-one-dimensional","angle-resolved polarized Raman spectroscopy","complex Raman tensor","density functional perturbation theory","phonon decay","four-phonon process"],"falsifier":"Measure angle-resolved polarized Raman spectra of the same Ag modes on monolayer and on flakes of several thicknesses: if the four-fold pattern is intrinsic anisotropic electron-phonon coupling, it should persist in the monolayer with the calculated C2v complex tensors; if birefringence were responsible, the fitted phase difference should grow linearly with thickness (Eq. 9) and vanish for the thinnest flakes.","tokens_in":34565,"feed_emoji":"🔬","tokens_out":5132,"duration_ms":47981,"temperature":0.7,"pith_summary":"This paper tries to establish that the unusual four-fold angular dependence of the Ag Raman modes in few-layer Ta2Ni3Te5 is not an artefact of classical Raman tensor fitting or sample birefringence but a fingerprint of intrinsically anisotropic electron-phonon interactions. The claim matters because low-symmetry quasi-one-dimensional materials are candidate platforms for excitonic, topological, and superconducting phases whose properties are set by how strongly phonons couple to electrons along different crystal directions, so a direct optical probe of that anisotropy would give experimenters a way to test those phases. The evidence combines angle-resolved polarized Raman spectra on 9- and 12-layer flakes with full quantum perturbation theory and density functional perturbation theory calculations that reproduce the measured polar plots only when anisotropic electron-phonon coupling is included. The paper further argues that temperature-dependent Raman shifts in several modes require four-phonon decay processes, and connects that higher-order anharmonicity to the same anisotropic coupling.","feed_headline":"Four-fold Raman pattern traced to anisotropic electron-phonon coupling","feed_subtitle":"Quantum DFPT calculations reproduce the unusual Ag-mode polar plots in quasi-1D Ta2Ni3Te5 flakes.","key_machinery":"The load-bearing object is the complex Raman tensor, whose off-diagonal phase differences control the Ag angular pattern. It is derived from a third-order perturbation-theory expression for Stokes Raman intensity (Eq. 13) in which the intermediate-state sum is constrained by symmetry-allowed optical transitions (dipole selection rules) and by the electron-phonon matrix element for emitting the phonon. The paper computes these complex tensors with density functional perturbation theory for the monolayer (point group C2v), relying on the statement that the monolayer A1 mode retains the same Raman tensor form as the bulk Ag mode of the measured few-layer flakes; those calculated tensors, with anisotropic electron-phonon coupling included, generate polar plots matching the measured four-fold Ag patterns.","core_discovery":"The central discovery is that the complex Raman tensor needed to fit the anomalous Ag-mode polar plots in Ta2Ni3Te5 has a microscopic origin in the anisotropy of electron-phonon interactions, rather than in classical light absorption phenomenology or linear birefringence. In the full quantum treatment of Stokes scattering, the intensity is built from two electron-photon matrix elements and one electron-phonon matrix element with explicit optical dipole selection rules; when the electron-phonon element is treated as isotropic, the calculation predicts two-fold Ag symmetry, contradicting the measured four-fold patterns, whereas including anisotropic electron-phonon interactions in DFPT reproduces the measured polar plots. The paper also shows that, in parallel polarization, absorptive complex tensors and birefringence produce mathematically identical Ag intensity formulas, and it argues on thickness and mode-by-mode phase-difference grounds that the operative mechanism is the intrinsic anisotropic interaction, not birefringence. Taken at face value, the result makes Ta2Ni3Te5 a concrete system in which anisotropic electron-phonon coupling is directly observable in an optical experiment.","pith_inferences":["If the monolayer-to-bulk Raman tensor correspondence holds, a thickness series of angle-resolved Raman measurements should show the four-fold Ag pattern persisting down to monolayer thickness, whereas a birefringence-dominated interpretation would predict the fitted phase difference to shrink linearly with thickness (Eq. 9).","The dominance of four-phonon decay in specific modes could be tested independently by measuring mode-resolved phonon lifetimes, for example with coherent phonon spectroscopy, and comparing them with the anharmonic constants extracted from the temperature fits.","Because the anisotropy is electronic in origin, electrostatic gating or doping should modulate the complex Raman tensors and thereby reshape the Ag polar plots, offering a tunable optical probe of the same electron-phonon coupling."],"forward_implications":["Angle-resolved polarized Raman spectroscopy combined with full quantum DFPT becomes a workflow for mapping anisotropic electron-phonon coupling in low-symmetry layered materials, not just in Ta2Ni3Te5.","The four-fold Ag patterns can serve as a symmetry-resolved fingerprint for identifying chain orientation in exfoliated flakes, since the lobe geometry is tied to the b-axis chain direction.","For several phonon modes the four-phonon decay channel dominates over the three-phonon channel, so models of phonon lifetimes and heat dissipation in this material must include quartic anharmonicity.","Because absorption and birefringence enter the same Ag intensity formula, experiments on other orthorhombic two-dimensional materials must separate the two before assigning unusual polar patterns to intrinsic anisotropic electron-phonon coupling."],"supporting_citations":[{"why":"Supplies the density-functional perturbation-theory framework used to compute resonance Raman intensities with explicit electron-phonon matrix elements.","marker":"[43]"},{"why":"Provides the classical Raman tensor formalism and D2h mode symmetries that define the Ag and B3g angular formulas.","marker":"[44]"},{"why":"Documents birefringence-directed Raman selection rules in black phosphorus, the orthorhombic comparison case that this paper extends.","marker":"[46]"},{"why":"Gives the Jones-matrix treatment of Raman scattering in optically anisotropic crystals that the paper adapts for the birefringence channel.","marker":"[47]"},{"why":"Supplies the full-quantum framework for anisotropic electron-photon and electron-phonon interactions in black phosphorus, the methodological template for Eq. 13.","marker":"[48]"},{"why":"Furnishes the three- and four-phonon anharmonic model used to fit the temperature-dependent Raman shifts.","marker":"[56]"},{"why":"Established the precedent of unusual angular dependence of Raman response in an anisotropic layered material, motivating the complex-tensor analysis.","marker":"[28]"}],"fun_headline_variants":["Anisotropic electron-phonon coupling explains four-fold Raman pattern","Ta2Ni3Te5 Raman anisotropy traced to electron-phonon interactions","Four-fold Raman symmetry from anisotropic electron-phonon coupling","Quantum theory ties Raman four-fold pattern to anisotropic EPI","Anisotropic interactions shape Raman polar plots in Ta2Ni3Te5"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculations are done for a monolayer with C2v symmetry, and the paper assumes that its A1 mode has the same Raman tensor form as the bulk Ag mode of the measured 9- and 12-layer flakes, without a bulk calculation or layer-by-layer validation.","fun_headline_variants_meta":{"raw":{"variants":["Anisotropic electron-phonon coupling explains four-fold Raman pattern","Ta2Ni3Te5 Raman anisotropy traced to electron-phonon interactions","Four-fold Raman symmetry from anisotropic electron-phonon coupling","Quantum theory ties Raman four-fold pattern to anisotropic EPI","Anisotropic interactions shape Raman polar plots in Ta2Ni3Te5"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000237,"raw_usage":{"total_tokens":1547,"prompt_tokens":1022,"completion_tokens":525,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":638,"completion_tokens_details":{"reasoning_tokens":434}},"tokens_in":638,"tokens_out":525,"duration_ms":4636,"temperature":1.0,"reasoning_tokens":434,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T10:13:01.367129+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure angle-resolved polarized Raman spectra of the same Ag modes on monolayer and on flakes of several thicknesses: if the four-fold pattern is intrinsic anisotropic electron-phonon coupling, it should persist in the monolayer with the calculated C2v complex tensors; if birefringence were responsible, the fitted phase difference should grow linearly with thickness (Eq. 9) and vanish for the thinnest flakes.","supporting_citations":[{"cited_title":"Nature Communications, 2024","cited_arxiv_id":null,"evidence_quote":"Supplies the density-functional perturbation-theory framework used to compute resonance Raman intensities with explicit electron-phonon matrix elements."},{"cited_title":"npj Computational Materials, 2023","cited_arxiv_id":null,"evidence_quote":"Provides the classical Raman tensor formalism and D2h mode symmetries that define the Ag and B3g angular formulas."},{"cited_title":"Small, 2016","cited_arxiv_id":null,"evidence_quote":"Supplies the full-quantum framework for anisotropic electron-photon and electron-phonon interactions in black phosphorus, the methodological template for Eq. 13."},{"cited_title":"Carbon, 2019","cited_arxiv_id":null,"evidence_quote":"Furnishes the three- and four-phonon anharmonic model used to fit the temperature-dependent Raman shifts."},{"cited_title":"Nature communications, 2015","cited_arxiv_id":null,"evidence_quote":"Established the precedent of unusual angular dependence of Raman response in an anisotropic layered material, motivating the complex-tensor analysis."}],"review_version":1}