Aylin Habibiyan†
a,
Shohei Yoshida†ab,
Rajas Sudhir Shaha and
Milana Trifkovic
*a
aDepartment of Chemical and Petroleum Engineering, University of Calgary, Calgary, T2N 1N4, Canada. E-mail: mtrifkov@ucalgary.ca
bFilms and Film Products Research Laboratories, Toray Industries Inc., Otsu, 520-8558, Japan
First published on 23rd September 2025
Hierarchical structures with bimodal porosity are crucial in diffusion and confinement-driven applications, such as catalysis and separation. This study introduces the first utilization of polymer blend nanocomposites as templates for isolating nanoparticle monoliths with bimodal porosity. We examined tri-continuous polymer blend nanocomposites of silica nanoparticles (SNPs) in polyethylene (PE), ethylene vinyl acetate (EVA), and polyethylene oxide (PEO) using three-channel confocal microscopy. This allowed visualization of their morphology and its evolution during quiescent annealing. The analysis extends to co-continuous polymer blend nanocomposites, with or without PEO. Our findings highlight the reinforcing effect of sequentially adding polymer phases in tri-continuous blends. This results in a refined morphology and strengthened three-dimensional particle network, as evidenced by a two-order-of-magnitude increase in the terminal modulus in frequency sweep rheometry. Conversely, co-continuous systems exhibit a significantly weaker particle network with a minimal increase in terminal storage modulus, making them prone to collapse during the polymer template removal. The interplay between domain size, nanoparticle jamming within one phase, and consequent particle network robustness enables the material to withstand deformation during polymer removal, facilitating the isolation of hierarchically structured monoliths. This novel templating method offers a scalable approach to fabricating hierarchically porous materials with potential applications in catalysis, energy storage, and gas separation.
New conceptsWe introduce a novel concept for the scalable fabrication of hierarchically porous monoliths using tri-continuous polymer blend nanocomposites as structural templates. Unlike most existing methods that rely on multi-step, solvent-intensive, and lab-scale synthesis routes, our continuous processing strategy integrates industrially viable extrusion and sintering methods to eliminate solvent use and simplify template removal, enabling a practical path toward scale-up. Our findings introduce a new paradigm where hierarchical porosity arises from a strong particle network formed via the unique confinement effect of nanoparticles within one of the tri-continuous phases. By comparing tri-continuous and co-continuous blends, we identify key factors critical to successful monolith isolation during polymer removal. Achieving high nanoparticle dispersion within the desired polymer phase, combined with a high jamming extent within small domains, is essential for preserving structure. When these conditions are not met—as observed in co-continuous systems—network collapse occurs during template removal. This work introduces new design parameters for creating robust porous materials and advances the understanding of how blend morphology influences particle reinforcement. The concept will be of broad interest to researchers in nanomaterials, polymer science, and porous materials for mass-transport limited applications such as catalysis, energy storage, gas separation, sensing, and drug delivery. |
Co-continuous and tri-continuous polymer blend nanocomposites are multi-phase systems. In these systems, two or more polymer phases form interconnected, continuous domains, which enable unique structural and functional properties.31 Tri-continuous morphology refers to a co-continuous polymer system separated by a third continuous phase, offering multiple interfaces required for broader applications.32 They are typically processed using scalable melt mixing techniques, which allow for precise control over phase morphology and the selective localization of nanoparticles.33,34 This tunable nanoparticle localization holds the potential to enable control over pore distribution and hierarchical structuring through the removal of sacrificial polymer components, making these systems highly promising templates for hierarchical porous materials. Despite this structural versatility and scalability, these systems have not yet been utilized for hierarchical porous material fabrication. Current applications primarily exploit their synergistic mechanical and functional properties in sectors such as automotive, aerospace, packaging, electronics, and biomedical devices.35–39 By broadening the range of polymer and nanoparticle combinations, these nanocomposites could provide a commercially viable platform for designing advanced porous materials tailored to catalysis, separation, energy storage, and other applications. The scalability of melt-blending methods further strengthens their potential for large-scale production.
Here, we report the first isolation of hierarchically porous materials using tri-continuous polymer blend nanocomposite as a sacrificial template. We establish a direct correlation between the spatial distribution of SNPs in PE/EVA/PEO tri-continuous blends and the morphological and rheological properties of the tri-continuous polymer blend nanocomposite system. Multi-channel LSCM observation with two distinct dyes enabled identification of each polymer phase, while the localization of untagged SNPs is detected using the reflection mode. We show that SNPs selectively localize and percolate within PEO, thereby inducing continuity in this particular phase. Furthermore, when PE, EVA and SNPs were blended prior to the addition of PEO, the tri-continuous structure undergoes further refinement. The obtained findings imply the existence of the unique confinement effect in tri-continuous polymer blend nanocomposites, which can be further exploited to achieve a strong particle network and consequently isolate porous materials. This method is more energy-intensive than soft gel-based methods due to the high processing temperatures and polymer template removal. However, its continuous processing nature and elimination of solvents make it a scalable and environmentally favorable approach for producing hierarchically structured materials. In particular, by combining tri-continuous polymer blend morphologies with nanoparticle jamming, our approach enables the fabrication of self-supporting monoliths that can be continuously manufactured through extrusion and shaped via an industrial process such as injection molding. This approach overcomes the scale-up limitations of batch-based methods such as bijels or emulsion templating, allowing efficient monolith production with tunable morphologies.
The wetting behavior of PE/PEO by EVA12 (or EVA25) were predicted by calculating the spreading coefficient based on Harkins theory, generalized by Torza and Mason.40,41 For an A/B/C ternary system, a set of three spreading coefficients can be calculated as follows:
λBAC = γBC − (γAB + γAC) | (1) |
λABC = γAC − (γAB + γBC) | (2) |
λACB = γAB − (γAC + γBC) | (3) |
During the initial stage of sequential blending approach, a co-continuous structure is established between PE and EVA, wherein SNPs are selectively localized within the EVA phase (Table S1). The disparity in viscosity between EVA12 and EVA25 (Fig. S1a) results in the generation of smaller domains in the PE/EVA12 system (1 ± 0.1 μm) compared to the PE/EVA25 system (1.7 ± 0.1 μm) where p-value ≪ 0.01. This, coupled with the migration of particles from the EVA phase to the PEO phase, results in the kinetic arrest of a finer microstructure in the PE/EVA12/PEO/SNP blend nanocomposite. Therefore, the average characteristic length became significantly smaller during sequential blending approach when PEO is added in the second step of blending. The comparison of characteristic size of the blends is illustrated in a histogram in Fig. S3. In the EVA12 system, the tri-continuous blend exhibits an increase in characteristic domain size compared to the neat blend, primarily due to the formation of large, interconnected PEO domains. In contrast, the sequentially processed sample yields the smallest domains. For EVA25, the domain size decreases gradually; however, it never reaches the reduced size observed in EVA12, since the system inherently begins with much larger PE and EVA domains, as shown in Fig. 1c.
Rheological investigation of the neat blends and PE/EVA/PEO/SNP composites as shown in Fig. 1e and f provides additional insights into the morphology and SNP packing within the PEO phase. The addition of SNP results in formation of the SNP particle network but its effect is highly dependent on the system as well as type of blending. The appearance of a plateau in viscoelastic moduli at low frequencies is a signature of the particle network.46 Simultaneous blending of all components in the PE/EVA25/PEO system yields a higher plateau modulus compared to the PE/EVA12/PEO system. This observation aligns with the morphological signature exhibited by these blends. The presence of smaller PEO domains in the PE/EVA25/PEO system contributes to higher jamming extent of SNP within these domains. This, in turn, enhances the solid-like behavior of the blends. In contrast, sequential blending leads to a significantly more pronounced plateau and an increase in G′ (storage modulus) compared to the simultaneous and pure blends. In jammed systems, the material exhibits solid-like behavior due to the constraints on particle movement. This is reflected in the high storage modulus (G′), which indicates the material's ability to store elastic energy. Despite the solid-like state, jammed systems exhibit a viscoelastic behavior characterized by energy dissipation (evidenced by a high loss modulus, G′′) resulting from particle rearrangement under applied shear stress. Under oscillatory shear, this viscoelasticity is more evident, as the material shows both solid-like and fluid-like responses depending on frequency and strain amplitude. The closeness of G′ and G′′ values (in Fig. 1e and f) in a frequency sweep suggests that the material's elastic and viscous properties are comparable. This balance is typical in viscoelastic systems near the jamming transition, where a material exhibits properties of both states.47
Fig. 2a and b illustrate the polymer template removal process with the specific thermal profile using a muffle furnace. Fig. 2c–e show the confocal images of the quaternary polymer blends after 30 min annealing and the SEM images after removal of polymer components through the heating process. Highly porous self-standing SNP structures were isolated only in the sequentially blended samples (Fig. 2d and e), while the larger PEO domains filled with SNPs were not preserved during the polymer removal process (Fig. 2c). The visual inspection of samples before and after the template removal process indicates retention of sample shape as a result of interconnected SNPs network during the polymer removal process of the sequentially blended sample. To confirm the jamming extent of nanoparticles, a water droplet was deposited on the PE/EVA12/PEO/SNP (sequential) surface for 30 s and then removed, selectively dissolving surface PEO and revealing the underlying particle network (Fig. 2f).
To assess the achievement of isolated hierarchical monoliths through the combination of high particle loading and small phase domains, co-continuous polymer blend nanocomposites of PE/EVA12/SNP and PP/EVA12/SNP were examined. PE/PEO/SNP and EVA/PEO/SNP blend nanocomposites were also prepared as shown in Fig. S4 and S5, but were not further investigated because of the lack of co-continuity in their morphology. Fig. 3 illustrates the co-continuous polymer blend nanocomposites that were prepared using both sequential and simultaneous blending approaches. The co-continuous polymer nanocomposite systems exhibited comparable domain sizes to the tri-continuous polymer nanocomposites, along with an SNP loading equivalent to a single polymer phase (30 wt% relative to EVA). While a porous structure of SNP monoliths was achieved, the post template removal samples exhibited structural collapse in each case. The calculated wetting parameters for PP/EVA12/SNP and PE/EVA12/SNP are −3.2 and −4.1, respectively (SI, Table S1), indicating that EVA is the thermodynamically preferred phase. However, high magnification LSCM images shown in the insets of 2D images of Fig. 3 indicate the presence of SNP aggregates at the blend interfaces. In the sequential blending process, pre-blending of the SNP particles with PP or PE phase, results in the formation of SNP aggregates. The subsequent addition of the thermodynamically favorable EVA12 phase induced a co-continuous morphology, creating Laplacian pressure at the highly curved blend interface. This pressure prevents the migration of the primary particles and small aggregates towards the EVA12 phase, while it is unfelt by the larger agglomerates,46 which cross the interface and localize within EVA12 phase. Although primary particles and small aggregates at the interface help suppress coarsening,43 micron-sized aggregates create localized stress concentrations and weak points within the interfacial network. Conversely, larger aggregates localized within EVA12 reduce dispersion extent and result in a less jammed SNP network, compromising mechanical integrity. During polymer removal, the capillary forces exerted on the rigid walls of these structures can cause deformation, making the interface susceptible to collapse. In contrast, in tri-continuous polymer blend nanocomposites, SNPs are well-dispersed within the PEO phase, and no particle aggregates are observed at the interface between the PEO and EVA phases (Fig. 1 and Fig. S2). The morphology of tri-continuous and co-continuous polymer nanocomposites is directly reflected in the inner pore structure of the remaining monoliths. This correlation is evident in the comparison between domain sizes of the annealed polymer nanocomposites and the pore sizes of the obtained monoliths. This dispersion extent and complete localization within the PEO phase enable interfacial rearrangement during polymer phase removal, thereby preserving structural integrity during the polymer removal process. Considering the interfacial tensions, SNPs exhibit the highest affinity for the PEO phase, with a low interfacial tension of 3.2 mN m−1. In contrast, the interfacial tension of SNP/EVA12 and SNP/EVA25 are 8.7 mN m−1 and 12.6 mN m−1, respectively, indicating a weaker interaction with EVA compared to PEO. This indicates that SNPs are more effectively distributed within the PEO phase in tri-continuous polymer blend nanocomposites, while in co-continuous blends, particles are more likely to aggregate at the blend interfaces. Consequently, the main difference between the co-continuous and tri-continuous systems examined here is the extent of SNP dispersion in the thermodynamically unfavorable phase during the first blending step and the differing affinities of SNPs for the thermodynamically preferred phase during the second blending step. Additionally, a key disadvantage of co-continuous systems is that processing at higher nanoparticle loadings (30 wt% relative to EVA12 phase) presents scalability challenges for our templating method compared to tri-continuous systems.
The rheological behaviour of the PP/EVA12/SNP and PE/EVA12/SNP blend nanocomposites confirm the significant difference in the SNP network formation between co-continuous and tri-continuous blend nanocomposites (Fig. 3e and f). The neat PP/EVA12 and PE/EVA12 blends exhibit slightly lower modulus than the PE/EVA12/PEO neat blend (Fig. 1e and f), indicating the addition of a third polymer phase and creation of an additional interface increases the modulus of pure blends. However, a significant difference is observed when comparing the co-continuous and tri-continuous blends that contain SNP within one confined phase. The storage and loss moduli of the tri-continuous blend nanocomposites are two orders of magnitude higher than those of the co-continuous blend nanocomposites (Fig. 1e, f and 3e, f). This increase in moduli demonstrates a higher extent of structural refinement in the tri-continuous blend nanocomposite, allowing isolation of SNP monoliths after the polymer template removal. In contrast, the co-continuous blend nanocomposites tend to collapse and deform after the template removal process. Fig. S6 also highlights the difference between co-continuous and tri-continuous polymer blend nanocomposites templates, as shown in the strain sweep data. The linear viscoelastic region is evident in all samples prior to the yield point. Beyond the yield strain, the tri-continuous blend nanocomposites exhibit a sharp drop in shear stress, indicating the de-jamming of the nanoparticles, where the applied shear breaks down the previously jammed nanoparticle network. Following this drop, the system stabilizes into a plateau, signifying the onset of a new steady flow regime. In contrast, this sharp drop is absent in both neat polymers and co-continuous blend nanocomposites, suggesting that the nanoparticles in the co-continuous system are not sufficiently jammed which confirms the structure collapse after polymer removal in co-continuous polymer blend nanocomposites.12,48
To demonstrate the complete removal of the polymer template, TGA was conducted on monoliths derived from co-continuous and tri-continuous blend nanocomposites and compared with the thermograms of pure polymers. The pure polymers exhibited a significant weight loss of up to 100% due to their degradation at temperatures above 400 °C (Fig. 4a). However, this characteristic weight loss due to polymer degradation was not observed in the monoliths, indicating the absence of polymer residue after template removal (Fig. 4b).
The thermogram for the pure silica nanoparticles used in this study shows the expected trend of two distinct weight loss stages. The first stage, up to 150 °C, is attributed to the release of adsorbed water on the hydrophilic silica nanoparticles. The second stage, occurring at temperatures above 250 °C, results from the removal of water produced during the non-catalyzed self-condensation reactions of silanol groups on the silica nanoparticle surfaces. For the monoliths, traces of adsorbed water are visible at temperatures below 150 °C, with a weight loss ranging from 1 to 3 wt%. However, the condensed water loss at temperatures above 250 °C is not observed in the monoliths, as the jammed silica nanoparticles have already undergone self-condensation reactions during the polymer removal process.49–52 The lower amount of absorbed water in the monoliths isolated from PE/EVA12/PEO/SNP and PE/EVA25/PEO/SNP samples is attributed to the highly jammed structure in these samples, reducing the penetration of water molecules into the jammed structure. however, in all cases, the weight loss difference remains within 2 wt%, indicating complete removal of the polymer phases.
The porosity of the PE/EVA12/PEO/SNP monolith prepared by the sequential method, which exhibited the most refined structure, was further evaluated by N2 adsorption–desorption analysis. The rapid uptake of N at p/p0 < 0.1 indicates the presence of micropores, while the Type IV isotherm with an H3 hysteresis loop is characteristic of mesopores larger than 4 nm (see Fig. 4c).53 The hysteresis loop, which does not exhibit limiting adsorption at high relative pressures, reflects the presence of meso/macropores formed as a result of particle jamming during fabrication.54 The nearly parallel adsorption and desorption branches further suggest that the mesoporous network is fully accessible and extends to the outer surface of the monolith, which is desirable for applications where enhanced mass transfer is required by shortening diffusion lengths in the micropore.55,56 The hierarchical porous structure was confirmed by combining density functional theory (DFT) and Barrett–Joyner–Halenda (BJH) analyses, used for the microporous and meso/macroporous regions, respectively. DFT analysis revealed micropores in the range of 1.2–1.7 nm, while BJH analysis identified mesopores with peaks at 21 nm and 40 nm, along with larger macropores centered around 93 nm (see Fig. 4d). It should be noted that our SNPs are ∼100 nm in diameter, and larger monolith pores present in the micron range are not detected by gas sorption analysis; however, their existence contributes to the relatively modest measured surface area (∼102 m2 g−1) and total pore volume (0.20 cm3 g−1), despite the highly porous nature of the monolith. This contrasts with silica monoliths produced by sol–gel processing,57 which typically have ∼15 nm silica nanoparticles,58 possess uniform small pores that yield very high surface areas, but lack hierarchical connectivity, thereby limiting pore accessibility. In comparison, the co-continuous and accessible hierarchical porosity of the PE/EVA12/PEO/SNP monoliths makes them promising candidates for catalysis, storage, and separation applications.
Since this is the first successful isolation of monoliths from the polymer blend nanocomposite systems, there is a significant potential for future research in this field. Further work should focus on modifying nanoparticles to achieve balanced dispersion in the thermodynamically unfavorable phase and strong affinity for the thermodynamically preferred phase across a range of practical polymer blend systems. While stable monoliths were not realized in the co-continuous systems studied, we do not rule out this possibility, recognizing that two-component polymer systems offer simplicity but face processing challenges with higher effective particle loadings. Additionally, refining monolith properties for targeted applications, such as cleantech and catalysis, will broaden their impact in high-performance materials. This approach not only tackles existing scalability challenges but also unlocks new opportunities in advanced material design, enabling the creation of multifunctional, high-performance monoliths with customized morphologies and enhanced properties for modern applications.
![]() | ||
Fig. 5 Visualization of separated channels of polymer blend images by LSCM, followed by the binarization and reconstruction using Avizo software. |
20× or 63× oil-immersion objective lens was used depending on the characteristic size of the derived blend microstructure. Every image was captured with 15 μm depth (z-step size 0.25 μm). Since the signal attenuation of SNP and Sudan Blue II channels decreased with increasing imaging depth, we compensated it by the laser gain adjustment to normalize the signal. The stack of 2D images from LSCM was binarized, reconstructed and analyzed in 3D using Thermo Scientific Avizo (CMC Microsystems) to calculate the characteristic length of PP/EVA blend interface. The software used a generalized marching cubes algorithm to generate a triangular mesh to the 3D interface and calculate the interfacial area, by adding up the area of the triangles. The average characteristic length, λ, was then determined by eqn (4),
![]() | (4) |
Footnote |
† These authors contributed equally to this work. |
This journal is © The Royal Society of Chemistry 2025 |