International Journal of Polymer Science
Volume 2018 (2018), Article ID 4903904, 9 pages
https://doi.org/10.1155/2018/4903904

Research Article

Experimental Study on the Transport of Light Gas Molecules through Low-Density Polyethylene Films

Riccardo Checchetto

Dipartimento di Fisica, Università di Trento, via Sommarive 14, Povo, 38123 Trento, Italy

Correspondence should be addressed to Riccardo Checchetto; riccardo.checchetto@unitn.it

Received 10 July 2018; Accepted 30 August 2018; Published 25 October 2018

Academic Editor: Miriam H. Rafailovich

Copyright © 2018 Riccardo Checchetto. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

An original experimental procedure for the study of gas permeation process through thin polymer films is presented. Employing mass spectroscopy techniques, this procedure allows the detection of the permeation flux with a signal-to-noise ratio large enough to obtain accurate measurements of the gas diffusivity also in processes with transient transport conditions lasting for short-interval times (~few seconds). The procedure is validated using as test material a thin low-density polyethylene (LDPE) film: the transport of four test gases with different molecular sizes and condensation properties (CO₂, N₂, D₂, and He) is studied in the 295 to 350 K temperature interval. The CO₂ diffusivity values well compare with values previously obtained studying the same LDPE film samples by integral permeation technique measuring the time lag value. Original data on the diffusivity of the He and D₂ penetrant molecules are reported: in the examined temperature range, the diffusivity values of these small-size penetrants are in the 10⁻⁶ cm²/s range and follow an Arrhenius behavior with temperature. The activation energy values for diffusion are 18.8 ± 0.4 and 10.0 ± 0.4 kJ/mol for D₂ and He, respectively.

1. Introduction

The study of mass transport process through polymer films is an important research area given its applicative interest for the development of innovative membranes for gas separation [1] and of gas barrier films to be used in the packaging technology [2]. This topic also presents fundamental interest: gas diffusion occurs by successive jumps of the penetrant molecules between cavities formed in the polymer matrix by the random thermal fluctuations of the macromolecular chains [3, 4]. These cavities constitute the free volume of the polymer which controls structural properties such as the glass transition temperature [5]. Penetrant molecules thus act as sensitive probe to explore the polymer-free volume and its variation, for example, with temperature and cross-linking density [6, 7] or with the addition of filler nanoparticles [8], and to analyze the structure of rigid cavities in biopolymer membranes [9].

The gas transport through polymer membranes obeys to the solution-diffusion mechanism [1] which assumes as transport parameters the molecular diffusivity (units: cm²/s) and solubility (units: mol·cm⁻³·Pa⁻¹) in the polymer layers (the gas permeability is given by the product of gas solubility and diffusivity, ). Measurements are carried out by gas-phase permeation technique often using planar membrane samples shaped in the form of a thin disc with a diameter much larger than the disc thickness . Keeping the sample temperature constant, at time (which is the trigger time of the experiment), one side of the planar membrane is exposed to the test gas kept at fixed pressure . This side of the membrane is called the high-pressure side (HPS). The other side is called the low-pressure side (LPS) and faces an evacuated analysis chamber kept in vacuum conditions. The gas permeation process consists in the transfer of gas molecules from the HPS of the membrane to the analysis chamber through the membrane layers and occurs in three microscopic steps:

(i)Absorption: penetrant molecules in the gas phase are absorbed by the HPS surface layers of the membrane. Here, their concentration immediately reaches its equilibrium value . Here, represents the coordinate of the HPS(ii)Diffusion step: molecules absorbed in the HPS of the membrane diffuse through the membrane layers to the LPS down to their concentration gradient according to the first Fick law, ; is the gas concentration at time in the membrane layers at depth below the HPS(iii)Desorption step: molecules reaching the LPS surface layers of the membrane desorb in the analysis chamber. When , the approximation holds ( is the coordinate of the LPS, that is, the membrane thickness)

The permeation flux can be obtained solving Fick’s second equation.

The previous diffusion problem admits for the following analytical solution:

Steady-state transport conditions are set when assumes the constant, time-stationary value .

Gas molecules permeate through the membrane of surface area at a rate . Molecules permeated in the analysis chamber form a rarefied gas at thermal equilibrium with the chamber walls at temperature .

The partial pressure of the test gas varies, as a function of time, according to the following relation:

where is the universal gas constant, is the volume of the analysis chamber, and is the pumping speed of the vacuum system.

