Yew-nique Geometry: A Mathematical Perspective of the Yew Tree

Abstract

In a world where most trees must choose between strength and flexibility, one species dares to have both.

Enter the yew: a slow growing conifer forged in harsh, dim, nutrient-poor environments, where every cell must pull its weight, literally. Beneath its poisonous bark, a secret arsenal is waiting in store: circle-packed cells, helical microfibrils, and a root network outsmarting the rules of biological connectivity. This essay undertakes a journey into the geometric universe inside yew wood, from honeycomb efficiency to spiral mechanics, where sound waves slow down and turn into warm resonance, while roots evolve from orderly trees to cyclic meshes worthy of graph theory legends. From microscopic helices through centuries-spanning topology, we show how the yew became nature's computational geometrician and why its design solutions continue to inspire engineers, luthiers, and mathematicians alike.

Introduction

Some Christians believe that Jesus was crucified on a yew tree. Formidably, it would be possible for one to walk by that same “tree of death” today (Wauters, 2022). Though this conifer is a common symbol of death, famously used for Voldemort's wand, it can live up to 4,000 years (Hageneder, 2013, p. 14). The tree’s dichotomous relationship between life and death is more than mythology: it is an extremely poisonous species that can also cure cancer.

Although the Pacific Yew (Taxus brevifolia), the Japanese Yew (Taxus cuspidata), and the English Yew (Taxus baccata) are found in different regions of the world, the genus Taxus is predominantly distributed across temperate areas of the Northern Hemisphere (Thakur & Kanwal, 2024).

It is especially puzzling to understand how a tree could survive for millennia under such unforgiving conditions. The “tree of life” first appeared on Pangea over 200 million years ago and has therefore had extraordinary time to fine-tune itself and adapt to its environment. To thrive in cool, shaded forests, nutrient-poor soils, and slow-changing climates, the yew took a lesson from the tortoise in Aesop’s fable: grow slowly and steadily (Moir et al., 2013). Constant exposure to low temperatures, as seen in Figure 1, exerted selective pressure on the tree, favoring traits that enhance tolerance to temperature fluctuations and snow accumulation (Kornienko, Shkirenko et al., 2025). This demands a paradoxical material strategy: be dense enough to resist compression, yet elastic enough to bend without breaking. Though the physical nature of this adaptation was previously discussed, the yew also solves these issues through mathematical design principles.

Fig. 1. Taxus baccata in deforming due to snow loads (Butler, 2022)

To understand how the yew behaves like a spring that resists the challenges of its environment, we must take a “magic school bus” journey into the woods itself. Like other softwoods, yew wood is built from cellulose microfibrils bundled into tracheids (Barnett & Bonham, 2004). However, it distinguishes itself by how it organizes them on the cellular scale. Its almost circular cells are tightly packed into hexagon-like arrangements, forming a natural honeycomb structure. It seems the yew has learned a lesson or two from bees: the honeycomb, with its unique hexagonal geometry, provides exceptional strength and resilience. This structure is so efficient that it is commonly used in 3D printing and construction (Malicse, 2025).

Furthermore, these cell walls are reinforced by unusually large microfibril angles (MFA). This makes the tree analogous to a spring, where helical spirals provide the ability to balance strength and flexibility by absorbing shock and torsional stress (Barnett & Bonham, 2004). Consequently, if the yew were buried under a heavy pile of snow, as seen in the figure above, it could deform under the load and subsequently return to its original shape once the stress was removed (Kornienko, Shkirenko et al., 2025).

The same helical architecture also shapes the conifer’s acoustic behavior. Its combination of low axial stiffness and high density slows the speed of sound waves passing through the wood, producing the warm, sustained resonance characteristic of historical yew lutes (Hartzell, 2008).

Finally, the tree of life’s longevity is tied not only to its regenerative properties, previously explored, but also to its mathematical root topology. The yew transitions from a simple, acyclic, tree-like root structure to a cyclic, mesh-like network depending on its developmental needs. The former is efficient but fragile, while the latter is robust against damage (Elshoff & Marcotty, 1978). This transformation mirrors fundamental principles in graph theory, where additional edges enhance redundancy and resilience.

