Stable polydisperse free-standing porous films made by mechanical deformation

Hsiao-Ping Hsu; Kurt Kremer

PMC · DOI:10.1039/d4sm00569d·August 2, 2024

Stable polydisperse free-standing porous films made by mechanical deformation

Hsiao-Ping Hsu, Kurt Kremer

PDF

Open Access

TL;DR

Researchers used simulations to show that stretching and cooling a mix of long and short polymer chains creates stable porous films.

Contribution

The study demonstrates that polydisperse polymer films can be made stable using mechanical deformation and cooling.

Findings

01

Long polymer chains stabilize the porous structure by spanning multiple pores.

02

Short chains aggregate inside the polymer matrix and avoid surfaces during cooling.

03

The process works similarly for polydisperse films as it does for monodisperse films.

Abstract

Using molecular dynamics simulations, we show that the methodology of making thin stable nanoporous monodisperse films by biaxial mechanical expansion and subsequent cooling into the glassy state, also works for polydisperse films. To test this, a bidisperse polymer system of an equal number of very long (≈72 entanglements) and short (≤4 entanglements) chains with a polydispersity index of 1.80 is considered. The void formation and the development of the local morphology upon expansion, relaxation, and cooling are investigated. As for the monodisperse case, long chains in thin porous polydisperse films extend over several pores, stabilizing the whole morphology. The short chains do not fill up the pores but tend to aggregate inside the polymer matrix and to avoid surface areas and reduce conformational constraints imposed by the surrounding, a scenario very similar to strain-induced…

Figures17

Click any figure to enlarge with its caption.

Fig. 1 — Snapshot configurations of six selected chains of N1 = 1900 out of the free-standing polydisperse film (upper) and hypothetical affine deformation (lower) of the initial conformation at λ = 1.0, 3.0, and 4.0, as indicated.

Fig. 2 — Two components of rescaled root-mean-square (rms) end-to-end distance, (〈Re,α2(N)〉λ/〈Re,α2(N)〉λ=1)1/2 (a) and rescaled rms radius of gyration, (〈Rg,α2(N)〉λ/〈Rg,α2(N)〉λ=1)1/2, bond orientation order parameter Qλ (c), and effective film thickness h (d), plotted versus the strain of λ. In (a) and (b), α = ‖ and ⊥ denote the components in the direction parallel and perpendicular to the expanding directions, respectively. In (a)–(c), data for two different chain lengths N1 = 1900 and N2 = 100 in expanded polydisperse film are shown, as indicated. Reference values are 〈Re,‖2(N1 = 1900)〉λ=1 ≈ 3567σ2, 〈Re,⊥2(N1 = 1900)〉λ=1 ≈ 1082σ2, 〈Re,‖2(N2 = 100)〉λ=1 ≈ 180σ2, 〈Re,⊥2(N2 = 100)〉λ=1 ≈ 85σ2 in (a), and 〈Rg,‖2(N1 = 1900)〉λ=1 ≈ 599σ2, 〈Rg,⊥2(N1 = 1900)〉λ=1 ≈ 211σ2, 〈Rg,‖2(N2 = 100)〉λ=1 ≈ 30σ2, 〈Rg,⊥2(N2 = 100)〉λ=1 ≈ 14σ2 in (b). In (a), (b) and (d), affine scaling laws are shown by straight lines for comparison.

Fig. 3 — Porosity ϕλ (a), pose size characterized by the mean pore size 〈Dpore〉λ and mean maximum pore size 〈D(max)pore〉λ (b), plotted as a function of strain λ for polydisperse film subject to expansion. Pore size distributions P(Dpore) plotted as a function of Dpore at several selected strain values of λ are shown in the inset of (b).

Fig. 4 — Two components of single chain structure factor Sc,α(qα,N) in the directions parallel (α = ‖) and perpendicular (α = ⊥) to the expanding directions for chains of N1 = 1900 and N2 = 100 at λ = 1.0 and 4.0, as indicated. The theoretically predicted scaling laws are shown by straight lines for comparison, cf. text.

Fig. 5 — Snapshot configurations of thin polydisperse porous films at λ ≈ 4.0 subject to relaxation at several selected relaxation times t and assumed thicknesses h, as indicated where all 1000 chains of chain length N2 = 100 are marked in orange colors (a), and six randomly selected chains of chain length N1 = 1900 are marked in different colors (b).

Fig. 6 — Porosity ϕ(t) (a), mean pore size 〈Dpore(t)〉 (b), and mean maximum pore size 〈D(max)pore(t)〉 (c), plotted versus relaxation time t for polydisperse (poly) and monodisperse (mono) porous films at λ ≈ 4.0, as indicated. In the inset of (a)–(c), we plot the same data for t > 0.25 × 106τ, respectively, versus 1/t. The straight lines indicate linear extrapolating of all data sets to t → ∞.

Fig. 7 — Histogram plot of pore size distribution P(Dpore) for thin polydisperse (a) and monodisperse (b) porous films at λ ≈ 4.0. Only data for thin porous films at the subsequent relaxation times t/τ = 0 and 1.2 × 106 are shown, as indicated.

Fig. 8 — Net stress σB(t) (a), and film thickness h(t) (b), plotted versus relaxation time t for polydisperse and monodisperse porous films at λ ≈ 4.0.