During the permeation experiments, the pressure in the analysis chamber is measured as a function of time: the analysis of the curve by (2) permits to measure the permeation flux and to obtain the values of the and transport parameters.

Gas-phase permeation experiments are carried out by integral methods [10, 11]: at , the pumping system is excluded (i.e., ) and permeated molecules accumulate in the analysis chamber. According to the solution-diffusion model, at time , the total number of permeated moles is given by [12]

and the partial pressure of the test gas in the analysis chamber increases as a function of time as .

In the integral permeation experiments, the pressure in the analysis chamber can be recorded by using low-pressure gauges (in fact, and low-vacuum conditions are rapidly set). Gas transport parameters are obtained by the analysis of the experimental permeation curves when steady-state transport conditions are set, that is, when the value linearly increases with time .

According to the model, in steady state conditions , thus, and

The measured rate of pressure increase provides the gas permeability by the relation . The intercept of the straight line fitting the measured pressure data with the time axis, according to (4), defines a parameter called time lag given by : this characteristic time is a measure of the lasting of transient transport conditions and permits the evaluation of the gas diffusivity .

The accuracy of this procedure in the evaluation of the gas transport parameter is limited by the background flow due to degassing from the chamber walls (H₂O and CO adsorbed during the venting procedure), to small vacuum leaks (H₂O vapor, N₂, and O₂), and to air permeation through Viton O-rings (H₂O vapor, N₂, and O₂). During the permeation test, increases the pressure of the gases forming the residual vacuum: the pressure measured in the analysis chamber is thus given by two contributions which cannot be resolved by the total pressure gauge.

For a given sample, the accuracy in the evaluation of the gas permeability rises by improving the signal-to-noise ratio , that is decreasing the sample thickness . The term affects the accuracy in the measure of and introduces large relative indetermination in the evaluation of the parameter which impedes to obtain accurate measures of the gas diffusivity studying the fast transport process as it occurs with thin film samples or analyzing the transport of penetrants with high values or membrane samples made of low-permeability materials [13–16]. Note also that in permeation processes with penetrant diffusivities in the 10⁻⁶ cm²/s order, using film samples with thickness in the 10 μm range, time lag values are in the few-second order in which the result is thus comparable to that of the response time of the low-pressure gauges.

In this paper, we present an original experimental procedure for the study of the gas permeation process which permits accurate evaluation of the gas diffusivity (and permeability) of thin polymer films and processes characterized by short transient times: this often occurs in the permeation of small-size penetrants and in high-temperature processes. The validation of the proposed experimental approach will be carried out using thin films made of low-density polyethylene (LDPE) which is a polymer used as gas barrier in packaging applications [17].

Measurements were carried out using test gases with different molecular sizes, molecular masses, and condensation properties and analyzing the transport process up to temperatures close to the LDPE melting temperature. Our measurements permitted us to obtain original data on the diffusivity of the light and small-size D₂ and He penetrant molecules, and analysis of the diffusivities permitted the evaluation of the D₂ and He activation energy for diffusion which are original data here presented.

2. Experimental Part

The study was carried out using membrane samples made of commercial semicrystalline low-density polyethylene (LDPE) kindly provided by ICBM-CNR (Naples, Italy). The X-ray diffraction pattern (XRD) of the LDPE films is presented in Figure 1. The spectrum shows sharp, narrow diffraction peaks due to the crystalline fraction of the material at 2θ = 21.1°, 2θ = 23.8°, and 2θ = 36° (see vertical arrows) while the amorphous component gives a very broad peak. The Bragg-Brentano analysis of the XRD peak at 2θ = 21.1° indicated that the average size of PE crystallites is about 12 nm. The crystalline fraction of the LDPE samples, as evaluated by DSC analysis, was 29 ± 1%. The melting temperature of the present LDPE samples was ~380 K [18].

Figure 1: XRD spectrum of the LDPE film (Cu Kα radiation,  Å). The arrow indicates the LDPE reflection peaks.

The study of the gas transport properties of the present samples was previously carried out by the procedure described in the introduction, and measurements were carried out in the 290 to 350 K temperature range using H₂, CO₂, and N₂ [18]. Reliable measurements of the permeability and diffusivity values were possible only with the CO₂ test gas. Diffusivity values ranged from ~1 × 10⁻⁷ cm²/s at 290 K to ~1 × 10⁻⁶ cm²/s at 350 K while the activation energy values for diffusion and permeation were 29.7 ± 0.4 kJ/mol and 36.4 ± 0.4 kJ/mol, respectively. Using the H₂ and N₂ test gases, it was possible to evaluate only the gas permeability because of the large relative indetermination in the short time lag values [18].

