{"id":"32ef605c-fd71-415e-b30c-365b9479b5e3","arxiv_id":"2507.01071","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":2,"one_line_summary":"Simulated creep of hierarchical Norway spruce models converges to a single proportionality factor across all directions and scales, so tissue topology alone cannot explain anisotropic creep; the authors argue non-linear material behavior is required.","lead":"This paper uses computer models of wood's cellular structure to test why Norway spruce creeps differently in different loading directions. The simulations collapse to a single elastic-to-viscoelastic ratio, so the authors conclude that tissue shape alone cannot explain the measured anisotropy and that non-linear cell wall responses are needed.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Identical KV spectra in every cell wall layer force the observed 'collapse'; the conclusion that topology cannot cause directional creep is an artifact of the shared-input assumption.","rationale":"The most load-bearing step in the argument is the assignment of identical KV parameters to every cell wall layer. The paper's strongest claim is a negative result about topology, but that negative result holds only within the class of models in which all constituents share a single time spectrum. Since this class is a set of measure zero within linear viscoelastic composites and the proportionality of effective compliance is guaranteed by linearity, the conclusion does not follow. The reader identified exactly this assumption; I agree. The proposed test—varying τ or γ in one layer while keeping everything else fixed—directly probes whether topology can generate directional creep when the input assumption is relaxed. If the test produces direction-dependent normalized creep, the paper's conclusion is falsified; if not, the concern is mitigated but the paper should still be revised to present the result as conditional on the shared-spectrum assumption. The rest of the paper—elastic validation, geometric modeling, and the inverse identification framework—is credible and useful, but it does not support the paper's strong conclusion.","tokens_in":15809,"tokens_out":4165,"duration_ms":48714,"concrete_test":"Rerun the Y-shaped RUC and growth-ring simulations with the S2 layer assigned a different KV spectrum (e.g., multiply the τ_i values by 10, or adopt a separate γ_i set for S2) while keeping all elastic constants, geometry, and other layers unchanged. Compute the normalized creep compliances Jc,AB(t)/C0,AB(t) under the nine load cases. If the relative deviations of γ_i,AB from their mean become comparable to the experimental spread in Fig. 8b, then the conclusion that topology cannot produce directional creep is invalidated. If the curves still collapse, the shared-spectrum assumption is less consequential.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central conclusion—that directional creep differences cannot emerge from tissue topology and therefore require nonlinear material behavior—is loaded into the model's input. In Sec. 2.1.2, all cell wall layers are assigned the same characteristic times τ_i = [0.1, 1, 10, 100] h, and Sec. 2.2.1 states that 'the ratios γve_i are identical for all layers.' Under this assumption each layer's creep compliance is a scalar multiple of its elastic compliance with the same time function. By the viscoelastic correspondence principle, any composite of such constituents has an effective creep compliance that is also a scalar multiple of its effective elastic compliance with that same time function, regardless of its topology or disorder. Consequently, the collapse of γve_i across all components and scales in Figs. 7–9 is a mathematical consequence, not an empirical finding about wood. The paper further claims this assumption 'was numerically verified before this study using RVEs with elastic fibers in a viscoelastic matrix,' but that verification does not test the load-bearing condition: that all layers share the same γ_i and τ_i. If layers had different KV spectra (e.g., different τ_i in the S2 layer vs. the compound middle lamella), the effective normalized creep could become direction-dependent even in linear viscoelasticity. The failure of the simulations to reproduce the experimental spread (Fig. 8b) therefore only demonstrates that the restricted model class cannot produce anisotropic creep; it does not support the conclusion that nonlinear material behavior is required. The strongest claim is thus an artifact of the identical-spectrum assumption.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a hierarchical multiscale finite-element model of Norway spruce, from individual cell wall layers (ML, P, S1, S2, S3) through regular Y-shaped unit cells and image-based disordered tissue RVEs up to a growth-ring RVE. Viscoelasticity is represented by a four-element Kelvin-Voigt chain in which each layer's creep compliance is proportional to its elastic compliance, with common proportionality factors gamma_i and fixed characteristic times tau_i = [0.1, 1, 10, 100] h. The gamma_i are identified inversely against macroscopic normalized creep measurements via a surrogate-based optimizer. The paper's central claim is that all simulated creep curves collapse onto a single universal set of proportionality factors regardless of tissue type, disorder, or geometric nonlinearity, and that therefore directional creep differences cannot originate from tissue topology but require nonlinear material behavior at stress concentration points.","tokens_in":16089,"tokens_out":4448,"duration_ms":55155,"significance":"If the central claim were established, the paper would justify a strong conclusion: that linear hierarchical models with a universal elastic-to-viscoelastic scaling cannot reproduce anisotropic creep of wood, so nonlinear cell wall behavior is necessary. The elastic part of the manuscript is a useful and largely convincing exercise: the predicted tissue and growth-ring elastic constants are compared with experimental and literature values, and the FEM framework with periodic boundary conditions is carefully described. However, the creep conclusion is not supported by the evidence. The collapse of the creep curves is a direct consequence of the constitutive assumption that every cell wall layer shares the same gamma_i and tau_i, and the paper does not test the alternative case in which different layers have different viscoelastic spectra. As submitted, the main scientific message about wood mechanics is therefore not established, although the multiscale modeling and identification pipeline could be valuable if the proportionality assumption is relaxed and the central question is revisited.","major_comments":[{"comment":"The model prescribes C_i^{-1} = C_0^{-1}/gamma_i and identical tau_i in every cell wall layer. Consequently each layer's compliance is C_0^{-1}(1 + sum_i gamma_i^{-1}(1 - exp(-t/tau_i))), i.e. the same scalar time function multiplies every layer's elastic compliance. For any linear elastic composite with fixed geometry and boundary conditions, scaling all constituent compliances by the same time function scales the effective compliance by that same function, independent of topology or disorder. The collapse shown in Figs. 7-9 is therefore a mathematical consequence of the input assumption, not an emergent numerical discovery. The statement in Sec. 2.2.1 that this assumption 'was numerically verified before this study using RVEs with elastic fibers in a viscoelastic matrix' refers to a different composite where fibers and matrix have different spectra and does not test the load-bearing condition that all layers share the same gamma_i and tau_i.","section":"Sec. 2.1.2, Eq. (4); Sec. 2.2.1"},{"comment":"The conclusion that 'directional differences do not emerge from topology, but must relate to non-linear material behavior at stress concentration points' is an overreach. The simulations only show that within the restricted model class in which every layer shares one gamma_i/tau_i spectrum, no directional creep emerges. If layers were assigned different spectra, for example different characteristic times for the S2 layer versus the compound middle lamella, topology could in principle produce direction-dependent effective creep even in linear viscoelasticity. The paper contains no simulations or analytical argument covering that case, so the central negative claim about topology is not supported by the presented evidence.","section":"Sec. 4, Conclusions"},{"comment":"The objective function in Eq. (9) compares normalized experimental creep curves for all AB components with the simulated normalized creep. Since the simulated curves collapse by construction to one master curve, the optimization cannot extract directional information from the experiments; it effectively fits a single curve. The resulting 'universal set' gamma_i is the same set that was prescribed as input for every layer in Sec. 2.2.1, so the agreement is circular rather than an independent identification. The paper should state clearly that the inverse procedure only calibrates the common spectrum, and cannot validate the proportionality assumption against the experimental spread shown in Fig. 8b.","section":"Eq. (9), Sec. 2.5"}],"minor_comments":[{"comment":"The positive-definiteness condition for the 3x3 subdeterminant is written as 'nu_12 nu_23 nu_31 < 0.5(1 - ...) < 0.5', which appears malformed; the intended inequalities should be corrected.","section":"Sec. 2.1.1"},{"comment":"The label '12mm' in the tissue-scale schematic is likely intended to be '12 um'; please verify the units.","section":"Fig. 1"},{"comment":"The text alternates between 'RUC' and 'RCU' (e.g., in Fig. 7 and Sec. 3.2.1); the terminology should be harmonized throughout.","section":"Sec. 3.2.1"},{"comment":"The phrase 'proving the analogy between elastic and viscoelastic deformation mechanisms' is stronger than justified; the simulations illustrate the analogy under the assumed proportionality but do not prove it for wood.","section":"Sec. 3.2.1"}],"recommendation":"reject","confidential_remarks":"The manuscript contains a potentially useful multiscale elastic modeling framework and an inverse identification pipeline, but the creep conclusion is predetermined by the constitutive assumption of identical Kelvin-Voigt spectra in all layers. The load-bearing issue is not a simple presentation fix: the central claim requires new simulations in which layers have different gamma_i or tau_i, and those simulations would likely change the conclusions. I would encourage the authors to rerun the analysis with layer-specific spectra and to frame the current result as a theorem about the restricted model class rather than a finding about wood. As it stands, the central claim is not publishable in its present form."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Hi,\n\nThe stress-test note is right, and it lands on reading the paper. The central conclusion—that tissue topology cannot produce directional creep, so nonlinear material behavior is required—is loaded into the model's input assumptions. In Sec. 2.1.2 all layers share the same τ_i = [0.1, 1, 10, 100] h, and Sec. 2.2.1 states that the ratios γ^ve_i are identical for all layers. Under that assumption, each layer's creep compliance is a scalar multiple of its elastic compliance with the same time function. The viscoelastic correspondence principle then forces any composite of such constituents to have effective creep compliance proportional to effective elastic compliance with the same time function, regardless of topology or disorder. So the universal collapse in Figs. 7–9 is a mathematical consequence, not an empirical discovery about wood.\n\nWhat the paper does well is the engineering framework. The hierarchical construction from cell wall layers to tissue RUCs/RVEs to a growth ring model is careful, the elastic predictions compare decently with meso- and macroscale measurements, and the surrogate-based inverse parameter identification is a practical contribution. The demonstration that the collapse holds across two scales is a useful consistency check for modelers who want to use the proportionality assumption.\n\nThe soft spot is that this consistency check is presented as a negative result about topology. The paper says the identical-spectrum assumption 'was numerically verified before this study using RVEs with elastic fibers in a viscoelastic matrix,' but that does not test the load-bearing condition: whether all layers share the same γ and τ. If the S2 layer and the compound middle lamella had different KV spectra, the normalized creep could become direction-dependent even in linear viscoelasticity. The failure of the simulations to reproduce the experimental spread in Fig. 8b therefore only shows that this restricted model class cannot produce anisotropy; it says nothing about wood itself. The final conclusion about nonlinear material behavior at cell corners is a reasonable hypothesis but is untested.\n\nThe paper is coherent on its own terms, and the authors are honest about several modeling simplifications, but they overstate the conclusion. I'd send it to a serious referee—the framework is worth publishing after the claim is properly conditionalized—but I would not cite the topology conclusion. If you need a reference for a hierarchical creep modeling workflow or the surrogate identification method, this is useful.","headline":"The paper's main conclusion is baked into its input assumptions; the useful part is the hierarchical modeling framework, not the negative result about topology.","tokens_in":16619,"tokens_out":2140,"would_cite":false,"duration_ms":24517,"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":"In hierarchical models of Norway spruce, the creep compliance tensor stays proportional to the elastic compliance tensor at every scale and in every tissue type, so the measured direction-dependence of creep must originate in nonlinear…","keywords":["orthotropic viscoelasticity","Norway spruce","wood creep","hierarchical cellular scaffold","Kelvin-Voigt chain","inverse parameter identification","surrogate optimization","microstructural disorder"],"falsifier":"Run the same hierarchical upscaling with layer-specific $\\gamma_i$ or $\\tau_i$ (or with a single layer whose creep is not proportional to its elasticity) and check whether the normalized creep curves $J_{c,AB}(t)/C^{-1}_{0,AB}$ separate by direction. A separation would show that topology alone can create orthotropic creep, contradicting the paper's conclusion that nonlinear stress-concentration behavior must be the cause.","tokens_in":15607,"feed_emoji":"🌲","tokens_out":10880,"duration_ms":110808,"temperature":0.7,"pith_summary":"The paper asks whether the direction-dependent creep measured in Norway spruce comes from its cellular architecture or from something in the cell-wall material itself. It builds hierarchical models that go from individual cell-wall layers to earlywood, transition wood, and latewood tissue and finally to a full growth ring, and it fits the unknown viscoelastic parameters of the cell-wall material against macroscopic creep curves. Under every topology tested, including geometrically nonlinear and image-based disordered models, the fitted creep-to-elasticity proportionality factors collapse onto one universal set for all components and all scales. Since real measurements show direction-dependent normalized creep, the paper concludes that tissue topology cannot produce this behavior and that nonlinear material response at stress concentration points such as cell corners must be the cause. If the conclusion is right, linear multi-scale models with a single scalar creep factor cannot predict experimentally observed orthotropic creep in wood.","feed_headline":"Cell shape alone cannot explain spruce's directional creep","feed_subtitle":"Models from cell wall to growth ring collapse onto one creep ratio; real direction-dependent creep must be nonlinear.","key_machinery":"The load-bearing object is a four-element Kelvin-Voigt chain whose creep compliance is $J_{c,AB}(t) = C^{-1}_{0,AB} \\sum_{i=1}^{4} (1/\\gamma_i^{\\mathrm{ve}})(1-e^{-t/\\tau_i})$, with fixed times $\\tau_i = [0.1, 1, 10, 100]$ hours and unknown dimensionless factors $\\gamma_i^{\\mathrm{ve}}$ linking each viscous element to the elastic compliance. The same $\\gamma_i^{\\mathrm{ve}}$ set is assigned to every cell-wall layer, so each upscaling step—from layer to Y-shaped cell to tissue to growth ring—preserves the proportionality between elastic and viscoelastic behavior. The paper then extracts component-wise $\\gamma_i$ from each model and shows that they collapse onto the input set, making the proportionality a property of the whole hierarchy rather than a single-scale assumption.","core_discovery":"The central result is that the effective viscoelastic compliance remains proportional to the effective elastic compliance at every hierarchical level of the model: cell wall layers, tissues, and the assembled growth ring. A single set of proportionality factors $\\gamma_i^{\\mathrm{ve}}$ with fixed characteristic times describes normalized creep in every direction, and this collapse is not broken by geometric nonlinearity or by replacing regular honeycomb cells with real disordered tissue images. The paper therefore asserts that the experimentally observed directional dependence of spruce creep does not emerge from tissue-scale topology. Instead, the missing directional dependence must come from nonlinear material behavior concentrated at stress concentration sites, where localized hinge-like deformation can develop.","pith_inferences":["If the collapse is generic rather than an artifact of the shared-spectrum assumption, then any measured directional spread in normalized creep is a direct experimental signal of nonlinearity, and the deviation of component-wise $\\gamma_i$ from the universal set could quantify the intensity of stress-concentration mechanisms.","The same proportionality-collapse should occur in other linear cellular composites provided all wall material has identical spectra; observing direction-dependent normalized creep in such a material would indicate nonlinearity.","Because the input assumption of identical $\\gamma_i$ and $\\tau_i$ in every layer pre-structures the result, a natural extension is to repeat the study with layer-specific spectra to establish whether compositional differences between S1, S2, and S3 layers change the conclusion.","A direct test is to measure creep of isolated earlywood tissue under radial versus tangential loading and compare component-wise normalized curves; non-coinciding curves would indicate that even a single tissue type carries nonlinear responses."],"forward_implications":["Linear multiscale wood models that assign one scalar creep-to-elasticity factor will not reproduce experimental orthotropic creep, so long-term deflection predictions for radial and tangential directions need nonlinear cell-wall input.","A single creep test in one direction can identify the whole linear viscoelastic spectrum of the growth ring, because all directions share the same proportionality factors in the model.","Adding geometric nonlinearity or realistic disordered cell geometry does not change the collapsed factors, so neither effect can rescue the topological explanation.","The next generation of wood creep models should place nonlinear constitutive laws at stress concentration sites such as cell corners, where hinge-like localized deformation can develop."],"supporting_citations":[{"why":"Supplies the macroscale normalized creep compliance measurements for all stress components that the inverse parameter identification targets and that show direction-dependent behavior.","marker":"Maas and Wittel (2025a)"},{"why":"Previous tissue-scale creep work that motivates the four-element Kelvin-Voigt chain and the fixed characteristic times used here.","marker":"Ferrara and Wittel (2025)"},{"why":"Provides the micromechanical cell-wall construction, layer properties, and tissue morphometry used to build the hierarchical models.","marker":"Ferrara et al. (2025)"},{"why":"Provides the Y-shaped repetitive unit cell geometry and periodic boundary condition framework for regular honeycomb tissue models.","marker":"Persson (2000)"},{"why":"Establishes the rheological wood model and the scalar proportionality assumption between elastic and viscoelastic compliances that the paper tests.","marker":"Hassani et al. (2015)"},{"why":"One of the standard references for assuming scalar proportionality between elastic and viscoelastic compliance in timber models.","marker":"Fortino et al. (2009)"},{"why":"Supplies the master-node technique used to impose periodic boundary conditions and extract the effective compliance tensors.","marker":"Rafsanjani et al. (2012)"}],"fun_headline_variants":["Spruce creep directionality resists tissue-shape explanation","Cell wall to ring: same creep ratio, but nonlinearity at stress sites matters","Why spruce creeps directionally: nonlinearity, not tissue shape","Creep in spruce: geometry collapses, nonlinearity drives direction"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Everything rests on assigning every cell-wall layer the same creep spectrum, with the same proportionality factors $\\gamma_i$ and the same characteristic times; if real layers differ in these parameters, tissue topology could in principle produce direction-dependent creep even in linear viscoelasticity.","fun_headline_variants_meta":{"raw":{"variants":["Spruce creep directionality resists tissue-shape explanation","Cell wall to ring: same creep ratio, but nonlinearity at stress sites matters","Why spruce creeps directionally: nonlinearity, not tissue shape","Creep in spruce: geometry collapses, nonlinearity drives direction"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000568,"raw_usage":{"total_tokens":2601,"prompt_tokens":769,"completion_tokens":1832,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":385,"completion_tokens_details":{"reasoning_tokens":1757}},"tokens_in":385,"tokens_out":1832,"duration_ms":14094,"temperature":1.0,"reasoning_tokens":1757,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T21:12:16.139208+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same hierarchical upscaling with layer-specific $\\gamma_i$ or $\\tau_i$ (or with a single layer whose creep is not proportional to its elasticity) and check whether the normalized creep curves $J_{c,AB}(t)/C^{-1}_{0,AB}$ separate by direction. A separation would show that topology alone can create orthotropic creep, contradicting the paper's conclusion that nonlinear stress-concentration behavior must be the cause.","supporting_citations":[{"cited_title":"PhD thesis","cited_arxiv_id":null,"evidence_quote":"Provides the Y-shaped repetitive unit cell geometry and periodic boundary condition framework for regular honeycomb tissue models."},{"cited_title":"Computer Methods in Applied Mechanics and Engineering 283, 1032–1060 (2015)","cited_arxiv_id":null,"evidence_quote":"Establishes the rheological wood model and the scalar proportionality assumption between elastic and viscoelastic compliances that the paper tests."},{"cited_title":"Mechanics of Time-Dependent Materials 13(4), 333–356 (2009)","cited_arxiv_id":null,"evidence_quote":"One of the standard references for assuming scalar proportionality between elastic and viscoelastic compliance in timber models."},{"cited_title":"Composites Science and Technology 72(6), 744–751 (2012) Villamizar Santamar´ ıa, S.: High-rise buildings made out of wood","cited_arxiv_id":null,"evidence_quote":"Supplies the master-node technique used to impose periodic boundary conditions and extract the effective compliance tensors."}],"review_version":1}