If trees had da Vincis, the yew would be the original one. It is a polymath of biological design that mastered geometry, physics, acoustics, and topology all at once.

Nature’s Computational Geometrician: Yew’s Solution to Packing

Yew, like many other trees, is composed of cellulose fibers which organize themselves into fiber aggregates, then into layers of material, then into cells. Wood cells are then able to formulate higher order structures which results in the wood that we can see with our naked eyes (Fig. 2). In general, wood cell arrangements are often modeled as irregular or regular hexagonal honeycombs, providing wood with its elasticity and multileveled structures (Rafsanjani et al., 2012).

Fig. 2. The hierarchical structure of wood tissue. While wood may seem like a homogenous structure at first glance, this material is made up of multiple hierarchical structures. Starting with organizations of growth rings made up of different cell alignments to fibril aggregates on the molecular scale, these hierarchical structures give wood unique properties that vary across different tree species (Rafsanjani et al., 2012).

However, in yew, wood cells can be modeled as circles due to their smaller cell size (Sun et al., 2025). While other wood cell types in both hardwoods and softwoods are highly variable, yew wood cells are organized and uniform, aligned in a linear fashion (Fig. 3). As previously discussed in our physics paper, the yew tree is particularly dense and elastic with high resistance to cracking and deformation. For example, in comparison to spruce trees with a density of around 570 kg/m³, yew has a higher density at around 700 kg/m³ (Keunecke et al., 2007). But these mechanical properties are not just owed to material chemistry, but also a result of the yew’s natural geometric optimization of cell shapes and organization.

Fig. 3. SEM (Scanning Electron Microscope) images of different hardwoods: (A) poplar wood, (B) olive wood, and (C) cocobolo wood. SEM images of different softwoods: (D) cedar wood, (E) pine wood, and (F) yew wood. Note that the scale bar in hardwoods (A-C) is 200 𝜇m while in softwoods (D-F) is 100 𝜇m (Adapted from Jia et al., 2017).

Tessellations and Circle Packing

To be dense, one must strive to pack as many things in a limited area/volume as possible. Imagine heading to McGill during morning rush hour in a crowded metro car: how do we put as many half-awake bioengineering students as possible into a rectangular space? This is the question that the yew tree must answer to achieve its density, and simultaneously, it is a question that mathematicians have attempted to answer through tessellations and circle packing problems. To begin, tessellation, or tiling, is a pattern of shapes that can fit together on a plane without gaps or overlapping which can extend infinitely in any direction (OpenStax, 2023). Shapes must meet at a vertex and must be combined using a transformation. Tessellation can be either regular, being made up of only one regular polygon like squares, triangles, or hexagons (Fig. 4A); semi-regular, made up of two or more regular polygons (Fig. 4B); or irregular, being made up of irregular polygons (Fig. 4C) (Aliu, 2024).

Fig. 4. Examples of (A) regular tessellation, (B) semi-regular tessellation, and (C) irregular tessellation (Adapted from Aliu, 2024).

Tessellations have long been observed in nature, like bee’s nests and turtle shells, and combined into human-made structures, like mosaics and floor tiles. In fact, an ancient conjecture called the “Honeycomb Conjecture”, that is now the Honeycomb Theorem, was proven in 1999 by Thomas C. Hales. It states that “any partition of the plane into regions of equal area has perimeter at least that of the regular hexagonal honeycomb tiling”, or in more simple words, hexagons can fill a space the best (Hales, 2001). This theorem also extends to circle packing, a mathematical study of the arrangement of circles on a plane with no overlapping. Axel Thue conjectured that regular hexagonal packing is the densest identical circle packing in the plane with a density of π / √12 = 0.90690 which is later proved by Hai-Chau Chang and Lih-Chung Wang (2010). Circle packing algorithms also aim to optimize the number of circles packed within various shapes in both 2D, like squares, rectangles, and circles, and in 3D, like cubes and spheres. Take a square region Ω, we want to pack n identical circles and search for the maximum radius r of n identical circles that fit into a square. The algorithm optimizes the variables r and (x_i, y_i), i ∈ I = {1,2, …, n}, and thus:

