Hollow-core polygonal pagodas, typically built of brick or stone in polygonal cross-sections such as hexagonal or octagonal forms, represent a distinctive structural typology that developed over more than a thousand years. Unlike solid masonry towers, these pagodas contain a central vertical void that reduces weight and, in many historical designs, housed stairways or reliquary chambers. The walls behave as thin, curved shells arranged in a polygonal geometry, and this configuration gives the structures a characteristic stiffness and stress distribution under lateral loading. When an earthquake strikes, the ground shaking induces bending moments that grow with height, and the masonry, which is strong in compression but weak in tension, can crack along horizontal or stepped joints long before any visible collapse occurs.

The core innovation of the new research lies in translating this complex three-dimensional problem into a tractable analytical formulation. Instead of modeling every brick and mortar joint, the method treats the pagoda cross-section as an equivalent polygonal hollow section whose bending capacity is governed by the masonry’s limited tensile strength. The authors derive closed-form expressions for the bending moment at which cracking initiates on the tension face and for the progressive reduction in stiffness that follows as the cracked zone spreads around the polygonal perimeter. Because the equations rely only on basic section geometry, material properties, and the distribution of self-weight, they can be evaluated with a spreadsheet or a short script, making the method accessible to practitioners who lack access to specialized finite-element software.

Central to the formulation is the recognition that bending failure in these structures is intimately linked to axial load. The weight of the overlying masonry above any given level compresses the walls and partially suppresses tensile cracking, so the bending resistance of the section increases with the magnitude of the compressive force, up to a limit set by crushing. Near the top of a pagoda, where the axial load is small, the section has little capacity to resist bending, which explains why upper tiers and the slender crowns of these towers are so frequently damaged in earthquakes. The simplified method captures this height-dependent vulnerability explicitly, allowing conservators to identify the critical stories of a specific pagoda where seismic demands are most likely to exceed capacity.

The polygonal geometry introduces further subtleties that the method addresses directly. In a circular hollow section, bending stress flows smoothly around the perimeter, but in a polygonal section the abrupt changes in wall direction at each corner concentrate stresses and alter the location of the neutral axis as cracking progresses. The analytical model accounts for the discrete geometry of the polygon by integrating the contribution of each flat wall segment to the section’s moment of inertia and to the cracking moment. This means the method can distinguish, for example, between a hexagonal and an octagonal plan of similar overall size, reflecting the historical observation that section shape influences both stiffness and the pattern of earthquake damage in surviving pagodas.

To validate the approach, the researchers compared their analytical predictions against established benchmarks for masonry behavior under combined axial load and bending. The simplified formulation reproduced the expected capacity trends, including the increase in cracking moment with axial compression and the sharp loss of flexural rigidity once the tensile strength of the masonry is exceeded. Because the method is deliberately conservative in its treatment of material strength, it tends to err on the side of safety, which is a desirable property for heritage assessment where the consequences of underestimating vulnerability are effectively irreversible. The authors emphasize that the model is intended as a first-line screening tool rather than a replacement for refined analysis of individual monuments, but that its transparency makes it uniquely valuable in that screening role.

The practical implications extend well beyond academic interest. Many of the most celebrated pagodas in China, Japan, and Korea are designated cultural properties subject to strict conservation requirements, and seismic safety evaluations are routinely demanded after significant earthquakes or as part of ongoing preservation planning. Full nonlinear finite-element modeling of a historic pagoda requires detailed surveys, material testing, and considerable expert effort, resources that are often unavailable for the hundreds of lesser-known structures scattered across seismically active regions. A simplified analytical method that needs only geometry, an estimate of masonry strength, and a design-level seismic demand offers a way to triage this large inventory, directing expensive detailed studies toward the structures that the screening analysis flags as most at risk.

The method also supports retrofit decision-making. Once the critical sections and stories of a pagoda are identified analytically, engineers can explore targeted interventions, such as confining bands of compatible reinforcement, grouting of degraded mortar joints, or减轻 of upper-tier mass, and then re-evaluate the bending capacity with the same closed-form equations to gauge the benefit of each measure. This iterative, low-cost loop between assessment and design is particularly important in heritage contexts, where interventions must be minimally invasive, reversible where possible, and justified with clear evidence of need. A tool that makes the mechanics of bending failure explicit, rather than hiding them inside a black-box simulation, is well suited to the interdisciplinary dialogue between engineers, archaeologists, and conservators that effective preservation demands.