Fig. 9 — Rescaled monomer density profiles for all monomers, ρ(z)/ρ0 (a), monomers in chains of N1 = 1900, ρ1(z)/ρ0 (b), and N2 = 100, ρ2(z)/ρ0 (c), plotted as a function of z at several selected subsequent relaxation times t, as indicated. The centers of films in the z-direction are matched at z = 0σ. The interfaces located at Z(lower)G and Z(upper)G determined from ρ(z) for films at t = 0τ and 1.2 × 106τ are indicated by dashed and solid arrows, respectively.

Fig. 10 — Collective structure factor S‖(q‖) in the expanding direction for expanded thin polydisperse (a) and monodisperse (b) films at λ ≈ 4.0, plotted versus the wave factor q‖ at several selected relaxation times t, as indicated. Theoretical predictions are shown by solid straight lines for comparison.

Fig. 11 — Two components of rescaled mean square radius of gyration 〈Rg,α2(N,t)〉/〈Rg,α2(N,t = 0)〉 (a) in the directions parallel (α = ‖) and perpendicular (α = ⊥) to the expending directions, and bond orientational order parameter Qλ(t) (b), plotted versus relaxation time t. Data are for two chain lengths N1 = 1900 and N2 = 100 in thin porous films at λ ≈ 4.0, as indicated.

Fig. 12 — Two components of single chain structure factor Sc,α(qα,N) in the directions parallel (α = ‖) and perpendicular (α = ⊥) to the expanding directions for two chain lengths N1 = 1900 and N2 = 100 in the expanded films. Data for several selected relaxation times t are shown, as indicated. Expected scaling laws are shown by straight lines for comparison, cf. text.

Fig. 13 — Snapshot configurations of free-standing thin polydisperse porous films at several selected temperatures T, as indicated where all 1000 chains of chain length N2 = 100 are marked in orange colors (a), and six randomly selected chains of chain length N1 = 1900 are marked in different colors (b).

Fig. 14 — Total potential energy per monomer U(T)/Ntot (a) and film thickness h(T) (b) plotted as a function of temperatures T for expanded polydisperse film at λ ≈ 4.0. The linear lines give the best fit of our data along the liquid and the glass branch. The corresponding Tg is marked by an arrow. The whole set of data of U(T) is shown in the inset of (a).

Fig. 15 — Porosity ϕ(T) (a), mean pore size 〈Dpore〉T and mean maximum pore size 〈D(max)pore〉T (b), plotted as a function of temperature T for expanded thin polydisperse films at λ ≈ 4.0. Estimates of Tg from h(T) and U(T) are indicated by solid black and dotted gray arrows, respectively, cf.Fig. 14.

Fig. 16 — Histogram plot of pore size distribution P(Dpore) at T = 1.0ε/kB (a) and (c) and 0.5ε/kB (b) and (d) right after relaxing for 30 000τ at each T for thin porous polydisperse (a) and (b) and monodisperse (c) and (d) films at λ ≈ 4.0 upon fast expansion subject to cooling.

Fig. 17 — Collective structure factor S‖(q‖) in the expanding direction for expanded thin polydisperse film at λ ≈ 4.0, plotted versus the wave factor q‖ at several selected temperatures T, as indicated. Theoretical predictions are shown by solid straight lines for comparison.

Funding1

—European Research Council10.13039/501100000781

Peer Reviews

No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.

Videos

No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.

Taxonomy

TopicsBlock Copolymer Self-Assembly · Material Dynamics and Properties · Pickering emulsions and particle stabilization

Full text

Introduction

I.

Porous polymer membranes are of high scientific and technological interest because of their potential applications in many fields such as water filtration, fuel cells or in pharmaceutical and biotechnology industries. This wide field of potential applications makes robust and simple preparation procedures very desirable.^1,2^ Typically the preparation of (nano)porous films takes advantage of incompatibility of components. Thus, a variety of block copolymers^3^ have been employed for designing nanostructured porous materials by changing block volume fractions, molecular weights, thermodynamic interactions between chemically different blocks (incompatibility), and solvent qualities.^4–9^ A collection of recently developed synthesis methods of generating micro-, meso-, and macro-porous materials with different properties is summarized in ref. 10. There exist also several reviews describing specific application driven requirements such as specific large absorption areas, high surface to volume ratios,^5,8,11^ and prospective developments.^12^ However, most methods of making porous polymer membranes are rather complex and often need very elaborate precision chemistry, and require processes controlled in detail.

In our recent work,^13–15^ we have demonstrated by simulation and by experiment that stable well controlled nanoporous films can be made just by mechanical deformation of highly entangled monodisperse polymer films and a subsequent quench of this nonequilibrium material into the glassy state. We have shown that the long chains extend over several bridges between pores. The relaxation is significantly slowed down by chain entanglements, which display a high density at regions of merging bridges between pores. This stabilizes the pores and prevents further growth and coalescence. In our previous study we have found typical pore diameters dp ≈ 5–10dT, dT being the reputation tube diameter. For the simulation part we considered a melt of chains of length N = 2000 ≈ 71Ne, Ne = 28 being the entanglement length while for experiment we used polystyrene (PS) of 1000 kDa ≈ 60Me with a polydispersity index of MW/MN = 1.03. Thus in both cases there were either no short chains or their volume fraction was negligible. Of course if one thinks of a broader application of such a process, which is just based on mechanical deformation of a pure polymer system without any additives or stabilizing chemical reactions, the demand of such a high monodispersity could be detrimental. Thus we investigate the influence of short chains on the properties of the mechanically deformed polymer films by studying a system containing an equal number of long and short chains. We show that the same process can be successfully applied to such polydisperse highly entangled polymer films. The short chains aggregate in the center of the polymer bridges, which form the pore walls and avoid contact with the surface. By that the average pore size is somewhat larger due to a slightly increased effective entanglement length, but the qualitative picture remains unchanged.