Maximize r:

Subject to: r < x_i < 1 - r AND r < y_i < 1 - r (i ∈ I) 

√((x_i - x_j)² + (y_i - y_j)²) ≥ 2r  (1 ≤ i < j ≤ n)

where (x_i', y_i'), i ∈ I denote the coordinates of the center of the circle i, and √((x_i - x_j)² + (y_i - y_j)²)  is the Euclidean distance separating the centers of circles i and j,   1 ≤ i < j ≤ n (Fig. 5A.) (Hifi & M’Hallah, 2009).

Fig. 5. (A) An example of N=924 circle packing optimized to fit within a square (Specht, 2025). (B) and (C) are examples of cross-sections of Taxus wood with different types of growth rings under a light microscope. The scale bars are 125 𝜇m. Note the rounded cells organized similarly to the idealized version in computed circle packing algorithms (Adapted from Ghimire et al., 2015).

To Be or Not to Be Dense: How Yew Has Its Cake and Eats It Too

Dense objects may be strong, but they are typically not very good at absorbing and redistributing energy: a piece of lead will not break if you were to punch it (however ill-advised), but it will also not return to its original form if it becomes dented. Yew hence must go a step beyond density and circle packing algorithms to organize its cells in such a way that its structural integrity is ensured even after impact. Its geometric solution is to create hierarchical structures of circles packed into rectangular or square cross-sections along growth rings, where cells are more compressed and have smaller radii. This not only increases the density of its cellsbut also increases its specific energy absorption (SEA), the amount of energy absorbed per unit mass of a material. A novel bio-inspired hierarchical circular and square cross-section honeycomb (BHSCH) was created by Sun and colleagues which mimics the cell structure of the yew tree with the goal of creating a crash-resistant material as seen in Figure 6 (2025).

Fig. 6. The organization of yew wood cells (A) and the subsequent modeling for hierarchical circular and square cross-section honeycomb (BHSCH) structures with (B) being the first order (n=1) structure, (C) being the second order (n=2) structure, and (D) being the third order (n=3) hierarchy structure. N refers to the number of layers, or hierarchies, the honeycomb possesses (Sun et al., 2025).

The specific energy absorption (SEA) is calculated with the following formula:

Where F(u) is the axial compression load, u_d is the displacement when reaching densification strain 𝜀_d (where the material is compressed so much it has no space to eliminate), and m is the total mass. When compared to other honeycomb structures, as seen in Figure 7, the third order yew honeycomb microstructure had a higher SEA at 3.23 kJ/kg, indicating that these structures are excellent as energy absorbers (Sun et al., 2025). The third order yew honeycomb structure also exhibits significant plastic deformation, allowing energy to dissipate across the material and cause less damage, a feature of yew previously discussed in our physics paper. The optimized geometrical structure of yew cells is thus another explanation of how the yew can achieve strain resistance far higher than other softwoods. For example, yew wood has been found to have a work of ultimate load, a measure of strength and resistance, that is more than twice as high as spruce at 173 kJ/m³ as compared to 83 kJ/m³ (Keunecke et al., 2007). Perhaps in organisms that do not need brains, being dense is in fact beneficial.

Fig. 7. Yew tree-based honeycomb structures compared to other bio-based honeycomb structures on an axes of specific energy absorption (SEA) and relative density (Sun et al., 2025).

Where Geometry Meets Wood: The Yew’s Spiral Design

While geometric packing explains how Yew wood achieves its exceptional density at the cellular scale, its mechanical behavior is governed by a second layer of geometry, the helical orientation of cellulose microfibrils, quantified by the microfibril angle. Wood may appear uniform at its surface, but at a microscopic scale; it is an intricately organized material. The mechanical behavior of the yew tree, whether it bends or snaps, emerges from geometric architectures invisible to the naked eye.

What is the Microfibril Angle?

The term microfibril angle (MFA) in wood science refers to the angle between the direction of the helical windings of cellulose microfibrils along the axis of the cell (Barnett & Bonham, 2004). A larger MFA would indicate a flat spiral, while a smaller one is a steep spiral. This concept parallels earlier discussions of helical structures in a physics context, where helices were analyzed for their ability to redistribute forces and energy (see Yew Physics). Here, the same mathematical principles apply to the biological structure of the wood cell wall. Mathematically, every tracheid cell wall is a helically reinforced cylinder, making the wood itself a geometric structure (Fig. 8).