Beyond its immediate engineering use, the study contributes to a broader scientific conversation about the seismic resilience of unreinforced masonry heritage. Researchers in earthquake engineering have increasingly recognized that simplified, mechanics-based models play an essential complementary role alongside high-fidelity simulation, because they reveal the governing physical parameters, support rapid parametric study, and facilitate the interpretation of observed damage patterns after real events. For hollow-core polygonal pagodas, the new work identifies axial load level, wall thickness, polygonal geometry, and masonry tensile strength as the dominant parameters controlling bending failure, providing a compact framework that future studies of aftershock vulnerability, soil-structure interaction, or cumulative damage can build upon.

As climate-driven hazard assessments and urban expansion place growing pressure on historic building stocks, tools that combine rigor with accessibility are likely to become increasingly important to the heritage sector. The simplified analytical method for earthquake-induced bending failure of hollow-core polygonal masonry pagodas offers exactly that combination: rooted in the mechanics of masonry, calibrated against recognized behavior, and simple enough for everyday professional use. For the custodians of these centuries-old towers, the study transforms an intimidating structural problem into a set of calculable quantities, bringing modern earthquake science to bear on some of humanity’s most enduring, and most endangered, architectural achievements.

Subject of Research: Simplified analytical prediction of earthquake-induced bending failure in hollow-core polygonal masonry pagodas

Article Title: A simplified analytical method for earthquake-induced bending failure of hollow-core polygonal masonry pagodas

Article References: Lu, W., Xiang, L., Wang, Y., & Li, D. (2026). A simplified analytical method for earthquake-induced bending failure of hollow-core polygonal masonry pagodas. npj Heritage Science. https://doi.org/10.1038/s40494-026-03000-w

Image Credits: AI Generated

DOI: 10.1038/s40494-026-03000-w

Keywords: masonry pagodas, seismic vulnerability, bending failure, heritage conservation, earthquake engineering, hollow-core sections, polygonal masonry, analytical modeling, structural assessment, unreinforced masonry, retrofit planning, cultural heritage

The work, carried out by Bijan Kumar Gangopadhyay, an independent researcher based in West Bengal, India, addresses a long-standing challenge in nanoscale thermodynamics. For decades, scientists have known that the properties of a material change dramatically when its dimensions shrink to the nanometer scale. The reason lies in simple arithmetic: as particles get smaller, an ever-larger fraction of their atoms sits on the surface rather than in the interior. Surface atoms are less tightly bound than their bulk counterparts because they have fewer neighbors, and this deficit of bonding partners destabilizes the crystal, allowing it to melt at lower temperatures and to contract under the pull of surface stress.

What has often been missing from earlier treatments, however, is a clean, explicit way to connect geometry to thermodynamics. Many existing models rely on average coordination numbers, empirical fitting parameters, or assumptions borrowed from macroscopic thermodynamics that become questionable at small sizes. The new model takes a different route. It begins with a direct, geometrical count of the fraction of atoms that reside on the surface of a nanoparticle of arbitrary shape, whether that particle is a sphere, a cube, a tetrahedron, or a thin film. From this surface-atom fraction, the model derives a relaxation factor that quantifies the effect of dangling bonds, the unsatisfied chemical bonds that terminate at any free surface.

The physical logic is straightforward. Every atom in the interior of a face-centered cubic metal such as gold or copper is surrounded by twelve nearest neighbors, giving it the full complement of bonding interactions that define the bulk cohesive energy. An atom on a flat surface, by contrast, may have only eight or nine neighbors, while an atom at a corner or edge of a faceted particle may have fewer still. These dangling bonds represent missing cohesive energy, and the more of them a particle has, relative to its total number of atoms, the more its average binding energy falls below the bulk value. Because melting occurs when thermal energy overcomes cohesive binding, a reduced average binding energy translates directly into a reduced melting temperature.

The same surface-atom fraction also governs the lattice parameter, the characteristic spacing between atoms in the crystal. Surface stress, arising from the imbalance of forces experienced by surface atoms, pulls the outer layers of the crystal inward, and this contraction propagates into the interior. Experimental measurements dating back to classic electron-diffraction studies of gold in the late 1960s and of copper and platinum in the early 1970s have confirmed that nanoscale metallic particles do indeed have smaller lattice constants than bulk crystals, with the deviation growing as particle size shrinks. The new analytical model reproduces this behavior by linking the relaxation factor, which describes how surface atoms adjust their positions and bonding, to the same geometric quantity that controls melting.