The outline of the paper is as follows: we first describe the system and investigate the morphological properties of polydisperse film subject to biaxial expansion in Section II. Then the developments of porous structures of thin expanded films upon subsequent relaxation in Section III and cooling in Section IV are analyzed. Finally, Section V contains our conclusions.

Morphology of polydisperse films upon biaxial deformation

II.

We consider a thick polymer film containing an equal number of chains of two lengths, N1 = 1900 and N2 = 100, respectively, which leads to a polydispersity index of MW/MN = 1.8. This system was prepared by starting from the fully equilibrated free-standing monodisperse polymer film containing 1000 weakly semiflexible chains of N = 2000 monomers of a simulation model based on a modification of the standard bead-spring model,^16–19^ which is well suited to study systems with free surfaces and the glass transition, at bulk melt (monomer number) density^14,20^ρ0 = 0.85σ^−3^. To prepare the bidisperse films a short piece of length N2 = 100 was cut off at the end of each chain of N = 2000. This film of now 2000 chains was equilibrated for about 3τR,N=N2, τR,N being the Rouse time of chains of size N, a time long enough for the short chains to move their own diameter. Alternatively, equilibrated polydisperse films can also be prepared via a soft-sphere approach and hierarchical backmapping.^20,21^ By this we arrived at eqilibrated films of effective thickness h ≈ 130σ ≈ 4.4R^(0)^g(N = N1) and two lateral dimensions with periodic boundary conditions, Lw = Lx = Ly ≈ 134σ. For comparison, R^(0)^g(N = N1 = 1900) ≈ 29.3σ and R^(0)^g(N = N2 = 100) ≈ 6.7σ are the root-mean square radii of gyration for bulk chains^22^ of size N, while we here have Rg(N1 = 1900) ≈ 28.5σ and Rg(N2 = 100) ≈ 6.7σ, showing that chain in thick films are very similar to those in bulk.^14,20^

This free standing film is subject to a simple “biaxial expansion” deformation. It is stretched into two lateral dimensions, i.e., equal-biaxial strain^23–25^ with periodic boundary conditions up to a maximum expansion of λ × λ ≈ 4 × 4 at T = 1.0ε/kB, while the thickness of the film is free to adjust, cf. Fig. S1 of ESI† (ref. 26). We follow the same protocol as in our previous work.^13,14^ In our earlier work we compared fast and slow deformation and found no significant differences in the results. Thus, we here apply the fast deformation rate, as shown in Fig. S1 and S2 (ESI†). The deformation rate is set that we can expect subchains of length of up to 0.6Ne ≈ 17 can equilibrate during the expansion, while the conformations of longer strands will be affected. Details are given in the appendix.

Snapshots of free-standing polydisperse polymer films for λ ≈ 1.0, 3.0, and 4.0 produced with an effective average strain rate *

τe ≈ 2.61 (see Fig. S2b, ESI†) are shown in Fig. 1. To demonstrate the conformational changes of long chains in the film, six randomly selected chains of N1 = 1900 are marked in different colors. Obviously, the overall shape of long chains follows an affine deformation while short chains do not (shown in Fig. S3 of ESI†). The overall morphologies are quite similar as they were observed for a monodisperse polymer films subject to a similar expansion rate τR,N=2000 ≈ 14 710, i.e., τe ≈ 2.88, also allowing for relaxation of subchains up to about 0.6N*e (see ref. 14).

Snapshot configurations of six selected chains of N1 = 1900 out of the free-standing polydisperse film (upper) and hypothetical affine deformation (lower) of the initial conformation at λ = 1.0, 3.0, and 4.0, as indicated.

Normalized strain-dependent chain extensions are shown in Fig. 2a and b. Affinely, the deformation follows λ for the parallel components and λ^−2^ for the perpendicular, respectively. For N1 = 1900 the in-plane expansion is affine all the way up to λ ≈ 4.0, while deviations for N2 = 100 develop around λ ≈ 2.6. In the perpendicular direction, long chains deform affinely (λ^−2^) only up to λ ≈ 2.4 in agreement with the adjustment of the film thickness, cf.Fig. 2d, while short chains deform nonaffinely almost from the very beginning. These global conformational changes also lead to characteristic changes in the bond orientational order parameter Qλ. Choosing the z-axis as a reference, where ϕz is the angle between any bond vector and the z-axis. For an isotropic distribution of bond directions Qλ = 0, while Qλ = −1/2, if all bonds would lie in the xy plane. The strain-dependent orientational order parameter Qλ of bond vectors along chains of N1 = 1900 shows that bond vectors in the extended films tend to lie randomly along the direction parallel to the interfaces. For short chains of N2 = 100, Qλ approaches a plateau value for λ > 3.5. The effective film thickness h determined from the monomer density profile (see Fig. S4a of ESI†) follows affine deformation up to λ ≈ 2.6 then approaches a plateau value approximately for λ ≥ 3 in the thin film regime, Fig. 2d.