Fig. 8. Cellulose microfibril in wood cell (Tabet & Aziz, 2013).

As outlined previously, wood cells packed together form a strong composite tube, whose mechanical properties depend strongly on the orientation of cellulose in their dominant S2 layer. Each cell has four distinct cell wall layers: Primary, S1, S2, and S3 (Fig. 9), but the S2 layer accounts for the variation in MFA across species (Tabet & Aziz, 2013). Differences in MFA values have a profound effect on the properties of wood, particularly in its stiffness.

Fig. 9. Layers of a mature cell (Tabet & Aziz, 2013).

While the microfibril angle describes orientation of cellulose within the S2 layer, the microfibrils themselves follow a precise geometric path: a helix. A helix is generated when a point rotates around a cylinder while simultaneously moving upward along its axis. In the case of a wood cell, the cylinder is the tracheid: the rotation occurs around its circumference, and the upward movement corresponds to the cell’s longitudinal direction. Thus, each cellulose microfibril can be described using the parametric equations of a helix (Weisstein, n.d.):

Where r is the radius of the cell wall layer, t is the angle around the cylinder, and c is the pitch per radian (constant vertical rise). This makes it clear that the tracheid wall is not merely biological tissue, but a regular, mathematically ordered spiral structure. This continues the idea introduced in our physics analysis, where the helical structure was treated as a spring capable of storing and releasing energy. Helical reinforcement is a recurring geometric strategy in nature, appearing in DNA, seashells, and keratin fibers. Microfibrils in yew follow the same helical winding characteristic of most softwoods.

Helix equations describe the path of a single microfibril, but the angle of the helix is what determines how the cell wall deforms under heavy loads. In a simple two-dimensional material, the maximum shear occurs when fibers are oriented at 45°, because planes inclined at 45° to the principal stress directions experience the greatest shear stress (Hellevik, 2018). Although wood cells are certainly not flat sheets, they are 3D cylindrical tubes, and fibers now wrapped around a cylinder rather than lying flat changes the mechanics. In a cylindrical wall, the load is shared between axial stretching and circumferential twisting. This behavior is analogous to a spring’s ability to deform by both extending and rotating, distributing strain more efficiently than a straight fiber. This shifts the optimal fiber orientation downward, into the 15-30° range (Fig. 10). Within that range, the cell wall doesn’t deform in just one direction: instead, the load is shared between stretching along the length of the cell and twisting around its circumference. This geometric reasoning provides an essential context for understanding the unusually large microfibril angle of yew wood.

Fig. 10. Parametric helix representing a cellulose microfibril oriented at the optimal 30° microfibril angle, modeled using (x, y, z) = (cost, sint, 0.577t). Figure created by author using Desmos 3D Calculator (Desmos, 2025).

Microfibril Angle in the Yew

The microfibril angle varies widely among different tree species, and this difference is crucial to understanding why the yew tree behaves so differently from most softwoods (Fig. 11). In typical conifers such as spruce and pine, the microfibrils in the S2 layer typically lie almost parallel to the longitudinal axis, producing small MFA values ranging from 0-5°. These small angles relate to high axial stiffness and limited elasticity, characteristics that define most softwoods (Barnett & Bonham, 2004). The yew, however, strays significantly from this model. Measurements of its S2 layer microfibrils show MFAs ranging from 16-31°, which are exceptionally high in comparison to spruce (Keunecke et al., 2007). Similarly, results reported mechanical tests place Yew’s MFA at 15-20° in transition wood, once again confirming that the species maintains usually high spiral orientations (Keunecke et al., 2008). Significantly, these values fall directly within the lower optimal range predicted for cylindrical fibre-reinforced structures, where load can be shared between axial stretching and circumferential twisting.

Fig. 11. SEM micrograph showing the microfibrils in the S2 layer of softwood (Tabet & Aziz, 2013).