In the procedure, we are here describing that the permeation experiment is carried out under continuous pumping of the analysis chamber which is done by turbomolecular or sputter-ion pumps: high vacuum conditions are thus maintained during the permeation test.

When the condition

is satisfied, then according to (2), the permeation flux is directly proportional to the instantaneous value, , of the test gas partial pressure

In the present vacuum system, and the condition are satisfied for each test gas in the examined temperature range.

The curves can be thus converted to experimental permeation flux curves through the relation

Following this procedure,

(i)the measured value of the permeation flux in stationary transport conditions permits the evaluation of the penetrant permeability by the relation (ii)the penetrant diffusivity can be obtained fitting the normalized permeation curves with the curve ; see (1)

Let us discuss now the advantages of this “constant-pumping” procedure. The continuous pumping of the analysis chamber during the experimental run maintains ultra high- (UHV) or high-vacuum (HV) conditions and, consequently, permits the use of a quadrupole mass spectrometer to detect the permeation signal that is the ionic current pertinent to the test gas. Knowing the QMS sensitivity for the test gas, , the signal provides a measure of the partial pressure of the permeating gas, The value is thus obtained, as a function of time, without noise contributions due to background flows of any other residual gas component.

This experimental procedure consequently allows

(i)very large the signal-to-noise ratio . In fact, where the noise is due to the fluctuations of the background partial pressure of the test gas, only. The is very large for all gases, in particular for those not forming the residual vacuum of the analysis chamber(ii)very low detection limit of the permeation flux. For each test gas, in fact, the term defines the detection limit (A/Pa). Consequently, the permeation signal can be efficiently detected also studying using low-permeability materials and the permeation fluxes can be detected also when (iii)that the high sampling rate of the QMS, adjustable from 0.5 ms/amu to 60 s/amu (amu: atomic mass unit), permits to record data with high time resolution

To test the effectiveness of this “constant-pumping” (CP) approach and for its validation, we used disc-shaped, planar samples with 13.5 ± 0.2 mm diameter and thickness  μm. Measurements were carried out in the 295 to 350 K temperature range with values of 40 kPa using test gas with different molecular sizes and condensation properties: carbon dioxide (CO₂), nitrogen (N₂), deuterium (D₂), and helium (He); see Table 1.

Table 1: Chemical-physical properties of the examined test gases [19].

Gas
Critical temperature (K)
Kinetic diameter (Å)
Molecular mass (amu)

D₂
38.2
2.89
4

He
5.19
2.6
4

CO₂
304.2
3.3
44

N₂
126.2
3.64
28

Experiments were carried out in a stainless steel analysis chamber pumped by turbomolecular pumps; the detection of the permeation signal is carried out by a QMS equipped with secondary electron multiplier (SEM) and grid-type ion source: this setup permits to detect partial pressures down to 10⁻¹² Pa.

A schematic diagram of the system is shown in Figure 2(a); for the sake of simplicity, in this figure, the tubular oven used for sample heating is not reported. In Figure 2(b), we present the results of the QMS calibration procedure to evaluate the sensitivity for the examined test gases. Starting from background pressure levels lower than 10⁻⁶ Pa, through a calibrated leak valve connected to the gas reservoir, the pressure of the test gas in the pumped analysis chamber is increased by a factor~10 and the corresponding QMS and IG signals were recorded. Repeating this procedure, couples of data were obtained and the slope of the line fitting them provides the value of the parameter. We obtained the following values  Pa/A,  Pa/A,  Pa/A, and  Pa/A.

(a)

(b)

Figure 2: (a) Schematic of the experimental apparatus. QMS: quadrupole mass spectrometer; HVS: ultra high-vacuum pumping system; VS1: medium vacuum pumping system (diaphragm pump and molecular drag pump); P1: ionization gauge; P2: low-pressure gauge; V1, V2, and V3: vacuum valves; LV: calibrated leak valve. The membrane sample (M) is vacuum sealed using two Viton O-rings. In this schematic diagram, the tubular oven for membrane sample heating is not reported. (b) QMS calibration curves. The QMS mass signal is reported as a function of the pressure value for the different test gases. Indetermination is inside the size of the symbol.