Two components of rescaled root-mean-square (rms) end-to-end distance, (〈Re,α2(N)〉λ/〈Re,α2(N)〉λ=1)1/2 (a) and rescaled rms radius of gyration, (〈Rg,α2(N)〉λ/〈Rg,α2(N)〉λ=1)1/2, bond orientation order parameter Qλ (c), and effective film thickness h (d), plotted versus the strain of λ. In (a) and (b), α = ‖ and ⊥ denote the components in the direction parallel and perpendicular to the expanding directions, respectively. In (a)–(c), data for two different chain lengths N1 = 1900 and N2 = 100 in expanded polydisperse film are shown, as indicated. Reference values are 〈Re,‖2(N1 = 1900)〉λ=1 ≈ 3567σ2, 〈Re,⊥2(N1 = 1900)〉λ=1 ≈ 1082σ2, 〈Re,‖2(N2 = 100)〉λ=1 ≈ 180σ2, 〈Re,⊥2(N2 = 100)〉λ=1 ≈ 85σ2 in (a), and 〈Rg,‖2(N1 = 1900)〉λ=1 ≈ 599σ2, 〈Rg,⊥2(N1 = 1900)〉λ=1 ≈ 211σ2, 〈Rg,‖2(N2 = 100)〉λ=1 ≈ 30σ2, 〈Rg,⊥2(N2 = 100)〉λ=1 ≈ 14σ2 in (b). In (a), (b) and (d), affine scaling laws are shown by straight lines for comparison.

The reduction of the normalized monomer density profile for λ > 2.3 shown in Fig. S4a (ESI†) indicates the onset of porosity ϕ (see Fig. 3a). Considering monomers in long and short chains separately, Fig. S4b and c show that ρ1(z) and ρ2(z) follow the same behavior as ρ(z). Here the porosity ϕ and pore size distribution P(Dpore) of pore diameter Dpore in expanded films is estimated following the definition given by Gubbins et al.,^27–30^ where ϕ and P(Dpore) depend on the accessible volume of a hard spherical test particle of size 1.0σ. I.e. the test particles explore all regions in the film, where the nearest monomer is at least a distance of 1σ away. This is a purely geometrical measure, as no interaction between test particles and monomers is considered. The porosity ϕ is then given by the percentage of void volume Vvoid compared to the total effective volume Vfilm = hLxLy of the films.^14^ Results of the strain-dependent porosity ϕλ, average pore size 〈Dpore〉λ, maximum pore size 〈D^(max)^pore〉λ, and probability distribution of pore size P(Dpore) are shown in Fig. 3. The porosity ϕ increases monotonically with the increase of strain λ. No significant change in 〈Dpore〉 while 〈D^(max)^pore〉 increases with the increase of strain λ slightly. P(Dpore) has a unimodal-like distribution and the distribution becomes broader with the increase of λ. Again this is very similar to the results of monodisperse films.

Porosity ϕλ (a), pose size characterized by the mean pore size 〈Dpore〉λ and mean maximum pore size 〈D(max)pore〉λ (b), plotted as a function of strain λ for polydisperse film subject to expansion. Pore size distributions P(Dpore) plotted as a function of Dpore at several selected strain values of λ are shown in the inset of (b).

To detect the anisotropy of the chain structure and to compare short and long chains the strain-dependent two components of single chain structure factor Sc,‖(q‖,N) and Sc,⊥(q⊥,N) are shown in Fig. 4. Initially, at λ = 1 both, Sc,‖(q‖,N) and Sc,⊥(q⊥,N) ∼ q⊥^−2^ follow the Gunier law decay Sc,α(qα,N) = N(1 − qα^2^Rg,α^2^(N)/3) for small qα as expected for ideal bulk chains. As the strain λ increases, Sc,⊥(q⊥,N) increases while Sc,‖(q‖,N) decreases which are consistent with the change in 〈Rg,⊥^2^(N)〉λ and 〈Rg,‖^2^(N)〉λ shown in Fig. 2b, respectively. In the thin film regime (λ ≥ 3.0), chains are highly stretched in the expansion plane. As for the monodisperse case chains Sc,‖(q‖,N) ∼ q‖^−4/3^ in an intermediate q‖ range. We relate this 2-d self-avoiding walks like structure to the pore structure in the film, as argued below. Again there is no difference between short and long chains in this regime. As for the monodisperse case the pores seem to introduce an effective excluded volume on shorter and intermediate length scales. Only on large length scales, q‖ ≤ 0.045σ^−1^ the ideal chain behavior Sc,‖(q‖,N) ∼ q‖^−2^ is recovered for N = N1 = 1900. In contrast, Sc,⊥(q⊥,N) ∼ q⊥^−2^, on short length scales (q⊥ > 2σ^−1^). On large length scales, a sharp interface described by a Porod law like scaling Sc,⊥(q⊥,N) ∼ q⊥^−4^ is observed for λ ≥ 3.0. However, the short chains seem to be a bit compressed in the perpendicular direction on shorter length scales.

Two components of single chain structure factor Sc,α(qα,N) in the directions parallel (α = ‖) and perpendicular (α = ⊥) to the expanding directions for chains of N1 = 1900 and N2 = 100 at λ = 1.0 and 4.0, as indicated. The theoretically predicted scaling laws are shown by straight lines for comparison, cf. text.

