pith. sign in

arxiv: 2606.04256 · v1 · pith:HIYDDOCZnew · submitted 2026-06-02 · ❄️ cond-mat.soft

Morphogenesis driven by nematic defects in active biological networks

Pith reviewed 2026-06-28 07:46 UTC · model grok-4.3

classification ❄️ cond-mat.soft
keywords morphogenesisnematic defectsactive biological networksHydra regenerationtopological defectsstress localizationtissue remodelinggrowth and relaxation
0
0 comments X

The pith

Topological defects in nematic tissues control where protrusions form and growth occurs during morphogenesis.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper builds a continuum model of tissues as active nematic polymer networks whose order parameter varies in space and drives the system through stress and growth. Morphogenesis is cast as a sequence of quasi-static states in which elasticity, stress-driven growth, and relaxation are coupled. Simulations prescribe defects to match an expected Hydra morphology and find that +1 defects localize stress to create protrusions while -1/2 defects produce little growth and hold structure steady. Varying the starting defect pattern then yields different outcomes such as tentacles or biaxial forms. A reader would care because the model supplies a direct mechanical route from defect placement to final shape.

Core claim

The central claim is that defect topology controls stress localization and shape evolution: +1 defects drive protrusion formation, while -1/2 defects act as structural stabilizers with minimal growth. By varying the initial defect configuration, the model produces diverse morphogenetic outcomes, including uniaxial regeneration, tentacle formation, and biaxial development.

What carries the argument

The nematic order parameter together with its topological defects inside a continuum description that couples order to elasticity, stress-driven growth, and adaptive relaxation.

If this is right

  • Prescribing +1 defects produces localized protrusions through stress concentration.
  • -1/2 defects produce minimal growth and maintain structural stability.
  • Different initial defect placements generate uniaxial, tentacled, or biaxial regeneration patterns.
  • Heterogeneous growth follows directly from spatial variations in the order parameter.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the mechanism holds, experimental placement of defects could be used to steer regeneration outcomes in living tissues.
  • The same stress-defect coupling might operate in other epithelial or cytoskeletal systems that exhibit nematic order.
  • A natural extension would test whether defects arise spontaneously or must be seeded to match the observed patterns.

Load-bearing premise

The model assumes that initial defect configurations can be prescribed to match the mature organism's expected morphology and that the resulting quasi-static equilibrium sequence will reproduce observed regeneration outcomes.

What would settle it

Direct imaging showing that +1 defects are absent from future protrusion sites in real Hydra, or that -1/2 defects do not correlate with stable regions, would falsify the claimed control by defect topology.

Figures

Figures reproduced from arXiv: 2606.04256 by Giulio G. Giusteri, L. Angela Mihai, Silvia Paparini.

