{"id":"5626cf8e-aec4-40bd-9155-e7035f791cd4","arxiv_id":"2502.01176","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"The nonlinear Murnaghan moduli l and m of a polystyrene-based material change with temperature about two orders of magnitude more strongly than the linear Lame moduli.","lead":"The authors measured how the nonlinear elastic moduli of a polystyrene-based material change with temperature between 25 and 65 degrees Celsius. They found these moduli are far more temperature-sensitive than the linear elastic moduli, which could matter for modeling and non-destructive testing of polymers.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Measurement protocol may couple pressure and effective path temperature, biasing the temperature slopes of l and m.","rationale":"The reader identified spatial temperature heterogeneity as the weakest assumption. I agree that heterogeneity matters, but the more load-bearing and more specific issue is that the measurement protocol introduces a systematic coupling between applied pressure and the effective temperature along the ultrasonic path, because pressure changes the thermal contact at the jaws. This coupling directly biases the pressure–velocity slopes that determine l and m, and it is acknowledged in the text. The numerical estimate using the reported dV/dT shows the bias could be comparable to the original signal, so the central claim of a two-orders-of-magnitude temperature susceptibility rests on an unverified assumption about this coupling. The proposed equilibrium control experiment is straightforward and would settle the question. Other aspects of the analysis — the linear fits, the frequency-shifted consistency check, and the agreement of zero-pressure velocities with literature — are internally coherent and support the manuscript. Therefore the verdict remains CONDITIONAL: the conclusion is plausible but requires an explicit test of this potential systematic error before it should be accepted without qualification.","tokens_in":10409,"tokens_out":10045,"duration_ms":111311,"concrete_test":"Perform control measurements at two or three temperatures (25 °C, 45 °C, 65 °C) in true thermal equilibrium: stabilize the sample for 40–50 minutes, hold T fixed, and step the static pressure over the same 3–15 MPa range, recording V(P) at each T. Compare the resulting α_j and derived l, m, n with those from the heating-protocol slopes. Agreement within quoted uncertainties rules out the pressure–temperature coupling bias; disagreement reveals its magnitude. A complementary check is to place internal thermocouples along the sample length to verify that the path-averaged temperature does not shift by more than ~0.5 °C across the pressure range.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The nonlinear moduli l, m, and n are derived from the pressure slopes α_j of M_j(P) at each nominal temperature. The protocol extracts these slopes by continuously heating the sample at each of five fixed initial pressures, rather than holding temperature constant and stepping pressure. The authors explicitly note in Section 2 that increasing pressure tightens the jaw contact, raises heat transfer, and can cause temperature fluctuations. Because the ultrasonic transit time is an integral over the 50 mm sample, a pressure-dependent spatial temperature distribution along the path would contaminate dV/dP with a spurious component (dV/dT)·(dT_eff/dP). The reported dV/dT is about −2 m/s/°C for longitudinal waves, and the genuine acousto-elastic velocity change over 16 MPa is only a few m/s. A pressure-induced shift of even 1 °C in the effective path temperature would therefore produce a slope bias comparable to the signal. The observed up-to-5 °C heterogeneity at 65 °C makes this mechanism plausible, and its magnitude is not propagated into the reported α_j uncertainties. If dT_eff/dP is non-negligible, the extracted l(T) and m(T) — and hence the two-orders-of-magnitude claim — could be systematically distorted.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports measurements of the Murnaghan third-order elastic moduli l, m, and n of a polystyrene-based glassy polymer as functions of temperature (25–65 °C) and ultrasonic frequency (0.7–3 MHz), using the acousto-elastic effect. Pressure slopes α_j of the effective moduli M_j(P,T) are extracted at each nominal temperature, and Eqs. (2)–(4) are used to obtain l, m, and n. The central claims are that the temperature susceptibilities of l and m are about two orders of magnitude larger than those of the linear Lamé moduli λ and μ (slopes b_l = -0.44 ± 0.06 GPa/°C, b_m = -0.12 ± 0.04 GPa/°C versus b_λ = -5.0 ± 0.5 MPa/°C, b_μ = -2.2 ± 0.5 MPa/°C), that n has no resolvable temperature dependence (b_n = 0.00 ± 0.03 GPa/°C), and that the temperature susceptibilities of l and m are nearly frequency-independent over the studied range, as demonstrated by the shifted-moduli analysis in Eqs. (7)–(9).","tokens_in":10607,"tokens_out":8046,"duration_ms":85670,"significance":"If the results are correct, they provide a rare quantitative data set on the temperature dependence of third-order elastic moduli of a glassy polymer and a striking contrast between linear and nonlinear thermal susceptibilities. The analysis chain is transparent: the moduli are obtained directly from measured acousto-elastic slopes with standard formulas, with no fitted theoretical model for the target quantities. The shifted-moduli consistency check is a useful way to display frequency independence of the temperature slopes. The paper also connects the observations to relaxation processes and time-temperature superposition, and it compares with prior results on vitreous silica and metal-matrix composites. The main risk is systematic temperature bias in the pressure slopes, which is not propagated into the reported uncertainties.","major_comments":[{"comment":"The pressure slopes α_j used in Eqs. (2)–(4) are not obtained under isothermal conditions. The protocol heats the sample from 25 °C to 65 °C at each fixed pressure, and the authors explicitly note that increasing pressure improves jaw contact and changes heat transfer, so the thermal state at different pressures is not identical. The reported surface-temperature heterogeneity reaches 5 °C at 65 °C (Fig. 1c), while the longitudinal velocity changes with temperature at roughly −2 m/s per °C (Table 1). Over the 16-MPa pressure span the genuine acousto-elastic velocity change is only of order a few m/s, so a pressure-induced effective-path-temperature shift of order 1 °C would bias α_j at the same level as the signal. This bias is not propagated into b_l and b_m or into the two-orders-of-magnitude claim. The authors should either quantify dT_eff/dP along the ultrasonic path from the thermal-imaging data, or provide a direct isothermal pressure-stepping validation at least at two temperatures.","section":"Section 2, protocol description and Fig. 2"},{"comment":"The conclusion that n is temperature-independent within error rests on b_n = 0.00 ± 0.03 GPa/°C, but the individual n(T) points are not shown with uncertainties, and n is obtained from the small difference α_y − α_z (Eq. (4)), where the two shear sensitivities are nearly equal. The reported slope uncertainty does not include correlated errors in α_y and α_z or the temperature-heterogeneity effect described in the previous comment. Please provide a table of individual l, m, and n values with uncertainties at each temperature and frequency, and state how the slope uncertainties were computed, including whether correlations between the α_j were accounted for.","section":"Section 3.2, Fig. 5 and slopes b_n"}],"minor_comments":[{"comment":"There is an apparent sign inconsistency: the text reports b_l = −0.44 ± 0.06 GPa/°C, while the label in Fig. 5(a) appears as '0.44 ± 0.06 GPa/°C'. In addition, the wording 'increased in their absolute values' should be reconciled with the sign of the plotted quantities, since the figures appear to show positive values of l and m that decrease with temperature.","section":"Section 3.2 and Fig. 5(a)"},{"comment":"The sentence 'This behavior is in good agreement with data previously published elsewhere and can be explained by rising viscosity of the material at elevating temperatures' appears to say the opposite of the intended meaning; viscosity of polymers typically decreases with increasing temperature.","section":"Section 3.1"},{"comment":"The text states that 'the correct value of the velocity is V_j = V_i(P)(1 + ε_x)' with ε_x = νP/E, and then says this correction is 'taken into account further when using Eqs. (2)–(4)'. It is unclear whether Eq. (1) already includes this correction or whether it is applied only later; please clarify to avoid double-counting or omission, especially because E is temperature-dependent and could slightly affect the temperature slopes.","section":"Section 2, Eq. (1)"},{"comment":"The individual data points in Figs. 4 and 5 are plotted without error bars. Adding error bars, at least at representative frequencies and temperatures, would help the reader assess the significance of the reported slopes and of the claimed frequency independence.","section":"Figures 4 and 5"},{"comment":"The phrase 'time-temperature superposition principle is reliable only for thermodynamically simple systems' would be more standard as 'thermorheologically simple systems'.","section":"Section 4"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the experimental approach is a direct extension of the authors' previous work. The principal concern is a possible systematic coupling between applied pressure and effective path temperature in the measurement protocol; this is a load-bearing issue for the central temperature-susceptibility claim, and it should be addressed quantitatively before publication. The circularity concern noted by the reader is not, in my view, serious: the moduli are extracted from measured slopes with standard formulas, and the shifted-moduli analysis is a consistency check rather than an independent determination."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nBottom line: this is a plausible and genuinely new dataset — first temperature dependence of Murnaghan moduli l, m, n for a polystyrene-based material over 25–65 °C — and the observation that the temperature susceptibility is roughly frequency independent is worth having. But the headline claim needs re-examination before it becomes citable, because the measurement protocol couples pressure and effective sample temperature, and the resulting systematic error is potentially the same size as the effect.\n\nWhat the paper does well: the analysis chain is transparent. Velocities come from phase shifts, M_j(P) is fitted linearly, and l, m, n follow from standard formulas. The authors are also honest about the temperature heterogeneity (up to 5 °C at 65 °C) and explain why they chose the heating-at-constant-pressure protocol. The shifted-moduli check (Eqs. 7–9) is a reasonable way to test frequency independence, and the reported collapse of those curves is decent evidence that the temperature slopes don't vary much across 0.7–3 MHz.\n\nSoft spots, in order of seriousness:\n\n1. Pressure–temperature coupling. The paper itself notes that increasing pressure tightens the jaw contact and raises heat transfer, which can change the temperature state. Since the measured velocity is a path average over the 50 mm sample, a pressure-dependent effective temperature would add a dV/dT term to the pressure slope. With dV/dT ≈ −2 m/s/°C and the genuine acoustoelastic change over 16 MPa only a few m/s, a 1 °C pressure-induced shift would be a large fraction of the signal. The authors report up to 5 °C spatial heterogeneity at the high end, but they do not estimate dT_eff/dP or propagate it into the α_j uncertainties. This is the main thing to fix.\n\n2. The “two orders of magnitude” wording is generous for m: b_m/b_µ ≈ 55, which is 1.7 orders. l is closer to 88, nearly two. Still a big effect, but it should be reported as roughly 1.5–2 orders.\n\n3. The plotted moduli in Figs. 4 and 5 have no error bars, which makes it hard to judge scatter, particularly for n. Adding them wouldn't be hard.\n\nMinor: the sample is a styrene copolymer with 10% EGDMA, not pure polystyrene; the title and abstract should carry that qualification.\n\nThe citation pattern looks fine; prior work is cited and the extension is clear.\n\nOverall, I'd send this to review. The central claim may survive or weaken after the temperature-coupling analysis, and the community would benefit from either outcome. I would not use the numbers in my own work until that is addressed.\n\nFor a reading group, this is a maybe — good case study in experimental systematic effects in acoustoelasticity.\n\nRecommendation: serious peer review, conditional on the authors quantifying the pressure–temperature coupling and adding error bars.","headline":"First temperature-dependent Murnaghan moduli for PS, though the headline effect may be partly biased by pressure–temperature coupling in the setup.","tokens_in":11149,"tokens_out":5171,"would_cite":false,"duration_ms":54407,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["43.35.Zc","62.20.de","81.05.Lg"],"model":"deepseek-v4-flash","headline":"The temperature susceptibility of the nonlinear elastic moduli l and m of polystyrene is two orders of magnitude larger than that of the linear moduli, while the n modulus is essentially temperature independent.","keywords":["temperature dependence","nonlinear elastic moduli","Murnaghan moduli","polystyrene","acousto-elastic effect","ultrasonic waves","glassy polymer","time-temperature superposition"],"falsifier":"A direct check would be to repeat the pressure-velocity sweeps at 55 and 65 °C with the sample held in a temperature bath or with local temperature sensors along the ultrasonic path, keeping heterogeneity below about 0.5 °C; if the extracted b_l and b_m change by more than their reported ±0.06 and ±0.04 GPa/°C uncertainties, the two-orders claim is compromised.","tokens_in":10202,"feed_emoji":"🔥","tokens_out":5971,"duration_ms":51827,"temperature":0.7,"pith_summary":"This paper measures how the nonlinear elastic moduli of polystyrene change with temperature, using ultrasonic wave velocities under applied pressure. It reports that the Murnaghan moduli l and m change with temperature about a hundred times faster than the linear Lamé moduli λ and µ in the 25–65 °C range, while the third nonlinear modulus n barely changes. The authors argue this matters because nonlinear elastic response governs how polymers behave under dynamic loads, and its temperature sensitivity is poorly known. If the measurements are right, temperature corrections to polymer nonlinearity are far larger than previously assumed, and the n modulus offers a convenient temperature-insensitive reference.","feed_headline":"Heat makes polystyrene's nonlinear stiffness change 100x faster","feed_subtitle":"Ultrasonic data show l and m moduli shift ~100x more than linear ones from 25 to 65 °C.","key_machinery":"The central object is the set of Murnaghan third-order elastic moduli l, m, n, which extend the linear Lamé description to leading nonlinear terms. They are extracted from the acousto-elastic effect: the slopes αx, αy, αz of the effective elastic moduli M_j(P,T) = ρ $V_j^{2}$ versus applied pressure for longitudinal and two shear wave polarizations. Equations (2)–(4) convert these slopes and the temperature-dependent Lamé moduli into l, m, n. The argument rests on the linearity of these pressure slopes and on the assumption that the measured wave velocities at each nominal temperature represent the sample's true state.","core_discovery":"Using the acousto-elastic effect, the authors measured velocities of longitudinal and shear ultrasonic waves in polystyrene under static pressures up to 16 MPa at temperatures from 25 to 65 °C and at four frequencies between 0.7 and 3 MHz. From the pressure slopes of the effective moduli they extracted the three Murnaghan third-order moduli l, m, and n via equations (2)–(4). The temperature slopes of the frequency-averaged moduli are b_l = −0.44 ± 0.06 GPa/°C, b_m = −0.12 ± 0.04 GPa/°C, and b_n = 0.00 ± 0.03 GPa/°C, compared with −5.0 ± 0.5 MPa/°C and −2.2 ± 0.5 MPa/°C for λ and µ. Thus l and m are roughly two orders of magnitude more temperature-sensitive than the linear moduli, while n shows no resolvable temperature dependence. The temperature susceptibility of l and m is essentially independent of frequency in the studied range.","pith_inferences":["If the two-orders ratio holds generally for glassy polymers, one would expect similar temperature susceptibility in other sub-Tg thermoplastics; a quick test is to repeat the same acousto-elastic measurement on PMMA or polycarbonate and compare b_l / b_λ.","The 5 °C spatial temperature spread reported at 65 °C could contribute systematic error; a direct check would be to measure the temperature profile along the ultrasonic path and propagate its uncertainty into the fitted slopes.","The frequency-independence of b_l and b_m combined with the strong frequency dependence of l and m suggests that relaxation processes shift in frequency with temperature; measuring over a wider frequency range or using broadband excitation could directly test the time-temperature superposition prediction.","The insensitivity of n to temperature and frequency may be characteristic of segmental-relaxation-driven coupling; this could be probed by testing polystyrene with different cross-linker content."],"forward_implications":["Temperature corrections to nonlinear elastic response of polystyrene (and likely similar glassy polymers) are dominated by the l and m terms; ignoring them in dynamic-load modeling could misestimate stress by a growing margin as temperature rises.","The near-zero temperature slope of n means shear-shear nonlinear coupling is stable across the studied range, potentially usable as a temperature-insensitive baseline in ultrasonic nonlinearity measurements.","Since the temperature susceptibility of l and m is frequency-independent while their absolute values depend strongly on frequency, models can separate a frequency-dependent modulus magnitude from a temperature-shift contribution.","The observed behavior is consistent with pressure shifting relaxation processes differently at different temperatures, and if the time-temperature superposition interpretation holds, measurements at higher temperatures can extend the effective frequency range of nonlinear modulus data."],"supporting_citations":[{"why":"Provides the acousto-elastic methodology and the equations relating pressure slopes of wave velocities to the Murnaghan moduli l, m, n.","marker":"[42]"},{"why":"Established the strong frequency dependence of l and m in polystyrene and the experimental setup that this paper extends to temperature.","marker":"[43]"},{"why":"Validates that reproducible ultrasonic velocity measurements in polystyrene can be made below 75 °C, defining the usable temperature range.","marker":"[44]"},{"why":"Reports a similar two-order ratio of temperature susceptibility between nonlinear and linear moduli in vitreous silica, used as a comparison for the new result.","marker":"[37]"},{"why":"Provides earlier measurements of elastic moduli and wave velocity versus pressure and temperature in plastics against which the linear temperature trends are checked.","marker":"[45]"}],"fun_headline_variants":["Polystyrene's nonlinear moduli heat up 100x faster","Temperature shifts nonlinear stiffness of polystyrene markedly","Nonlinear elastic moduli of polystyrene heat sensitivity 100x","Polystyrene nonlinear moduli are 100x more heat-sensitive"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The 5 °C spatial temperature variation inside the sample at the high end of the range is not propagated into the slopes, so the nominal temperature may differ from the actual temperature along the ultrasonic path by several degrees; if that bias is systematic, the reported temperature susceptibilities could be distorted.","fun_headline_variants_meta":{"raw":{"variants":["Polystyrene's nonlinear moduli heat up 100x faster","Temperature shifts nonlinear stiffness of polystyrene markedly","Nonlinear elastic moduli of polystyrene heat sensitivity 100x","Polystyrene nonlinear moduli are 100x more heat-sensitive"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000476,"raw_usage":{"total_tokens":2381,"prompt_tokens":988,"completion_tokens":1393,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":604,"completion_tokens_details":{"reasoning_tokens":1336}},"tokens_in":604,"tokens_out":1393,"duration_ms":10810,"temperature":1.0,"reasoning_tokens":1336,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T16:17:33.820411+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct check would be to repeat the pressure-velocity sweeps at 55 and 65 °C with the sample held in a temperature bath or with local temperature sensors along the ultrasonic path, keeping heterogeneity below about 0.5 °C; if the extracted b_l and b_m change by more than their reported ±0.06 and ±0.04 GPa/°C uncertainties, the two-orders claim is compromised.","supporting_citations":[{"cited_title":"Relative Variations of Nonlinear Elastic Moduli in Polystyrene-Based Nanocomposites","cited_arxiv_id":null,"evidence_quote":"Provides the acousto-elastic methodology and the equations relating pressure slopes of wave velocities to the Murnaghan moduli l, m, n."},{"cited_title":"Frequency Dependence of Nonlinear Elastic Moduli of Polystyrene","cited_arxiv_id":null,"evidence_quote":"Established the strong frequency dependence of l and m in polystyrene and the experimental setup that this paper extends to temperature."},{"cited_title":"Equation of State of Polystyrene and Polymethyl- methacrylate from Ultrasonic Measurements at Moderate Pressures","cited_arxiv_id":null,"evidence_quote":"Validates that reproducible ultrasonic velocity measurements in polystyrene can be made below 75 °C, defining the usable temperature range."},{"cited_title":"Temperature Dependences of the Third-Order Elastic Constants and Acoustic Mode Vibrational Anharmonicity of Vitreous Silica","cited_arxiv_id":null,"evidence_quote":"Reports a similar two-order ratio of temperature susceptibility between nonlinear and linear moduli in vitreous silica, used as a comparison for the new result."},{"cited_title":"Variation of Elastic Moduli and Wave Velocity with Pressure and Temperature in Plastics","cited_arxiv_id":null,"evidence_quote":"Provides earlier measurements of elastic moduli and wave velocity versus pressure and temperature in plastics against which the linear temperature trends are checked."}],"review_version":1}