We also calculate the two components of the collective scattering function of the films, S‖(q‖) and S⊥(q⊥), respectively (see Fig. S5 of ESI†). The intensity of S‖(q‖) increases with λ on large and intermediate length scales while on short length scales (q‖ > 2σ^−1^), it remains unchanged, showing that the local monomer packing is still conserved. The conserved peak at shows that the inter-monomer packing distance of still remains the same. On large and intermediate length scales, S‖(q‖) ∼ q‖^−2^ is observed at λ ≥ 3.0. The sharp local minima at nq = hq⊥/(2π), nq = 1, 2,… in the curves of S⊥(q⊥) confirm that the film thickness estimate h shown in Fig. 2d is consistent with the estimate from S⊥(q⊥). Even the q⊥^−4^ envelope for λ = 4 is reasonably well displayed, which is an indication of a rather uniform film thickness. Altogether, the data are almost indistinguishable from the monodisperse case. There is, however one interesting difference for λ = 1.0. The low q‖ regime indicates a slightly smaller compressibility or spatial density inhomogeneity for the polydisperse film compared to the monodisperse one. This indicates that short chain additions seem to level out fluctuations more effectively.

Relaxation of expanded polydisperse film at T = 1.0ε/kB

III.

So far we have studied the initial structure of the expanded films and found that the differences between mono- and polydisperse films are only marginal and that the dominance of the very long chains to determine the film properties remains unchanged. We now turn to the relaxation of the expanded film at the process temperature, kept at λ ≈ 4.0. Our previous study on monodisperse films at λ ≈ 4.0 has shown a strong retardation of conformational relaxation. This eventually stabilized nanoscopic pores, which then slowly grew under fixed strain. The overall pore diameter was limited by entanglement constraints. The central question here is, whether and how the shorter chains influence the stability and growth of pores. It should be kept in mind that T = 1.0ε/kB ≈ 1.5T^(0)^g, T^(0)^g ≈ 0.67ε/kB being the glass transition temperature of a bulk polymer melt.^31,32^ To investigate the influence of polydispersity we follow the same relaxation protocol as in ref. 13 and 14 and observe the relaxation of the films by MD simulations in the NVT ensemble over a time of up to t = 1.2 × 10^6^τ. Fig. 5 shows the morphological changes of expanded polydisperse polymer film, see Fig. S6 in the ESI,† for more details. We observe the same nucleation and slowed coarsening of the hole structures as for the monodisperse case with, however, somewhat larger holes. There is one important difference between a typical nucleation and the growth mechanism in our present system. Small pores, which disappear by shrinking and not by merging with other pores, since the long polymers extend over several holes merging is prevented. Fig. 5 already points to three important features. The first, just as for the monodisperse films, long chains extend over several pores, stabilizing the whole morphology. Second, the short chains seem to tend to aggregate inside the polymer matrix and avoid the surface regions. This is even better seen, when we look at thin slices of the film. Slices of thickness 3σ are shown in Fig. S7 in the ESI.† And third, this aggregation scenario is very similar to strain-induced segregation between the strained long and relaxed short chains is observed due to the reduction of conformational constrains of short chains imposed by the surrounding.

Snapshot configurations of thin polydisperse porous films at λ ≈ 4.0 subject to relaxation at several selected relaxation times t and assumed thicknesses h, as indicated where all 1000 chains of chain length N2 = 100 are marked in orange colors (a), and six randomly selected chains of chain length N1 = 1900 are marked in different colors (b).

Times covered range up to t = 1.2 × 10^6^ ≈ 530τe, corresponding to the Rouse time of subchains of length of Ns ≈ 644. Similarly, as observed for monodisperse films upon fast expansion and in contrast to slow expansion, initially no well defined pore structure is seen. After a short time well defined pores nucleate. The pore sizes increase accompanied by some drop in the number of pores. The growth of the pores slows down significantly after about half of the relaxation time, which corresponds to the Rouse time of subchains of about 16Ne. Pore sizes in the fast expanded monodisperse film appear only weakly smaller than in the polydisperse film. This relaxation retardation is also demonstrated by the six marked chains whose conformations only marginally change, obviously due to the topological constraints these highly entangled chains encounter.^33–35^ Moreover, in all cases the shape of pores tends to become spherical to minimize surface tension.

The small difference in pore size and porosity is confirmed by the direct measurement of the porosity ϕ(t), the pore size Dpore(t), and the maximum pore size D^(max)^pore(t) which all increase upon time (see Fig. 6). For both films the average pore sizes and maximum pore sizes continue to grow only very slowly with increasing time. The average pore diameter at t = 1.2 × 10^6^τ corresponds to about 46σ ≈ 9dT for monodisperse film. Here we observe 56σ for the polydisperse film. Assuming that this size is directly determined by the tube diameter, we estimate the shift in dT upon cutting off the short chains, since they have ample time to relax completely in the course of the relaxation. With N1 = 1900 for the polydisperse system the packing length p = N1/(ρ1R^2^(N1)) is increased by a factor of 1/0.95, ρ1 being the density of just the long chains. Taking the relation between plateau modulus and p and Ne, respectively, one would expect a shift of dT of about 5%, which is relatively close to the shift in the pore diameter.^36^ Assuming a linear extrapolation to 1/t→0, ϕ(t) ≈ 61% and 〈D^(max)^pore(t)〉 ≈ 120σ converge within error to the same values for both cases while 〈Dpore(t)〉 is larger for expanded polydiseperse film, 〈Dpore(t)〉 ≈ 75σ ≈ 14dT, and 64σ ≈ 13dT for the monodisperse film, respectively. Again this shift is roughly within the expected range. A more detailed theoretical account, however, is needed to make this argument more quantitative. The probability distributions of pore size Dpore, P(Dpore) are shown in Fig. 7. At t = 0τ, P(Dpore) has a unimodal distribution for both cases, and it becomes a much broader multimodal distribution at t = 1.2 × 10^6^τ. For polydisperse film, results of P(Dpore) show that the probability of finding larger pore size Dpore increases while it decreases for small pore sizes, as illustrated in Fig. 5, similar as for the expanded monodisperse film upon fast expansion.