A recurring challenge in engineering is achieving a balance between strength and flexibility, as often most structural materials excel in one domain at the expense of the other. Highly stiff materials will fail when bent, while flexible materials lose their ability to bear heavy loads. The yew presents a natural solution to this problem through the geometry of its cell wall microstructure. The steep helical trajectory of its cellulose microfibrils around the tracheid wall, producing such high microfibril angles, alters how the cell wall redistributes strain under load.

Stress-Strain Curves

Stress-strain curves describe how a material deforms under load by plotting stress (σ) against strain (ε), giving the relationship between the applied pressure and amount of deformation (Roylance, 2001). In mathematical terms, the slope of this curve, E = dσ/dε, represents the instantaneous stiffness, or Young’s modulus, a measure of a solid materials stiffness or resistance to compressive deformation (Vincent & Engler, 2015). Mechanical tests on isolated yew fibers reveal that they exhibit a biphasic stress-strain curve: an initial low-stiffness region followed by a steeper second phase (Keunecke et al., 2008). This behavior contrasts with spruce fibers, with an MFA of 0-5º, which display a nearly single-phase linear-elastic response (Fig. 12).

Fig. 12. Mechanical comparison of Yew and Spruce Fibers. This figure shows the distinct biphasic stress-strain curve of yew fibers, in contrast to spruce, which exhibits a linear response (Keunecke et al., 2008).

For the Yew, the two-stage response reflects the mechanics of its large microfibril angle. At low strains, the helix can rotate and accommodate deformation, reducing the initial modulus. As strain increases, load shifts to the cellulose backbone, increasing stiffness. The microfibril helix allows the material to absorb deformation first, then stiffen, preventing sudden failure (Keunecke et al., 2008).

Helices in Harmony: Why Yew Sounds Better

Beyond structural loading, the yew’s spiral microstructure also plays a critical role in its acoustic performance. The Yew wood is historically known for lute construction, a stringed musical instrument (Fig. 13) (Hartzell, 2008).

Fig. 13. Yew-wood body lute (Titmuss, 2020).

Sound propagation in wood depends on both the stiffness and mass density, described mathematically, by:

where v is the wave speed, E the Young’s modulus, and ρ the density (OpenStax, 2016). The large microfibril angle of the Yew wood reduces its longitudinal stiffness, thereby lowering the value of E in the wave speed equation and lowering the overall wave speed, v. This causes vibrations to travel more slowly through the instrument (Keunecke et al., 2006). Such properties explain why Renaissance luthiers, those who make stringed instruments, selected yew for lute bowls and ribs. Its geometric fiber architecture produced sound that was both warm and sustained, while its stiffness allowed thin, curved shells that resisted warping. Historical records show that yew-bodied instruments were so valued that workshops were only limited to a few per year (Hartzell, 2008).

Beyond the yew’s musical and acoustic properties, its complex and intricate root system can also be described mathematically.

Topology of the Yew

We must first understand how the yew tree’s root system is organized. The yew possesses a highly flexible and interconnected root structure that can be regarded as a biological network. In addition to their primary roots, yews often produce supplemental roots-secondary structures that arise when branches or trunks come into contact with the ground. These roots may merge with preexisting ones, forming a network that enhances nutrient transport and overall stability. This self-sustaining framework enables ancient yews to live for centuries, even as portions of the trunk deteriorate. The root pattern of a mature tree differs from that of a young one. The former exhibits a cyclic topology, while the latter remains acyclic. In its early developmental stages, the root system displays a simple, hierarchical structure similar to a branching tree or a directed acyclic graph (DAG). A DAG is a type of graph in which nodes are connected by one-way links that do not form cycles (Fig. 14). DAGs are commonly used to illustrate dependencies, causal relationships, or sequences of events that occur in a specific order.

Fig. 14. Example of a directed acyclic graph (DAG) (Directed Acyclic Graph (DAG), 2025).

Each root segment extends outward from the main root, producing lateral branches in a non-repetitive, unidirectional pattern. In this configuration, the flow of nutrients and water moves directionally from distal root tips toward the plant’s vascular center and upward to the shoots without feedback loops or circular connections. This acyclic structure facilitates efficient resource transport and growth during the early stages of development, before the root system evolves into a more interconnected and complex form.

