{"id":"dc4f36de-fe8c-4456-b175-ffcc24b3d773","arxiv_id":"2607.21682","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":3.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A symmetry-based analysis of PZT shows that its tetragonal-to-rhombohedral transition appears in powder Raman spectra as continuous intensity redistribution rather than new modes.","lead":"This paper uses standard crystallographic symmetry rules to predict which physical properties and vibrations are allowed in different phases of the piezoelectric ceramic PZT. It then re-fits previously measured Raman spectra to argue that the material's phase change across the morphotropic boundary is a smooth transfer of intensity, not the appearance of new vibration modes.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Fixed 13-component basis makes 'no new modes' a built-in assumption; CMCD check only on x=0.44 leaves central claim of continuous redistribution unsupported across MPB.","rationale":"The reader's weakest_assumption identified the fixed 13-component decomposition as the soft spot, and my analysis agrees: the fitting strategy enforces the conclusion of 'no new modes' by construction. The paper's group-theoretic tensor analysis is standard and correct, but the experimental re-analysis is the part that carries the headline claim. The concern is not that the authors are wrong—their interpretation is consistent with prior literature (Buixaderas et al.) and may well be correct—but that the evidence presented does not uniquely support the claim. The CMCD cross-check at x=0.44 is helpful but does not validate the composition range where the MPB and possible monoclinic phase would matter. Lack of raw data and parameter uncertainties further weaken falsifiability. These issues are addressable with a focused re-analysis, so the appropriate verdict remains CONDITIONAL, not rejection. My read therefore does not change the reader's verdict, and agreement is 'agree'.","tokens_in":23146,"tokens_out":4536,"duration_ms":47280,"concrete_test":"Digitize the published Raman spectra (Fig. 4, from Ref. [13]) and refit each composition independently using a Bayesian mixture model with the number of components a free parameter (e.g., reversible-jump MCMC or BIC over 10–16 components), using the same pseudo-Voigt lineshape and baseline. If the Bayes-optimal number of components is not 13 for all x, or if adding a monoclinic A'/A''-type component near the E(TO2)/E(TO3) region yields a significant improvement (e.g., ΔBIC > 10) for any x≥0.50, the central claim of no emergent modes and continuous redistribution loses support. If 13 components remain Bayes-optimal at every composition, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (Sec. IX) that the tetragonal-to-rhombohedral transition appears as a continuous redistribution among 13 components without new modes rests on the fitting protocol described in Sec. VIII A-B. In that protocol, the number of components is fixed a priori at 13 for every composition, and the mixing parameter η_j is fixed from unconstrained fits to x=0.44. Consequently, the fit cannot detect any new mode: a spectral feature that might signal phase-specific or monoclinic splitting is automatically absorbed into changes in the positions, widths, and amplitudes of the existing 13 components. The reported RMS residuals of 0.53–0.69% are therefore not evidence against missing modes, since broad pseudo-Voigt components with free centers, widths, and amplitudes are flexible enough to accommodate additional spectral weight within these residual levels. The only independent validation of the number of components, CMCD analysis (Fig. S1 and Sec. VIII B), was applied solely to the reference x=0.44 composition; it does not constrain x=0.50–0.58, where the MPB coexistence/monoclinic question is most acute. No uncertainty estimates are given for the fitted parameters, so the 'continuous' evolution of amplitudes and widths could be consistent with noise or with a different decomposition. Thus the absence of new modes and the attribution of subpeaks to local disorder are to a significant degree consequences of the model choice rather than independent empirical findings.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a symmetry-based workflow for ferroelectric perovskites, deriving allowed pyroelectric, dielectric, and piezoelectric tensor forms and Raman-active mode classifications for the cubic, tetragonal, rhombohedral, and monoclinic phases of PZT. It applies this framework to reanalyze previously reported room-temperature powder Raman spectra (x = 0.44–0.58) across the morphotropic phase boundary, using a fixed 13-component pseudo-Voigt decomposition. The central claim is that the tetragonal-to-rhombohedral transition appears as a continuous redistribution of spectral weight among these 13 components, with no new modes emerging and no long-range monoclinic signature; the additional subpeaks are attributed to local disorder and anharmonicity.","tokens_in":23511,"tokens_out":9194,"duration_ms":87532,"significance":"The group-theoretic portion is a clear and internally consistent presentation: the Bhagavantam–Venkatarayudu and correlation methods agree (Eqs. 27–28), and the tensor forms in Table I match standard textbook results. The proposed workflow and the suggestion to track the tetragonal B1 mode in oriented samples are useful methodological contributions. If the experimental conclusion were robust, the paper would provide a valuable reference baseline for interpreting disorder-broadened powder Raman spectra of PZT. However, the central experimental claim is not uniquely supported by the analysis as presented: the fixed 13-component basis, the lack of composition-resolved validation, and the absence of uncertainty estimates mean that the 'no new modes' conclusion is partly a consequence of the model choice rather than an independent empirical finding. The manuscript is therefore of interest, but the experimental section needs substantial strengthening before the central claim can be accepted.","major_comments":[{"comment":"The central conclusion in Sec. IX—that the transition manifests as continuous redistribution 'without the emergence of new modes'—is not supported by the fitting protocol. The number of components is fixed a priori to 13 for every composition, so the model cannot represent a new branch; free centers, widths, and amplitudes of broad pseudo-Voigt peaks can absorb extra spectral weight. The reported RMS residuals of 0.53–0.69% therefore do not rule out additional peaks. The authors should provide a model comparison that allows phase-specific or extra components (e.g., BIC/AIC or an F-test) or apply an independent peak-count method to all compositions.","section":"Sec. VIII A–B; Fig. 5; Eq. (35)"},{"comment":"The CMCD analysis lists 10 'reliable' curvature maxima, 5 'additional' features, and 2 LO3 maxima—17 features in total—yet concludes that the results support 'approximately 13 underlying spectral contributions'. This numerical inconsistency undermines the claimed independent consistency check. Moreover, CMCD is performed only on the x = 0.44 reference; no such validation is provided for x = 0.50–0.58, where the MPB and monoclinic questions are most acute.","section":"Sec. VIII B"},{"comment":"The justification for using a single decomposition basis across the phase diagram is that tetragonal and rhombohedral PZT have 'the same number ... of Raman-active phonon branches'. This is inaccurate: P4mm gives 3A1 + B1 + 4E, all Raman-active, while R3m gives 3A1 + A2 + 4E with A2 silent. Thus the number of Raman-active branches differs (8 vs 7). This weakens the symmetry argument for fixing the same 13-component basis for all compositions up to x = 0.58.","section":"Sec. VI and Sec. VIII B"},{"comment":"No uncertainty estimates are reported for the fitted parameters (ω0j, Γj, Aj). The claims of 'continuous evolution' and of specific amplitude transfers (e.g., E(TO2) from 66 to 20 arb. units) cannot be distinguished from noise or from an alternative decomposition. In addition, the Lorentzian–Gaussian mixing parameter ηj is fixed from unconstrained fits to x = 0.44; a sensitivity analysis allowing ηj to float or testing the dependence of the conclusions on this choice is needed.","section":"Sec. VIII A and VIII D"},{"comment":"The attribution of the extra split components (E(TO'2), E(TO'3), A1(TO'3)) to local disorder is partly circular: these components are first introduced to remove structured residuals, and their persistence is then interpreted as evidence of local symmetry breaking. An independent falsifiable test is required—for example, temperature dependence of the subpeaks, comparison with polarized single-crystal or epitaxial-film data, or a quantitative comparison with a phase-mixture or monoclinic basis—to support the disorder interpretation rather than treat it as an assumption of the fitting procedure.","section":"Sec. VIII C and IX"}],"minor_comments":[{"comment":"Residual plots are not shown in the main text. Including them would support the claim that the fits are free of systematic structure.","section":"Fig. 4"},{"comment":"The Savitzky–Golay smoothing parameters (window size, polynomial order) used for CMCD should be stated, since derivative-based peak counting is sensitive to smoothing.","section":"Sec. VIII B"},{"comment":"The fitted peak positions, widths, and amplitudes should be plotted with error bars or confidence intervals; without them, the 'continuous evolution' claim is hard to evaluate.","section":"Sec. VIII D / Fig. S2"},{"comment":"The phrase 'revealing the tetragonal-to-rhombohedral transition as a continuous redistribution of spectral intensity rather than emergence of new Raman modes' should be softened to reflect that this conclusion depends on the chosen decomposition basis.","section":"Abstract"},{"comment":"The text says 'the ideal structure supports up to 13 Raman-active components' but then states 'only approximately ten should be experimentally resolvable'; the relationship between these numbers and the final 13-component fit should be clarified.","section":"Sec. VIII B"}],"recommendation":"major_revision","confidential_remarks":"The symmetry part of the manuscript is solid and could be published as a methods-oriented contribution. The experimental section, however, overclaims: the fixed 13-component basis and the single-composition CMCD check mean that the 'no new modes' and 'no monoclinic signature' conclusions are not independently established. I would send the paper back for major revision, asking for composition-resolved peak-count validation, uncertainty quantification, and an explicit alternative-basis model comparison."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThe one thing to know: this is not a paper with a new physics claim; it is a careful, mostly correct worked example of standard symmetry methods applied to PZT, followed by a plausible but under-supported reanalysis of powder Raman data. The tensor tables and mode classifications are textbook results, and the Raman conclusion repeats what Buixaderas et al. already reported. If you work on PZT Raman, it is a useful compendium; if you are looking for new results, there are not many.\n\nWhat it does well: the group theory branch is sound. The BV and correlation decompositions agree, the Raman tensors match Nye/Rousseau, and the authors are transparent about the ideal-mode framework. The fitting workflow in Fig. 5 is a sensible statement of practice, and the CMCD check is the right kind of independent sanity check, even though the paper admits it only \"supports, though does not uniquely determine\" the 13-component basis.\n\nThe soft spots are exactly where the reader and the stress-test note point. The 13-component basis is fixed a priori for all compositions, and the mixing parameter is pinned from x=0.44. That choice makes the \"no new modes\" conclusion partly a built-in assumption: any extra spectral weight at other compositions is absorbed by the free positions, widths, and amplitudes. The CMCD validation was only performed at x=0.44, and no parameter uncertainties are reported, so the \"continuous evolution\" claim is not quantitatively grounded. The paper does partly concede this: it calls the Axe-band amplitudes and widths qualitative. But the abstract and conclusions state the continuous redistribution more strongly than the evidence supports. Raw spectra and a test with phase-specific mode counts would cure most of this.\n\nAlso, the authors cite and agree with Buixaderas et al., which is honest and actually weakens the novelty claim. The citation pattern itself looks fine.\n\nWho is this for? A practitioner who wants a self-contained symmetry reference for PZT tensor properties and Raman mode assignments, and a cautionary example of how fitting basis choice can drive interpretation. It deserves a serious referee because the worked derivations are useful and the data-analysis critique matters, but it needs revision: report uncertainties, provide raw data, and soften the central claim to what the fitting actually proves. I would send it out, expecting major revision.","headline":"A sound, standard symmetry reference plus a plausible-but-under-supported Raman reanalysis; the main conclusion is partly built into the fixed 13-component fit.","tokens_in":23960,"tokens_out":2342,"would_cite":false,"duration_ms":24844,"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 single symmetry-based 13-phonon description accounts for PZT's room-temperature Raman spectra across the morphotropic phase boundary, with the tetragonal-to-rhombohedral transition appearing as a continuous redistribution of spectral inte","keywords":["ferroelectric perovskites","PZT","morphotropic phase boundary","Raman spectroscopy","point-group symmetry","physical-property tensors","vibrational mode classification","local disorder"],"falsifier":"A reproducible monoclinic A'/A'' doublet near 280 cm^-1, or a new Raman branch appearing as composition crosses the morphotropic phase boundary in polarization-resolved or low-temperature spectra, would contradict the no-new-modes, no-long-range-monoclinic claim; so would a fit using a monoclinic symmetry basis that achieves systematically lower residuals than the 13-component tetragonal basis.","tokens_in":23010,"feed_emoji":"🔬","tokens_out":5115,"duration_ms":51598,"temperature":0.7,"pith_summary":"The paper tries to establish that crystallographic point-group symmetry alone provides a consistent baseline for both macroscopic property tensors and microscopic lattice vibrations, and that this baseline resolves how the ferroelectric perovskite PZT evolves across its morphotropic phase boundary. Applying the symmetry-derived mode classification to previously reported room-temperature powder Raman spectra for compositions x = 0.44 to 0.58, the paper argues that a single set of 13 phonon components fits all compositions without adding new modes. The tetragonal-to-rhombohedral transition is therefore a continuous transfer of spectral weight among existing components, and the extra subpeaks observed in some modes reflect local disorder and anharmonicity rather than a distinct crystallographic phase. A sympathetic reader would care because this gives a physically constrained way to interpret disorder-broadened powder Raman data and to distinguish long-range symmetry changes from local structural heterogeneity in a benchmark piezoelectric material.","feed_headline":"No new Raman modes across PZT's morphotropic phase boundary","feed_subtitle":"Thirteen symmetry-allowed phonon components fit every composition; the phase change is just spectral-weight transfer.","key_machinery":"The machinery is a two-branch symmetry workflow. The macroscopic branch applies the invariance condition that property tensors must be unchanged by every point-group operation, yielding the independent components of pyroelectric, dielectric, and piezoelectric tensors. The microscopic branch uses factor-group analysis and site-symmetry correlation of atomic displacements to decompose the vibrational representation into irreducible representations and construct Raman tensors. For the spectra, the load-bearing tool is the 13-component pseudo-Voigt decomposition with a composition-independent Lorentzian–Gaussian mixing parameter fixed from the x = 0.44 reference; a derivative-based curvature ana","core_discovery":"From the point groups of PZT's cubic, tetragonal, rhombohedral, and monoclinic phases, the paper derives the allowed pyroelectric, dielectric, and piezoelectric tensor forms and the Raman-active phonon irreps, then uses the predictions to decompose prior room-temperature powder Raman spectra (x = 0.44 to 0.58). A single 13-component pseudo-Voigt basis, the count expected from tetragonal symmetry with LO–TO splitting, fits every composition without new peaks. Fitted positions, widths, and amplitudes evolve smoothly; spectral weight transfers from E(TO2) to its shoulder E(TO'2) as rhombohedral character grows. The paper concludes that the tetragonal-to-rhombohedral transition is a continuous i","pith_inferences":["If the no-new-modes claim is right, the continuous spectral-weight transfer strengthens the polarization-rotation and adaptive-nanodomain picture of the morphotropic phase boundary: the macroscopic phase change is a gradual repopulation of local vibrational environments rather than a first-order switch between two rigid lattices.","A testable extension is to fit the same spectra with a monoclinic symmetry basis and with composition-dependent mode counts; if a monoclinic basis does not reduce residuals, the no-long-range-monoclinic conclusion holds, but if it does, the 13-component basis may be underfitting.","Polarization-resolved Raman on oriented PZT films across the morphotropic phase boundary should show the tetragonal B1 mode intensity declining continuously with Zr content, providing a quantitative probe of the same transition that powder data only imply.","Low-temperature Raman on the same compositions could separate thermal broadening from genuine symmetry change: if monoclinic A'/A'' splitting appears on cooling, the room-temperature conclusion is a broadening artifact; if it does not, monoclinic order is absent even at low temperature."],"forward_implications":["Across x = 0.44 to 0.58, the same 13 phonon components describe room-temperature powder Raman spectra, so composition change is tracked by continuous peak-parameter shifts rather than by mode creation or loss.","Because tetragonal and rhombohedral PZT have the same zone-center Raman mode count, Raman spectra alone cannot distinguish the phases; identifying the transition requires following spectral-weight transfer or a symmetry marker such as the tetragonal B1 mode.","The absence of resolvable monoclinic A'/A'' splitting at room temperature indicates that any monoclinic distortion near the morphotropic phase boundary is short-range and does not appear as a distinct long-range phase in powder Raman data.","Persistent subpeak structure in E and A1 modes across all compositions implies that local B-site disorder and anharmonicity are intrinsic to the morphotropic phase boundary region, not anomalies of one composition.","The symmetry workflow itself is transferable: the same point-group operations can generate tensor forms and vibrational classifications for other ferroelectric perovskites and complex functional materials."],"fun_headline_variants":["PZT phase transition is spectral-weight redistribution, not new modes","Single 13-mode fit explains all PZT spectra across boundary","PZT's morphotropic transition is continuous intensity transfer","Symmetry-based fit reveals PZT transition as peak shift, not new modes"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The 13-component pseudo-Voigt basis, with its mixing parameter locked to the x = 0.44 spectrum, is assumed to stay adequate for every composition through x = 0.58; if a different number of peaks or a phase-specific basis fits the data equally well, the conclusion that no new modes emerge is not uniquely supported.","fun_headline_variants_meta":{"raw":{"variants":["PZT phase transition is spectral-weight redistribution, not new modes","Single 13-mode fit explains all PZT spectra across boundary","PZT's morphotropic transition is continuous intensity transfer","Symmetry-based fit reveals PZT transition as peak shift, not new modes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00109,"raw_usage":{"total_tokens":4390,"prompt_tokens":741,"completion_tokens":3649,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":485,"completion_tokens_details":{"reasoning_tokens":3577}},"tokens_in":485,"tokens_out":3649,"duration_ms":26695,"temperature":1.0,"reasoning_tokens":3577,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T07:44:15.433006+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A reproducible monoclinic A'/A'' doublet near 280 cm^-1, or a new Raman branch appearing as composition crosses the morphotropic phase boundary in polarization-resolved or low-temperature spectra, would contradict the no-new-modes, no-long-range-monoclinic claim; so would a fit using a monoclinic symmetry basis that achieves systematically lower residuals than the 13-component tetragonal basis.","supporting_citations":[],"review_version":1}