Porosity ϕ(t) (a), mean pore size 〈Dpore(t)〉 (b), and mean maximum pore size 〈D(max)pore(t)〉 (c), plotted versus relaxation time t for polydisperse (poly) and monodisperse (mono) porous films at λ ≈ 4.0, as indicated. In the inset of (a)–(c), we plot the same data for t > 0.25 × 106τ, respectively, versus 1/t. The straight lines indicate linear extrapolating of all data sets to t → ∞.

Histogram plot of pore size distribution P(Dpore) for thin polydisperse (a) and monodisperse (b) porous films at λ ≈ 4.0. Only data for thin porous films at the subsequent relaxation times t/τ = 0 and 1.2 × 106 are shown, as indicated.

The finding that the system relaxation slows down significantly also is supported by the reduction of restoring forces per unit area, σB(t), shown in Fig. 8. We observe a dramatic reduction in net restoring stress σB(t) = |Pzz(t) − (Pxx(t) + Pyy(t))/2| after a very short initial time of about (0.2–0.3) × 10^6^τ. At the same time Pzz(t) remains at 0.0ε/σ^3^. Furthermore, after that the time-dependent stress is almost indistinguishable between these two cases, again in accord with the visual inspections of the membranes. The results of h(t) in Fig. 8b show that the initial reduction of film thickness for both cases follows the same curve. This reverse behavior for the films upon fast expansion also appears in the change of monomer density profile ρ(z) (see Fig. 9). Eventually the expanded polydisperse film upon relaxation remains somewhat thicker in agreement with the apparently slightly larger porosity and thus larger pores.

Net stress σB(t) (a), and film thickness h(t) (b), plotted versus relaxation time t for polydisperse and monodisperse porous films at λ ≈ 4.0.

Rescaled monomer density profiles for all monomers, ρ(z)/ρ0 (a), monomers in chains of N1 = 1900, ρ1(z)/ρ0 (b), and N2 = 100, ρ2(z)/ρ0 (c), plotted as a function of z at several selected subsequent relaxation times t, as indicated. The centers of films in the z-direction are matched at z = 0σ. The interfaces located at Z(lower)G and Z(upper)G determined from ρ(z) for films at t = 0τ and 1.2 × 106τ are indicated by dashed and solid arrows, respectively.

The above described scenario is well supported by the in expansion plane collective structure factor S‖(q‖) in Fig. 10 at several relaxation times t. The region around the amorphous halo at remains unchanged upon relaxation for all times. Thus the local bead packing is not affected by our processes. With the increase of time, the signature of sharp pore surfaces, Porod scaling S‖(q‖) ∼ q‖^−4^, is stabilized and extended a little further to larger length scales, as expected by the slow increase of porosity^37^ϕ(t) (see Fig. 6). For both cases, the initially fuzzy interfaces sharpen and already after short relaxation time the same scaling is observed. The initially large q‖^−2^ regime narrows down to a small region. In all cases, S‖(q‖) reaches a shallow maximum/plateau at low q‖, roughly corresponding to distances around 100σ, corresponding to 2–3 average pore diameters and reminding of a semidilute 2-d liquid of hard disks (i.e. the pores). It should be noted that the late time data are almost indistinguishable from the same data obtained from films, which were expanded in a much slower process, see ref. 14.

Collective structure factor S‖(q‖) in the expanding direction for expanded thin polydisperse (a) and monodisperse (b) films at λ ≈ 4.0, plotted versus the wave factor q‖ at several selected relaxation times t, as indicated. Theoretical predictions are shown by solid straight lines for comparison.