In Figure 3(a), we present the mass spectrum of the background vacuum at  1.2 × 10⁻⁶ Pa, after one night of turbomolecular pumping: the spectrum reports the QMS mass signal pertinent to the gases forming the residual vacuum. It can be seen that the background vacuum is mostly composed of water vapor (H₂O⁺, m/e = 18, fragmentation products OH⁺ at m/e = 17, and O⁺ at m/e = 16), carbon monoxide (CO⁺, m/e = 28, fragmentation product C⁺ at m/e = 12), and molecular hydrogen (H₂⁺, m/e = 2, fragmentation product H⁺ at m/e = 1) [19]. In Figure 3(b), we report the QMS mass signal for N₂ (upper panel) and D₂ (lower panel) recorded in the permeation test carried out at 295 K with the described “constant-pumping” procedure: vertical arrows indicate the trigger time .

(a)

(b)

Figure 3: (a) Mass spectrum of the background vacuum at  Pa after one night of turbomolecular pumping. The dominant contributions are related to the signals pertinent to H₂O, CO, and H₂ (m/e = 18, 28, and 2, respectively). Their fragmentation products are also reported; see text. (b) QMS mass signals and (m/e = 4 and 28, respectively) recorded in the permeation test. Arrows indicate the trigger time of the experiment.

Looking at the curves for , we observe that the QMS background signals are in the 10⁻¹² and 10⁻¹⁰ A order for D₂ and N₂ while the fluctuations are ~2 × 10⁻¹² A for both gases: the is thus larger than 10 already in the initial instants of the permeation transient () and reaches values in the 10³ order in stationary transport conditions.

The fluctuations in the 10⁻¹² A order permit to evaluate the detection limit of the present experimental approach by the relation . Given the effective values in the 50 L/s order, QMS sensitivity , see Figure 1(b), and membrane surface area , then is in the 10⁻¹⁴ mol·cm⁻²·s⁻¹ order (;  L·kPa·mol−1·K⁻¹): well compare with the upper value of the oxygen transmission rate (OTR) for gas barrier films to be used in the organic light-emitting device (OLED) technology, ~3 × 10⁻¹⁴ mol·cm⁻²·s⁻¹ [20].

The permeation flux is obtained from the curves in Figure 3(a) subtracting the background QMS signal from the detected QMS signal and using the relation

3. Results and Discussion

In Figure 4(a), we present the curves as obtained at 295 K exposing the HPS at 40 kPa pressure. In the figure, symbols are the experimental data (relative indeterminations are inside the size of the symbol) obtained by (8) while continuous lines are their numerical fit. All permeation curves are characterized by an initial interval time where penetrant occurs in transient conditions; the permeation flux then reaches a constant value indicating that stationary transport conditions are set. Fits were obtained first measuring the value of the permeation flux in stationary transport condition relation to evaluate the gas permeability value and then fitting the normalized curve with the function to evaluate the gas diffusivity .

(a)

(b)

Figure 4: (a) Permeation flux curves (open symbols) pertinent to CO₂, N₂, D₂, and He obtained at  K and 40 kPa. Lines are fitting of the experimental curves; see text. (b) Normalized curves (open symbols) pertinent to D₂ and He obtained at  K and  kPa. Lines are the curve fittings by the function; see text.

Comparing the curves in Figure 4(a), we observe that transient transport conditions last for longer interval times with the CO₂ and N₂ penetrants than with D₂ and He which have smaller molecular size. The difference between the permeation curves in transient transport conditions pertinent to the small-size test molecules, namely, D₂ and He, can be appreciated by their plots in Figure 4(b) where the normalized permeation flux curves are presented.

In Figures 5(a)–5(d), we present the experimental permeation curves obtained at temperatures between 295 and 350 K for CO₂, N₂, D₂, and He with 40 kPa (symbols): the continuous lines are their fitting by the function. Note that the high signal-to-noise ratio of the present procedure permits to clearly distinguish the permeation curves obtained in the 295 to 350 K temperature range also using the small-size penetrants.

(a)

(b)

(c)

(d)

Figure 5: (a) Experimental curves (symbols) pertinent to CO₂ obtained at temperatures between  K and  K with  kPa. Lines are fitting of the experimental curves; see text. (b) Experimental curves (symbols) pertinent to N₂ obtained at temperatures between  K and  K with  kPa. Lines are fitting of the experimental curves; see text. (c) Experimental curves (symbols) pertinent to D₂ obtained at temperatures between  K and  K with  kPa. Lines are fitting of the experimental curves; see text. (d) Experimental curves (symbols) pertinent to He. To avoid confusion, in the figure, we only report curves obtained at 295, 330, and 351 K with  kPa. Lines are fitting of the experimental curves; see text.

