{"id":"051d1a3f-8894-49b3-91b9-2d90d739a441","arxiv_id":"1908.08935","paper_version":4,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":1,"one_line_summary":"Glass bead packs show clean, low-noise slow dynamic recovery after impulsive, harmonic, and quasi-static conditioning, with quasi-static conditioning reported for the first time.","lead":"This paper shows that unconsolidated glass bead packs, which are simpler than rocks or cement, display slow dynamic nonlinearity: a mechanical jolt softens the material and it recovers logarithmically in time. The authors demonstrate the effect with three conditioning methods, including one never reported before, using ultrasonic waves to track tiny stiffness changes.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (6) strain calibration is internally inconsistent: Table 1 and Fig. 8 use ε/F values that differ by factors of 10, so the strain-scaling and threshold conclusions in Sec. V rest on an unvalidated x-axis.","rationale":"The reader's conditional verdict is appropriate: the qualitative slow-dynamics observation is well supported by clean drops, log-time recoveries, and the foam/shim control experiments, but the quantitative strain-scaling and threshold claims rest on fragile inputs. The reader identified the post hoc T0 choice and the model-based strain estimates as the fragile inputs. I agree, and add a more specific and checkable problem: the numbers in Eq. (6), Table 1, and Fig. 8 do not self-consistently follow from the stated physical parameters. This is not a question of a slightly uncertain model parameter; it is an apparent order-of-magnitude inconsistency in the calibration of the independent variable of the paper's main quantitative figure. Resolving it is necessary before the linear-regime and threshold conclusions in Sec. V are taken as established. The central qualitative claim, however, does not depend on the absolute strain calibration, so the verdict should remain conditional rather than being upgraded or downgraded. The recommended concrete test is a direct recomputation plus an independent strain measurement, which would settle whether the factor-of-ten issue is a typographical error or a substantive calibration flaw.","tokens_in":13712,"tokens_out":17051,"duration_ms":168862,"concrete_test":"Recompute every entry in Table 1 and Fig. 8 directly from Eq. (6) with the stated parameters L = 33 mm, m_pack = 221 g, c = 825 m/s and the reported forces (F_imp = 128 N, F_harm = 2.8 N, F_QS = 9.81 N). If the resulting strains differ from the plotted values by a factor of 10, remeasure the pump strain independently—for example with a strain gauge on the bead pack or a displacement sensor on the load—and refit the slope-versus-strain relation in Fig. 8b using corrected x-axis values and full fit uncertainties, including the sensitivity to T0 choice.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central qualitative observation—log-time recovery after impulsive, harmonic, and quasi-static conditioning—is credible and survives this concern. The weakest load-bearing point is the calibration of the conditioning strain used for the quantitative claims in Sec. V and Fig. 8. Eq. (6) states ε = L F / (m_pack c^2) with L = 33 mm, m_pack = 221 g, and c = 825 m/s; this gives ε/F = 0.033 / (0.221 × 825^2) = 2.19×10^-7 N^-1. Using that coefficient with the reported forces yields: quasi-static ε = 2.15×10^-6 for F = 9.81 N, harmonic ε = 6.15×10^-7 for F = 2.8 N, and impulsive ε = 2.8×10^-5 for F = 128 N. The values quoted in Table 1 and Fig. 8 instead correspond to coefficients of roughly 2.2×10^-6 N^-1 (harmonic), 2.2×10^-7 N^-1 (quasi-static), and 2.2×10^-8 N^-1 (impulsive), i.e. mutually inconsistent by factors of 10 depending on which entry is read at face value. Because the x-axis of Fig. 8b is the rms pump strain, an order-of-magnitude error in this axis would directly shift the claimed linear-above-one-microstrain regime and the proposed sub-6×10^-7 threshold. The secondary conclusion therefore cannot be accepted until the strain calibration is independently validated and the apparent arithmetic inconsistency in Eq. (6) is resolved.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents an experimental study of slow dynamic nonlinearity in unconsolidated glass bead packs, using ultrasonic pulses and coda wave interferometry (CWI) to measure changes in wave speed (stretch). The authors apply three low-frequency conditioning methods: impulsive (dropped ball), harmonic (60 Hz shaker), and quasi-static (adding/subtracting a 1 kg mass). They report a drop in stretch followed by logarithmic recovery after all three pump types, with the quasi-static method identified as novel. They also report a preliminary dependence of the recovery slope on the estimated rms pump strain, finding an apparent linear regime above about one microstrain and a leveling off below it, and use this to suggest that any threshold strain for slow dynamics in this system is below 6e-7 if it exists at all. Control experiments with brass side walls and without rubber shims indicate that the observed effects originate in the bead pack rather than in the support apparatus.","tokens_in":14010,"tokens_out":7042,"duration_ms":71932,"significance":"If the central qualitative claim holds, the paper provides a clean, low-noise venue for studying slow dynamics, with the advantage of a simpler and better-characterized microstructure than rocks or cement. The demonstration of slow dynamics after quasi-static conditioning is a new experimental result. The CWI-based stretch measurements achieve precision better than about 1e-6, and the control experiments for the foam and rubber shims strengthen the attribution of the effect to the bead pack. These are genuine strengths and make the system potentially useful for future parametric studies. The quantitative strain-threshold conclusions in Sec. V are more fragile: they depend on estimates of the conditioning strain and on slopes fitted without reported uncertainties, so they should be treated as preliminary unless the calibration is independently validated.","major_comments":[{"comment":"The quantitative claims in Sec. V—especially the linear-above-one-microstrain behavior and the 'below 6e-7' threshold statement—rest on the absolute values of the estimated pump strain. The strain values in Table 1 are actually consistent with the coefficient in Eq. (6): using epsilon/F = 2.19e-7 N^-1 gives 2.8e-5 for F = 128 N, 6.15e-7 for F = 2.8 N, and 2.15e-6 for F = 9.81 N, so there is no internal factor-of-ten inconsistency. However, Eq. (6) assumes a uniform uniaxial strain over the whole pack with a single wave speed c = 825 m/s, and Eq. (7) models the 87 kg load plus pack as a single-degree-of-freedom spring-mass system with a resonance at 20 Hz. These are reasonable order-of-magnitude estimates, but they are not independently validated. Because the x-axis of Fig. 8b and the numerical value of the putative threshold depend directly on this calibration, the paper should either provide an independent check (e.g., a direct displacement measurement of the load or pack boundary, or a comparison with static deflection under the added mass) or explicitly frame the strain axis and the threshold as order-of-magnitude estimates with a stated uncertainty. As written, the threshold conclusion in Sec. V is stronger than the calibration supports.","section":"Sec. IV, Eq. (6), Eq. (7); Table 1; Fig. 8"},{"comment":"The zero time T0 for each recovery is adjusted post hoc to make the log-time plots linear, as acknowledged in footnote 5. The recovery slopes used in Table 1 and Fig. 8b are obtained by fitting from 15 s to 3 min, but no error bars or sensitivity analysis are reported. The slope values are the central quantitative output of the paper, and the apparent linear-versus-plateau structure in Fig. 8b is the basis for the threshold discussion. The authors should report at least the variation of the fitted slopes as T0 is varied within the plausible three-second interval, and ideally provide confidence intervals from the fits. Without this, it is difficult to assess whether the apparent leveling off at low strain in Fig. 8b is significant or within the fitting uncertainty.","section":"Sec. IV.A, footnote 5; Table 1; Fig. 8b"},{"comment":"The low-strain plateau in Fig. 8b is inferred from very few data points, and footnote 7 states that measurements below the lowest strain were contaminated by drift. If the plateau is defined by only one or two points, the claim that the slope 'levels off' below about one microstrain should be worded more cautiously. Additional data points at intermediate strains, or a clear statement of how many independent measurements define each point, would make the threshold discussion more convincing. As it stands, the distinction between a genuine plateau and a scatter-limited flattening is not established.","section":"Sec. V, Fig. 8"}],"minor_comments":[{"comment":"The phrase 'a creep-like phenomena' should be 'a creep-like phenomenon' (singular) or 'creep-like phenomena' without the indefinite article.","section":"Abstract"},{"comment":"Some of the superscript and subscript characters are garbled in the transcript (e.g., the density unit 'g/cm1' should be 'g/cm^3', and several strain exponents such as '10#$' appear as raw symbols). Please ensure the final typeset version uses correct superscripts; this is a readability issue that should not affect the physics.","section":"Sec. II, Eq. (2) and surrounding text"},{"comment":"The table lists slopes for the harmonic conditioning with 'shaker off' and 'shaker on'. It would help to state explicitly in the caption whether the 'shaker on' slope refers to the conditioning phase or to the recovery after the shaker is turned off, and to define the sign convention used for both.","section":"Table 1"},{"comment":"Panel (b) would benefit from error bars or at least a statement of the estimated uncertainty in each slope value. As currently plotted, the visual impression of a linear regime and a plateau may overstate the precision of the measurements.","section":"Fig. 8"},{"comment":"Reference [24] appears to be a thesis citation without a title or institution; if it is publicly available, adding the full reference would improve completeness.","section":"References"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper's core observation is solid: unconsolidated glass bead packs show clean log-time recovery in ultrasonic stretch after impulsive, harmonic, and quasi-static conditioning, and the quasi-static protocol is genuinely new. The experimental care is real—floating walls, dead-weight load, foam controls, brass-cylinder and rubber-shim checks—and the CWI stretch data look far less irregular than earlier bead-pack slow-dynamics work. I buy the qualitative result.\n\nThe quantitative side is weaker, exactly as the reader's report says. The recovery slopes have no error bars, T0 is adjusted post hoc to force log-linearity (footnote 5), and the strain-dependence curve in Fig. 8 rests on five points, with the lowest strain point near the noise floor. The claim that the threshold is below 6e-7 if it exists at all is plausible but not established. The strain estimates themselves come from a simplified model (Eq. 6), and an absolute error there would shift the x-axis of Fig. 8b, though the relative ordering of points would survive.\n\nI checked the stress-test note's arithmetic inconsistency. It does not hold up. Eq. (6) with L=0.033 m, m=0.221 kg, c=825 m/s gives epsilon/F = 2.19e-7 N^-1. Plugging in the reported forces yields exactly the Table 1 strains: 2.8e-5 for 128 N, 6.15e-7 for 2.8 N, 2.15e-6 for 9.81 N. All three entries are mutually consistent. The stress-test's claim of factors-of-ten discrepancy comes from dividing the other way or misreading superscripts. So that specific concern is a red herring.\n\nThe citation pattern is fine; the Hertz cutoff and quasi-static stretch comparisons are genuine theory-data checks, and the CWI method is externally validated. The paper does not overclaim: it explicitly flags T0 as post hoc and the temperature-drift caveat on the lowest strain point is honest.\n\nWho gets value from this: anyone working on nonlinear mesoscopic elasticity, granular acoustics, or seismic wave speed recovery. It deserves a serious referee. The referee should push for uncertainty quantification on the slopes, an independent T0 measurement, and more strain points, but the central observation is a useful contribution. I'd cite it for the quasi-static result and the clean bead-pack venue.","headline":"Careful bead-pack study convincingly shows slow dynamics after three pump types, including a new quasi-static protocol; the quantitative strain-dependence claims are shakier, and the stress-test's arithmetic inconsistency does not hold up.","tokens_in":14586,"tokens_out":2720,"would_cite":true,"duration_ms":25769,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Unconsolidated glass bead packs exhibit slow dynamic nonlinearity after impulsive, harmonic, and quasi-static conditioning, with quasi-static pumping demonstrated for the first time.","keywords":["slow dynamics","glass bead packs","coda wave interferometry","nonlinear elasticity","logarithmic recovery","strain threshold","granular media"],"falsifier":"Record the pump onset with an independent trigger and refit the recoveries; if the slopes change beyond stated precision, the log-time fits are artifacts of the chosen $T_0$. Mount a strain gauge in or on the bead pack to measure conditioning strain directly rather than through Eqs. (6) and (7); if the leveling-off strain is not near $10^{-6}$, the threshold claim fails.","tokens_in":13468,"feed_emoji":"🔬","tokens_out":7213,"duration_ms":71258,"temperature":0.7,"pith_summary":"This paper argues that a simple laboratory material—loose glass beads packed under a dead-weight load—can serve as a clean test bed for slow dynamics, the poorly understood loss of stiffness followed by logarithmic recovery seen in rocks and concrete. Using ultrasound and coda wave interferometry to detect tiny wave-speed changes, the authors show that the bead pack recovers logarithmically after three kinds of low-frequency conditioning: a dropped ball, a 60 Hz shaker, and the quasi-static addition or removal of a 1 kg mass. The quasi-static conditioning result, they report, has not appeared in the literature before. The paper also finds that the recovery rate scales linearly with the root-mean-square pump strain above about one microstrain and levels off below it, suggesting that any threshold strain for slow dynamics in this system lies below $6\\times10^{-7}$, if a sharp threshold exists at all. If these claims hold, bead packs give theorists a simpler, better-characterized microstructure in which to test mechanisms of slow dynamics.","feed_headline":"Loose glass beads heal logarithmically like rocks","feed_subtitle":"Slow, log-time recovery appears even under quasi-static load changes, giving theorists a simpler material to study.","key_machinery":"The load-bearing measurement is coda wave interferometry (CWI) applied to multiply scattered ultrasonic waves. A broadband pulse is sent through the bead pack every 0.01 s, 100 received signals are averaged every few seconds, and each averaged signal is cross-correlated with a reference signal in short signal-time windows. The time delay at maximum correlation is plotted against signal-time and fit to a line $\\mathcal{T}_{ni} = \\sigma_n t_i$; the slope, called stretch, equals the relative change in wave speed and therefore the relative change in modulus. Because the ultrasonic coda has explored many bead contacts, tiny stiffness changes show up as a measurable stretch. The three conditioning methods supply the perturbation, and the stretch recovery curves provide the slow-dynamics signature.","core_discovery":"The central claim is that unconsolidated glass bead packs exhibit reproducible slow dynamic nonlinearity when probed by diffuse ultrasound, with a stiffness drop during low-frequency mechanical conditioning and a recovery that is linear in the logarithm of time. After impulsive, harmonic, and quasi-static pumping, the measured coda-wave stretch $\\sigma$ follows $\\sigma = m \\ln(T - T_0) + \\text{const.}$, with slopes of order $10^{-6}$ per log cycle. The quasi-static pump—adding or removing a 1 kg weight—produces slow logarithmic relaxation after both loading and unloading, showing that the system relaxes toward equilibrium regardless of the sign of the imposed strain. In addition, the recovery slope $m$ is linear in the rms pump strain for strains at and above roughly $10^{-6}$, and it flattens at lower strains, which the paper reads as evidence that the threshold for slow dynamics is much lower than previously estimated for bead packs, or absent.","pith_inferences":["One testable extension is to vary humidity or apply hydrophobic bead coatings and watch whether the log-time recovery rate shifts; that would discriminate moisture-assisted contact healing from purely mechanical creep.","A longer-duration, temperature-controlled quasi-static recovery measurement could test whether the astronomically long extrapolated recovery reflects an unusually flat activation-energy distribution or a genuinely different relaxation process.","If the absence of a sharp strain threshold survives direct strain measurement, then small ambient vibrations may continuously excite slow dynamics in granular media, meaning even quiescent measurements of modulus in such materials are already affected by the measurement process."],"forward_implications":["A simpler material is now available for testing theories of slow dynamics, with controllable grain size, surface treatment, humidity, and saturation.","Quasi-static conditioning extends slow dynamics to a pump mode relevant to tectonic loading, since adding and removing a load both produce logarithmic recovery.","The linear dependence of recovery slope on rms strain above about one microstrain and the leveling off below it indicate that the threshold strain in bead packs, if any, is below $6\\times10^{-7}$.","Ultrasonic probing via CWI can detect slow dynamics even when resonant frequencies do not shift, suggesting that reported threshold strains in sandstones should be re-examined at lower strains.","The paper notes that extrapolated full-recovery times differ enormously by pump type—about 50 hours for the dropped ball versus effectively infinite for the quasi-static case—and leaves this discrepancy as an open puzzle."],"supporting_citations":[{"why":"Establishes the prior result that recovery slope scales linearly with rms pump strain and levels off at low strain, the comparison for this paper's Fig. 8.","marker":"[2]"},{"why":"Shows impulsive conditioning and early-time log recovery in sandstone and cement paste, motivating the ball-drop pump.","marker":"[11]"},{"why":"Uses impact conditioning with coda wave interferometry in concrete, the proximate template for the pump-probe scheme.","marker":"[12]"},{"why":"Reports prior slow dynamics in glass bead packs by resonant spectroscopy and estimates a threshold strain near one microstrain, the claim this paper challenges.","marker":"[31]"},{"why":"Introduces coda wave interferometry and the stretch parameter used throughout the paper.","marker":"[33]"},{"why":"Provides comparable diffuse-wave transport measurements in bead packs, grounding the ultrasonic characterization.","marker":"[39]"},{"why":"Reports a threshold strain for Berea sandstone, used as reference for the threshold discussion.","marker":"[47]"},{"why":"Reports a threshold strain for Fontainebleau sandstone, used as reference for the threshold discussion.","marker":"[48]"}],"fun_headline_variants":["Bead packs show log-time recovery after stress","Ultrasound tracks slow healing in glass beads","Simple beads reveal slow dynamics like rocks","Log-time stiffness recovery in bead packs","Quasi-static loads trigger slow bead relaxation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper's quantitative conclusions would collapse if direct strain measurement showed that the conditioning strains estimated from the simplified force-to-strain relations are systematically wrong, since the linear-above-one-microstrain and no-sharp-threshold claims are built on those estimates.","fun_headline_variants_meta":{"raw":{"variants":["Bead packs show log-time recovery after stress","Ultrasound tracks slow healing in glass beads","Simple beads reveal slow dynamics like rocks","Log-time stiffness recovery in bead packs","Quasi-static loads trigger slow bead relaxation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.0002,"raw_usage":{"total_tokens":1359,"prompt_tokens":913,"completion_tokens":446,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":529,"completion_tokens_details":{"reasoning_tokens":380}},"tokens_in":529,"tokens_out":446,"duration_ms":5221,"temperature":1.0,"reasoning_tokens":380,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:24:41.901222+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Record the pump onset with an independent trigger and refit the recoveries; if the slopes change beyond stated precision, the log-time fits are artifacts of the chosen $T_0$. Mount a strain gauge in or on the bead pack to measure conditioning strain directly rather than through Eqs. (6) and (7); if the leveling-off strain is not near $10^{-6}$, the threshold claim fails.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the prior result that recovery slope scales linearly with rms pump strain and levels off at low strain, the comparison for this paper's Fig. 8."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Shows impulsive conditioning and early-time log recovery in sandstone and cement paste, motivating the ball-drop pump."},{"cited_title":"Tremblay, E","cited_arxiv_id":null,"evidence_quote":"Uses impact conditioning with coda wave interferometry in concrete, the proximate template for the pump-probe scheme."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports prior slow dynamics in glass bead packs by resonant spectroscopy and estimates a threshold strain near one microstrain, the claim this paper challenges."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces coda wave interferometry and the stretch parameter used throughout the paper."},{"cited_title":"Jia, Phys","cited_arxiv_id":null,"evidence_quote":"Provides comparable diffuse-wave transport measurements in bead packs, grounding the ultrasonic characterization."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports a threshold strain for Berea sandstone, used as reference for the threshold discussion."},{"cited_title":"Pasqualini, K","cited_arxiv_id":null,"evidence_quote":"Reports a threshold strain for Fontainebleau sandstone, used as reference for the threshold discussion."}],"review_version":1}