Figure 1
Figure 1. Figure 1: Image of nematic actin fiber organization in the ectoderm of a small, mature Hydra [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Schematic of the nematic organization of actin fibers at successive stages of regenera [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Schematics of the multiplicative decomposition of [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Figures (a)-(c) reprinted or adapted from [19] and [22]. [PITH_FULL_IMAGE:figures/full_fig_p010_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: N = 25, s = 0.4, κ = 0.01. Successive equilibrium configurations illustrat￾ing the quasi-static morphological deformation of the Hydra tissue fragment (from left to right) for steps k ∈ {0, 3, 5, 8, 15, 25}. The corresponding parameters take the values s (k) 0 ∈ {0, 0.04, 0.1, 0.14, 0.25, 0.4} and ε (k) ∈ {π/2, 1.34, 0.578, 0.121, 0, 0}, respectively. The sequences show the stepwise expansion of the active… view at source ↗
Figure 6
Figure 6. Figure 6: N = 25, s = 0.4, κ = 0.01. Apical view (prospective mouth). At each simulation step k, the colormap shows the spatial distribution of local mechanical stress. A pronounced high-stress halo is observed to propagate toward the apical pole, closely following the outward expansion of the active nematic band. In our stress-driven growth framework, new material is preferentially incorporated in regions of elevat… view at source ↗
Figure 7
Figure 7. Figure 7: N = 25, s = 0.4, κ = 0.01. Basal view (prospective foot). At each simulation step k, the colormap shows the distribution of local mechanical stress. A pronounced high-stress halo propagates toward the basal pole, closely following the outward expansion of the active nematic band. Within the stress-driven growth framework, new cells are preferentially incorporated in regions of elevated mechanical stress. H… view at source ↗
Figure 8
Figure 8. Figure 8: Hydra with two heads configuration: (a) stereographic projection showing the planar [PITH_FULL_IMAGE:figures/full_fig_p024_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: N = 50, s = 0.6, κ = 0.01. Successive equilibrium configurations illustrating the quasi-static morphological deformation of the Hydra bicephalic configuration (from left to right) for steps k ∈ {0, 10, 20, 30, 40, 50}. The corresponding parameters take the values s (k) 0 ∈ {0, 0.12, 0.24, 0.36, 0.48, 0.6} and δ (k) ∈ {5, 30.4, 60.3, 90.2, 120, 150}, respectively. The color map displays the logarithm of the… view at source ↗
Figure 10
Figure 10. Figure 10: N = 50, s = 0.6. Apical view. At each simulation step k, the color map shows the spatial distribution of local mechanical stress Ψ(k) e . High-stress regions locate around the heads and the foot. At the apical poles, the stress generated by nematic distortions is progressively relaxed, leading to the gradual enlargement of the oral apertures, while the foot remains the most stressed region. The simulation… view at source ↗
Figure 11
Figure 11. Figure 11: N = 50, s = 0.6, k = 20. Magnified view of the actin fiber organization at the branching junction between the two heads and the main trunk. 26 [PITH_FULL_IMAGE:figures/full_fig_p026_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Hydra with tentacles configuration: (a) stereographic projection showing planar [PITH_FULL_IMAGE:figures/full_fig_p028_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: N = 80, s = 0.6, κ = 0.01. Successive equilibrium configurations illustrat￾ing the quasi-static morphological deformation of the Hydra tissue fragment (from left to right) for steps k ∈ {0, 16, 32, 48, 64, 80}. The corresponding parameters take the values s (k) 0 ∈ {0, 0.12, 0.24, 0.36, 0.48, 0.6} and δ (k) ∈ {5, 34, 63, 92, 121, 150}, respectively. The sequence shows successive equilibrium configurations… view at source ↗
Figure 14
Figure 14. Figure 14: N = 80, s = 0.6. Apical view. The right half of the regenerating tissue displays the director field tangent to the inner and outer surfaces of the shell. The color map represents the elastic energy density Ψ(k) e , providing a direct measure of the local stress distribution. During the initial steps, high-stress regions are observed, localized around the head and the tentacles. As the simulation progresse… view at source ↗
Figure 15
Figure 15. Figure 15: N = 80, s = 0.6. The sequence shows a detailed view of the director field around the frontal −1/2 defects, color-coded by the local elastic stress distribution Ψ(k) e . As regeneration proceeds, the characteristic three-fold nematic configuration becomes progressively stretched along the direction of the head and the adjacent tentacle. The localization of elastic energy around the −1/2 defects at advanced… view at source ↗
read the original abstract

Cellular morphogenesis, the process by which biological tissues acquire shape and structure, remains a fundamental challenge in understanding pattern formation and the coordinated remodeling of cellular assemblies. Under appropriate conditions, cytoskeletal filaments can organize into a nematic phase exhibiting partial orientational order. Topological defects within this nematic organization generate localized mechanical stresses that destabilize the tissue and promote deformation and structural rearrangements to relieve internal stresses. We develop a continuum framework that models living tissues as active biological networks represented as nematic polymer networks capable of heterogeneous growth and remodeling. The model captures macroscopic effects through spatial variations in the fiber order parameter which drives the system away from equilibrium. Morphogenesis is described as a sequence of quasi-static equilibrium states governed by the coupling between nematic order, elasticity, stress-driven growth, and adaptive relaxation. Finite element simulations illustrate Hydra regeneration and development when topological defects are prescribed according to the mature organism's expected morphology. The results show that defect topology controls stress localization and shape evolution: $+1$ defects drive protrusion formation, while $-1/2$ defects act as structural stabilizers with minimal growth. By varying the initial defect configuration, we model diverse morphogenetic outcomes, including uniaxial regeneration, tentacle formation, and biaxial development.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 1 minor

Summary. The manuscript develops a continuum model treating living tissues as active nematic polymer networks with heterogeneous growth and remodeling. Morphogenesis is represented as a sequence of quasi-static equilibria arising from the coupling of nematic order, elasticity, stress-driven growth, and adaptive relaxation. Finite-element simulations of Hydra regeneration and development are performed with topological defects (+1 and -1/2) prescribed in the initial order-parameter field to match the mature organism's expected morphology; the results indicate that +1 defects localize stress to drive protrusions while -1/2 defects stabilize regions with minimal growth, and that varying the prescribed defect pattern reproduces uniaxial regeneration, tentacle formation, and biaxial outcomes.

Significance. The framework supplies a concrete mechanical link between nematic defect topology and localized growth in a specific biological system. The finite-element implementation permits systematic exploration of different initial defect configurations and their effect on shape evolution. If the defects were instead shown to emerge dynamically from the active-network equations and if the model outputs were compared quantitatively to experimental regeneration data, the approach could strengthen mechanistic understanding of defect-driven morphogenesis.

major comments (3)
  1. [Abstract] Abstract: the claim that 'defect topology controls stress localization and shape evolution' and that '+1 defects drive protrusion formation' rests on simulations in which initial defect positions are prescribed by hand to match the target mature morphology. This setup demonstrates consistency (prescribed defects produce the expected shapes) but does not establish that defects arise spontaneously from the active-network dynamics without foreknowledge of the final pattern.
  2. [Abstract] Abstract (model description): the quasi-static equilibrium sequence explicitly decouples defect motion from growth and remodeling. Because defect positions remain fixed inputs rather than evolving quantities, the reported stress localization cannot be attributed to a self-consistent dynamical feedback between activity, order, and shape change.
  3. [Abstract] Abstract: no validation data, mesh-convergence tests, parameter-sensitivity analysis, or direct quantitative comparison to experimental Hydra regeneration timelines or morphologies are supplied. The central illustrative claim therefore depends on unexamined modeling choices whose robustness is not demonstrated.
minor comments (1)
  1. [Abstract] Abstract: the statement that the model 'captures macroscopic effects through spatial variations in the fiber order parameter' would benefit from an explicit statement of the free-energy or constitutive relation used for the nematic-elastic coupling.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive and detailed comments. We agree that the model prescribes defect positions and employs a quasi-static approximation without dynamic defect evolution or quantitative experimental validation. We will revise the abstract and add supporting numerical tests to clarify the model's scope and assumptions.

read point-by-point responses
  1. Referee: [Abstract] Abstract: the claim that 'defect topology controls stress localization and shape evolution' and that '+1 defects drive protrusion formation' rests on simulations in which initial defect positions are prescribed by hand to match the target mature morphology. This setup demonstrates consistency (prescribed defects produce the expected shapes) but does not establish that defects arise spontaneously from the active-network dynamics without foreknowledge of the final pattern.

    Authors: We agree that defects are prescribed in the initial order-parameter field to match expected morphology, as explicitly stated in the abstract. The simulations demonstrate mechanical consistency between prescribed defect topology and resulting shapes but do not address spontaneous emergence. We will revise the abstract to replace 'controls' with 'influences' and to state explicitly that defects are prescribed inputs. revision: yes

  2. Referee: [Abstract] Abstract (model description): the quasi-static equilibrium sequence explicitly decouples defect motion from growth and remodeling. Because defect positions remain fixed inputs rather than evolving quantities, the reported stress localization cannot be attributed to a self-consistent dynamical feedback between activity, order, and shape change.

    Authors: The referee correctly notes that defect positions are held fixed during the sequence of quasi-static equilibria. This is a deliberate modeling choice to isolate the mechanical consequences of defect topology. We will revise the abstract and model description to emphasize the fixed-defect, quasi-static approximation and its limitations. revision: yes

  3. Referee: [Abstract] Abstract: no validation data, mesh-convergence tests, parameter-sensitivity analysis, or direct quantitative comparison to experimental Hydra regeneration timelines or morphologies are supplied. The central illustrative claim therefore depends on unexamined modeling choices whose robustness is not demonstrated.

    Authors: We will add mesh-convergence tests and a parameter-sensitivity study for activity strength and relaxation timescales in the revised manuscript. Direct quantitative comparison to experimental regeneration timelines is not performed in the current work, which remains qualitative; we will note this as a limitation and direction for future study. revision: partial

Circularity Check

1 steps flagged

Defect positions prescribed to match target morphology make outcomes by construction

specific steps
  1. fitted input called prediction [Abstract]
    "Finite element simulations illustrate Hydra regeneration and development when topological defects are prescribed according to the mature organism's expected morphology. The results show that defect topology controls stress localization and shape evolution: +1 defects drive protrusion formation, while -1/2 defects act as structural stabilizers with minimal growth. By varying the initial defect configuration, we model diverse morphogenetic outcomes, including uniaxial regeneration, tentacle formation, and biaxial development."

    Defect configurations are chosen to match the target morphology in advance; the reported shape evolution is then the direct output of evolving equilibria under that fixed prescription, so the claimed control by defect topology reduces to a consistency check on the inputs rather than an independent derivation.

full rationale

The central claim that defect topology controls morphogenesis is demonstrated only after the model explicitly sets initial defect locations by hand to reproduce the mature organism's expected morphology. The subsequent quasi-static evolution then necessarily produces the corresponding shapes (protrusions at +1 defects, stabilization at -1/2), rendering the reported control a direct consequence of the prescribed inputs rather than an emergent prediction from the active network dynamics. This matches the fitted-input-called-prediction pattern and accounts for the observed partial circularity; the remainder of the continuum framework appears self-contained.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract supplies no information on free parameters, axioms, or invented entities; all fields left empty.

pith-pipeline@v0.9.1-grok · 5747 in / 1022 out tokens · 18579 ms · 2026-06-28T07:46:13.980670+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

41 extracted references · 33 canonical work pages

  1. [1]

    Bladon P, Terentjev EM, Warner M. 1994. Deformation-induced orientational tran- sitions in liquid crystals elastomer, Journal of Physics II France 4(1), 75–91 (doi:10.1051/jp2:1994100)

  2. [2]

    Bailles A, Serafini G, Andreas H, Zechner C, Modes CD, Tomancak P. 2025. Anisotropic stretch biases the self-organization of actin fibers in multicellular Hydra ag- gregates, Proceedings of the National Academy of Sciences 122(32), e2423437122 (doi: 10.1073/pnas.2423437122)

  3. [3]

    Blanch-Mercader C, Yashunsky V, Garcia S, Duclos G, Giomi L, Silberzan P. 2018. Tur- bulent dynamics of epithelial cell cultures, Physical Review Letters 120, 208101 (doi: 10.1103/PhysRevLett.120.208101)

  4. [4]

    Doostmohammadi A, Ign´ es-Mullol J, Yeomans JM, Sagu’es F. 2018. Active nematics, Na- ture Communication 9, 3246 (doi: 10.1038/s41467-018-05666-8)

  5. [5]

    Ericksen JL. 1991. Liquid crystals with variable degree of orientations, Archive for Rational Mechanics and Analysis 113, 97-120 (doi:10.1007/BF00380413)

  6. [6]

    Goodbrake C, Goriely A, Yavari A. 2021. The mathematical foundations of anelasticity: existence of smooth global intermediate configurations, Proceedings of the Royal Society A 477, 20200462 (doi: 10.1098/rspa.2020.0462)

  7. [7]

    Goriely A. 2017. The Mathematics and Mechanics of Biological Growth, Springer-Verlag, New York. 36

  8. [8]

    Guillamat P, Mirza W, Bal PK, G´ omez-Gomz´ ales M, Roca-Cusachs P, Arroyo M, Trepat X. 2026. Guidance of cellular nematic elastomers into shape-programmable living surfaces, Science 392, 317-323 (doi: 10.1126/science.adz9174)

  9. [9]

    Hernandez A, Zhu C, Giomi L. 2026. Mechanics of axis formation in Hydra, arXiv a2601.05220 (doi:arxiv.org/abs/2601.05220)

  10. [10]

    Hoffmann LA, Carenza LN, Eckert J, Giomi L. 2022. Theory of defect-mediated morpho- genesis, Science Advances 8, eabk2712 (doi:10.1126/sciadv.abk2712)

  11. [11]

    Hoffmann LA, Carenza LN, Giomi L. 2023. Tuneable defect-curvature coupling and topo- logical transitions in active shells, Soft Matter 19, 3423–3435 (doi:10.1039/D3SM00071K)

  12. [12]

    Hopf H. 1983. Differential Geometry in the Large, Springer-Verlag, Berlin, Heidelberg

  13. [13]

    Huang R, Ogden RW, Penta R. 2025. Modelling volumetric growth in soft solids via residual stress, Journal of Elasticity 157(4), 73 (doi: 10.1007/s10659-025-10164-2)

  14. [14]

    Kawaguchi K, Kageyama R, Sano M. 2017. Topological defects control collective dynamics in neural progenitor cell cultures, Nature 545, 212 (doi: 10.1038/nature22321)

  15. [15]

    Copenhagen K, Alert R, Wingreen N, Shaevitz J. 2021. Topological defects promote layer formation in myxococcusxanthus colonies, Nature Physics 17, 211 (doi: 10.1038/s41567- 020-01056-4)

  16. [16]

    Keber FC, Loiseau E, Sanchez T, DeCamp SJ, Giomi L, Bowick MJ, Marchetti MC, Dogic Z, Bausch AR. 2014. Topology and dynamics of active nematic vesicles, Science 345 (6201), 1135-1139 (doi: 10.1126/science.1254784)

  17. [17]

    Kumar N, Zhang R, de Pablo JJ, Gardel ML. 2018. Tunable structure and dynamics of active liquid crystals, Science Advances 4, eaat7779 (doi: 10.1126/sciadv.aat7779)

  18. [18]

    Lee P, Wolgemuth CW. 2011. Crawling cells can close wounds without purse strings or sig- naling, PLOS Computational Biology 7(3), e1002007 (doi: 10.1371/journal.pcbi.1002007)

  19. [19]

    Livshits A, Shani-Zerbib L, Maroudas-Sacks Y. 2017. Structural inheritance of the actin cytoskeletal organization determines the body axis in regenerating Hydra, Cell Reports 18, 1410-1421 (doi: 10.1016/j.celrep.2017.01.036)

  20. [20]

    Lubarda VA. 2004. Constitutive theories based on the multiplicative decomposition of de- formation gradient: thermoelasticity, elastoplasticity and biomechanics, Applied Mechanics Reviews 57(2), 95-108 (doi: 10.1115/1.1591000)

  21. [21]

    Marchetti MC, Joanny JF, Ramaswamy S, Liverpool TB, Prost J, Rao M, Simha RA

  22. [22]

    Hydrodynamics of soft active matter, Reviews of Modern Physics 85, 1143 (doi: 10.1103/RevModPhys.85.1143)

  23. [23]

    Maroudas-Sacks Y, Garion L, Shani-Zerbib L, Livshits A, Braun E, Keren K. 2020. Topo- logical defects in the nematic order of actin fibers as organization centers of Hydra mor- phogenesis, Nature Physics 17(2), 1-9 (doi: 10.1101/2020.03.02.972539)

  24. [24]

    Maroudas-Sacks Y, Suganthan S, Garion L, Ascoli-Abbina Y, Westfried A, Dori N, Pasv- inter I, Popovi´ c M, Keren K. 2025. Mechanical strain focusing at topological defect sites in regenerating Hydra, Development 152(4) (doi: 10.1242/dev.204514). 37

  25. [25]

    Mihai LA. 2022. Stochastic Elasticity: A Nondeterministic Approach to the Nonlinear Field Theory, Springer Cham, Switzerland (doi: 10.1007/978-3-031-06692-4)

  26. [26]

    Napoli G, Paparini S. 2016. Nematic bubbles and the breaking of spherical symmetry, arXiv 2603.07082 (doi: arxiv.org/abs/2603.07082)

  27. [27]

    Narayan V, Ramaswamy S, Menon N. 2007. Long-lived giant number fluctuations in a swarming granular nematic, Science 317, 105 (doi: 10.1126/science.1140414)

  28. [28]

    Paparini S, Giusteri GG, Mihai LA. 2025. Shape instabilities driven by topological defects in nematic polymer networks, Journal of Elasticity 157, 69 (doi: 10.1007/s10659-025-10160-6)

  29. [29]

    Pearce DJG, Thibault C, Chaboche Q, Blanch-Mercader C. 2025. Passive defect driven morphogenesis in nematic membranes, Physical Review Letters 134(1), 018402 (doi: 10.1103/PhysRevLett.134.018402)

  30. [30]

    P´ erez-Gonz´ alez C, Alert R, Blanch-Mercader C, G´ omez-Gonz´ alez M, Kolodziej T, Bazel- lieres E, Casademunt J, Trepat X. 2019. Active wetting of epithelial tissues, Nature Physics 15(1), 79-88 (doi: 10.1038/s41567-018-0279-5)

  31. [31]

    Perros T, Biquet-Bisquert A, Ben Meriem Z, Delarue M, Joseph P, Marcq P, Cochet- Escartin O. 2025. Mechanical characterization of regenerating Hydra tissue spheres, Bio- physical Journal 123, 1792-1803 (doi: 10.1016/j.bpj.2024.05.022)

  32. [32]

    Poincar´ e, H. 1886. Sur les courbes d´ efinies par les ´ equations diff´ erentielles (quatri` eme par- tie), Journal de Math´ ematiques Pures et Appliqu´ ees, 4e s´ erie, 2, 151–217

  33. [33]

    Ravichandran Y, Vogg M, Kruse K, Pearce DJG, Roux A. 2025. Topology changes of Hydra define actin orientation defects as organizers of morphogenesis, Science Advances 11(3), eadr9855 (doi: 10.1126/sciadv.adr98)

  34. [34]

    Rosso R, Virga EG, Kralj S. 2012. Parallel transport and defects on nematic shells, Con- tinuum Mechanics and Thermodynamics 24, 643-664 (doi:10.1007/s00161-012-0259-4)

  35. [35]

    Sadik S, Yavari A. 2017. On the origins of the idea of the multiplicative decomposition of the deformation gradient, Mathematics and Mechanics of Solids 22(4), 771-772 (doi: 10.1177/1081286515612280)

  36. [36]

    Saw TB, Doostmohammadi A, Nier V, Kocgozlu L, Thampi S, Toyama Y, Marcq P, Lim C, Yeomans J, Ladoux B. 2017. Topological defects in epithelia govern cell death and extrusion, Nature 544, 212 (doi: 10.1038/nature21718)

  37. [37]

    Stoker JJ. 1989. Differential Geometry, Wiley, Wiley Classics Library

  38. [38]

    Wang Z, Marchetti MC, Brauns F. 2023. Patterning of morphogenetic anisotropy fields. Proceedings of the National Academy of Sciences 120, e2220167120 (doi: 10.1073/pnas.2220167120)

  39. [39]

    Warner M, Gelling KP, Vilgis TA. 1988. Theory of nematic networks, The Journal of Chemical Physics 88(6), 4008–4013 (doi:10.1063/1.453852)

  40. [40]

    Warner M, Terentjev EM. 2007. Liquid Crystal Elastomers, paper back, Oxford University Press, Oxford, UK

  41. [41]

    Warner M, Wang X. 1991. Elasticity and phase behavior of nematic elastomers, Macro- molecules 24, 4932-4941 (doi: 10.1021/ma0 0 017a033). 38