The Arrhenius plots of the permeability and diffusivity values are reported in Figures 6(a) and 6(b), respectively. Permeability values are reported in practical Barrer units

(a)

(b)

Figure 6: (a) Arrhenius plot of the gas permeability as obtained by analysis of the permeation flux value in stationary transport conditions. Permeability values are presented in practical units (Barrer). (b) Arrhenius plot of the gas diffusivity as obtained by the fitting of curves; see Figures 5(a)–5(d).

The activation energy value for permeation and the apparent value of the activation energy for diffusion were obtained fitting the permeability and diffusivity data by the following relations for the gas diffusivity and permeability

and values of the examined test gases are reported in Table 2.

Table 2: Activation energy for permeation and activation energy for diffusion as obtained by the Arrhenius plot of the permeability and diffusivity data in Figures 6(a) and 6(b).

Test gas
(kJ/mol)
(kJ/mol)

CO₂
36.5 ± 0.5
34.0 ± 0.3

N₂
43.2 ± 0.8
38.2 ± 0.5

D₂
33.2 ± 0.8
18.8 ± 0.4

He
37.2 ± 0.7
10.0 ± 0.4

It is worthy to observe that the diffusivity values have been obtained by the analysis of the whole permeation curve which is composed, given the high sampling rate of the QMS, by a large ensemble of experimental data. This procedure thus allows more accurate evaluation of the gas diffusivity than using the value of a single parameter, as in the time lag method.

The obtained diffusivity and permeability values for CO₂ as well as the pertinent activation energies are in well agreement with values previously obtained by the “constant-volume” approach [18] and with other literature data; see references therein. The same consideration holds for the permeability values and activation energy values for H₂ and N₂.

To conclude this discussion, it is worthy to comment the obtained results considering the physical-chemical properties of the employed test gases. Considering the diffusion data in Figure 6(b) and the penetrant properties in Table 1, we conclude that the penetrant diffusivity through the present LDPE samples decreases by increasing the molecular size. Moreover, comparing the permeability data in Figure 6(a) and the molecular physical-chemical properties in Table 1, we conclude that the gas permeability of the present LDPE samples increases with the critical temperature of the test gas evidencing the important role of the gas condensation properties. In fact, CO₂ and N₂ have similar molecular size (see Table 1) and diffusivity (see Figure 6(b)) but the CO₂ permeability is larger than that of N₂. Note that the same considerations hold also comparing the H₂ and He transport parameters.

This scaling can be discussed considering the theoretical results presented in a recent paper of Dutta and Bhatia on the free volume structure and gas transport properties of polyethylene carried out by equilibrium molecular dynamic simulations [21]. The CO₂ and N₂ self-diffusivity values indicated by the authors are in the 10⁻⁶ cm²/s range, in accord with the present experimental evaluation. They also suggested that the temperature dependence of the penetrant self-diffusivity follows the Arrhenius behavior, as also observed in the present study. Simulations also revealed that the free volume of PE is composed by voids with size ranging from 1.5 to 3 Å. The larger D₂ and He diffusivity compared with CO₂ and N₂ can be explained as consequence of the availability of diffusion sites with a diameter close to the D₂ and He molecular sizes and thus is more efficient in accommodating these small-size penetrants.

4. Conclusions

We have described a procedure for the experimental study of gas permeation through polymer films. The procedure is based on the measure of the permeation signal in a continuously pumped analysis chamber by a calibrated quadrupole mass spectrometer. This procedure allows (i) very low detection limit of the permeation flux, ~10⁻¹⁴ mol·cm⁻²·s⁻¹, and (ii) high signal-to-noise ratio in transient transport conditions also when studying polymer films having gas barrier properties. The improved performances permit to resolve the kinetics of the He and D₂ transport processes, to obtain their diffusivity values up to temperatures close to those of the LDPE melting and the pertinent activation energies of diffusion. This procedure can be applied for the study of the gas transport processes with the recently developed polymer materials for packaging applications to be applied in fields such as OLED technology. It is also of interest for the community of polymer scientists interested to the study of the correlation between the polymer-free volume structure and transport properties.

Data Availability

Experimental data supporting the described results can be obtained upon request to the corresponding author.

Conflicts of Interest

The author declares that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The help of M. Zanolli in the gas transport experiments is gratefully acknowledged.

References

