1. Introduction
One of the actual directions in creating materials with new functional properties is the introduction of nanometer-sized particles with various geometric parameters and physical properties into polymer materials. Most often, various ceramic and carbon-containing particles are used for implementation, which allows researchers to vary the properties of the selected matrix in a wide range [
1,
2]. It is known that the transport and selective characteristics of MMMs vary in a wide range depending on the type of polymer and particles [
3,
4,
5,
6,
7,
8]. In such structures, the properties of MMMs can change according to different scenarios and there is still no universal model that describes the observed effects. It is known that the characteristics of MMMs depend not only on the type of embedded particles, but also on the parameters of the structures formed by the particles in the matrix [
9,
10,
11,
12,
13].
Experiments show that the properties of MMMs change in a threshold manner when a certain concentration of embedded particles is reached. Such a sharp change in properties can be explained only within the framework of the percolation approach, assuming that at a critical concentration a connected structure (percolation cluster) is formed, which determines the additional flow of the target component in the MMMs. According to the classical percolation theory, when the critical concentration is exceeded, percolation structures become homogeneous and a monotonous change in the properties of such materials should be observed. However, as experimental studies show, fluctuations in properties are observed in real MMMs and the results obtained may differ significantly from forecasts. At the same time, the same material may show an improvement in properties for some gases and a deterioration for others. Despite the twenty-year history of research on this problem, these effects are regularly observed at present.
In work [
14], in PVA/PEG membranes filled with multilayer carbon nanotubes with a modified surface (SE-MWCNTs), with a CNT concentration of 0.5%, a maximum permeability to CO
2 and a selectivity of CO
2/N
2 and CO
2/CH
4 are observed. The permeability to other gases first increased and then decreased, i.e., it varied non-linearly in a narrow range.
When adding functionalized poly(ethyleneglycol) (PEG) carbon nanotubes, it was found that the permeability to CO
2 first increases and then decreases with an increase in the content of CNT-PEG in the dry state [
15]. Water vapor in the gas stream further improves the permeability to CO
2. At a relative humidity of 100% for a hybrid membrane containing 3 wt.% CNT-PEG, a permeability to CO
2 of 369.1 barrels was obtained with a CO
2/N
2 selectivity of 110.8, which exceeds the upper limit of Robson. At the same time, the permeability of CNT-PEG to CO
2 varies non-linearly depending on the concentration of CNTs. An increase in the mass concentration of CNTs from 0 to 3% leads to an increase in the permeability to CO
2 from 90 to 260 barrels, and an increase in the concentration of CNTs from 3 to 20% leads to a decrease in the permeability to CO
2 from 260 to 110 barrels.
Carbon nanotubes were wrapped by poly(styrenesulfonate) (PSSA) and poly(vinylpyrrolidone) (PVP). Wrapped nanotubes were incorporated into poly (vinyl alcohol) matrix [
16]. Depending on the MMM manufacturing method, two different effects were obtained from the introduction of CNTs. In the first case, an increase in the concentration of CNTs from 0 to 3% within the margin of error did not affect the flow of the azeotropic IPA/H
2O feed mixture at 50 °C, and an increase in the concentration from 3% to 5% led to a doubling of the flow. In the second, an increase in the CNT concentration from 0 to 3% led to a 3.5-fold increase in the flux, and a further increase in the flux led to a decrease in the flux from 0.14 to 0.08 kg/(m
2∙h).
The incompatibility between filler and polymer chains accompanied by particle agglomeration has a detrimental effect on the performance of mixed-matrix membranes (MMMs) [
17]. To obviate intermolecular forces of different additives, functional groups (-COOH, -NCO, and -NH
2) were grafted on the surface of multi-walled carbon nanotubes (MWCNTs) which were then incorporated as fillers in the poly(ether-block-amide) (PEBA) polymeric matrix in the range of 0.1–1 wt% loading. The permeability of MMMs to all the studied gases first increased with an increase in the concentration of CNTs, and then decreased.
In work [
18], MMMs were fabricated from a highly microporous Tröger’s base (TB) ladder polymer matrix (ITTB) containing triptycene as contortion center, and COOH-functionalized single-wall carbon nanotubes (CNTs) as filler ranging from 0.6 to 2.0 wt%. By the strong acid–base interaction between COOH (from the CNTs) and tertiary-amine (from the TB) as well as π-π interaction between the triptycene unit and the CNT skeleton, the resulting MMMs had good compatibility. The permeability of the obtained MMMs to all gases (O
2, H
2, CO
2) varied non-linearly with the CNT concentration. An increase in the concentration of CNTs from 0 to 1% led to an increase in the permeability to all gases under study, and a further increase in the concentration from 1 to 2% led to a decrease in the permeability to all gases except hydrogen; for hydrogen, the permeability did not change within the error. In recent studies on the introduction of CNTs into polymers, it is reported that the upper limit of Robson has been overcome [
14,
15,
16,
17,
18].
In two main approaches, such materials are described within the framework of self-consistent models (for example, Maxwell, Kang–Jones–Nair, etc. [
19]) in which the interaction of particles is not taken into account. However, these models can only explain monotonous changes in the properties of MMMs. Some of the observed threshold effects are described within the framework of the percolation theory approach (for example, the models of Shen, Kirkpatrick, etc. [
20]), which allows us to take into account cooperative effects in such structures. This approach is applicable to both convective and diffusive transport in MMMs and allows us to take into account the influence of morphology and the interaction of particles and the matrix.
In work [
21], the permeation process of gas molecules through MMMs embedding impermeable spherical and cuboid nanoparticles of different aspect ratios was simulated using a three-dimensional finite-difference solution of Fick’s second law of diffusion. The authors showed that, for the same volume fraction, the relative permeability of MMMs with cuboid nanoparticles with an aspect ratio larger than unity was lower than the relative permeability of membranes with spherical nanoparticles. This led to an analytical model, which involves only two geometrical parameters, the ratio of the projected area available for diffusion and the maximum surface area for diffusion, which depends on the aspect ratio of the nanoparticle, and the relative thickness of the nanoparticle. A model for the prediction of the mass transport through ideal MMMs for pervaporation and gas separation processes has been introduced [
22]. A resistance-based model was used in conjunction with a three-dimensional finite difference (FD) numerical solution to derive a semi-empirical model for calculating the effective permeability of ideal mixed-matrix membranes. The authors claim that the extended RB model was theoretically able to accurately predict the effective permeability of ideal MMMs over a large range of filler volume fraction. However, in this work, the authors also consider the ideal morphology, which is difficult to achieve in polymers with embedded nanoparticles due to the physical and chemical characteristics of the nanoparticles–matrix interface [
23]. However, these models and parameters of percolation structures are proposed for systems in which the dimensions of the matrix exceed the dimensions of embedded particles, which is not true for real selective layers. This fundamental difference, in our opinion, determines all the observed deviations of model predictions from experimental results. As will be shown below, the surface effects (the proportion of the modified material) and the commensurability of the matrix and particle sizes critically affect the conditions for the formation of the percolation structure.
The modeling of percolation structures in isotropic 2D and 3D systems has shown that an increase in the aspect number of particles (transition to tubular or planar objects) leads to a significant decrease in percolation concentration (
Figure 1) [
24]. A non-zero probability of physical interaction (contact) with embedded 3D particles exists in an effective volume near the particle, which can be described by a sphere with a diameter equal to the largest particle size. However, the proportion of space that this object actually occupies in the matrix will be determined by all three of its dimensions. Therefore, with an increase in the aspect number of particles (the ratio of two characteristic linear ones), the fraction of the particle volume in the effective volume decreases. At the same time, the probability of formation of a percolation cluster (interaction of individual particles) increases with an increase in the effective volume, or, accordingly, with an increase in the aspect number. The extremum of the critical concentration observed in
Figure 1 corresponds to particles with an aspect number close to one, at which the physical and effective volume are equivalent (without taking into account the effects of particle interaction through the matrix), and corresponds to the theoretical limit values of the percolation threshold in 3D isotropic space. At the same time, both for tubular particles with a large aspect number (CNTs) and for flat particles with a small aspect number (graphene), the effective volume is significantly larger than the real one and the percolation transition is observed at lower concentrations. Experimentally, this effect was confirmed for various MMMs with CNTs and other particles [
3,
9,
10,
11,
25,
26]. However, in other studies, there is a significant deviation of properties similar to MMMs [
27,
28,
29,
30,
31,
32,
33,
34,
35,
36]. It is precisely this difference in properties for identical materials that confirms our thesis about fluctuations in the parameters of the percolation transition, which cannot be explained only by particle sizes.
We expect that, in addition to the properties of the particles, the percolation parameters are controlled by the geometry of the membranes. Complex experimental and model studies of the conditions for the formation of a percolation cluster and its characteristics make it possible to build a consistent model of such structures in MMMs, which is necessary for the further development of the technology of synthesis and production of MMMs.
2. Calculation
The threshold change in the properties of materials with embedded nanoparticles is considered in the theory of percolation, in which it is believed that an abrupt change in the properties of materials with embedded nanoparticles occurs as a result of the formation of interconnected new structures, which are called a percolation cluster. This approach is often used for descriptions of changes in the electrical and thermal conductivity of polymers when nanoparticles are introduced into them.
Mass transfer in a polymer with carbon nanotubes embedded in its structure is possible:
- (1)
Through the polymer;
- (2)
Through the polymer/CNT interfacial layer, which results from poor adhesion between polymer chains and CNTs;
- (3)
Through the internal cylindrical channels of the CNTs.
And if we consider the flow through regions (1)–(3) as independent, then we can write the rule of mixtures, in which the fraction of the volume of the membrane with high permeability should be related to the volume of the percolation cluster, as
where C
x is the fraction of the polymer volume occupied by the percolation cluster, K
x is the permeability coefficient of the percolation cluster, K
p is the polymer permeability coefficient, and dR = R – R
0 is the thickness of the interfacial layer. The thickness of polymeric membranes with embedded CNTs differs slightly from the length of the CNTs; such a system cannot be considered infinitely large when calculating the parameters of a CNT percolation cluster, and surface effects must be taken into account. Expression 1 has a general form for a composite material parameter and illustrates the linear correlation of properties in a two-component material based on the rule of mixtures. It shows that the contribution of each component is proportional to its concentration and the observed nonlinear changes in properties can be associated with a nonlinear change in the number of particles affecting the properties of the matrix.
To determine the parameters of a percolation cluster of CNTs in systems of finite size, a software package in C# was developed, which allows researchers to determine the probability of formation and the strength of a percolation cluster for various geometric objects (sphere, sphere with an impenetrable core, spherical cylinder, spherical cylinder with an impenetrable core, segment, circle, and rectangle), and set the distribution function of geometric objects by size.
Due to the fact that carbon nanotubes are the object of study in this work, the algorithm for calculating the parameters of a percolation cluster was considered using the example of spherocylinders with impenetrable cores in the form of a spherocylinder.
In our model, each spherical cylinder is surrounded by a shell with a thickness of dR = R − R0, which forms regions of the polymer, the permeability of which increases due to the interaction with the surface of the nanotubes, and is determined by the interaction between the CNTs and the polymer. Spherocylinders form a cluster when shells overlap, and if a cluster connects opposite faces of a 3D matrix (upper and lower), then it is considered percolative. To achieve this, a three-dimensional matrix is randomly filled with spherical cylinders, the distance between them is determined, a list of formed clusters is formed, and it is checked whether there is a percolative one among them.
To determine the probability of formation of a percolation cluster, a given number of iterations is carried out with the same number of spherocylinders added to the system.
To determine the strength of the percolation cluster, the ratio of the number of spherical cylinders that make up the percolation cluster to the number of spherical cylinders added to the system is calculated as
where N
∞ is the number of spherical cylinders that make up the percolation cluster, and N is the number of spherocylinders added to the system.
The volume concentration of spherocylinders was determined as
where R
0 is the radius and L is the length of the spherocylinders, N is the number of spherocylinders, and V is the volume of the matrix.
The average value of the volume fraction occupied by high-permeability regions of the polymer is
Here,
is the number of spherocylinders in the percolation cluster. To verify the algorithm and program, the percolation threshold was calculated for various geometric objects [
37].
Since a spherical cylinder with an impenetrable core in the form of a spherical cylinder is a CNT model in our work, we will use CNT from now on in the text. We will give the dimensionless value of the geometric dimensions of the matrix and spherocylinders with units of measurement tied to the dimensions of polymer membranes and CNTs.
4. Conclusions
The purpose of our study was to model the conditions for the formation of a percolation structure in MMMs and to study the “stability” of such a transition when the geometric characteristics of embedded particles change. Since the transport properties of the MMMs are associated with the formation of a percolation cluster, it becomes crucial to understand the real parameters of such a cluster in the MMMs—the concentration of formation, the cluster’s strength, and its stability—when changing the geometric characteristics of the structure. The relationship between the parameters of percolation structures and the properties of MMMs is not the subject of this study, since there are no representative comprehensive experimental studies of both functional and structural properties of MMMs. At the same time, as has been repeatedly shown, the characteristics of the percolation cluster and particle sizes affect various macroscopic properties of the MMMs and the prediction of the parameters of such structures allows researchers to vary the functional properties of the MMMs. Let us consider several scenarios in which an erroneous assessment of real geometry leads to the absence of the predicted properties of MMMs.
Strong affinity. Let us suppose that at a certain concentration of particles with an assumed radius and length (R
0 = 25 nm, L = 4 µm) in a geometrically isotropic matrix, we obtained a percolation structure of a certain strength. As can be seen from
Figure 4, with an assumed layer thickness of dR = 100 nm, we can expect the formation of a percolation cluster with a capacity of about 80% at a particle concentration of more than 0.35%. However, if the actual thickness of the modified layer is 65 nm, then the cluster strength at the same concentration will be about 20%. That is, the proportion of MMMs in which specific transport associated with embedded particles is implemented will decrease by four times, which will significantly worsen the transport properties of such MMMs.
Geometric anisotropy of the matrix. Aiming to obtain a good percolation, for certain particles (R
0 = 25 nm, L = 5 µm, and dR = 100 nm) in a geometrically isotropic matrix, a CNT concentration of more than 0.2% was chosen, which guaranteed a 100% probability of forming a percolation cluster with a strength of about 50% (
Figure 8). However, with a decrease in the actual thickness of the film, the cluster strength is reduced by more than two times already with the ratio of lateral and normal matrix sizes at 5 µm (
Figure 8b). A further increase, despite the fulfillment of the concentration conditions of percolation (
Figure 8b), may lead to zero strength of the percolation cluster, and in such systems there will be a strong heterogeneity of the transport properties of the MMMs.
Size of CNTs. Too high a concentration of CNTs with inaccurate determination of their characteristics can also lead to negative changes in the percolation structure. For example, an excessive amount (0.5%) of CNTs (R
0 = 25 nm, L = 2 µm, and dR = 100 nm) has been added to the isotropic matrix. With such particle sizes, a percolation cluster is formed at concentrations of about 0.4% (
Figure 6); however, an error in determining the actual radius of the particles (for example, R
0 = 30 nm) leads to a situation where the concentration is insufficient to guarantee the formation of a percolation cluster, and the strength of such a percolation cluster is less than 20%. Thus, most of the formed MMMs will not contain a percolation structure. Errors in the actual CNT length lead to similar consequences. For example, for particles of a given size (R
0 = 25 nm, dR = 100 nm, and length L = 5 µm), we determine the required concentration of 0.2%, with some excess. However, if the actual length of the CNTs is equal to 4 µm, the percolation cluster is practically not formed, since the probability is less than 10% (
Figure 2) and its capacity is about 10%. In such a system, the change in transport properties will be significantly lower than expected and in a small number of experiments. Thus, an error in the geometric dimensions of the particles by 20–25% leads to an erroneous choice of the required concentration and, as a consequence, the absence of the necessary percolation structure in the MMMs.
A simultaneous implementation of the scenarios of variation and inaccuracy described above in the evaluation of geometric factors of MMM elements will lead to the instability of the observed properties of MMMs and the absence of percolation even at particle concentrations significantly exceeding the threshold value. From the point of view of the classical percolation theory, this result can be associated with the variance of probabilistic characteristics of percolation structures in systems of finite size. The transition of the system to a percolation state is possible both as a result of a change in the concentration of particles and by reducing the size of the system. For example, as shown in
Figure 9, by cutting out cells (matrices) of various sizes from an infinite percolation cluster, we can fix the presence of a percolation cluster with a given frequency. At the same time, the characteristic scale of clusters and the strength of the resulting clusters depends both on the proximity to the critical concentration and on the size of the system.
In scaling theory, the cluster structure near the transition is related to the size of the problem space [
38], and below the percolation threshold, the average cluster size determines the correlation length (
Figure 9). At the percolation threshold, the correlation length diverges. Since the correlation length above the percolation threshold has a finite value, the infinite cluster is self-similar. One can consider the correlation length as the most probable distance up to which an infinite cluster is self-similar and can be considered as a fractal. At large length scales, the cluster structure is no longer self-similar, and the cluster can be considered as homogeneous.
Thus, when describing the formation of percolation structures in films (systems of finite size), it is necessary to take into account multidirectional trends when changing the size and concentration of particles in the matrix. On the one hand, an increase in the concentration of particles and the correlation radius leads to an increase in the size of agglomerates, the formation of a percolation cluster, and the transition of the system to a homogeneous state. On the other hand, a decrease in the size of the matrix (system) leads to an increase in the variance of the probability of formation of a percolation cluster at a given concentration and the distribution of clusters by size.
Using the characteristics of the percolation cluster makes it possible to model the functional properties and parameters of the MMMs with a known geometry of all elements. At the same time, the influence of the percolation structures/geometric parameters on properties (for example, the permeability) depends on the characteristics of the process (phenomenon). For example, the permeability of the matrix and the structure of the percolation cluster determine the magnitude of the permeability of the MMMs and the area of observed changes. That is, nonlinear changes in various properties of MMMs are associated with the appearance of a percolation structure in the matrix when the contribution from individual particles becomes macroscopic (Expression 1). At the same time, the destruction of the cluster or a decrease in its capacity leads to a decrease in the effect of the cluster and the contribution of the matrix becomes decisive.
These processes lead to the formation of unstable conditions at certain points of the matrix and the percolation cluster becomes less “powerful”. That is, by reducing the thick-ness of the matrix, we can observe a deterioration in transport properties at a fixed concentration of particles, compared with the percolation cluster in thicker films. It is also possible to expect constant permeability with a change in concentration, in the case of a transition to a homogeneous region of particle concentrations. At the same time, at high concentrations of nanoparticles, a process of increasing the size of particle agglomerates can be observed, which leads to the degradation of the percolation structure.