As the tree matures, its topology becomes cyclic. Supplemental and adventitious roots develop and reconnect with primary and lateral roots (Fig. 15), forming loops within the network.

Fig. 15. The yew root system (Yew - Taxus Baccata | Plants | Kew, n.d.).

These new connections introduce redundancy, ensuring that water, nutrients, and signaling molecules can still reach their destinations even if a root segment is damaged or obstructed. The transition from a strictly acyclic to a partially cyclic structure greatly enhances the tree’s robustness and flexibility. In essence, the mature root network begins to resemble a partially meshed network found in engineering, where multiple alternative pathways ensure continued resource transport and resilience during environmental stress or physical damage (Fig. 16).

Fig. 16. Example of a partially mesh network (Gillis, 2021).

Euler’s Königsberg bridge problem can help us visualize cyclic topologies. In his problem, the city’s landmasses were represented as vertices (V) and the bridges connecting them as edges (E). Euler proved that it was impossible to cross each bridge exactly once and return to the starting point because more than two vertices had an odd degree (an odd number of edges connected to them). This insight founded graph theory and established that a network has a Eulerian circuit only if all vertices have an even degree, and a Eulerian trail if exactly two have an odd degree (Holloway, 1974). Applying this to root topology, an early root system behaves like an acyclic tree graph, where:

and the cyclomatic number (a metric that quantifies the number of linearly independent cycles in a graph, representing the minimum number of edges to remove to make the graph acyclic):

indicating no loops (in other words each edge is a “bridge” essential for connectivity). As the tree matures, new root connections increase E, producing β>0 and creating cycles that reduce dependence on any single path (Elshoff & Marcotty, 1978). This transforms the structure into a partially connected mesh, conveying to a resilient network in engineering, where redundancy ensures continuous nutrient and water flow even if one pathway fails. If β=0, it implies a young perfectly acyclic yew, whereas a large β shows many reconnections and loops that provide mechanical stability. In biological language, a rise in E relative to V indicates the formation of more cross-links or root fusions. An increased cyclomatic number thus measures not just topological complexity (Elshoff & Marcotty, 1978), but also the tree’s ability to redirect resources and uphold structural stability as portions of the root system deteriorate. Additional graph parameters, like edge connectivity (λ) and vertex connectivity (κ), can indicate the number of root connections or junctions that need to be removed for the network to become disconnected. Mature yews probably exhibit higher λ and κ values than younger ones, indicating their exceptional resilience to environmental disturbances.

Conclusion

Across its structure, the yew tree demonstrates how geometry becomes a survival strategy. To survive for millennia under harsh climatic conditions, the yew learns to be densely packed, filling its available space so that its wood remains resilient and strong. Its small circular cells align in an optimized fashion consistent with circle-packing principles, allowing it to achieve remarkable density. Yet it does not sacrifice elasticity, enabling the tree to return to its original shape after heavy snow loads. This geometric packing produces hierarchical honeycomb structures in which circles arranged within square-like frameworks optimize specific energy absorption.

Furthermore, its unusually steep microfibril angles apply helical mechanics to balance bending flexibility with load-bearing capacity. This geometric design arranges the MFA between 15–30°, forming a fiber-reinforced cylindrical structure capable of resolving applied loads into both axial and circumferential components. This allows the helical microfibrils to rotate and redistribute strain rather than concentrating it. By reducing effective longitudinal modulus and altering wave-propagation speed, the yew not only produces tougher wood but also creates acoustic properties valued for centuries in musical instrument design.

This resilience is further enhanced by its root system, which evolves from a simple acyclic graph into a redundant, mesh-like network that mirrors engineered systems built for durability and fault tolerance.

Each of these design solutions explains the tree’s longevity: how it endures an unforgiving climate, remains both dense and elastic enough to withstand snow loads, and stabilizes itself with a sturdy, self-reinforcing base. Seen through this lens, the yew tree is more than a botanical curiosity; it is a natural demonstration of how living systems can harness mathematics to solve structural, mechanical, and ecological challenges.