S. Matteucci, Y. Yampolskii, B. D. Freeman, and I. Pinneau, Materials Science of Membranes for Gas and Vapor Separation, John Wiley and Sons, Chichester, 2006.

A. Lopez-Rubio, E. Almenar, P. Hernandez-Munoz, R. Catala, and R. Gavara, “Overview of active polymer-based packaging technologies for food applications,” Food Reviews International, vol. 20, no. 4, pp. 357–387, 2004.

J. S. Vrentas and J. M. Vrentas, “Evaluation of the free volume theory of diffusion,” Journal of Polymer Science Part B: Polymer Physics, vol. 41, no. 5, pp. 501–507, 2003.

H. Fujiita, “Notes on free volume theories,” Polymer Journal, vol. 23, no. 12, pp. 1499–1506, 1991.

R. P. White and J. E. G. Lipson, “Polymer free volume and its connection to the glass transition,” Macromolecules, vol. 49, no. 11, pp. 3987–4007, 2016.

P. N. Patil, D. Roilo, R. S. Brusa, A. Miotello, and R. Checchetto, “Influence of nano-level molecular packing on the gas transport properties in amine-modified epoxy resins,” Polymer, vol. 58, pp. 130–138, 2015.

P. N. Patil, D. Roilo, R. S. Brusa et al., “Free volumes and gas transport in polymers: amine-modified epoxy resins as a case study,” Physical Chemistry Chemical Physics, vol. 18, no. 5, pp. 3817–3824, 2016.

D. Roilo, P. N. Patil, R. S. Brusa et al., “Polymer rigidification in graphene based nanocomposites: gas barrier effects and free volume reduction,” Polymer, vol. 121, pp. 17–25, 2017.

D. Roilo, C. A. Maestri, M. Scarpa et al., “Cellulose nanofibrils films: molecular diffusion through elongated nanocavities,” Journal of Physical Chemistry C, vol. 121, no. 28, pp. 15437–15447, 2017.

B. Flaconnèche, J. Martin, and M. H. Klopffer, “Transport properties of gases in polymers: experimental methods,” Oil & Gas Science and Technology, vol. 56, no. 3, pp. 245–259, 2001.

R. M. Felder, “Estimation of gas transport coefficients from differential permeation, integral permeation, and sorption rate data,” Journal of Membrane Science, vol. 3, no. 1, pp. 15–27, 1978.

J. Crank, The Mathematics of Diffusion, Claredon, Oxford, 2nd edition, 1999.

W. H. Lin and T. S. Chung, “Gas permeability, diffusivity, solubility, and aging characteristic of 6FDA-durene polyimide membranes,” Journal of Membrane Science, vol. 186, no. 2, pp. 183–193, 2011.

J. Ahn, W. J. Chung, I. Pinneau, and M. D. Guiver, “Polysulfone/silica nanoparticle mixed-matrix membranes for gas separation,” Journal of Membrane Science, vol. 314, no. 1-2, pp. 123–133, 2008.

J. P. G. Villaluenga and B. Seoane, “Experimental estimation of gas-transport properties of linear low-density polyethylene membranes by an integral permeation method,” Journal of Applied Polymer Science, vol. 82, no. 12, pp. 3013–3021, 2001.

W. J. Koros, A. H. Chan, and D. R. Paul, “Sorption and transport of various gases in polycarbonate,” Journal of Membrane Science, vol. 2, pp. 165–190, 1977.

K. Marsch and B. Bugusu, “Food packaging roles, materials and environmental issues,” Journal of Food Science, vol. 72, no. 3, pp. R39–R55, 2007.

R. Checchetto, A. Miotello, L. Nicolais, and G. Carotenuto, “Gas transport through nanocomposite membrane composed by polyethylene with dispersed graphite nanoplatelets,” Journal of Membrane Science, vol. 463, pp. 196–204, 2014.

P. A. Redhead, J. P. Hobson, and E. V. Kornelsen, The Physical Basis of Ultrahigh Vacuum, American Institute of Physics, New York, NY, USA, 1993.

J. S. Lewis and M. S. Weaver, “Thin-film permeation-barrier technology for flexible organic light-emitting devices,” IEEE Journal of Selected Topics in Quantum Electronics, vol. 10, no. 1, pp. 45–57, 2004.

R. C. Dutta and S. K. Bhatia, “Transport diffusion of light gases in polyethylene using atomistic simulations,” Langmuir, vol. 33, no. 4, pp. 936–946, 2017.