{"id":"0209de96-24de-486a-affc-68682eeea1b3","arxiv_id":"2507.01179","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Applying the energy fluctuation model to Raman linewidth data, the authors extract activation energies and relaxation times for molecular modes in the ferroelastic transition of TMACd(N3)3.","lead":"This paper uses an energy fluctuation model to explain how the widths of certain Raman peaks change as a hybrid perovskite material transforms from one crystal shape to another. It extracts how long molecular motions take to settle and the energy barriers they cross, which could help design better cooling materials.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Relaxation times are computed from Eq. (11) as τ = Γ/Q², whose units are cm⁻¹, not seconds; the reported ns values and the claimed hierarchy of slow torsional/librational renormalization are therefore unsupported.","rationale":"The reader's weakest_assumption focused on whether the strain-based order parameter Q is a valid proxy over the fitted range. That is a legitimate concern, but the more load-bearing and directly checkable problem is dimensional: Eq. (11) cannot produce a relaxation time from Γ/Q² without a missing conversion, and Eq. (4) contains a unit inconsistency that affects every fit. The reader's rationale did flag both of these, so there is substantial agreement, but the primary single concern here is not the same as the stated weakest assumption. The linewidth fits themselves may still support the claim that the EF model captures the Raman broadening, and the activation energies could survive an Arrhenius re-analysis. For that reason, the appropriate verdict remains CONDITIONAL rather than a full rejection: the central qualitative claim is plausible, but the reported ns relaxation times and the accompanying molecular-dynamics interpretation are not supported by the equations as written. A focused unit audit and refit would settle whether the quantitative claims survive.","tokens_in":17323,"tokens_out":4346,"duration_ms":55610,"concrete_test":"Perform a unit and re-derivation audit of Eqs. (4) and (11) against the original references (Schaack and Winterfeldt; Laulicht; and the source of Eq. 11). Concretely, recompute the Fig. 3 relaxation curves using the correct damping relation τ = 1/(2π c Γ) and compare with the reported ns values; then refit the Table 1 linewidths with the standard numerator T(T−TcQ²) and compare the fitted A′ values and resulting Table 2 activation energies. If the corrected relaxation times are picoseconds rather than nanoseconds, or if the corrected fit parameters shift by more than the stated uncertainties, the paper's quantitative conclusions require revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's most distinctive quantitative outputs are the relaxation times in nanoseconds and the mode-dependent activation energies. Both rest on equations with dimensional inconsistencies. Eq. (11) defines τ = Γ/Q², where Γ is a Raman linewidth in cm⁻¹ and Q is dimensionless, so τ has units of cm⁻¹; no stated conversion factor turns this into seconds. Physically, a damped oscillator has a lifetime that decreases as the linewidth increases, whereas Eq. (11) makes τ grow with Γ and diverge as Q→0, so the reported 'longer relaxation times near Tc' behavior is an artifact of the formula rather than a measured dynamical feature. Eq. (4) has a separate problem: its numerator is written T(T−Q²), mixing units of K with the dimensionless Q², whereas the standard EF expression requires T(T−TcQ²) (or the equivalent); the printed form changes the temperature dependence used in every fit. Because Table 2 activation energies come from Arrhenius fits of the fitted Γ over two arbitrary temperature intervals, and no raw data or code are provided, the fitted parameters cannot be independently cross-checked. The EF model may still describe the linewidth anomaly, but the derived relaxation times are not physically meaningful as reported.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript applies the energy-fluctuation (EF) model to the temperature-dependent Raman linewidths of six modes in TMACd(N3)3 across the ferroelastic gamma-to-delta transition at T_C = 322 K. The order parameter is obtained from the symmetry-adapted spontaneous strain combination (e0^2 + et^2), fitted to a first-order Landau 2-4-6 potential, and is then inserted into the EF linewidth expression. From the fitted linewidths the authors extract activation energies via Arrhenius plots and relaxation times via Eq. (11), reporting nanosecond-scale values and a mode-dependent hierarchy that is compared with DMACd(N3)3. The paper concludes that methyl torsion and azide librational modes couple strongly to the symmetry-breaking strain and renormalize more slowly after the transition.","tokens_in":17463,"tokens_out":6058,"duration_ms":180242,"significance":"If quantitatively sound, this work would extend the EF-model methodology to a ferroelastic hybrid organic-inorganic perovskite using a strain-derived order parameter, and would provide activation energies and relaxation times relevant to barocaloric and order-disorder materials. The qualitative picture that torsional and librational modes are particularly coupled to the ferroelastic distortion is plausible and consistent with the structural discussion. The explicit use of published lattice-parameter data to construct a first-order Landau order parameter is a concrete strength. However, the quantitative outputs currently rest on two equations with dimensional or consistency problems, so the numerical values, including the claimed comparison with DMACd(N3)3 and the relaxation-time hierarchy, are not established by the manuscript as written. No raw data or fitting code are provided, which further limits independent verification.","major_comments":[{"comment":"The energy-fluctuation linewidth expression is written as Gamma_EF = Gamma0' + A' [T(T - Q^2)/(T - T_C(1 - Q^2))]^(1/2). Since T is in kelvin and Q is normalized and dimensionless, the factor T(T - Q^2) is dimensionally inconsistent; the standard EF expression contains T(T - T_C Q^2) or the equivalent. This changes the temperature dependence used in every fit in Section 3 and therefore affects the Gamma_EF values that are the basis of Tables 1 and 2. The fits and all derived activation energies need to be redone with the corrected expression before the quantitative conclusions can be accepted.","section":"2, Eq. (4)"},{"comment":"Eq. (11) defines tau = Gamma/Q^2. With Gamma in cm^-1 and Q dimensionless, tau has units of cm^-1, not seconds. No conversion factor or physical derivation is supplied, so the reported nanosecond relaxation times in Fig. 3 and in the text (for example, 10.2 ns for L(N3-) and 7.7 ns for tau(CH3)) are unsupported. In addition, because this definition makes tau grow with Gamma and diverge as Q approaches zero, the claimed 'longer relaxation times near T_C' behavior is an artifact of the formula unless a proper model-based derivation is provided. The relaxation-time analysis must be rederived with consistent units or removed.","section":"2, Eq. (11)"},{"comment":"The text states that the EF model was fitted to linewidth data for T < T_C, but Table 1 lists temperature intervals extending to 343 K, which is above T_C = 322 K where the order parameter is set to zero. The manuscript should clarify which data points were included and how Eq. (4) is evaluated for T > T_C. If high-temperature points were included, the fits and the extracted parameters need to be revised, and the statement about fitting below T_C corrected.","section":"3, Table 1"}],"minor_comments":[{"comment":"Please write explicitly how the normalized order parameter Q entering Eq. (4) is obtained from (e0^2 + et^2); the text mentions Q/Q_max but does not give the exact relation between Q and the strain components used in the fits.","section":"2, after Eq. (8)"},{"comment":"The caption for Fig. 3(f) labels the 3032 cm^-1 mode as nu_s(CH3), but Table 1 and the text identify this mode as nu_as(CH3).","section":"3, Fig. 3"},{"comment":"Reference [47] (Montgomery et al., Human Pathology) appears to be an unrelated medical citation and is not appropriate as a source for the energy-fluctuation model; the relevant literature for TGS/TGSe should be cited instead.","section":"References"},{"comment":"Eq. (5a) writes '0 < (T_C - T) < T', which is not a meaningful temperature range as printed; this is presumably a typographical issue since Eq. (5) is not used in the subsequent first-order analysis.","section":"2, Eq. (5)"},{"comment":"No raw FWHM data or fitting code are provided; since the linewidth data are reused from Ref. [24], a table of the digitized FWHM values or a clear data-availability statement would allow the fits and extracted parameters to be checked independently.","section":"Reproducibility"}],"recommendation":"major_revision","confidential_remarks":"The two dimensional and consistency problems in Eqs. (4) and (11) are load-bearing because the paper's central quantitative outputs are activation energies and nanosecond relaxation times. The linewidth fits themselves look plausible, and the strain-based order parameter is a reasonable idea, so a corrected version with refits and a proper derivation or removal of the relaxation-time formula could become publishable. I recommend major revision rather than rejection, but the revision must be substantive and not merely cosmetic."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague, the short version: this is a reasonable reanalysis of Raman linewidths for the ferroelastic transition in TMACd(N3)3 using the energy fluctuation model, but the two headline outputs—relaxation times in ns and activation energies—rest on equations with dimensional errors. The paper is worth a second look after correction, not as is.\n\nWhat's actually new: first application of the EF model to this compound, with the order parameter built from symmetry-adapted spontaneous strain (e0^2+et^2) from a Landau 2-4-6 potential. That's a sensible way to handle a first-order ferroelastic transition, and the fit giving Ttr - Tc* = 8 K is consistent with weakly first-order behavior. The linewidth fits to six modes look reasonable in the figures, and the discussion connecting activation energies to hydrogen bonding (or its absence) in TMA vs DMA is plausible, though it inherits the parameter problems.\n\nWhere it falls down: Eq. (4) as printed is dimensionally inconsistent—the numerator T(T-Q^2) mixes kelvin with dimensionless Q^2. The standard expression has T_C in the numerator, T(T - T_C Q^2). This changes the temperature dependence in every fit, so the extracted Γ_EF values and the activation energies built from them are questionable. Eq. (11) is worse: τ = Γ/Q^2 has units of cm^-1, not seconds. Reporting ns values requires a conversion factor that never appears, and the formula itself makes τ grow with Γ and diverge as Q→0, which is the opposite of physical lifetime behavior. The claimed slow renormalization of torsional/librational modes near Tc is an artifact of the formula, not a measured dynamical feature. The Arrhenius activation energies are also split into two arbitrary temperature intervals with no error propagation, so the reported U values need independent validation. No code or raw data are provided, so the fits cannot be cross-checked.\n\nOn balance, the central idea—strain-coupled EF model for this compound—is plausible, and the linewidth anomaly itself is real. But the quantitative results as presented are unsupported. I'd send it to a referee who can demand correction of the units and refitting, but I would not cite the relaxation times in the current form. If the authors fix Eqs. (4) and (11) and redo the fits, the paper could be a solid contribution to the barocaloric HOIP literature. For now, treat it as a promising draft with a load-bearing units error.","headline":"A plausible strain-coupled EF model for TMACd(N3)3 whose reported relaxation times and activation energies rest on dimensional errors; the paper needs correction before its quantitative claims can be trusted.","tokens_in":18135,"tokens_out":2796,"would_cite":false,"duration_ms":28258,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["78.30.-j","64.70.K-","63.20.-e"],"model":"deepseek-v4-flash","headline":"The energy fluctuation model, fed a strain-derived order parameter, reproduces the Raman linewidths across the 322 K ferroelastic transition in TMACd(N3)3 and yields mode-specific activation energies and relaxation times.","keywords":["hybrid organic-inorganic perovskite","ferroelastic phase transition","Raman linewidth","energy fluctuation model","order-disorder transition","spontaneous strain","activation energy","relaxation time"],"falsifier":"Measure the integrated intensity of a superlattice reflection across the $\\gamma\\to\\delta$ transition in TMACd(N3)3: if the normalized direct order parameter does not track the normalized $e_0^2+e_t^2$ curve used in Eq. (4), the extracted activation energies and relaxation times are artifacts of a wrong $Q$.","tokens_in":16992,"feed_emoji":"🔬","tokens_out":8925,"duration_ms":90314,"temperature":0.7,"pith_summary":"This paper asks whether the energy fluctuation (EF) model can describe how Raman mode widths change as the hybrid perovskite TMACd(N3)3 crosses its 322 K first-order ferroelastic transition. The authors feed the model a strain-derived order parameter rather than the usual pseudospin expression, because the transition is first-order and discontinuous. They find that the EF model tracks the critical broadening of six internal and lattice modes, and from the fits they extract activation energies and relaxation times for each mode. The values put the methyl torsion and azide libration near or above the thermal energy $k_B T_C$, and those modes renormalize more slowly than stretching modes, suggesting that they are the motions that carry the order-disorder mechanism.","feed_headline":"Strain model reproduces Raman linewidths at a 322 K phase transition","feed_subtitle":"Torsional modes renormalize slowly and activation energies reach kB Tc, exposing what drives the switch.","key_machinery":"The central object is the EF linewidth equation $\\Gamma_{EF}=\\Gamma_0' + A'[T(T-Q^2)/(T-T_C(1-Q^2))]^{1/2}$, where $Q$ is the order parameter; the authors supply $Q$ not from the usual pseudospin expression but from the normalized symmetry-adapted strain $\\sqrt{e_0^2+e_t^2}$, with $e_0=e_1-e_2$ and $e_t=(2e_3-e_1-e_2)/\\sqrt3$, whose temperature dependence is fixed by the standard first-order solution of a 2-4-6 Landau potential. This strain-fed $Q$ is the bridge that lets static lattice-parameter data enter a dynamical linewidth formula, and the same $Q$ converts fitted linewidths into relaxation times through $\\tau=\\Gamma/Q^2$.","core_discovery":"The paper's central claim is that the energy fluctuation (EF) model, supplied with an order parameter derived from symmetry-adapted spontaneous strain, quantitatively reproduces the critical broadening of Raman linewidths in TMACd(N3)3 across its 322 K $\\gamma \\to \\delta$ ferroelastic transition. The linewidths of six modes $L(\\mathrm{N_3^-})$, $\\tau(\\mathrm{CH_3})$, $\\nu_s(\\mathrm{NC_4})$, $\\nu_s(\\mathrm{N_3})$, $\\nu_s(\\mathrm{CH_3})$, and $\\nu_{as}(\\mathrm{CH_3})$ follow $\\Gamma_{EF} = \\Gamma_0' + A' [T(T-Q^2)/(T-T_C(1-Q^2))]^{1/2}$ below $T_C$, with $Q$ normalized from $(e_0^2+e_t^2)$ via a first-order Landau 2-4-6 solution. From these fits the authors extract activation energies up to roughly 24.7 meV and relaxation times up to roughly 10.2 ns, with the methyl torsion and azide libration showing the slowest renormalization. They interpret the near-$k_B T_C$ barriers and the absence of hydrogen bonds around the TMA cation as evidence that these torsional and librational motions, coupled to symmetry-breaking strain, drive the ferroelastic order-disorder mechanism.","pith_inferences":["If strain is a faithful proxy for the pseudospin order parameter, the same workflow could predict relaxation hierarchies in other first-order hybrid perovskite transitions from static lattice-parameter data alone, without dynamical measurements.","A sharper test would compare the EF-derived relaxation times with independent values from NMR or inelastic neutron scattering; agreement would confirm that the linewidth broadening is dominated by order-parameter fluctuations rather than ordinary anharmonic decay.","The near-$k_B T_C$ activation-energy pattern suggests a criterion for identifying the driving mode in a hybrid perovskite: the mode whose barrier matches the thermal energy at $T_C$ is the one that must freeze or reorient for the transition to occur.","Because the paper ties larger activation barriers to the absence of hydrogen bonds, one can test this by synthesizing TMA analogues with different B-site metals or with partial deuteration and checking whether the methyl-torsion barrier shifts as predicted."],"forward_implications":["The same strain-fed EF fitting procedure can be reused on any first-order ferroelastic transition where lattice parameters are known, turning Raman linewidths into a probe of mode-specific activation energies and relaxation times.","The near-$k_B T_C$ activation energies for the methyl torsion and azide libration imply that these reorientational motions are thermally active at the transition and participate in the symmetry breaking.","Longer relaxation times for $\\tau(\\mathrm{CH_3})$ and $L(\\mathrm{N_3^-})$, about 7 to 10 ns, identify the slow renormalizing degrees of freedom, which should be the modes most sensitive to external stimuli such as pressure in barocaloric applications.","Because the same motions in the DMA analogue cost less energy, the absence of N-H...N hydrogen bonds and the larger TMA cation raise the potential barrier, giving a design rule for tuning transition dynamics in azide perovskites."],"supporting_citations":[{"why":"It supplies the experimental Raman linewidths and the transition temperature $T_C=322$ K that are fitted throughout the paper.","marker":"[24]"},{"why":"It establishes the above-room-temperature ferroelastic transition in this compound and identifies azide swing and TMA rotation as the driving motions.","marker":"[38]"},{"why":"It provides the Ising pseudospin-phonon Hamiltonian from which the critical linewidth formula descends.","marker":"[42]"},{"why":"It extends the pseudospin-phonon model and gives the phonon frequency and width expressions that the EF linewidth equation approximates.","marker":"[43]"},{"why":"It derives the simplified EF and pseudospin-phonon linewidth formulas, Eqs. (3) and (4), used to fit the measured widths.","marker":"[45]"},{"why":"It is the original source of the energy fluctuation model expression for the temperature-dependent critical linewidth.","marker":"[49]"},{"why":"It supplies the symmetry-adapted strain formalism and the 2-4-6 Landau solution used to build the order parameter Q from lattice parameters.","marker":"[55]"},{"why":"It provides the DMACd(N3)3 activation energies and relaxation times that serve as the comparison benchmark.","marker":"[35]"}],"fun_headline_variants":["Strain model matches Raman widths at ferroelastic transition","Energy fluctuation model fits Raman linewidths near 322 K","Why torsional modes lag at the 322 K ferroelastic switch","Activation energies reach kB Tc in strain-coupled transition"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The strain-derived order parameter from the Landau fit is assumed to be a faithful proxy for the pseudospin order parameter in the EF equation across the whole fitted temperature range; if strain decouples from the true order parameter near $T_C$, the extracted activation energies and relaxation times are not physically meaningful.","fun_headline_variants_meta":{"raw":{"variants":["Strain model matches Raman widths at ferroelastic transition","Energy fluctuation model fits Raman linewidths near 322 K","Why torsional modes lag at the 322 K ferroelastic switch","Activation energies reach kB Tc in strain-coupled transition"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00079,"raw_usage":{"total_tokens":3614,"prompt_tokens":1209,"completion_tokens":2405,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":825,"completion_tokens_details":{"reasoning_tokens":2337}},"tokens_in":825,"tokens_out":2405,"duration_ms":18285,"temperature":1.0,"reasoning_tokens":2337,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T20:58:54.605412+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the integrated intensity of a superlattice reflection across the $\\gamma\\to\\delta$ transition in TMACd(N3)3: if the normalized direct order parameter does not track the normalized $e_0^2+e_t^2$ curve used in Eq. (4), the extracted activation energies and relaxation times are artifacts of a wrong $Q$.","supporting_citations":[{"cited_title":"Levola and R","cited_arxiv_id":null,"evidence_quote":"It supplies the experimental Raman linewidths and the transition temperature $T_C=322$ K that are fitted throughout the paper."},{"cited_title":"Zhang, Y","cited_arxiv_id":null,"evidence_quote":"It establishes the above-room-temperature ferroelastic transition in this compound and identifies azide swing and TMA rotation as the driving motions."},{"cited_title":"Trzebiatowska, M","cited_arxiv_id":null,"evidence_quote":"It provides the Ising pseudospin-phonon Hamiltonian from which the critical linewidth formula descends."},{"cited_title":"Kurt, Calculation Of The Damping Constant (FWHM), The Relaxation Time, And The Activation Energy As A Function Of Temperature For DmaCd(N3)3, J Mol Struct 1244, 130901 (2021)","cited_arxiv_id":null,"evidence_quote":"It extends the pseudospin-phonon model and gives the phonon frequency and width expressions that the EF linewidth equation approximates."},{"cited_title":"Yurtseven and A","cited_arxiv_id":null,"evidence_quote":"It derives the simplified EF and pseudospin-phonon linewidth formulas, Eqs. (3) and (4), used to fit the measured widths."},{"cited_title":"Aizu, Possible Species of Ferromagnetic, Ferroelectric, and Ferroelastic Crystals, Phys Rev B 2, 754 (1970)","cited_arxiv_id":null,"evidence_quote":"It is the original source of the energy fluctuation model expression for the temperature-dependent critical linewidth."},{"cited_title":"Montgomery et al., Reproducibility of the diagnosis of dysplasia in Barrett esophagus: A reaffirmation, Hum Pathol 32, 368 (2001)","cited_arxiv_id":null,"evidence_quote":"It supplies the symmetry-adapted strain formalism and the 2-4-6 Landau solution used to build the order parameter Q from lattice parameters."},{"cited_title":"Asaji, Y","cited_arxiv_id":null,"evidence_quote":"It provides the DMACd(N3)3 activation energies and relaxation times that serve as the comparison benchmark."}],"review_version":1}