This comes along with the restricted conformational relaxation of the individual chains, as characterized by their linear dimensions and structure factor. The time-dependent two components of 〈Rg,α^2^(N,t)〉 parallel (α = ‖) and perpendicular (α = ⊥), and the bond orientational order parameter Qλ(t) are shown in Fig. 11. For long chains (N = N1) 〈Rg,‖^2^(N,t)〉 only decreases marginally during initial relaxation and then remains almost unchanged throughout the whole relaxation time. 〈Rg,⊥^2^(N,t)〉 increases slightly with time t after a very short time initial decrease, but then gets stuck at a value compatible with the thickness of the thin film. Note that marginal relaxation indicates that the chain retraction inside the tube as predicted by the Doi-Edwards and GLaMM tube models^38,39^ here is strongly retarded. For short chains (N = N2), chain retraction is not observed since chains are only weakly entangled.^34^ 〈Rg,‖^2^(N,t)〉 decreases while 〈Rg,⊥^2^(N,t)〉 increases during initial relaxation and eventually move towards the equilibrium value. Especially 〈Rg,‖^2^(N,t)〉 reaches the unperturbed bulk value. Similarly the bond orientational order parameter Q(t) displays a significant relaxation delay towards the isotropic phase for long chains while an isotropic distribution of bond directions is found for short chains. The two components of the single chain structure factor, Sc,‖(q‖,N) and Sc,⊥(q⊥,N), only change slightly with time t for N = N1 while for N = N2, Sc,‖(q‖,N) ≈ Sc,⊥(q⊥,N) ∼ q⊥,‖^−2^ at t = 1.2 × 10^6^ ≈ 42τR,N1, cf.Fig. 12, as expected for ideal chains. The instantaneously observed crossover for N = N1 from a two dimensional self-avoiding walk like structure (Sc,‖(q‖,N) ∼ q‖^−4/3^) to ideal random walk like structure (Sc,‖(q‖,N) ∼ q‖^−2^) at low q‖ is even shifted to lower q‖ and hardly visible anymore. However, the perpendicular component still displays a pronounced Porod power law (Sc,⊥(q⊥,N) ∼ q⊥^−4^) indicating the sharp surface.

Two components of rescaled mean square radius of gyration 〈Rg,α2(N,t)〉/〈Rg,α2(N,t = 0)〉 (a) in the directions parallel (α = ‖) and perpendicular (α = ⊥) to the expending directions, and bond orientational order parameter Qλ(t) (b), plotted versus relaxation time t. Data are for two chain lengths N1 = 1900 and N2 = 100 in thin porous films at λ ≈ 4.0, as indicated.

Two components of single chain structure factor Sc,α(qα,N) in the directions parallel (α = ‖) and perpendicular (α = ⊥) to the expanding directions for two chain lengths N1 = 1900 and N2 = 100 in the expanded films. Data for several selected relaxation times t are shown, as indicated. Expected scaling laws are shown by straight lines for comparison, cf. text.

Stabilization by cooling

IV.

Above we have shown, that the conformational and morphological relaxation of the expanded films is significantly retarded. The same and ultimate stabilization can be achieved by cooling the systems down towards the glass transition temperature Tg.^32,40^ For this we follow a stepwise cooling protocol with the fixed cooling rate Γ = 8.3 × 10^−7^ε/(kBτ) = ΔT/Δt (ΔT = 0.025ε/kB, Δt = 30 000τ) by NVT MD simulations with Langevin thermostat, just as also applied for the monodisperse case.^13,14^ We apply this to the polydisperse film of film thickness h ≈ 16.0σ at λ ≈ 4.0, right after deformation. For this cooling rate chains of length O(100) can full relax around and weakly below T = 1.0ε/kB. Thus the short chains of N2 = 100 can equilibrate completely while for N1 = 1900 only subchains of similar lengths can relax. Fig. 13 illustrates this. Similar structures are observed for the thin monodisperse polymer film upon fast expansion.^14^ However, the pore sizes are slightly larger in the polydisperse system. The clustering of the short chain indicates some strain-induced segregation between the strained long chains and relaxed short chains. At the same time short chains seem to avoid the surfaces. However, this does not affect the suitability of the resulting porous polydisperse film, since long chains extend over several pore envelopes.

Snapshot configurations of free-standing thin polydisperse porous films at several selected temperatures T, as indicated where all 1000 chains of chain length N2 = 100 are marked in orange colors (a), and six randomly selected chains of chain length N1 = 1900 are marked in different colors (b).

To determine the apparent glass transition temperature Tg we resort to the total potential energy U(T) shown in Fig. 14, which gives T^U^g = 0.69(5)ε/kB. A good measure of the density shift upon cooling is difficult due to the pore structure. Another estimate of Tg from the change in film thickness h(T) based from the monomer density profile (see Fig. S8, ESI†) gives T^h^g = 0.74(4)ε/kB, respectively. The latter value is a bit too high, compared to independent studies of thin films.^32^ However, there is no influence of polydispersity on Tg. The glass transition is around 0.70ε/kB as found for thin monodisperse nanoporous films at ref. 14λ ≈ 4.0. The porous structures characterized by porosity ϕ(T) and by pore size Dpore(T) are presented in Fig. 15. As T decreases ϕ(T) is increasing, approaching a plateau value in a temperature region, where also the monomer density profile ρ(z) converges. Both 〈Dpore(T)〉 and 〈D^(max)^pore(T)〉 behave similar as ϕ. They all first increase with the decrease of T for T > Tg, and then tend to reach a plateau approximately around Tg. The resulting porosity ϕ ≈ 33%, 〈Dpore〉 ≈ 25σ, and 〈D^(max)^pore〉 ≈ 50σ for T < Tg are slightly larger than the expanded monodisperse polymer film upon fast expansion (ϕ ≈ 29%, 〈Dpore〉 ≈ 22σ, 〈Dmax〉 ≈ 42σ). The pore size distribution P(Dpore) presented in Fig. 16 at T = 1.0ε/kB is no longer a unimodal-like distribution after short chains are relaxed. At T = 0.5ε/kB, the multmodal distribution is slightly broader comparing to the expanded monodisperse film upon fast expansion.