Applying the framework to gold and copper nanoparticles across a wide range of sizes, the study finds good agreement with available experimental measurements of both melting temperature and lattice parameter. The comparison covers spherical particles, cubes, tetrahedra, and nanofilms, demonstrating that a single set of analytical expressions can handle geometries that differ radically in their surface-to-volume ratios. A thin film, with two dominant surfaces and a thickness of only a few nanometers, has a far larger fraction of surface atoms than a sphere of comparable characteristic dimension, and the model captures the consequences: stronger melting-point depression and more pronounced lattice contraction.

One of the study’s clearest findings concerns what the author calls the shape factor. As the shape factor increases, reflecting a geometry with a larger surface-to-volume ratio, both melting-point depression and lattice contraction intensify. This provides a practical design rule for experimentalists: if you want to tune the thermal behavior of a metallic nanostructure, changing its shape can be as consequential as changing its size. A tetrahedral particle and a spherical particle containing the same number of atoms will not melt at the same temperature, because their surface atoms carry different weights in the overall energy balance.

The implications extend beyond gold and copper. The model is formulated for face-centered cubic metals in general, and its analytical simplicity means it can be evaluated with pencil and paper rather than computationally expensive simulations. Molecular dynamics and Monte Carlo approaches remain indispensable for capturing the full atomistic detail of nanoscale systems, including surface reconstructions, facet-specific chemistry, and thermal fluctuations, but they are costly and often difficult to interpret. An analytical expression that captures the leading-size and shape effects gives researchers a fast screening tool and a physical baseline against which simulations and experiments can be compared. The author suggests the framework could be extended to other thermodynamic properties of metallic nanomaterials, such as Debye temperature, specific heat, and thermal expansion, which previous studies have shown follow related size-dependent trends.

The scientific pedigree of the problem is long. Researchers have proposed liquid-drop models, coordination-number models, and various semi-empirical relations to explain melting-point depression since the phenomenon was first systematically studied. What distinguishes the present contribution is its explicit geometrical foundation: rather than treating the surface-atom fraction as an adjustable parameter, the model computes it directly from particle shape, and then ties it transparently to the physics of dangling bonds. This makes the model physically transparent in a way that purely fitted formulas are not, and it explains why different shapes produce different depressions of the melting point without requiring shape-specific calibration.

For technologists, the stakes are real. Gold nanoparticles are workhorses of catalysis, plasmonics, biomedical imaging, and drug delivery, and copper nanoparticles are increasingly important in electronics, antimicrobial coatings, and thermal interface materials. In all of these applications, the particles are processed, annealed, and operated at temperatures where their reduced melting points matter. Sintering, coalescence, and shape changes during manufacturing are governed by the same surface thermodynamics that the new model describes. A reliable analytical prediction of when a given nanoparticle will begin to soften and rearrange could help engineers choose processing windows that preserve the carefully engineered shapes on which device performance depends. Conversely, controlled melting could be exploited to fuse particles into desired architectures. As nanomaterials continue to move from laboratory curiosities to manufactured components, compact predictive tools of this kind are likely to become standard equipment in the nanoscale designer’s toolkit.

Subject of Research: Size- and shape-dependent melting temperature and lattice parameter behavior of gold and copper metallic nanoparticles

Article Title: Unified analytical model for size- and shape-dependent melting temperature and lattice parameter of gold and copper nanoparticles

Article References: Gangopadhyay, B. K. (2026). Unified analytical model for size- and shape-dependent melting temperature and lattice parameter of gold and copper nanoparticles. Journal of Nanoparticle Research, 28(10), Article 249. https://doi.org/10.1007/s11051-026-06770-3

Image Credits: AI Generated

DOI: 10.1007/s11051-026-06770-3

Keywords: metallic nanoparticles, melting-point depression, lattice contraction, surface-atom fraction, surface relaxation, gold nanoparticles, copper nanoparticles, nanoscale thermodynamics, analytical modeling, face-centered cubic metals, nanofilms, surface stress