Total potential energy per monomer U(T)/Ntot (a) and film thickness h(T) (b) plotted as a function of temperatures T for expanded polydisperse film at λ ≈ 4.0. The linear lines give the best fit of our data along the liquid and the glass branch. The corresponding Tg is marked by an arrow. The whole set of data of U(T) is shown in the inset of (a).

Porosity ϕ(T) (a), mean pore size 〈Dpore〉T and mean maximum pore size 〈D(max)pore〉T (b), plotted as a function of temperature T for expanded thin polydisperse films at λ ≈ 4.0. Estimates of Tg from h(T) and U(T) are indicated by solid black and dotted gray arrows, respectively, cf.Fig. 14.

Histogram plot of pore size distribution P(Dpore) at T = 1.0ε/kB (a) and (c) and 0.5ε/kB (b) and (d) right after relaxing for 30 000τ at each T for thin porous polydisperse (a) and (b) and monodisperse (c) and (d) films at λ ≈ 4.0 upon fast expansion subject to cooling.

The observed very strong similarity to the properties of monodisperse films also hold for the collective in-plane scattering function, S‖(q‖) as shown in Fig. 17. With the decrease of T, S‖(q‖) increases on large length scale, and levels off in a broad maximum/shoulder below q ≈ 0.2σ^−1^ related to the average pore size. More locally, surfaces become flat and sharp for T < Tg, i.e., S‖(q‖) ∼ q‖^−4^. The changes in local wall structure of pores and microscopic monomer packing also remains unchanged compared to the monodisperse example.

Collective structure factor S‖(q‖) in the expanding direction for expanded thin polydisperse film at λ ≈ 4.0, plotted versus the wave factor q‖ at several selected temperatures T, as indicated. Theoretical predictions are shown by solid straight lines for comparison.

Conclusions

V.

In summary, we have extended the concept of employing entanglements to create and stabilize nanoporous polymer films^13,14^ to polydisperse melts. By choosing a bidisperse system of an equal number of very long (N1 = 1900) and short (N2 = 100) chains, i.e. MW/MN = 1.8, the system contains short chains which can fully equilibrate throughout the film expansion process while the long chains cannot relax, leading to essentially conserved entanglements and by this to a stabilization of the pores. As in our first study, no additional chemical processing or stabilization beyond quenching the system into the glassy state is needed. Despite the addition of short chains even well above Tg the relaxation of the expanded film is significantly slowing down, making the whole process rather robust. The glass transition temperature itself for films at strain of λ ≈ 4.0 remains within our error bars at the bulk value. The temperature dependence of the film thickness, which displays a kink at a slightly higher temperature, points to additional stabilization mechanisms, which await a more detailed investigation. As in our previous study the pore properties are determined by the entanglement structure of the underlying polymer melt. However, the polydispersity leads to slightly larger pores of an extrapolated 〈Dpore〉 ≈ 14dT, which is within the expected regime, considering faster relaxation and reduced entanglement density due to the shorter chains. Importantly, the short chains do not diffuse to the surfaces of the pores and thus do not fill up smaller pores. They stay within the bulk of the polymer matrix and by that can assume equilibrium conformations. The pores themselves act like an excluded volume interaction on the conformation of the chains winding around them, which is seen in the intermediate exponent in the scattering function. Altogether this study demonstrates that stable nanoporous polymer films can be made by mechanical deformation also for polydisperse systems as long as there are enough long chains which extend over several pores and by that stabilize the overall structure.

Author contributions

H.-P. H. designed research, performed the molecular dynamics simulations, analyzed the data, and wrote the paper. K. K. designed research, analyzed the data and wrote the paper.

Data availability

Data for this article are available on https://github.com/dehphsu/Publication_data/tree/b0d0357891fb1a7b0696163018f787479d28dac6/SM_2024.

Conflicts of interest

There are no conflicts to declare.

Supplementary Material

SM-020-D4SM00569D-s001

Bibliography40

The reference list from the paper itself. Each links out to its DOI / PubMed record.

1Lin H. Ding Y. J. Polym. Sci.2020582433243410.1002/pol.20200622 · doi ↗
2Abetz V. Brinkmann T. Sözbilir M. Chem. Teach. Int.2021314115410.1515/cti-2020-0023 · doi ↗
3Bates F. S. Fredrickson G. H. Phys. Today 199952323810.1063/1.882522 · doi ↗
4Abetz V. Macromol. Rapid Commun.201536102210.1002/marc.20140055625451792 · doi ↗ · pubmed ↗
5Abetz V. Kremer K. Mueller M. Reiter G. Macromol. Chem. Phys.2019220180033410.1002/macp.201800334 · doi ↗
6Sankhala K. Koll J. Abetz V. Membranes 2020108310.3390/membranes 1005008332353997 PMC 7281245 · doi ↗ · pubmed ↗
7Bang J. Kim S. H. Drockenmuller E. Misner M. J. Russell T. P. Hawker C. J. J. Am. Chem. Soc.20061287622762910.1021/ja 060814116756319 · doi ↗ · pubmed ↗
8Ariga K. Nishikawa M. Mori T. Takeya J. Shrestha L. K. Hill J. P. Sci. Technol. Adv. Mater.201920519510.1080/14686996.2018.155310830787960 PMC 6374972 · doi ↗ · pubmed ↗